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

Modeling the Homologous Recombination Process: Methods, Successes and Challenges

Homologous recombination (HR) is a fundamental process common to all species. HR aims to faithfully repair DNA double strand breaks. HR involves the formation of nucleoprotein filaments on DNA single strands (ssDNA) resected from the break. The nucleoprotein filaments search for homologous regions in the genome and promote strand exchange with the ssDNA homologous region in an unbroken copy of the genome. HR has been the object of intensive studies for decades. Because multi-scale dynamics is a fundamental aspect of this process, studying HR is highly challenging, both experimentally and using computational approaches. Nevertheless, knowledge has built up over the years and has recently progressed at an accelerated pace, borne by increasingly focused investigations using new techniques such as single molecule approaches. Linking this knowledge to the atomic structure of the nucleoprotein filament systems and the succession of unstable, transient intermediate steps that takes place during the HR process remains a challenge; modeling retains a very strong role in bridging the gap between structures that are stable enough to be observed and in exploring transition paths between these structures. However, working on ever-changing long filament systems submitted to kinetic processes is full of pitfalls. This review presents the modeling tools that are used in such studies, their possibilities and limitations, and reviews the advances in the knowledge of the HR process that have been obtained through modeling. Notably, we will emphasize how cooperative behavior in the HR nucleoprotein filament enables modeling to produce reliable information.

UNILIPID, a Methodology for Energetically Accurate Prediction of Protein Insertion into Implicit Membranes of Arbitrary Shape

The insertion of proteins into membranes is crucial for understanding their function in many biological processes. In this work, we present UNILIPID, a universal implicit lipid-protein description as a methodology for dealing with implicit membranes. UNILIPID is independent of the scale of representation and can be applied at the level of all atoms, coarse-grained particles down to the level of a single bead per amino acid. We provide example implementations for these scales and demonstrate the versatility of our approach by accurately reflecting the free energy of transfer for each amino acid. In addition to single membranes, we describe the analytical implementation of double membranes and show that UNILIPID is well suited for modeling at multiple scales. We generalize to membranes of arbitrary shape. With UNILIPID, we provide a methodological framework for a simple and general parameterization tuned to reproduce a selected reference hydrophobicity scale. The software we provide along with the methodological description is optimized for specific user features such as real-time response, live visual analysis, and virtual reality experiences.

Primordial bioenergy sources: The two facets of adenosine triphosphate

Life requires energy to exist, to reproduce and to survive. Two major hypotheses have been put forward concerning the source of this energy at the very early stages of life evolution: (i) abiotic organics either brought to Earth by comets and/or meteorites, or produced at its atmosphere, and (ii) mineral surface-dependent bioinorganic catalytic reactions. Considering the latter possibility, I propose that, besides being a precursor of nucleic acids, adenosine triphosphate (ATP), which probably was used very early to improve the fidelity of nucleic acid polymerization, played an essential role in the transition between mineral-bound protocells and their free counterparts. Indeed, phosphorylation by ATP renders carboxylate groups electrophilic enough to react with nucleophiles such as amines, an effect that, thanks to their Lewis acid character, also have dehydrated metal ions on mineral surfaces. Early ATP synthesis for metabolic processes most likely depended on substrate level phosphorylation. However, the exaptation of a hexameric helicase-like ATPase and a transmembrane H⁺ pump (which evolved to counteract the acidity caused by fermentation reactions within the protocell) generated a much more efficient membrane-bound ATP synthase that uses chemiosmosis to make ATP.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

A stretched conformation of DNA with a biological role?

We have discovered a well-defined extended conformation of double-stranded DNA, which we call Σ-DNA, using laser-tweezers force-spectroscopy experiments. At a transition force corresponding to free energy change ΔG = 1·57 ± 0·12 kcal (mol base pair)-1 60 or 122 base-pair long synthetic GC-rich sequences, when pulled by the 3'-3' strands, undergo a sharp transition to the 1·52 ± 0·04 times longer Σ-DNA. Intriguingly, the same degree of extension is also found in DNA complexes with recombinase proteins, such as bacterial RecA and eukaryotic Rad51. Despite vital importance to all biological organisms for survival, genome maintenance and evolution, the recombination reaction is not yet understood at atomic level. We here propose that the structural distortion represented by Σ-DNA, which is thus physically inherent to the nucleic acid, is related to how recombination proteins mediate recognition of sequence homology and execute strand exchange. Our hypothesis is that a homogeneously stretched DNA undergoes a 'disproportionation' into an inhomogeneous Σ-form consisting of triplets of locally B-like perpendicularly stacked bases. This structure may ensure improved fidelity of base-pair recognition and promote rejection in case of mismatch during homologous recombination reaction. Because a triplet is the length of a gene codon, we speculate that the structural physics of nucleic acids may have biased the evolution of recombinase proteins to exploit triplet base stacks and also the genetic code.

Direct Single-Molecule Observation of Mode and Geometry of RecA-Mediated Homology Search

Genomic integrity, when compromised by accrued DNA lesions, is maintained through efficient repair via homologous recombination. For this process the ubiquitous recombinase A (RecA), and its homologues such as the human Rad51, are of central importance, able to align and exchange homologous sequences within single-stranded and double-stranded DNA in order to swap out defective regions. Here, we directly observe the widely debated mechanism of RecA homology searching at a single-molecule level using high-speed atomic force microscopy (HS-AFM) in combination with tailored DNA origami frames to present the reaction targets in a way suitable for AFM-imaging. We show that RecA nucleoprotein filaments move along DNA substrates via short-distance facilitated diffusions, or slides, interspersed with longer-distance random moves, or hops. Importantly, from the specific interaction geometry, we find that the double-stranded substrate DNA resides in the secondary DNA binding-site within the RecA nucleoprotein filament helical groove during the homology search. This work demonstrates that tailored DNA origami, in conjunction with HS-AFM, can be employed to reveal directly conformational and geometrical information on dynamic protein-DNA interactions which was previously inaccessible at an individual single-molecule level.

Cooperative RecA clustering: the key to efficient homology searching

The mechanism by which pre-synaptic RecA nucleoprotein filaments efficiently locate sequence homology across genomic DNA remains unclear. Here, using atomic force microscopy, we directly investigate the intermediates of the RecA-mediated homologous recombination process and find it to be highly cooperative, involving multiple phases. Initially, the process is dominated by a rapid ‘association’ phase, where multiple filaments interact on the same dsDNA simultaneously. This cooperative nature is reconciled by the observation of localized dense clusters of pre-synaptic filaments interacting with the observed dsDNA molecules. This confinement of reactive species within the vicinity of the dsDNA, is likely to play an important role in ensuring that a high interaction rate between the nucleoprotein filaments and the dsDNA can be achieved. This is followed by a slower ‘resolution’ phase, where the synaptic joints either locate sequence homology and progress to a post-synaptic joint, or dissociate from the dsDNA. Surprisingly, the number of simultaneous synaptic joints decreases rapidly after saturation of the dsDNA population, suggesting a reduction in interaction activity of the RecA filaments. We find that the time-scale of this decay is in line with the time-scale of the dispersion of the RecA filament clusters, further emphasising the important role this cooperative phenomena may play in the RecA-facilitated homology search.

Mechanism of the formation of the RecA-ssDNA nucleoprotein filament structure: a coarse-grained approach

In prokaryotes, the RecA protein catalyzes the repair and strand exchange of double-stranded DNA. RecA binds to single-stranded DNA (ssDNA) and forms a presynaptic complex in which the protein polymerizes around the ssDNA to form a right-handed helical nucleoprotein filament structure. In the present work, the mechanism for the formation of the RecA-ssDNA filament structure is modeled using coarse-grained molecular dynamics simulations. Information from the X-ray structure was used to model the protein itself but not its interactions; the interactions between the protein and the ssDNA were modeled solely by electrostatic, aromatic, and repulsive energies. For the present study, the monomeric, dimeric, and trimeric units of RecA and 4, 8, and 11 NT-long ssDNA, respectively, were studied. Our results indicate that monomeric RecA is not sufficient for nucleoprotein filament formation; rather, dimeric RecA is the elementary binding unit, with higher multimeric units of RecA facilitating filament formation. Our results reveal that loop region flexibility at the primary binding site of RecA is essential for it to bind the incoming ssDNA, that the aromatic residues present in the loop region play an important role in ssDNA binding, and that ATP may play a role in guiding the ssDNA by changing the electrostatic potential of the RecA protein.

Mechanics and Single-Molecule Interrogation of DNA Recombination

The repair of DNA by homologous recombination is an essential, efficient, and high-fidelity process that mends DNA lesions formed during cellular metabolism; these lesions include double-stranded DNA breaks, daughter-strand gaps, and DNA cross-links. Genetic defects in the homologous recombination pathway undermine genomic integrity and cause the accumulation of gross chromosomal abnormalities-including rearrangements, deletions, and aneuploidy-that contribute to cancer formation. Recombination proceeds through the formation of joint DNA molecules-homologously paired but metastable DNA intermediates that are processed by several alternative subpathways-making recombination a versatile and robust mechanism to repair damaged chromosomes. Modern biophysical methods make it possible to visualize, probe, and manipulate the individual molecules participating in the intermediate steps of recombination, revealing new details about the mechanics of genetic recombination. We review and discuss the individual stages of homologous recombination, focusing on common pathways in bacteria, yeast, and humans, and place particular emphasis on the molecular mechanisms illuminated by single-molecule methods.

Integrating multi-scale data on homologous recombination into a new recognition mechanism based on simulations of the RecA-ssDNA/dsDNA structure

RecA protein is the prototypical recombinase. Members of the recombinase family can accurately repair double strand breaks in DNA. They also provide crucial links between pairs of sister chromatids in eukaryotic meiosis. A very broad outline of how these proteins align homologous sequences and promote DNA strand exchange has long been known, as are the crystal structures of the RecA-DNA pre- and postsynaptic complexes; however, little is known about the homology searching conformations and the details of how DNA in bacterial genomes is rapidly searched until homologous alignment is achieved. By integrating a physical model of recognition to new modeling work based on docking exploration and molecular dynamics simulation, we present a detailed structure/function model of homology recognition that reconciles extremely quick searching with the efficient and stringent formation of stable strand exchange products and which is consistent with a vast body of previously unexplained experimental results.

The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems

The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.

Innovative interactive flexible docking method for multi-scale reconstruction elucidates dystrophin molecular assembly

At present, our molecular knowledge of dystrophin, the protein encoded by the DMD gene and mutated in myopathy patients, remains limited. To get around the absence of its atomic structure, we have developed an innovative interactive docking method based on the BioSpring software in combination with Small-angle X-ray Scattering (SAXS) data. BioSpring allows interactive handling of biological macromolecules thanks to an augmented Elastic Network Model (aENM) that combines the spring network with non-bonded terms between atoms or pseudo-atoms. This approach can be used for building molecular assemblies even on a desktop or a laptop computer thanks to code optimizations including parallel computing and GPU programming. By combining atomistic and coarse-grained models, the approach significantly simplifies the set-up of multi-scale scenarios. BioSpring is remarkably efficient for the preparation of numeric simulations or for the design of biomolecular models integrating qualitative experimental data restraints. The combination of this program and SAXS allowed us to propose the first high-resolution models of the filamentous central domain of dystrophin, covering repeats 11 to 17. Low-resolution interactive docking experiments driven by a potential grid enabled us to propose how dystrophin may associate with F-actin and nNOS. This information provides an insight into medically relevant discoveries to come.

Model of a DNA-protein complex of the architectural monomeric protein MC1 from Euryarchaea

In Archaea the two major modes of DNA packaging are wrapping by histone proteins or bending by architectural non-histone proteins. To supplement our knowledge about the binding mode of the different DNA-bending proteins observed across the three domains of life, we present here the first model of a complex in which the monomeric Methanogen Chromosomal protein 1 (MC1) from Euryarchaea binds to the concave side of a strongly bent DNA. In laboratory growth conditions MC1 is the most abundant architectural protein present in Methanosarcina thermophila CHTI55. Like most proteins that strongly bend DNA, MC1 is known to bind in the minor groove. Interaction areas for MC1 and DNA were mapped by Nuclear Magnetic Resonance (NMR) data. The polarity of protein binding was determined using paramagnetic probes attached to the DNA. The first structural model of the DNA-MC1 complex we propose here was obtained by two complementary docking approaches and is in good agreement with the experimental data previously provided by electron microscopy and biochemistry. Residues essential to DNA-binding and -bending were highlighted and confirmed by site-directed mutagenesis. It was found that the Arg25 side-chain was essential to neutralize the negative charge of two phosphates that come very close in response to a dramatic curvature of the DNA.

Coarse-grain modelling of protein-protein interactions

Here, we review recent advances towards the modelling of protein-protein interactions (PPI) at the coarse-grained (CG) level, a technique that is now widely used to understand protein affinity, aggregation and self-assembly behaviour. PPI models of soluble proteins and membrane proteins are separately described, but we note the parallel development that is present in both research fields with three important themes: firstly, combining CG modelling with knowledge-based approaches to predict and refine protein-protein complexes; secondly, using physics-based CG models for de novo prediction of protein-protein complexes; and thirdly modelling of large scale protein aggregates.

Simplified biased random walk model for RecA-protein-mediated homology recognition offers rapid and accurate self-assembly of long linear arrays of binding sites

Inspired by RecA-protein-based homology recognition, we consider the pairing of two long linear arrays of binding sites. We propose a fully reversible, physically realizable biased random walk model for rapid and accurate self-assembly due to the spontaneous pairing of matching binding sites, where the statistics of the searched sample are included. In the model, there are two bound conformations, and the free energy for each conformation is a weakly nonlinear function of the number of contiguous matched bound sites.

Haptic-driven applications to molecular modeling: state-of-the-art and perspectives

Drug design is a creative process that combines different scientific expertise. With the development of increasingly powerful computers, disciplines such as molecular modeling and, in particular, drug design, are becoming an important component of drug discovery. However, modern software often limits the user interaction with the computer calculation, reducing the potential for researchers to use their knowledge in the design process. For this reason, interactive methodologies have been investigated in recent years. In particular, haptic-driven simulators offer the possibility for users to drive and control the modeling simulations, efficiently combining human knowledge and computational power. In this article, we will discuss the state-of-the-art and future perspectives of such methodologies.

Modelling Nucleic Acid Structure and Flexibility: From Atomic to Mesoscopic Scale

This chapter surveys some of the recent developments in coarse-grained modelling of nucleic acids. We first discuss models based on pseudoatoms, effective spherical particles representing groups of atoms. A major part of the chapter is devoted to models in which bases or base pairs are represented as independent, interacting rigid bodies. Two popular definitions of internal coordinates, as implemented in the programs 3DNA and Curves+, are outlined from a common perspective. Recently developed rigid base and basepair models with nonlocal quadratic interactions are presented. A statistical mechanical description of the models on their full phase space yields exact relations between model parameters and expected values of some state functions. We estimated shape and stiffness parameters for nonlocal rigid base and basepair models of a DNA oligomer containing A-tract. The parameterization is based on atomic-resolution molecular dynamics simulation data. We found that the rigid base model is consistent with a local interaction pattern, while interactions in the rigid basepair model are visibly non-local, in agreement with earlier findings. Differences in shape and stiffness parameters obtained using Curves+ and 3DNA coordinates are found to be small for structures within the B-DNA family. Anharmonic effects, coarser models, and other approaches to describe nucleic acid structure and flexibility are discussed.

Ca2+ improves organization of single-stranded DNA bases in human Rad51 filament, explaining stimulatory effect on gene recombination

Human RAD51 protein (HsRad51) catalyses the DNA strand exchange reaction for homologous recombination. To clarify the molecular mechanism of the reaction in vitro being more effective in the presence of Ca²⁺ than of Mg²⁺, we have investigated the effect of these ions on the structure of HsRad51 filament complexes with single- and double-stranded DNA, the reaction intermediates. Flow linear dichroism spectroscopy shows that the two ionic conditions induce significantly different structures in the HsRad51/single-stranded DNA complex, while the HsRad51/double-stranded DNA complex does not demonstrate this ionic dependence. In the HsRad51/single-stranded DNA filament, the primary intermediate of the strand exchange reaction, ATP/Ca²⁺ induces an ordered conformation of DNA, with preferentially perpendicular orientation of nucleobases relative to the filament axis, while the presence of ATP/Mg²⁺, ADP/Mg²⁺ or ADP/Ca²⁺ does not. A high strand exchange activity is observed for the filament formed with ATP/Ca²⁺, whereas the other filaments exhibit lower activity. Molecular modelling suggests that the structural variation is caused by the divalent cation interfering with the L2 loop close to the DNA-binding site. It is proposed that the larger Ca²⁺ stabilizes the loop conformation and thereby the protein–DNA interaction. A tight binding of DNA, with bases perpendicularly oriented, could facilitate strand exchange.

ATTRACT and PTools: open source programs for protein-protein docking

The prediction of the structure of protein-protein complexes based on structures or structural models of isolated partners is of increasing importance for structural biology and bioinformatics. The ATTRACT program can be used to perform systematic docking searches based on docking energy minimization. It is part of the object-oriented PTools library written in Python and C++. The library contains various routines to manipulate protein structures, to prepare and perform docking searches as well as analyzing docking results. It also intended to facilitate further methodological developments in the area of macromolecular docking that can be easily integrated. Here, we describe the application of PTools to perform systematic docking searches and to analyze the results. In addition, the possibility to perform multi-component docking will also be presented.

Accounting for large amplitude protein deformation during in silico macromolecular docking

Rapid progress of theoretical methods and computer calculation resources has turned in silico methods into a conceivable tool to predict the 3D structure of macromolecular assemblages, starting from the structure of their separate elements. Still, some classes of complexes represent a real challenge for macromolecular docking methods. In these complexes, protein parts like loops or domains undergo large amplitude deformations upon association, thus remodeling the surface accessible to the partner protein or DNA. We discuss the problems linked with managing such rearrangements in docking methods and we review strategies that are presently being explored, as well as their limitations and success.