It Takes Tau to Tango: Investigating the Fuzzy Interaction between the R2-Repeat Domain and Tubulin C-Terminal Tails.

The microtubule-associated protein (MAP) tau plays a key role in the regulation of microtubule assembly and spatial organization. Tau hyperphosphorylation affects its binding on the tubulin surface and has been shown to be involved in several pathologies such as Alzheimer's disease. As the tau binding site on the microtubule lays close to the disordered and highly flexible tubulin C-terminal tails (CTTs), these are likely to impact the tau-tubulin interaction. Since the disordered tubulin CTTs are missing from the available experimental structures, we used homology modeling to build two complete models of tubulin heterotrimers with different isotypes for the ß-tubulin subunit (ßI/αI/ßI and ßIII/αI/ßIII). We then performed long timescale classical Molecular Dynamics simulations for the tau-R2/tubulin assembly (in systems with and without CTTs) and analyzed the resulting trajectories to obtain a detailed view of the protein interface in the complex and the impact of the CTTs on the stability of this assembly. Additional analyses of the CTT mobility in the presence, or in the absence, of tau also highlight how tau might modulate the CTT activity as hooks that are involved in the recruitment of several MAPs. In particular, we observe a wrapping phenomenon, where the ß-tubulin CTTs form a loop over tau-R2, thus stabilizing its interaction with the tubulin surface and simultaneously reducing the CTT availability for interactions with other MAPs.

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.

Modeling the Dynamics of Protein-Protein Interfaces, How and Why?

Protein-protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein-protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein-protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein-protein interface.

Moving pictures: Reassessing docking experiments with a dynamic view of protein interfaces.

The modeling of protein assemblies at the atomic level remains a central issue in structural biology, as protein interactions play a key role in numerous cellular processes. This problem is traditionally addressed using docking tools, where the quality of the models is based on their similarity to a single reference experimental structure. However, using a static reference does not take into account the dynamic quality of the protein interface. Here, we used all-atom classical Molecular Dynamics simulations to investigate the stability of the reference interface for three complexes that previously served as targets in the CAPRI competition. For each one of these targets, we also ran MD simulations for ten models that are distributed over the High, Medium and Acceptable accuracy categories. To assess the quality of these models from a dynamic perspective, we set up new criteria which take into account the stability of the reference experimental protein interface. We show that, when the protein interfaces are allowed to evolve along time, the original ranking based on the static CAPRI criteria no longer holds as over 50% of the docking models undergo a category change (which can be either toward a better or a lower accuracy group) when reassessing their quality using dynamic information.

Weaving DNA strands: structural insight on ATP hydrolysis in RecA-induced homologous recombination.

Homologous recombination is a fundamental process in all living organisms that allows the faithful repair of DNA double strand breaks, through the exchange of DNA strands between homologous regions of the genome. Results of three decades of investigation and recent fruitful observations have unveiled key elements of the reaction mechanism, which proceeds along nucleofilaments of recombinase proteins of the RecA family. Yet, one essential aspect of homologous recombination has largely been overlooked when deciphering the mechanism: while ATP is hydrolyzed in large quantity during the process, how exactly hydrolysis influences the DNA strand exchange reaction at the structural level remains to be elucidated. In this study, we build on a previous geometrical approach that studied the RecA filament variability without bound DNA to examine the putative implication of ATP hydrolysis on the structure, position, and interactions of up to three DNA strands within the RecA nucleofilament. Simulation results on modeled intermediates in the ATP cycle bring important clues about how local distortions in the DNA strand geometries resulting from ATP hydrolysis can aid sequence recognition by promoting local melting of already formed DNA heteroduplex and transient reverse strand exchange in a weaving type of mechanism.

The positioning of Chi sites allows the RecBCD pathway to suppress some genomic rearrangements.

Bacterial recombinational repair of double-strand breaks often begins with creation of initiating 3' single-stranded DNA (ssDNA) tails on each side of a double-strand break (DSB). Importantly, if the RecBCD pathway is followed, RecBCD creates a gap between the sequences at 3' ends of the initiating strands. The gap flanks the DSB and extends at least to the nearest Chi site on each strand. Once the initiating strands form ssDNA-RecA filaments, each ssDNA-RecA filament searches for homologous double-stranded DNA (dsDNA) to use as a template for the DNA synthesis needed to fill the gap created by RecBCD. Our experimental results show that the DNA synthesis requires formation of a heteroduplex dsDNA that pairs >20 contiguous bases in the initiating strand with sequence matched bases in a strand from the original dsDNA. To trigger synthesis, the heteroduplex must be near the 3' end of the initiating strand. Those experimentally determined requirements for synthesis combined with the Chi site dependence of the function of RecBCD and the distribution of Chi sites in bacterial genomes could allow the RecBCD pathway to avoid some genomic rearrangements arising from directly induced DSBs; however, the same three factors could promote other rearrangements.

Slow extension of the invading DNA strand in a D-loop formed by RecA-mediated homologous recombination may enhance recognition of DNA homology.

DNA recombination resulting from RecA-mediated strand exchange aided by RecBCD proteins often enables accurate repair of DNA double-strand breaks. However, the process of recombinational repair between short DNA regions of accidental similarity can lead to fatal genomic rearrangements. Previous studies have probed how effectively RecA discriminates against interactions involving a short similar sequence that is embedded in otherwise dissimilar sequences but have not yielded fully conclusive results. Here, we present results of in vitro experiments with fluorescent probes strategically located on the interacting DNA fragments used for recombination. Our findings suggest that DNA synthesis increases the stability of the recombination products. Fluorescence measurements can also probe the homology dependence of the extension of invading DNA strands in D-loops formed by RecA-mediated strand exchange. We examined the slow extension of the invading strand in a D-loop by DNA polymerase (Pol) IV and the more rapid extension by DNA polymerase LF-Bsu We found that when DNA Pol IV extends the invading strand in a D-loop formed by RecA-mediated strand exchange, the extension afforded by 82 bp of homology is significantly longer than the extension on 50 bp of homology. In contrast, the extension of the invading strand in D-loops by DNA LF-Bsu Pol is similar for intermediates with ≥50 bp of homology. These results suggest that fatal genomic rearrangements due to the recombination of small regions of accidental homology may be reduced if RecA-mediated strand exchange is immediately followed by DNA synthesis by a slow polymerase.

Residues in the fingers domain of the translesion DNA polymerase DinB enable its unique participation in error-prone double-strand break repair.

The evolutionarily conserved Escherichia coli translesion DNA polymerase IV (DinB) is one of three enzymes that can bypass potentially deadly DNA lesions on the template strand during DNA replication. Remarkably, however, DinB is the only known translesion DNA polymerase active in RecA-mediated strand exchange during error-prone double-strand break repair. In this process, a single-stranded DNA (ssDNA)-RecA nucleoprotein filament invades homologous dsDNA, pairing the ssDNA with the complementary strand in the dsDNA. When exchange reaches the 3' end of the ssDNA, a DNA polymerase can add nucleotides onto the end, using one strand of dsDNA as a template and displacing the other. It is unknown what makes DinB uniquely capable of participating in this reaction. To explore this topic, we performed molecular modeling of DinB's interactions with the RecA filament during strand exchange, identifying key contacts made with residues in the DinB fingers domain. These residues are highly conserved in DinB, but not in other translesion DNA polymerases. Using a novel FRET-based assay, we found that DinB variants with mutations in these conserved residues are less effective at stabilizing RecA-mediated strand exchange than native DinB. Furthermore, these variants are specifically deficient in strand displacement in the absence of RecA filament. We propose that the amino acid patch of highly conserved residues in DinB-like proteins provides a mechanistic explanation for DinB's function in strand exchange and improves our understanding of recombination by providing evidence that RecA plays a role in facilitating DinB's activity during strand exchange.

RecA requires two molecules of Mg2+ ions for its optimal strand exchange activity in vitro.

Mg2+ ion stimulates the DNA strand exchange reaction catalyzed by RecA, a key step in homologous recombination. To elucidate the molecular mechanisms underlying the role of Mg2+ and the strand exchange reaction itself, we investigated the interaction of RecA with Mg2+ and sought to determine which step of the reaction is affected. Thermal stability, intrinsic fluorescence, and native mass spectrometric analyses of RecA revealed that RecA binds at least two Mg2+ ions with KD ≈ 2 mM and 5 mM. Deletion of the C-terminal acidic tail of RecA made its thermal stability and fluorescence characteristics insensitive to Mg2+ and similar to those of full-length RecA in the presence of saturating Mg2+. These observations, together with the results of a molecular dynamics simulation, support the idea that the acidic tail hampers the strand exchange reaction by interacting with other parts of RecA, and that binding of Mg2+ to the tail prevents these interactions and releases RecA from inhibition. We observed that binding of the first Mg2+ stimulated joint molecule formation, whereas binding of the second stimulated progression of the reaction. Thus, RecA is actively involved in the strand exchange step as well as bringing the two DNAs close to each other.

ATP hydrolysis provides functions that promote rejection of pairings between different copies of long repeated sequences.

During DNA recombination and repair, RecA family proteins must promote rapid joining of homologous DNA. Repeated sequences with >100 base pair lengths occupy more than 1% of bacterial genomes; however, commitment to strand exchange was believed to occur after testing ∼20-30 bp. If that were true, pairings between different copies of long repeated sequences would usually become irreversible. Our experiments reveal that in the presence of ATP hydrolysis even 75 bp sequence-matched strand exchange products remain quite reversible. Experiments also indicate that when ATP hydrolysis is present, flanking heterologous dsDNA regions increase the reversibility of sequence matched strand exchange products with lengths up to ∼75 bp. Results of molecular dynamics simulations provide insight into how ATP hydrolysis destabilizes strand exchange products. These results inspired a model that shows how pairings between long repeated sequences could be efficiently rejected even though most homologous pairings form irreversible products.

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

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.

Mobility and Core-Protein Binding Patterns of Disordered C-Terminal Tails in ß-Tubulin Isotypes.

Although they play a significant part in the regulation of microtubule structure, dynamics, and function, the disordered C-terminal tails of tubulin remain invisible to experimental structural methods and do not appear in the crystallographic structures that are currently available in the Protein Data Bank. Interestingly, these tails concentrate most of the sequence variability between tubulin isotypes and are the sites of the principal post-translational modifications undergone by this protein. Using homology modeling, we developed two complete models for the human αI/ßI- and αI/ßIII-tubulin isotypes that include their C-terminal tails. We then investigated the conformational variability of the two ß-tails using long time-scale classical molecular dynamics simulations that revealed similar features, notably the unexpected presence of common anchoring regions on the surface of the tuulin dimer, but also distinctive mobility or interaction patterns, some of which could be related to the tail lengths and charge distributions. We also observed in our simulations that the C-terminal tail from the ßI isotype, but not the ßIII isotype, formed contacts in the putative binding site of a recently discovered peptide that disrupts microtubule formation in glioma cells. Hindering the binding site in the ßI isotype would be consistent with this peptide's preferential disruption of microtubule formation in glioma, whose cells overexpress ßIII, compared to normal glial cells. While these observations need to be confirmed with more intensive sampling, our study opens new perspectives for the development of isotype-specific chemotherapy drugs.

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 poor homology stringency in the heteroduplex allows strand exchange to incorporate desirable mismatches without sacrificing recognition in vivo.

RecA family proteins are responsible for homology search and strand exchange. In bacteria, homology search begins after RecA binds an initiating single-stranded DNA (ssDNA) in the primary DNA-binding site, forming the presynaptic filament. Once the filament is formed, it interrogates double-stranded DNA (dsDNA). During the interrogation, bases in the dsDNA attempt to form Watson-Crick bonds with the corresponding bases in the initiating strand. Mismatch dependent instability in the base pairing in the heteroduplex strand exchange product could provide stringent recognition; however, we present experimental and theoretical results suggesting that the heteroduplex stability is insensitive to mismatches. We also present data suggesting that an initial homology test of 8 contiguous bases rejects most interactions containing more than 1/8 mismatches without forming a detectable 20 bp product. We propose that, in vivo, the sparsity of accidental sequence matches allows an initial 8 bp test to rapidly reject almost all non-homologous sequences. We speculate that once the initial test is passed, the mismatch insensitive binding in the heteroduplex allows short mismatched regions to be incorporated in otherwise homologous strand exchange products even though sequences with less homology are eventually rejected.

Investigating the Structural Variability and Binding Modes of the Glioma Targeting NFL-TBS.40-63 Peptide on Tubulin.

NFL-TBS.40-63 is a 24 amino acid peptide corresponding to the tubulin-binding site located on the light neurofilament subunit, which selectively enters glioblastoma cells, where it disrupts their microtubule network and inhibits their proliferation. We investigated its structural variability and binding modes on a tubulin heterodimer using a combination of NMR experiments, docking, and molecular dynamics (MD) simulations. Our results show that, while lacking a stable structure, the peptide preferentially binds on a specific single site located near the ß-tubulin C-terminal end, thus giving us precious hints regarding the mechanism of action of the NFL-TBS.40-63 peptide's antimitotic activity at the molecular level.

Highly mismatch-tolerant homology testing by RecA could explain how homology length affects recombination.

In E. coli, double strand breaks (DSBs) are resected and loaded with RecA protein. The genome is then rapidly searched for a sequence that is homologous to the DNA flanking the DSB. Mismatches in homologous partners are rare, suggesting that RecA should rapidly reject mismatched recombination products; however, this is not the case. Decades of work have shown that long lasting recombination products can include many mismatches. In this work, we show that in vitro RecA forms readily observable recombination products when 16% of the bases in the product are mismatched. We also consider various theoretical models of mismatch-tolerant homology testing. The models test homology by comparing the sequences of Ltest bases in two single-stranded DNAs (ssDNA) from the same genome. If the two sequences pass the homology test, the pairing between the two ssDNA becomes permanent. Stringency is the fraction of permanent pairings that join ssDNA from the same positions in the genome. We applied the models to both randomly generated genomes and bacterial genomes. For both randomly generated genomes and bacterial genomes, the models show that if no mismatches are accepted stringency is ∼ 99% when Ltest = 14 bp. For randomly generated genomes, stringency decreases with increasing mismatch tolerance, and stringency improves with increasing Ltest. In contrast, in bacterial genomes when Ltest ∼ 75 bp, stringency is ∼ 99% for both mismatch-intolerant and mismatch-tolerant homology testing. Furthermore, increasing Ltest does not improve stringency because most incorrect pairings join different copies of repeats. In sum, for bacterial genomes highly mismatch tolerant homology testing of 75 bp provides the same stringency as homology testing that rejects all mismatches and testing more than ∼75 base pairs is not useful. Interestingly, in vivo commitment to recombination typically requires homology testing of ∼ 75 bp, consistent with highly mismatch intolerant testing.

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.

Building Biological Relevance Into Integrative Modelling of Macromolecular Assemblies.

Recent advances in structural biophysics and integrative modelling methods now allow us to decipher the structures of large macromolecular assemblies. Understanding the dynamics and mechanisms involved in their biological function requires rigorous integration of all available data. We have developed a complete modelling pipeline that includes analyses to extract biologically significant information by consistently combining automated and interactive human-guided steps. We illustrate this idea with two examples. First, we describe the ryanodine receptor, an ion channel that controls ion flux across the cell membrane through transitions between open and closed states. The conformational changes associated with the transitions are small compared to the considerable system size of the receptor; it is challenging to consistently track these states with the available cryo-EM structures. The second example involves homologous recombination, in which long filaments of a recombinase protein and DNA catalyse the exchange of homologous DNA strands to reliably repair DNA double-strand breaks. The nucleoprotein filament reaction intermediates in this process are short-lived and heterogeneous, making their structures particularly elusive. The pipeline we describe, which incorporates experimental and theoretical knowledge combined with state-of-the-art interactive and immersive modelling tools, can help overcome these challenges. In both examples, we point to new insights into biological processes that arise from such interdisciplinary approaches.

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.

When Order Meets Disorder: Modeling and Function of the Protein Interface in Fuzzy Complexes.

The degree of proteins structural organization ranges from highly structured, compact folding to intrinsic disorder, where each degree of self-organization corresponds to specific functions: well-organized structural motifs in enzymes offer a proper environment for precisely positioned functional groups to participate in catalytic reactions; at the other end of the self-organization spectrum, intrinsically disordered proteins act as binding hubs via the formation of multiple, transient and often non-specific interactions. This review focusses on cases where structurally organized proteins or domains associate with highly disordered protein chains, leading to the formation of interfaces with varying degrees of fuzziness. We present a review of the computational methods developed to provide us with information on such fuzzy interfaces, and how they integrate experimental information. The discussion focusses on two specific cases, microtubules and homologous recombination nucleoprotein filaments, where a network of intrinsically disordered tails exerts regulatory function in recruiting partner macromolecules, proteins or DNA and tuning the atomic level association. Notably, we show how computational approaches such as molecular dynamics simulations can bring new knowledge to help bridging the gap between experimental analysis, that mostly concerns ensemble properties, and the behavior of individual disordered protein chains that contribute to regulation functions.