Machine learning to estimate the local quality of protein crystal structures.

Low-resolution electron density maps can pose a major obstacle in the determination and use of protein structures. Herein, we describe a novel method, called quality assessment based on an electron density map (QAEmap), which evaluates local protein structures determined by X-ray crystallography and could be applied to correct structural errors using low-resolution maps. QAEmap uses a three-dimensional deep convolutional neural network with electron density maps and their corresponding coordinates as input and predicts the correlation between the local structure and putative high-resolution experimental electron density map. This correlation could be used as a metric to modify the structure. Further, we propose that this method may be applied to evaluate ligand binding, which can be difficult to determine at low resolution.

Real-time tracking reveals catalytic roles for the two DNA binding sites of Rad51.

During homologous recombination, Rad51 forms a nucleoprotein filament on single-stranded DNA to promote DNA strand exchange. This filament binds to double-stranded DNA (dsDNA), searches for homology, and promotes transfer of the complementary strand, producing a new heteroduplex. Strand exchange proceeds via two distinct three-strand intermediates, C1 and C2. C1 contains the intact donor dsDNA whereas C2 contains newly formed heteroduplex DNA. Here, we show that the conserved DNA binding motifs, loop 1 (L1) and loop 2 (L2) in site I of Rad51, play distinct roles in this process. L1 is involved in formation of the C1 complex whereas L2 mediates the C1-C2 transition, producing the heteroduplex. Another DNA binding motif, site II, serves as the DNA entry position for initial Rad51 filament formation, as well as for donor dsDNA incorporation. Our study provides a comprehensive molecular model for the catalytic process of strand exchange mediated by eukaryotic RecA-family recombinases.

Combination of coarse-grained molecular dynamics simulations and small-angle X-ray scattering experiments.

The combination of molecular dynamics (MD) simulations and small-angle X-ray scattering (SAXS), called the MD-SAXS method, is efficient for investigating protein dynamics. To overcome the time-scale limitation of all-atom MD simulations, coarse-grained (CG) representations are often utilized for biomolecular simulations. In this study, we propose a method to combine CG MD simulations with SAXS, termed the CG-MD-SAXS method. In the CG-MD-SAXS method, the scattering factors of CG particles for proteins and nucleic acids are evaluated using high-resolution structural data in the Protein Data Bank, and the excluded volume and the hydration shell are modeled using two adjustable parameters to incorporate solvent effects. To avoid overfitting, only the two parameters are adjusted for an entire structure ensemble. To verify the developed method, theoretical SAXS profiles for various proteins, DNA/RNA, and a protein-RNA complex are compared with both experimental profiles and theoretical profiles obtained by the all-atom representation. In the present study, we applied the CG-MD-SAXS method to the Swi5-Sfr1 complex and three types of nucleosomes to obtain reliable ensemble models consistent with the experimental SAXS data.

Rotation Mechanism of Molecular Motor V1-ATPase Studied by Multiscale Molecular Dynamics Simulation.

Enterococcus hirae V1-ATPase is a molecular motor composed of the A3B3 hexamer ring and the central stalk. In association with ATP hydrolysis, three catalytic AB pairs in the A3B3 ring undergo conformational changes, which lead to a 120° rotation of the central stalk. To understand how the conformational changes of three catalytic pairs induce the 120° rotation of the central stalk, we performed multiscale molecular dynamics (MD) simulations in which coarse-grained and all-atom MD simulations were combined using a fluctuation matching methodology. During the rotation, a catalytic AB pair spontaneously adopted an intermediate conformation, which was not included in the initial inputs of the simulations and was essentially close to the "bindable-like" structure observed in a recently solved crystal structure. Furthermore, the creation of a space between the bindable-like and tight pairs was required for the central stalk to rotate without steric hindrance. These cooperative rearrangements of the three catalytic pairs are crucial for the rotation of the central stalk.

Molecular modeling and molecular dynamics simulations of recombinase Rad51.

The Rad51 ATPase plays central roles in DNA homologous recombination. Yeast Rad51 dimer structure in the active form of the filament was constructed using homology modeling techniques, and all-atom molecular dynamics (MD) simulations were performed using the modeled structure. We found two crucial interaction networks involving ATP: one is among the γ-phosphate of ATP, K(+) ions, H352, and D374; the other is among the adenine ring of ATP, R228, and P379. Multiple MD simulations were performed in which the number of bound K(+) ions was changed. The simulated structures suggested that K(+) ions are indispensable for the stabilization of the active dimer and resemble the arginine and lysine fingers of other P-loop containing ATPases and GTPases. MD simulations also showed that the adenine ring of ATP mediates interactions between adjacent protomers. Furthermore, in MD simulations starting from a structure just after ATP hydrolysis, the opening motion corresponding to dissociation from DNA was observed. These results support the hypothesis that ATP and K(+) ions function as glue between protomers.

Characterisation of an intrinsically disordered protein complex of Swi5-Sfr1 by ion mobility mass spectrometry and small-angle X-ray scattering.

It is now recognized that intrinsically disordered proteins (IDPs) play important roles as hubs in intracellular networks, and their structural characterisation is of significance. However, due to their highly dynamic features, it is challenging to investigate the structures of IDPs solely by conventional methods. In the present study, we demonstrate a novel method to characterise protein complexes using electrospray ionization ion mobility mass spectrometry (ESI-IM-MS) in combination with small-angle X-ray scattering (SAXS). This method enables structural characterisation even of proteins that have difficulties in crystallisation. With this method, we have characterised the Schizosaccharomyces pombe Swi5-Sfr1 complex, which is expected to have a long disordered region at the N-terminal portion of Sfr1. ESI-IM-MS analysis of the Swi5-Sfr1 complex revealed that its experimental collision cross-section (CCS) had a wide distribution, and the CCS values of the most dominant ions were ∼56% of the theoretically calculated value based on the SAXS low-resolution model, suggesting a significant size reduction in the gas phase. The present study demonstrates that the newly developed method for calculation of the theoretical CCSs of the SAXS low-resolution models of proteins allows accurate evaluation of the experimental CCS values of IDPs provided by ESI-IM-MS by comparing with the low-resolution solution structures. Furthermore, it was revealed that the combination of ESI-IM-MS and SAXS is a promising method for structural characterisation of protein complexes that are unable to crystallise.

Fission yeast Swi5-Sfr1 protein complex, an activator of Rad51 recombinase, forms an extremely elongated dogleg-shaped structure.

In eukaryotes, DNA strand exchange is the central reaction of homologous recombination, which is promoted by Rad51 recombinases forming a right-handed nucleoprotein filament on single-stranded DNA, also known as a presynaptic filament. Accessory proteins known as recombination mediators are required for the formation of the active presynaptic filament. One such mediator in the fission yeast Schizosaccharomyces pombe is the Swi5-Sfr1 complex, which has been identified as an activator of Rad51 that assists in presynaptic filament formation and stimulates its strand exchange reaction. Here, we determined the 1:1 binding stoichiometry between the two subunits of the Swi5-Sfr1 complex using analytical ultracentrifugation and electrospray ionization mass spectrometry. Small-angle x-ray scattering experiments revealed that the Swi5-Sfr1 complex displays an extremely elongated dogleg-shaped structure in solution, which is consistent with its exceptionally high frictional ratio (f/f(0)) of 2.0 ± 0.2 obtained by analytical ultracentrifugation. Furthermore, we determined a rough topology of the complex by comparing the small-angle x-ray scattering-based structures of the Swi5-Sfr1 complex and four Swi5-Sfr1-Fab complexes, in which the Fab fragments of monoclonal antibodies were specifically bound to experimentally determined sites of Sfr1. We propose a model for how the Swi5-Sfr1 complex binds to the Rad51 filament, in which the Swi5-Sfr1 complex fits into the groove of the Rad51 filament, leading to an active and stable presynaptic filament.

Homology-modelled structure of the ßB2B3-crystallin heterodimer studied by ion mobility and radical probe MS.

Ion mobility MS was employed to study the structure of the ßB2B3-crystallin heterodimer following its detection by ESI-TOF MS. The results demonstrate that the heterodimer has a similar cross-section (3 165 Å(2)) and structure to the ßB2B2-crystallin homodimer. Several homology-modelled structures for the ßB2B3 heterodimer were constructed and assessed in terms of their calculated collision cross-sections and whether the solvent accessibilities of reactive amino acid side chains throughout the ßB3 subunit are in accord with measured oxidation levels in radical probe MS protein footprinting experiments. The ßB2B3 heterodimer AD model provides the best representation of the heterodimer's structure overall following a consideration of both the ion mobility and radical probe MS data.