Chantal Prévost

Accounting for loop flexibility during protein-protein docking.

Although reliable docking can now be achieved for systems that do not undergo important induced conformational change upon association, the presence of flexible surface loops, which must adapt to the steric and electrostatic properties of a partner, generally presents a major obstacle. We report here the first docking method that allows large loop movements during a systematic exploration of the possible arrangements of the two partners in terms of position and rotation. Our strategy consists in taking into account an ensemble of possible loop conformations by a multi‐copy representation within a reduced protein model. The docking process starts from regularly distributed positions and orientations of the ligand around the whole receptor. Each starting configuration is submitted to energy minimization during which the best‐fitting loop conformation is selected based on the mean‐field theory. Trials were carried out on proteins with significant differences in the main‐chain conformation of the binding loop between isolated form and complexed form, which were docked to their partner considered in their bound form. The method is able to predict complexes very close to the crystal complex both in terms of relative position of the two partners and of the geometry of the flexible loop. We also show that introducing loop flexibility on the isolated protein form during systematic docking largely improves the predictions of relative position of the partners in comparison with rigid‐body docking. Proteins 2006.

Docking macromolecules with flexible segments

We address a major obstacle to macromolecular docking algorithms by presenting a new method that takes into account the induced conformational adjustment of flexible loops situated at a protein/macromolecule interface. The method, MC2, is based on a multiple copy representation of the loops, coupled with a Monte Carlo conformational search of the relative position of the macromolecules and their side chain conformations. The selection of optimal loop conformations takes place during Monte Carlo cycling by the iterative adjustment of the weight of each copy. We describe here the parameterization of the method and trials on a protein‐DNA complex of known 3‐D structure, involving the Drosophila prd paired domain protein and its target oligonucleotide Wenqing, X. et al., Cell 1995, 80, 639 . We demonstrate that our algorithm can correctly configure and position this protein, despite its relatively complex interactions with both grooves of DNA.

Geometry of the DNA strands within the RecA nucleofilament: role in homologous recombination

Homologous recombination consists of exchanging DNA strands of identical or almost identical sequence. This process is important for both DNA repair and DNA segregation. In prokaryotes, it involves the formation of long helical filaments of the RecA protein on DNA. These filaments incorporate double-stranded DNA from the cells genetic material, recognize sequence homology and promote strand exchange between the two DNA segments. DNA processing by these nucleofilaments is characterized by large amplitude deformations of the double helix, which is stretched by 50% and unwound by 40% with respect to B-DNA. In this article, information concerning the structure and interactions of the RecA, DNA and ATP molecules involved in DNA strand exchange is gathered and analyzed to present a view of their possible arrangement within the filament, their behavior during strand exchange and during ATP hydrolysis, the mechanism of RecA-promoted DNA deformation and the role of DNA deformation in the process of homologous recombination. In particular, the unusual characteristics of DNA within the RecA filament are compared to the DNA deformations locally induced by architectural proteins which bind in the DNA minor groove. The possible role and location of two flexible loops of RecA are discussed.

Deforming DNA: from physics to biology.

The DNA double helix has become a modern icon which symbolizes our understanding of the molecular basis of life. It is less widely recognized that the double helix proposed by Watson and Crick more than half a century ago is a remarkably adaptable molecule that can undergo major conformational rearrangements without being irreversibly damaged. Indeed, DNA deformation is an intrinsic feature of many of the biological processes in which it is involved. Over the last two decades, single-molecule experiments coupled with molecular modeling have transformed our understanding of DNA flexibility, while the accumulation of high-resolution structures of DNA-protein complexes have demonstrated how organisms can exploit this property as a useful feature for preserving, reading, replicating, and packaging the genetic message. In this Minireview we summarize the information now available on the extreme--and the less extreme--deformations of the double helix.

PTools: an opensource molecular docking library.

BackgroundMacromolecular docking is a challenging field of bioinformatics. Developing new algorithms is a slow process generally involving routine tasks that should be found in a robust library and not programmed from scratch for every new software application.ResultsWe present an object-oriented Python/C++ library to help the development of new docking methods. This library contains low-level routines like PDB-format manipulation functions as well as high-level tools for docking and analyzing results. We also illustrate the ease of use of this library with the detailed implementation of a 3-body docking procedure.ConclusionThe PTools library can handle molecules at coarse-grained or atomic resolution and allows users to rapidly develop new software. The library is already in use for protein-protein and protein-DNA docking with the ATTRACT program and for simulation analysis. This library is freely available under the GNU GPL license, together with detailed documentation.

Insights on protein-DNA recognition by coarse grain modelling

Coarse grain modelling of macromolecules is a new approach, potentially well adapted to answer numerous issues, ranging from physics to biology. We propose here an original DNA coarse grain model specifically dedicated to protein‐DNA docking, a crucial, but still largely unresolved, question in molecular biology. Using a representative set of protein‐DNA complexes, we first show that our model is able to predict the interaction surface between the macromolecular partners taken in their bound form. In a second part, the impact of the DNA sequence and electrostatics, together with the DNA and protein conformations on docking is investigated. Our results strongly suggest that the overall DNA structure mainly contributes in discriminating the interaction site on cognate proteins. Direct electrostatic interactions between phosphate groups and amino acid side chains strengthen the binding. Overall, this work demonstrates that coarse grain modeling can reveal itself a precious auxiliary for a general and complete description and understanding of protein‐DNA association mechanisms.

A Molecular Model for RecA-Promoted Strand Exchange via Parallel Triple-Stranded Helices

A number of studies have concluded that strand exchange between a RecA-complexed DNA single strand and a homologous DNA duplex occurs via a single-strand invasion of the minor groove of the duplex. Using molecular modeling, we have previously demonstrated the possibility of forming a parallel triple helix in which the single strand interacts with the intact duplex in the minor groove, via novel base interactions (Bertucat et al., J. Biomol. Struct. Dynam. 16:535-546). This triplex is stabilized by the stretching and unwinding imposed by RecA. In the present study, we show that the bases within this triplex are appropriately placed to undergo strand exchange. Strand exchange is found to be exothermic and to result in a triple helix in which the new single strand occupies the major groove. This structure, which can be equated to so-called R-form DNA, can be further stabilized by compression and rewinding. We are consequently able to propose a detailed, atomic-scale model of RecA-promoted strand exchange. This model, which is supported by a variety of experimental data, suggests that the role of RecA is principally to prepare the single strand for its future interactions, to guide a minor groove attack on duplex DNA, and to stabilize the resulting, stretched triplex, which intrinsically favors strand exchange. We also discuss how this mechanism can incorporate homologous recognition.

Modeling the early stage of DNA sequence recognition within RecA nucleoprotein filaments

Homologous recombination is a fundamental process enabling the repair of double-strand breaks with a high degree of fidelity. In prokaryotes, it is carried out by RecA nucleofilaments formed on single-stranded DNA (ssDNA). These filaments incorporate genomic sequences that are homologous to the ssDNA and exchange the homologous strands. Due to the highly dynamic character of this process and its rapid propagation along the filament, the sequence recognition and strand exchange mechanism remains unknown at the structural level. The recently published structure of the RecA/DNA filament active for recombination (Chen et al., Mechanism of homologous recombination from the RecA-ssDNA/dsDNA structure, Nature 2008, 453, 489) provides a starting point for new exploration of the system. Here, we investigate the possible geometries of association of the early encounter complex between RecA/ssDNA filament and double-stranded DNA (dsDNA). Due to the huge size of the system and its dense packing, we use a reduced representation for protein and DNA together with state-of-the-art molecular modeling methods, including systematic docking and virtual reality simulations. The results indicate that it is possible for the double-stranded DNA to access the RecA-bound ssDNA while initially retaining its Watson–Crick pairing. They emphasize the importance of RecA L2 loop mobility for both recognition and strand exchange.

Structure/function relationships in RecA protein-mediated homology recognition and strand exchange.

Abstract RecA family proteins include RecA, Rad51, and Dmc1. These recombinases are responsible for homology search and strand exchange. Homology search and strand exchange occur during double-strand break repair and in eukaryotes during meiotic recombination. In bacteria, homology search begins when RecA binds an initiating single-stranded DNA (ssDNA) in the primary DNA-binding site to form the presynaptic filament. The filament is a right-handed helix, where the initiating strand is bound deep within the filament. Once the presynaptic filament is formed, it interrogates nearby double-stranded DNA (dsDNA) to find a homologous sequence; therefore, we provide a detailed discussion of structural features of the presynaptic filament that play important functional roles. The discussion includes many diagrams showing multiple filament turns. These diagrams illustrate interactions that are not evident in single turn structures. The first dsDNA interactions with the presynaptic filament are insensitive to mismatches. The mismatch insensitive interactions lead to dsDNA deformation that triggers a homology testing process governed by kinetics. The first homology test involves ∼8 bases. Almost all interactions are rejected by this initial rapid test, leading to a new cycle of homology testing. Interactions that pass the initial rapid test proceed to a slower testing stage. That slower stage induces nonhomologous dsDNA to reverse strand exchange and begin a new cycle of homology testing. In contrast, homologous dsDNA continues to extend the heteroduplex strand-exchange product until ATP hydrolysis makes strand exchange irreversible.

A model for parallel triple helix formation by RecA: single-single association with a homologous duplex via the minor groove.

The nucleoproteic filaments of RecA polymerized on single stranded DNA are able to integrate double stranded DNA in a coaxial arrangement (with DNA stretched by a factor 1.5), to recognize homologous sequences in the duplex and to perform strand exchange between the single stranded and double stranded molecules. While experimental results favor the hypothesis of an invasion of the minor groove of the duplex by the single strand, parallel minor groove triple helices have never been isolated or even modeled, the minor groove offering little space for a third strand to interact. Based on an internal coordinate modeling study, we show here that such a structure is perfectly conceivable when the two interacting oligomers are stretched by a factor 1.5, in order to open the minor groove of the duplex. The model helix presents characteristics that coincide with known experimental data on unwinding, base pair inclination and inter-proton distances. Moreover, we show that extension and unwinding stabilize the triple helix. New patterns of triplet interaction via the minor groove are presented.