NMR and molecular dynamics studies of the hydration of a zinc finger-DNA complex

Molecular dynamics modelling of the interaction of a synthetic zinc-finger miniprotein with DNA

We report the modelling of the DNA complex of an artificial miniprotein composed of two zinc finger modules and an AT-hook linking peptide. The computational study provides for the first time a structural view of these types of complexes, dissecting interactions that are key to modulate their stability. The relevance of these interactions was validated experimentally. These results confirm the potential of this type of computational approach for studying peptide-DNA complexes and suggest that they could be very useful for the rational design of non-natural, DNA binding miniproteins.

Synthesis and evaluation of 1,2,3-dithiazole inhibitors of the nucleocapsid protein of feline immunodeficiency virus (FIV) as a model for HIV infection

We disclose a series of potent anti-viral 1,2,3-dithiazoles, accessed through a succinct synthetic approach from 4,5-dichloro-1,2,3-dithiazolium chloride (Appel's salt). A series of small libraries of compounds were screened against feline immunodeficiency virus (FIV) infected cells as a model for HIV. This approach highlighted new structure activity relationship understanding and led to the development of sub-micro molar anti-viral compounds with reduced toxicity. In addition, insight into the mechanistic progress of this system is provided via advanced QM-MM modelling. The 1,2,3-dithiazole represents a versatile scaffold with potential for further development to treat both FIV and HIV.

Molecular Dynamics Simulations of Reduced and Oxidized TFIIIA Zinc Fingers Free and Interacting with 5S RNA

Interactions of zinc finger (ZF) proteins with nucleic acids and proteins play an important role in DNA transcription and repair, biochemical recognition, and protein regulation. The release of Zn²⁺ through oxidation of cysteine thiolates is associated with disruption of gene expression and DNA repair, preventing tumor growth. Multi-microsecond molecular dynamics (MD) simulations were carried out to examine the effect of Cys oxidation on the ZF456 fragment of transcription factor III A (TFIIIA) and its complex with 5S RNA. In the absence of 5S RNA, the reduced ZF456 peptide undergoes conformational changes in the secondary structure due to the reorientation of the intact ZF domains. Upon oxidation, the individual ZF domains unfold to various degrees, yielding a globular ZF456 peptide with ZF4 and ZF6, responsible for base-specific hydrogen bonds with 5S RNA, losing their ββα-folds. ZF5, on the other hand, participates in nonspecific interactions through its α-helix that conditionally unravels early in the simulation. In the presence of RNA, oxidation of the ZF456 peptide disrupts the key hydrogen bonding interactions between ZF5/ZF6 and 5S RNA. However, interactions with ZF4 are dependent on the protonation state of His119.

Specifically bound BZIP transcription factors modulate DNA supercoiling transitions

Torsional stress on DNA, introduced by molecular motors, constitutes an important regulatory mechanism of transcriptional control. Torsional stress can modulate specific binding of transcription factors to DNA and introduce local conformational changes that facilitate the opening of promoters and nucleosome remodelling. Using all-atom microsecond scale molecular dynamics simulations together with a torsional restraint that controls the total twist of a DNA fragment, we address the impact of torsional stress on DNA complexation with a human BZIP transcription factor, MafB. We gradually over- and underwind DNA alone and in complex with MafB by 0.5° per dinucleotide step, starting from the relaxed state to a maximum of 5° per dinucleotide step, monitoring the evolution of the protein-DNA contacts at different degrees of torsional strain. Our computations show that MafB changes the DNA sequence-specific response to torsional stress. The dinucleotide steps that are susceptible to absorbing most of the torsional stress become more torsionally rigid, as they are involved in protein-DNA contacts. Also, the protein undergoes substantial conformational changes to follow the stress-induced DNA deformation, but mostly maintains the specific contacts with DNA. This results in a significant asymmetric increase of free energy of DNA twisting transitions, relative to free DNA, where overtwisting is more energetically unfavourable. Our data suggest that specifically bound BZIP factors could act as torsional stress insulators, modulating the propagation of torsional stress along the chromatin fibre, which might promote cooperative binding of collaborative DNA-binding factors.

Molecular Dynamics Simulations of Dna and Its Complexes

This article describes how classical molecular simulation methods are being used to gain a molecular-level understanding of the interaction mechanisms responsible for DNA–ligand recognition, and that govern the response of DNA to ligand binding. Case studies using a variety of different ligands—including small pharmaceutical drugs, proteins and lipids—are used to illustrate the power of modern molecular dynamics simulation methods for understanding how we may control the function and structure of DNA.

Modulation of Solvation and Molecular Recognition of a Lipid Bilayer under Dynamical Phase Transition

It is well accepted in contemporary biology that an ∼30 Å thick lipid bilayer film around living cells is a matter of life and death as the film typically delimits the environments that serve as a crucial margin. The dynamic organization of lipid molecules both across the lipid bilayer and in the lateral dimension are known to be crucial for cellular transport and molecular recognition by important biological macromolecules. Here, we study dilute (20 mM) Dioctadecyldimethylammonium bromide (DODAB) vesicles at different temperatures in aqueous dispersion with well-defined phases namely liquid crystalline, gel and subgel. The spectroscopic studies on two fluorescent probes 8-anilino-1-naphthalene sulfonic acid ammonium salt (ANS) and Coumarin 500 (C500), former in the head group region of the lipid-water interface and later located deeper in the lipid bilayer follow dynamics (solvation and fluidity) of their local environments in the vesicles. Binding of an anti-tuberculosis drug rifampicin has also been studied employing Förster resonance energy transfer (FRET) technique. The molecular insight concerning the effect of dynamical organization of the lipid molecules on the local dynamics of aqueous environments in different phases leading to molecular recognition becomes evident in our study.

Lysine Side-Chain Dynamics in the Binding Site of Homeodomain/DNA Complexes As Observed by NMR Relaxation Experiments and Molecular Dynamics Simulations

An important but poorly characterized contribution to the thermodynamics of protein-DNA interactions is the loss of entropy that occurs from restricting the conformational freedom of amino acid side chains. The effect of restricting the flexibility of several side chains at a protein-DNA interface may be comparable in many cases to the other factors that determine the binding thermodynamics and may, therefore, play a key role in dictating the binding affinity and/or specificity. Because the entropic contributions, including the presence and influence of side-chain dynamics, are especially difficult to estimate based on structural information, it is important to pursue experimental and theoretical studies that can provide direct information regarding these issues. We report on studies of a model system, the homeodomain/DNA complex, focusing on the Lys50 class of homeodomains where a key lysine residue in position 50 was shown previously to be critical for binding site specificity. NMR methodology was employed for determining the dynamics of lysine side-chain amino groups via 15N relaxation measurements in the Lys50-class homeodomains from the Drosophila protein Bicoid and the human protein Pitx2. In the case of Pitx2, complexes with both a consensus and a nonconsensus DNA binding site were examined. NMR-derived order parameters indicated moderate to substantial conformational freedom for the lysine NH3+ group in the complexes studied. To complement the experimental NMR measurements, molecular dynamics simulations were performed for the consensus complexes to gain further, detailed insights regarding the dynamics of the Lys50 side chain and other important residues in the protein-DNA interface.

Zinc finger protein binding to DNA: an energy perspective using molecular dynamics simulation and free energy calculations on mutants of both zinc finger domains and their specific DNA bases

Energy calculations based on MM-GBSA were employed to study various zinc finger protein (ZF) motifs binding to DNA. Mutants of both the DNA bound to their specific amino acids were studied. Calculated energies gave evidence for a relationship between binding energy and affinity of ZF motifs to their sites on DNA. ΔG values were -15.82(12), -3.66(12), and -12.14(11.6) kcal/mol for finger one, finger two, and finger three, respectively. The mutations in the DNA bases reduced the value of the negative energies of binding (maximum value for ΔΔG = 42Kcal/mol for F1 when GCG mutated to GGG, and ΔΔG = 22 kcal/mol for F2, the loss in total energy of binding originated in the loss in electrostatic energies upon mutation (r = .98). The mutations in key amino acids in the ZF motif in positions-1, 2, 3, and 6 showed reduced binding energies to DNA with correlation coefficients between total free energy and electrostatic was .99 and with Van der Waal was .93. Results agree with experimentally found selectivity which showed that Arginine in position-1 is specific to G, while Aspartic acid (D) in position 2 plays a complicated role in binding. There is a correlation between the MD calculated free energies of binding and those obtained experimentally for prepared ZF motifs bound to triplet bases in other reports (), our results may help in the design of ZF motifs based on the established recognition codes based on energies and contributing energies to the total energy.

Direct probing of solvent accessibility and mobility at the binding interface of polymerase (Dpo4)-DNA complex

Water plays essential structural and dynamical roles in protein-DNA recognition through contributing to enthalpic or entropic stabilization of binding complex and by mediating intermolecular interactions and fluctuations for biological function. These interfacial water molecules are confined by the binding partners in nanospace, but in many cases they are highly mobile and exchange with outside bulk solution. Here, we report our studies of the interfacial water dynamics in the binary and ternary complexes of a polymerase (Dpo4) with DNA and an incoming nucleotide using a site-specific tryptophan probe with femtosecond resolution. By systematic comparison of the interfacial water motions and local side chain fluctuations in the apo, binary, and ternary states of Dpo4, we observed that the DNA binding interface and active site are dynamically solvent accessible and the interfacial water dynamics are similar to the surface hydration water fluctuations on picosecond time scales. Our molecular dynamics simulations also show the binding interface full of water molecules and nonspecific weak interactions. Such a fluid binding interface facilitates the polymerase sliding on DNA for fast translocation whereas the spacious and mobile hydrated active site contributes to the low fidelity of the lesion-bypass Y-family DNA polymerase.

Localized frustration and binding-induced conformational change in recognition of 5S RNA by TFIIIA zinc finger

Protein TFIIIA is composed of nine tandemly arranged Cys2His2 zinc fingers. It can bind either to the 5S RNA gene as a transcription factor or to the 5S RNA transcript as a chaperone. Although structural and biochemical data provided valuable information on the recognition between the TFIIIIA and the 5S DNA/RNA, the involved conformational motions and energetic factors contributing to the binding affinity and specificity remain unclear. In this work, we conducted MD simulations and MM/GBSA calculations to investigate the binding-induced conformational changes in the recognition of the 5S RNA by the central three zinc fingers of TFIIIA and the energetic factors that influence the binding affinity and specificity at an atomistic level. Our results revealed drastic interdomain conformational changes between these three zinc fingers, involving the exposure/burial of several crucial DNA/RNA binding residues, which can be related to the competition between DNA and RNA for the binding of TFIIIA. We also showed that the specific recognition between finger 4/finger 6 and the 5S RNA introduces frustrations to the nonspecific interactions between finger 5 and the 5S RNA, which may be important to achieve optimal binding affinity and specificity.

Study of base pair mutations in proline-rich homeodomain (PRH)-DNA complexes using molecular dynamics

Proline-rich homeodomain (PRH) is a regulatory protein controlling transcription and gene expression processes by binding to the specific sequence of DNA, especially to the sequence 5'-TAATNN-3'. The impact of base pair mutations on the binding between the PRH protein and DNA is investigated using molecular dynamics and free energy simulations to identify DNA sequences that form stable complexes with PRH. Three 20-ns molecular dynamics simulations (PRH-TAATTG, PRH-TAATTA and PRH-TAATGG complexes) in explicit solvent water were performed to investigate three complexes structurally. Structural analysis shows that the native TAATTG sequence forms a complex that is more stable than complexes with base pair mutations. It is also observed that upon mutation, the number and occupancy of the direct and water-mediated hydrogen bonds decrease. Free energy calculations performed with the thermodynamic integration method predict relative binding free energies of 0.64 and 2 kcal/mol for GC to AT and TA to GC mutations, respectively, suggesting that among the three DNA sequences, the PRH-TAATTG complex is more stable than the two mutated complexes. In addition, it is demonstrated that the stability of the PRH-TAATTA complex is greater than that of the PRH-TAATGG complex.

Structural analysis of zinc-finger (TTK) + [Cu(BPA)]2+ /[Cu(IDB)]2+ + DNA complexes: an investigation by molecular dynamics simulation

In the present study, the molecular dynamics simulation technique is employed to investigate the hydrogen abstraction possibility from sugar of DNA in two designed complexes of copper-based chemical nuclease [Cu(BPA)](2+) bis(2-pyridylmethyl) amine (BPA) or [Cu(IDB)](2+) N,N-bis(2-benzimidazolylmethyl) amine (IDB) bound to the zinc finger protein Tramtrack (TTK). The simulated results show that each of the designed complexes can form a stable conformation within 30 ns of simulation time with the substrate OOH(-) and an 18-base pair (bp) DNA segment and is located in the major groove of the DNA segment. The active terminal O atom of the OOH(-) substrate is found in close proximity to the target C2'H, C3'H, C4'H or C5'H proton of the DNA in TTK + [Cu(BPA)OOH](+) + DNA or TTK + [Cu(IDB)OOH](+) + DNA complex, which is crucial to propose the hydrogen abstraction possibility that is responsible for the DNA cleavage. The positions of copper-based chemical nucleases bound to TTK may substantially influence the hydrogen abstraction possibility. The structures and sizes of ligands in copper-based nucleases are also found to have influence on the order of difficulty of the hydrogen abstraction from the sugars of DNA.

Exploring translocation of proteins on DNA by NMR

While an extensive body of knowledge has accumulated on the structures of transcription factors, DNA and their complexes from both NMR and crystallography, much less is known at a molecular level regarding the mechanisms whereby transcription factors locate their specific DNA target site within an overwhelming sea of non-specific DNA sites. Indirect kinetic data suggested that three processes are involved in the search procedure: jumping by dissociation of the protein from the DNA followed by re-association at another site, direct transfer from one DNA molecule or segment to another, and one-dimensional sliding. In this brief perspective I summarize recent NMR developments from our laboratory that have permitted direct characterization of the species and molecular mechanisms involved in the target search process, including the detection of highly transient sparsely-populated states. The main tool in these studies involves the application of paramagnetic relaxation enhancement, supplemented by z-exchange spectroscopy, lineshape analysis and residual dipolar couplings. These studies led to the first direct demonstration of rotation-coupled sliding of a protein along the DNA and the direct transfer of a protein from one DNA molecule to another without dissociating into free solution.

Biomolecularmodeling and simulation: a field coming of age

We assess the progress in biomolecular modeling and simulation, focusing on structure prediction and dynamics, by presenting the field’s history, metrics for its rise in popularity, early expressed expectations, and current significant applications. The increases in computational power combined with improvements in algorithms and force fields have led to considerable success, especially in protein folding, specificity of ligand/biomolecule interactions, and interpretation of complex experimental phenomena (e.g. NMR relaxation, protein-folding kinetics and multiple conformational states) through the generation of structural hypotheses and pathway mechanisms. Although far from a general automated tool, structure prediction is notable for proteins and RNA that preceded the experiment, especially by knowledge-based approaches. Thus, despite early unrealistic expectations and the realization that computer technology alone will not quickly bridge the gap between experimental and theoretical time frames, ongoing improvements to enhance the accuracy and scope of modeling and simulation are propelling the field onto a productive trajectory to become full partner with experiment and a field on its own right.

Interaction identification of Zif268 and TATA(ZF) proteins with GC-/AT-rich DNA sequence: A theoretical study

Molecular dynamics (MD) simulations for Zif268 (a zinc-finger-protein binding specifically to the GC-rich DNA)-d(A(1) G(2) C(3) G(4) T(5) G(6) G(7) G(8) C(9) A(10) C(11) )(2) and TATA(ZF) (a zinc-finger-protein recognizing the AT-rich DNA)-d(A(1) C(2) G(3) C(4) T(5) A(6) T(7) A(8) A(9) A(10) A(11) G(12) G(13) )(2) complexes have been performed for investigating the DNA binding affinities and specific recognitions of zinc fingers to GC-rich and AT-rich DNA sequences. The binding free energies for the two systems have been further analyzed by using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. The calculations of the binding free energies reveal that the affinity energy of Zif268-DNA complex is larger than that of TATA(ZF) -DNA one. The affinity between the zinc-finger-protein and DNA is mainly driven by more favorable van-der-Waals and nonpolar/solvation interactions in both complexes. However, the affinity energy difference of the two binding systems is mainly caused by the difference of van-der-Waals interactions and entropy components. The decomposition analysis of MM-PBSA free energies on each residue of the proteins predicts that the interactions between the residues with the positive charges and DNA favor the binding process; while the interactions between the residues with the negative charges and DNA behave in the opposite way. The interhydrogen-bonds at the protein-DNA interface and the induced intrafinger hydrogen bonds between the residues of protein for the Zif268-DNA complex have been identified at some key contact sites. However, only the interhydrogen-bonds between the residues of protein and DNA for TATA(ZF) -DNA complex have been found. The interactions of hydrogen-bonds, electrostatistics and van-der-Waals type at some new contact sites have been identified. Moreover, the recognition characteristics of the two studied zinc-finger-proteins have also been discussed.

Modeling the hydration layer around proteins: HyPred

Protein hydration plays an integral role in determining protein function and stability. We develop a simple method with atomic level precision for predicting the solvent density near the surface of a protein. A set of proximal radial distribution functions are defined and calculated for a series of different atom types in proteins using all-atom, explicit solvent molecular dynamic simulations for three globular proteins. A major improvement in predicting the hydration layer is found when the protein is held immobile during the simulations. The distribution functions are used to develop a model for predicting the hydration layer with sub-1-Angstrom resolution without the need for additional simulations. The model and the distribution functions for a given protein are tested in their ability to reproduce the hydration layer from the simulations for that protein, as well as those for other proteins and for simulations in which the protein atoms are mobile. Predictions for the density of water in the hydration shells are then compared with high occupancy sites observed in crystal structures. The accuracy of both tests demonstrates that the solvation model provides a basis for quantitatively understanding protein solvation and thereby predicting the hydration layer without additional simulations.

Investigating dynamic and energetic determinants of protein nucleic acid recognition: analysis of the zinc finger zif268-DNA complexes

BACKGROUND

Protein-DNA recognition underlies fundamental biological processes ranging from transcription to replication and modification. Herein, we present a computational study of the sequence modulation of internal dynamic properties and of intraprotein networks of aminoacid interactions that determine the stability and specificity of protein-DNA complexes.

RESULTS

To this aim, we apply novel theoretical approaches to analyze the dynamics and energetics of biological systems starting from MD trajectories. As model system, we chose different sequences of Zinc Fingers (ZF) of the Zif268 family bound with different sequences of DNA. The complexes differ for their experimental stability properties, but share the same overall 3 D structure and do not undergo structural modifications during the simulations. The results of our analysis suggest that the energy landscape for DNA binding may be populated by dynamically different states, even in the absence of major conformational changes. Energetic couplings between residues change in response to protein and/or DNA sequence variations thus modulating the selectivity of recognition and the relative importance of different regions for binding.

CONCLUSIONS

The results show differences in the organization of the intra-protein energy-networks responsible for the stabilization of the protein conformations recognizing and binding DNA. These, in turn, are reflected into different modulation of the ZF's internal dynamics. The results also show a correlation between energetic and dynamic properties of the different proteins and their specificity/selectivity for DNA sequences. Finally, a dynamic and energetic model for the recognition of DNA by Zinc Fingers is proposed.

Electrostatic hot spot on DNA-binding domains mediates phosphate desolvation and the pre-organization of specificity determinant side chains

A major obstacle towards elucidating the molecular basis of transcriptional regulation is the lack of a detailed understanding of the interplay between non-specific and specific protein–DNA interactions. Based on molecular dynamics simulations of C₂H₂ zinc fingers (ZFs) and engrailed homeodomain transcription factors (TFs), we show that each of the studied DNA-binding domains has a set of highly constrained side chains in preset configurations ready to form hydrogen bonds with the DNA backbone. Interestingly, those domains that bury their recognition helix into the major groove are found to have an electrostatic hot spot for Cl⁻ ions located on the same binding cavity as the most buried DNA phosphate. The spot is characterized by three protein hydrogen bond donors, often including two basic side chains. If bound, Cl⁻ ions, likely mimicking phosphates, steer side chains that end up forming specific contacts with bases into bound-like conformations. These findings are consistent with a multi-step DNA-binding mechanism in which a pre-organized set of TF side chains assist in the desolvation of phosphates into well defined sites, prompting the re-organization of specificity determining side chains into conformations suitable for the recognition of their cognate sequence.

Physical basis of metal-binding specificity in Escherichia coli NikR

In Escherichia coli and other bacteria, nickel uptake is regulated by the transcription factor NikR. Nickel binding at high-affinity sites in E. coli NikR (EcNikR) facilitates EcNikR binding to the nik operon, where it then suppresses transcription of genes encoding the nickel uptake transporter, NikABCDE. A structure of the EcNikR-DNA complex suggests that a second metal-binding site is also present when NikR binds to the nik operon. Moreover, this co-crystal structure raises the question of what metal occupies the second site under physiological conditions: K(+), which is present in the crystal structure, or Ni(2+), which has been proposed to bind to low- as well as high-affinity sites on EcNikR. To determine which ion is preferred at the second metal-binding site and the physical basis for any preference of one ion over another in both the second metal-binding site and the high-affinity sites, we conducted a series of detailed molecular simulations on the EcNikR structure. Simulations that place Ni(2+) at high-affinity sites lead to stable trajectories with realistic ion-ligand distances and geometries, while simulations that place K(+) at these sites lead to conformational changes in the protein that are likely unfavorable for ion binding. By contrast, simulations on the second metal site in the EcNikR-DNA complex lead to stable trajectories with realistic geometries regardless of whether K(+) or Ni(2+) occupies this site. Electrostatic binding free energy calculations, however, suggest that EcNikR binding to DNA is more favorable when the second metal-binding site contains K(+). An analysis of the energetic contributions to the electrostatic binding free energy suggests that, while the interaction between EcNikR and DNA is more favorable when the second site contains Ni(2+), the large desolvation penalty associated with moving Ni(2+) from solution to the relatively buried second site offsets this favorable interaction term. Additional free energy simulations that account for both electrostatic and non-electrostatic effects argue that EcNikR binding to DNA is most favorable when the second site contains a monovalent ion the size of K(+). Taken together, these data suggest that the EcNikR structure is most stable when Ni(2+) occupies high-affinity sites and that EcNikR binding to DNA is more favorable when the second site contains K(+).

Statistical thermodynamics of biomembranes

An overview of the major issues involved in the statistical thermodynamic treatment of phospholipid membranes at the atomistic level is summarized: thermodynamic ensembles, initial configuration (or the physical system being modeled), force field representation as well as the representation of long-range interactions. This is followed by a description of the various ways that the simulated ensembles can be analyzed: area of the lipid, mass density profiles, radial distribution functions (RDFs), water orientation profile, deuterium order parameter, free energy profiles and void (pore) formation; with particular focus on the results obtained from our recent molecular dynamic (MD) simulations of phospholipids interacting with dimethylsulfoxide (Me(2)SO), a commonly used cryoprotective agent (CPA).

Assessing the performance of implicit solvation models at a nucleic acid surface

Implicit solvation models are popular alternatives to explicit solvent methods due to their ability to "pre-average" solvent behavior and thus reduce the need for computationally-expensive sampling. Previously, we have demonstrated that Poisson-Boltzmann models for polar solvation and integral-based models for nonpolar solvation can reproduce explicit solvation forces in a low-charge density protein system. In the present work, we examine the ability of these continuum models to describe solvation forces at the surface of a RNA hairpin. While these models do not completely describe all of the details of solvent behavior at this highly-charged biomolecular interface, they do provide a reasonable description of average solvation forces and therefore show significant promise for developing more robust implicit descriptions of solvent around nucleic acid systems for use in biomolecular simulation and modeling. Additionally, we observe fairly good transferability in the nonpolar model parameters optimized for protein systems, suggesting its robustness for modeling general nonpolar solvation phenomena in biomolecular systems.

Water-mediated interactions between DNA and PhoB DNA-binding/transactivation domain: NMR-restrained molecular dynamics in explicit water environment

The solution structure of the complex between the transcription factor PhoB DNA-binding/transactivation domain and DNA was determined by NMR spectroscopy and simulated annealing in a periodic boundary box of explicit water with the particle mesh Ewald method. The refined structures provided better convergence and better local geometry compared with the structures determined in vacuum. The hydrogen bond interactions between the PhoB domain and DNA in the aqueous environment were fully formed. The complex structure was found to be very similar to the crystal structure, particularly at the PhoB-DNA interface, much more so than expected from the vacuum structure. These results indicate the importance of the proper treatment of electrostatic and hydration influences in describing protein-DNA interactions. The hydration structures observed for the refined structures contained most of the crystal waters as a subset. We observed that various water-mediated PhoB-DNA interactions contributed to the molecular recognition between PhoB and DNA.

WATGEN: an algorithm for modeling water networks at protein-protein interfaces

Water molecules at protein-protein interfaces contribute to the close packing of atoms and ensure complementarity between the protein surfaces, as well as mediating polar interactions. Therefore, modeling of interface water is of importance in understanding the structural basis of biomolecular association. We present an algorithm, WATGEN, which predicts locations for water molecules at a protein-protein or protein-peptide interface, given the atomic coordinates of the protein and peptide. A key element of the WATGEN algorithm is the prediction of water sites that can form multiple hydrogen bonds that bridge the binding interface. Trial calculations were performed on water networks predicted by WATGEN at 126 protein-peptide interfaces (X-ray resolutions <or= 2.0 A), using different criteria for water placement. The energies of the predicted water networks were evaluated in AMBER8 and used in the choice of parameters for WATGEN. The 126 interfaces include 1264 experimentally determined bridging water sites, and the WATGEN algorithm predicts 72 and 88% of these sites within 1.5 and 2.0 A, respectively. The predicted number of water molecules at each interface was much higher than the number of water molecules identified experimentally. Therefore, random placement of the same number of water molecules as that predicted at each interface was performed as a control, and resulted in only 22 and 40% of water sites placed within 1.5 and 2.0 A of experimental sites, respectively. Based on these data, we conclude that WATGEN can accurately predict the location of water molecules at a protein-peptide interface, and this may be of value for understanding the energetics and specificity of biomolecular association.

Is There a Dynamic DNA-Protein Interface in the Transcription Factor IIIA-5 S rRNA Gene Complex?

Others have proposed that several amino acid side chains exhibit considerable conformational mobility at the DNA-protein interface in the transcription factor IIIA.5 S rRNA gene complex and that the rapid movements of these side chains permit them to make fluctuating contacts with adjacent bp in the DNA target site. This "dynamic interface" model makes biochemical predictions concerning the consequences of truncating specific amino acid side chains and the effects of these truncations on sequence selectivity in DNA binding. The model also makes predictions concerning the effects of DNA sequence context on the apparent energetic contributions to binding made by individual bp. We have tested these predictions, and our results are inconsistent with any significant energetic role being played by the contact of multiple bp by conformationally mobile amino acid side chains. They do, however, show that some individual amino acids affect the recognition of multiple bp through mechanisms other than direct interaction.

Role of the linker domain and the 203-214 N-terminal residues in the human topoisomerase I DNA complex dynamics

The influence of the N-terminal residues 203-214 and the linker domain on motions in the human topoisomerase I-DNA complex has been investigated by comparing the molecular dynamics simulations of the system with (topo70) or without (topo58/6.3) these regions. Topo58/6.3 is found to fluctuate more than topo70, indicating that the presence of the N-terminal residues and the linker domain dampen the core and C-terminal fluctuations. The simulations also show that residues 203-207 and the linker domain participate in a network of correlated movements with key regions of the enzyme, involved in the human topoisomerase I catalytic cycle, providing a structural-dynamical explanation for the better DNA relaxation activity of topo70 when compared to topo58/6.3. The data have been examined in relation to a wealth of biochemical, site-directed mutagenesis and crystallographic data on human topoisomerase I. The simulations finally show the occurrence of a network of direct and water mediated hydrogen bonds in the proximity of the active site, and the presence of a water molecule in the appropriate position to accept a proton from the catalytic Tyr-723 residue, suggesting that water molecules have an important role in the stabilization and function of this enzyme.

The role of water in protein-DNA recognition

Is it by design or by default that water molecules are observed at the interfaces of some protein-DNA complexes? Both experimental and theoretical studies on the thermodynamics of protein-DNA binding overwhelmingly support the extended hydrophobic view that water release from interfaces favors binding. Structural and energy analyses indicate that the waters that remain at the interfaces of protein-DNA complexes ensure liquid-state packing densities, screen the electrostatic repulsions between like charges (which seems to be by design), and in a few cases act as linkers between complementary charges on the biomolecules (which may well be by default). This review presents a survey of the current literature on water in protein-DNA complexes and a critique of various interpretations of the data in the context of the role of water in protein-DNA binding and principles of protein-DNA recognition in general.

From sequence to structure and back again: approaches for predicting protein-DNA binding

Gene regulation in higher organisms is achieved by a complex network of transcription factors (TFs). Modulating gene expression and exploring gene function are major aims in molecular biology. Furthermore, the identification of putative target genes for a certain TF serve as powerful tools for specific targeting of rational drugs.

Detecting the short and variable transcription factor binding sites (TFBSs) in genomic DNA is an intriguing challenge for computational and structural biologists. Fast and reliable computational methods for predicting TFBSs on a whole-genome scale offer several advantages compared to the current experimental methods that are rather laborious and slow. Two main approaches are being explored, advanced sequence-based algorithms and structure-based methods.

The aim of this review is to outline the computational and experimental methods currently being applied in the field of protein-DNA interactions. With a focus on the former, the current state of the art in modeling these interactions is discussed. Surveying sequence and structure-based methods for predicting TFBSs, we conclude that in order to achieve a sound and specific method applicable on genomic sequences it is desirable and important to bring these two approaches together.

Molecular dynamics simulations of metalloproteins

Molecular dynamics simulations are now commonly applied to metalloproteins, despite the challenges introduced by the presence of metal ions. Force field parameters are nowadays available also for these 'exotic' atoms and several biological systems have been successfully studied. Some of the most relevant results and methodological advancements are reviewed.

Molecular dynamics simulation of the RNA complex of a double-stranded RNA-binding domain reveals dynamic features of the intermolecular interface and its hydration

The interaction between double-stranded RNA (dsRNA) and the third double-stranded domain (dsRBD) from Drosophila Staufen protein represents a paradigm to understand how the dsRBD protein family, one of the most common RNA-binding protein units, binds dsRNA. The nuclear magnetic resonance (NMR) structure of this complex and the x-ray structure of another family member revealed the stereochemical basis for recognition, but also raised new questions. Although the crystallographic studies revealed a highly ordered interface containing numerous water-mediated contacts, NMR suggested extensive residual motion at the interface. To address how interfacial motion contributes to molecular recognition in the dsRBD-dsRNA system, we conducted a 2-ns molecular dynamics simulation of the complex derived from Staufen protein and of the separate protein and RNA components. The results support the observation that a high degree of conformational flexibility is retained upon complex formation and that this involves interfacial residues that are critical for dsRBD-dsRNA binding. The structural origin of this residual flexibility is revealed by the analysis of the trajectory of motion. Individual basic side chains switch continuously from one RNA polar group to another with a residence time seldom exceeding 100 ps, while retaining favorable interaction with RNA throughout much of the simulation. Short-lived water molecules mediate some of these interactions for a large fraction of the trajectory studied here. This result indicates that water molecules are not statically associated with the interface, but continuously exchange with the bulk solvent on a 1-10-ps time scale. This work provides new insight into dsRBD-dsRNA recognition and builds upon a growing body of evidence, suggesting that short-lived dynamic interactions play important roles in protein-nucleic acid interactions.

PvuII-endonuclease induces structural alterations at the scissile phosphate group of its cognate DNA

We investigated the PvuII endonuclease with its cognate DNA by means of molecular dynamics simulations. Comparing the complexed DNA with a reference simulation of free DNA, we saw structural changes at the scissile phosphodiester bond. At this GpC step, the enzyme induces the highest twist and axial rise, inclination is increased and the minor groove widened. The distance between the scissile phosphate group and the phosphate group of the following thymine base is shortened significantly, indicating a substrate-assisted catalysis. A feasible reason for this vicinity is the catalytically important amino acid residue lysine 70, which bridges the free oxygen atoms of the successive phosphate groups. Due to this geometry, a compact reaction pocket is formed where a water molecule can be held, thus bringing the reaction partners for hydrolysis into contact. The O1-P-O2 angle of the scissile nucleotide is decreased, probably due to a complexation of the negative oxygen atoms through protein and solvent contacts.

Simulations of nucleic acids and their complexes

Recent years have seen considerable progress in simulations of nucleic acids. Improvements in force fields, simulation techniques and protocols, and increasing computer power have all contributed to making nanosecond-scale simulations of both DNA and RNA commonplace. The results are already helping to explain how nucleic acids respond to their environment and to their base sequence and to reveal the factors underlying recognition processes by probing biologically important nucleic acid-protein interactions and medically important nucleic acid-drug complexation. This Account summarizes methodological progress and applications of molecular dynamics to nucleic acids over the past few years and tries to identify remaining challenges.

RNA solvation: a molecular dynamics simulation perspective

With the availability of accurate methods to treat the electrostatic long-range interactions, molecular dynamics simulations have resulted in refined dynamical models of the structure of the hydration shell around RNA motifs. The models reviewed here range from basic Watson-Crick to more specific noncanonical base pairs, from "simple" double helices to RNA molecules displaying more complex tertiary folds, and from DNA/RNA hybrid double helices to RNA hybrids formed with a chemically modified strand.

Dynamics of hydration in hen egg white lysozyme

We investigate the hydration dynamics of a small globular protein, hen egg-white lysozyme. Extensive simulations (two trajectories of 9 ns each) were carried out to identify the time-scales and mechanism of water attachment to this protein. The location of the surface and integral water molecules in lysozyme was also investigated. Three peculiar temporal scales of the hydration dynamics can be discerned: two among these, with sub-nanosecond mean residence time, tau(w), are characteristic of surface hydration water; the slower time-scale (tau(w) approximately 2/3 ns) is associated with buried water molecules in hydrophilic pores and in superficial clefts. The computed tau(w) values in the two independent runs fall in a similar range and are consistent with each other, thus adding extra weight to our result. The tau(w) of surface water obtained from the two independent trajectories is 20 and 24 ps. In both simulations only three water molecules are bound to lysozyme for the entire length of the trajectories, in agreement with nuclear magnetic relaxation dispersion estimates. Locations other than those identified in the protein crystal are found to be possible for these long-residing water molecules. The dynamics of the hydration water molecules observed in our simulations implies that each water molecule visits a multitude of residues during the lifetime of its bound with the protein. The number of residues seen by a single water molecule increases with the time-scale of its residence time and, on average, is equal to one only for the water molecules with shorter residence time. Thus, tau(w) values obtained from inelastic neutron scattering and based on jump-diffusion models are likely not to account for the contribution of water molecules with longer residence time.

Structure and hydration of the DNA-human topoisomerase I covalent complex

The structure and hydration of reconstituted human topoisomerase I comprising the core and the carboxyl-terminal domains in covalent complex with 22-basepair DNA duplex has been investigated by molecular dynamics simulation. The structure and the intermolecular interactions were found to be well maintained over the simulation. The complex displays a high degree of flexibility of the contact area, confirmed by the presence of numerous water-mediated protein-DNA hydrogen bonds comparable in quantity and distribution to the direct ones. The interaction between the enzyme and the solvent also provides the key for interpreting the experimental reduction of activity or affinity observed upon single residue mutation. Finally, four long lasting water molecules are observed in the proximity of the active site, one of which in the appropriate position to accept a proton from the active Tyr723.