Multiple protein-DNA interfaces unravelled by evolutionary information, physico-chemical and geometrical properties.

Interactions between proteins and nucleic acids are at the heart of many essential biological processes. Despite increasing structural information about how these interactions may take place, our understanding of the usage made of protein surfaces by nucleic acids is still very limited. This is in part due to the inherent complexity associated to protein surface deformability and evolution. In this work, we present a method that contributes to decipher such complexity by predicting protein-DNA interfaces and characterizing their properties. It relies on three biologically and physically meaningful descriptors, namely evolutionary conservation, physico-chemical properties and surface geometry. We carefully assessed its performance on several hundreds of protein structures and compared it to several machine-learning state-of-the-art methods. Our approach achieves a higher sensitivity compared to the other methods, with a similar precision. Importantly, we show that it is able to unravel 'hidden' binding sites by applying it to unbound protein structures and to proteins binding to DNA via multiple sites and in different conformations. It is also applicable to the detection of RNA-binding sites, without significant loss of performance. This confirms that DNA and RNA-binding sites share similar properties. Our method is implemented as a fully automated tool, [Formula: see text], freely accessible at: http://www.lcqb.upmc.fr/JET2DNA. We also provide a new dataset of 187 protein-DNA complex structures, along with a subset of 82 associated unbound structures. The set represents the largest body of high-resolution crystallographic structures of protein-DNA complexes, use biological protein assemblies as DNA-binding units, and covers all major types of protein-DNA interactions. It is available at: http://www.lcqb.upmc.fr/PDNAbenchmarks.

The static and dynamic structural heterogeneities of B-DNA: extending Calladine-Dickerson rules.

We present a multi-laboratory effort to describe the structural and dynamical properties of duplex B-DNA under physiological conditions. By processing a large amount of atomistic molecular dynamics simulations, we determine the sequence-dependent structural properties of DNA as expressed in the equilibrium distribution of its stochastic dynamics. Our analysis includes a study of first and second moments of the equilibrium distribution, which can be accurately captured by a harmonic model, but with nonlocal sequence-dependence. We characterize the sequence-dependent choreography of backbone and base movements modulating the non-Gaussian or anharmonic effects manifested in the higher moments of the dynamics of the duplex when sampling the equilibrium distribution. Contrary to prior assumptions, such anharmonic deformations are not rare in DNA and can play a significant role in determining DNA conformation within complexes. Polymorphisms in helical geometries are particularly prevalent for certain tetranucleotide sequence contexts and are always coupled to a complex network of coordinated changes in the backbone. The analysis of our simulations, which contain instances of all tetranucleotide sequences, allow us to extend Calladine-Dickerson rules used for decades to interpret the average geometry of DNA, leading to a set of rules with quantitative predictive power that encompass nonlocal sequence-dependence and anharmonic fluctuations.

Internal Coordinate Normal Mode Analysis: A Strategy To Predict Protein Conformational Transitions.

We analyze the capacity of normal mode analysis in internal coordinates' space to generate large-amplitude structural deformations that can describe the conformational changes occurring upon the binding of proteins to other species. We also analyze how many modes need to be studied to capture a given transition and whether a combination of two modes is better than using a single mode. The technique is tested on known unbound-to-bound structural transitions for a set of single- or multidomain proteins. The results suggest that this approach is a promising way to generate structures for protein docking or for more refined molecular dynamics simulations.

A molecular dynamics study of adenylyl cyclase: The impact of ATP and G-protein binding.

Adenylyl cyclases (ACs) catalyze the biosynthesis of cyclic adenosine monophosphate (cAMP) from adenosine triphosphate (ATP) and play an important role in many signal transduction pathways. The enzymatic activity of ACs is carefully controlled by a variety of molecules, including G-protein subunits that can both stimulate and inhibit cAMP production. Using homology models developed from existing structural data, we have carried out all-atom, microsecond-scale molecular dynamics simulations on the AC5 isoform of adenylyl cyclase and on its complexes with ATP and with the stimulatory G-protein subunit Gsα. The results show that both ATP and Gsα binding have significant effects on the structure and flexibility of adenylyl cyclase. New data on ATP bound to AC5 in the absence of Gsα notably help to explain how Gsα binding enhances enzyme activity and could aid product release. Simulations also suggest a possible coupling between ATP binding and interactions with the inhibitory G-protein subunit Gαi.

Sequence-dependent response of DNA to torsional stress: a potential biological regulation mechanism.

Torsional restraints on DNA change in time and space during the life of the cell and are an integral part of processes such as gene expression, DNA repair and packaging. The mechanical behavior of DNA under torsional stress has been studied on a mesoscopic scale, but little is known concerning its response at the level of individual base pairs and the effects of base pair composition. To answer this question, we have developed a geometrical restraint that can accurately control the total twist of a DNA segment during all-atom molecular dynamics simulations. By applying this restraint to four different DNA oligomers, we are able to show that DNA responds to both under- and overtwisting in a very heterogeneous manner. Certain base pair steps, in specific sequence environments, are able to absorb most of the torsional stress, leaving other steps close to their relaxed conformation. This heterogeneity also affects the local torsional modulus of DNA. These findings suggest that modifying torsional stress on DNA could act as a modulator for protein binding via the heterogeneous changes in local DNA structure.

Decomposing protein-DNA binding and recognition using simplified protein models.

We analyze the role of different physicochemical factors in protein/DNA binding and recognition by comparing the results from all-atom molecular dynamics simulations with simulations using simplified protein models. These models enable us to separate the role of specific amino acid side chains, formal amino acid charges and hydrogen bonding from the effects of the low-dielectric volume occupied by the protein. Comparisons are made on the basis of the conformation of DNA after protein binding, the ionic distribution around the complex and the sequence specificity. The results for four transcription factors, binding in either the minor or major grooves of DNA, show that the protein volume and formal charges, with one exception, play a predominant role in binding. Adding hydrogen bonding and a very small number of key amino acid side chains at the all-atom level yields results in DNA conformations and sequence recognition close to those seen in the reference all-atom simulations.

Structure and Dynamics of a 197 bp Nucleosome in Complex with Linker Histone H1.

Linker histones associate with nucleosomes to promote the formation of higher-order chromatin structure, but the underlying molecular details are unclear. We investigated the structure of a 197 bp nucleosome bearing symmetric 25 bp linker DNA arms in complex with vertebrate linker histone H1. We determined electron cryo-microscopy (cryo-EM) and crystal structures of unbound and H1-bound nucleosomes and validated these structures by site-directed protein cross-linking and hydroxyl radical footprinting experiments. Histone H1 shifts the conformational landscape of the nucleosome by drawing the two linkers together and reducing their flexibility. The H1 C-terminal domain (CTD) localizes primarily to a single linker, while the H1 globular domain contacts the nucleosome dyad and both linkers, associating more closely with the CTD-distal linker. These findings reveal that H1 imparts a strong degree of asymmetry to the nucleosome, which is likely to influence the assembly and architecture of higher-order structures.

Analyzing DNA curvature and its impact on the ionic environment: application to molecular dynamics simulations of minicircles.

We propose a method for analyzing the magnitude and direction of curvature within nucleic acids, based on the curvilinear helical axis calculated by Curves+. The method is applied to analyzing curvature within minicircles constructed with varying degrees of over- or under-twisting. Using the molecular dynamics trajectories of three different minicircles, we are able to quantify how curvature varies locally both in space and in time. We also analyze how curvature influences the local environment of the minicircles, notably via increased heterogeneity in the ionic distributions surrounding the double helix. The approach we propose has been integrated into Curves+ and the utilities Canal (time trajectory analysis) and Canion (environmental analysis) and can be used to study a wide variety of static and dynamic structural data on nucleic acids.

The Role of Unconventional Hydrogen Bonds in Determining BII Propensities in B-DNA.

An accurate understanding of DNA backbone transitions is likely to be the key for elucidating the puzzle of the intricate sequence-dependent mechanical properties that govern most of the biologically relevant functions of the double helix. One factor believed to be important in indirect recognition within protein-DNA complexes is the combined effect of two DNA backbone torsions (ε and ζ) which give rise to the well-known BI/BII conformational equilibrium. In this work we explain the sequence-dependent BII propensity observed in RpY steps (R = purine; Y = pyrimidine) at the tetranucleotide level with the help of a previously undetected C-H···O contact between atoms belonging to adjacent bases. Our results are supported by extensive multimicrosecond molecular dynamics simulations from the Ascona B-DNA Consortium, high-level quantum mechanical calculations, and data mining of the experimental structures deposited in the Protein Data Bank.

Extra views on structure and dynamics of DNA loops on nucleosomes studied with molecular simulations.

It has been shown experimentally that the action of the RSC chromatin remodeler leads to the formation of an irregular, partially remodeled nucleosome, termed a remosome. The remosome contains an extra 30-40 base pairs of DNA compared to a canonical nucleosome. Large-scale molecular simulations have provided information on the probable structure of remosomes and have explained why they remain stable in the absence of RSC. Here we explain how these simulations were carried out and what the resulting remosome models imply in terms of the mechanism of action of RSC. We notably show that local kinks within DNA are key in explaining how extra DNA can be in added to nucleosomes without unduly disturbing DNA-histone binding.

Protein-DNA interfaces: a molecular dynamics analysis of time-dependent recognition processes for three transcription factors.

We have studied the dynamics of three transcription factor-DNA complexes using all-atom, microsecond-scale MD simulations. In each case, the salt bridges and hydrogen bond interactions formed at the protein-DNA interface are found to be dynamic, with lifetimes typically in the range of tens to hundreds of picoseconds, although some interactions, notably those involving specific binding to DNA bases, can be a hundred times longer lived. Depending on the complex studied, this dynamics may or may not lead to the existence of distinct conformational substates. Using a sequence threading technique, it has been possible to determine whether DNA sequence recognition is sensitive or not to such conformational changes, and, in one case, to show that recognition appears to be locally dependent on protein-mediated cation distributions.

DNA minicircles clarify the specific role of DNA structure on retroviral integration.

Chromatin regulates the selectivity of retroviral integration into the genome of infected cells. At the nucleosome level, both histones and DNA structure are involved in this regulation. We propose a strategy that allows to specifically study a single factor: the DNA distortion induced by the nucleosome. This strategy relies on mimicking this distortion using DNA minicircles (MCs) having a fixed rotational orientation of DNA curvature, coupled with atomic-resolution modeling. Contrasting MCs with linear DNA fragments having identical sequences enabled us to analyze the impact of DNA distortion on the efficiency and selectivity of integration. We observed a global enhancement of HIV-1 integration in MCs and an enrichment of integration sites in the outward-facing DNA major grooves. Both of these changes are favored by LEDGF/p75, revealing a new, histone-independent role of this integration cofactor. PFV integration is also enhanced in MCs, but is not associated with a periodic redistribution of integration sites, thus highlighting its distinct catalytic properties. MCs help to separate the roles of target DNA structure, histone modifications and integrase (IN) cofactors during retroviral integration and to reveal IN-specific regulation mechanisms.

Structure and dynamics of DNA loops on nucleosomes studied with atomistic, microsecond-scale molecular dynamics.

DNA loop formation on nucleosomes is strongly implicated in chromatin remodeling and occurs spontaneously in nucleosomes subjected to superhelical stress. The nature of such loops depends crucially on the balance between DNA deformation and DNA interaction with the nucleosome core. Currently, no high-resolution structural data on these loops exist. Although uniform rod models have been used to study loop size and shape, these models make assumptions concerning DNA mechanics and DNA-core binding. We present here atomic-scale molecular dynamics simulations for two different loop sizes. The results point to the key role of localized DNA kinking within the loops. Kinks enable the relaxation of DNA bending strain to be coupled with improved DNA-core interactions. Kinks lead to small, irregularly shaped loops that are asymmetrically positioned with respect to the nucleosome core. We also find that loop position can influence the dynamics of the DNA segments at the extremities of the nucleosome.

Long-timescale dynamics of the Drew-Dickerson dodecamer.

We present a systematic study of the long-timescale dynamics of the Drew-Dickerson dodecamer (DDD: d(CGCGAATTGCGC)2) a prototypical B-DNA duplex. Using our newly parameterized PARMBSC1 force field, we describe the conformational landscape of DDD in a variety of ionic environments from minimal salt to 2 M Na(+)Cl(-) or K(+)Cl(-) The sensitivity of the simulations to the use of different solvent and ion models is analyzed in detail using multi-microsecond simulations. Finally, an extended (10 µs) simulation is used to characterize slow and infrequent conformational changes in DDD, leading to the identification of previously uncharacterized conformational states of this duplex which can explain biologically relevant conformational transitions. With a total of more than 43 µs of unrestrained molecular dynamics simulation, this study is the most extensive investigation of the dynamics of the most prototypical DNA duplex.

Dynamics and recognition within a protein-DNA complex: a molecular dynamics study of the SKN-1/DNA interaction.

Molecular dynamics simulations of the Caenorhabditis elegans transcription factor SKN-1 bound to its cognate DNA site show that the protein-DNA interface undergoes significant dynamics on the microsecond timescale. A detailed analysis of the simulation shows that movements of two key arginine side chains between the major groove and the backbone of DNA generate distinct conformational sub-states that each recognize only part of the consensus binding sequence of SKN-1, while the experimentally observed binding specificity results from a time-averaged view of the dynamic recognition occurring within this complex.

Internal Normal Mode Analysis (iNMA) Applied to Protein Conformational Flexibility.

We analyze the capacity of normal modes to predict observed protein conformational changes, and, notably, those induced by the formation of protein-protein complexes. We show that normal modes calculated in internal coordinate space (ICS) provide better predictions. For a large test set, using the ICS approach describes the conformational changes more completely, and with fewer low-frequency modes than the equivalent Cartesian coordinate modes, despite the fact that the internal coordinate calculations were restricted to torsional angles. This can be attributed to the fact that the use of ICS extends the range over which movements along the corresponding eigenvectors remain close to the true conformational energy hypersurface. We also show that the PaLaCe coarse-grain protein model performs better than a simple elastic network model. We apply ICS normal-mode analysis to protein complexes and, by extending the approach of Sunada and Go̅, [Sunada, S.; Go̅, N. J. Comput. Chem. 1995, 16, 328-336], we show that we can couple an accurate view of the Cartesian coordinate movements induced by ICS modes with the detection of the key residues responsible for the movements.

Analyzing ion distributions around DNA: sequence-dependence of potassium ion distributions from microsecond molecular dynamics.

Microsecond molecular dynamics simulations of B-DNA oligomers carried out in an aqueous environment with a physiological salt concentration enable us to perform a detailed analysis of how potassium ions interact with the double helix. The oligomers studied contain all 136 distinct tetranucleotides and we are thus able to make a comprehensive analysis of base sequence effects. Using a recently developed curvilinear helicoidal coordinate method we are able to analyze the details of ion populations and densities within the major and minor grooves and in the space surrounding DNA. The results show higher ion populations than have typically been observed in earlier studies and sequence effects that go beyond the nature of individual base pairs or base pair steps. We also show that, in some special cases, ion distributions converge very slowly and, on a microsecond timescale, do not reflect the symmetry of the corresponding base sequence.

µABC: a systematic microsecond molecular dynamics study of tetranucleotide sequence effects in B-DNA.

We present the results of microsecond molecular dynamics simulations carried out by the ABC group of laboratories on a set of B-DNA oligomers containing the 136 distinct tetranucleotide base sequences. We demonstrate that the resulting trajectories have extensively sampled the conformational space accessible to B-DNA at room temperature. We confirm that base sequence effects depend strongly not only on the specific base pair step, but also on the specific base pairs that flank each step. Beyond sequence effects on average helical parameters and conformational fluctuations, we also identify tetranucleotide sequences that oscillate between several distinct conformational substates. By analyzing the conformation of the phosphodiester backbones, it is possible to understand for which sequences these substates will arise, and what impact they will have on specific helical parameters.

Unraveling the sequence-dependent polymorphic behavior of d(CpG) steps in B-DNA.

We have made a detailed study of one of the most surprising sources of polymorphism in B-DNA: the high twist/low twist (HT/LT) conformational change in the d(CpG) base pair step. Using extensive computations, complemented with database analysis, we were able to characterize the twist polymorphism in the d(CpG) step in all the possible tetranucleotide environment. We found that twist polymorphism is coupled with BI/BII transitions, and, quite surprisingly, with slide polymorphism in the neighboring step. Unexpectedly, the penetration of cations into the minor groove of the d(CpG) step seems to be the key element in promoting twist transitions. The tetranucleotide environment also plays an important role in the sequence-dependent d(CpG) polymorphism. In this connection, we have detected a previously unexplored intramolecular C-H···O hydrogen bond interaction that stabilizes the low twist state when 3'-purines flank the d(CpG) step. This work explains a coupled mechanism involving several apparently uncorrelated conformational transitions that has only been partially inferred by earlier experimental or theoretical studies. Our results provide a complete description of twist polymorphism in d(CpG) steps and a detailed picture of the molecular choreography associated with this conformational change.