Cryo-EM structure and B-factor refinement with ensemble representation.

Cryo-EM experiments produce images of macromolecular assemblies that are combined to produce three-dimensional density maps. Typically, atomic models of the constituent molecules are fitted into these maps, followed by a density-guided refinement. We introduce TEMPy-ReFF, a method for atomic structure refinement in cryo-EM density maps. Our method represents atomic positions as components of a Gaussian mixture model, utilising their variances as B-factors, which are used to derive an ensemble description. Extensively tested on a substantial dataset of 229 cryo-EM maps from EMDB ranging in resolution from 2.1-4.9 Å with corresponding PDB and CERES atomic models, our results demonstrate that TEMPy-ReFF ensembles provide a superior representation of cryo-EM maps. On a single-model basis, it performs similarly to the CERES re-refinement protocol, although there are cases where it provides a better fit to the map. Furthermore, our method enables the creation of composite maps free of boundary artefacts. TEMPy-ReFF is useful for better interpretation of flexible structures, such as those involving RNA, DNA or ligands.

CASP15 cryo-EM protein and RNA targets: Refinement and analysis using experimental maps.

CASP assessments primarily rely on comparing predicted coordinates with experimental reference structures. However, experimental structures by their nature are only models themselves-their construction involves a certain degree of subjectivity in interpreting density maps and translating them to atomic coordinates. Here, we directly utilized density maps to evaluate the predictions by employing a method for ranking the quality of protein chain predictions based on their fit into the experimental density. The fit-based ranking was found to correlate well with the CASP assessment scores. Overall, the evaluation against the density map indicated that the models are of high accuracy, and occasionally even better than the reference structure in some regions of the model. Local assessment of predicted side chains in a 1.52 Å resolution map showed that side-chains are sometimes poorly positioned. Additionally, the top 118 predictions associated with 9 protein target reference structures were selected for automated refinement, in addition to the top 40 predictions for 11 RNA targets. For both proteins and RNA, the refinement of CASP15 predictions resulted in structures that are close to the reference target structure. This refinement was successful despite large conformational changes often being required, showing that predictions from CASP-assessed methods could serve as a good starting point for building atomic models in cryo-EM maps for both proteins and RNA. Loop modeling continued to pose a challenge for predictors, and together with the lack of consensus amongst models in these regions suggests that modeling, in combination with model-fit to the density, holds the potential for identifying more flexible regions within the structure.

RIBFIND2: Identifying rigid bodies in protein and nucleic acid structures.

Molecular structures are often fitted into cryo-EM maps by flexible fitting. When this requires large conformational changes, identifying rigid bodies can help optimize the model-map fit. Tools for identifying rigid bodies in protein structures exist, however an equivalent for nucleic acid structures is lacking. With the increase in cryo-EM maps containing RNA and progress in RNA structure prediction, there is a need for such tools. We previously developed RIBFIND, a program for clustering protein secondary structures into rigid bodies. In RIBFIND2, this approach is extended to nucleic acid structures. RIBFIND2 can identify biologically relevant rigid bodies in important groups of complex RNA structures, capturing a wide range of dynamics, including large rigid-body movements. The usefulness of RIBFIND2-assigned rigid bodies in cryo-EM model refinement was demonstrated on three examples, with two conformations each: Group II Intron complexed IEP, Internal Ribosome Entry Site and the Processome, using cryo-EM maps at 2.7-5 Å resolution. A hierarchical refinement approach, performed on progressively smaller sets of RIBFIND2 rigid bodies, was clearly shown to have an advantage over classical all-atom refinement. RIBFIND2 is available via a web server with structure visualization and as a standalone tool.

Cooperative amyloid fibre binding and disassembly by the Hsp70 disaggregase.

Although amyloid fibres are highly stable protein aggregates, a specific combination of human Hsp70 system chaperones can disassemble them, including fibres formed of α-synuclein, huntingtin, or Tau. Disaggregation requires the ATPase activity of the constitutively expressed Hsp70 family member, Hsc70, together with the J domain protein DNAJB1 and the nucleotide exchange factor Apg2. Clustering of Hsc70 on the fibrils appears to be necessary for disassembly. Here we use atomic force microscopy to show that segments of in vitro assembled α-synuclein fibrils are first coated with chaperones and then undergo bursts of rapid, unidirectional disassembly. Cryo-electron tomography and total internal reflection fluorescence microscopy reveal fibrils with regions of densely bound chaperones, preferentially at one end of the fibre. Sub-stoichiometric amounts of Apg2 relative to Hsc70 dramatically increase recruitment of Hsc70 to the fibres, creating localised active zones that then undergo rapid disassembly at a rate of ~ 4 subunits per second. The observed unidirectional bursts of Hsc70 loading and unravelling may be explained by differences between the two ends of the polar fibre structure.

Base-pair resolution analysis of the effect of supercoiling on DNA flexibility and major groove recognition by triplex-forming oligonucleotides.

In the cell, DNA is arranged into highly-organised and topologically-constrained (supercoiled) structures. It remains unclear how this supercoiling affects the detailed double-helical structure of DNA, largely because of limitations in spatial resolution of the available biophysical tools. Here, we overcome these limitations, by a combination of atomic force microscopy (AFM) and atomistic molecular dynamics (MD) simulations, to resolve structures of negatively-supercoiled DNA minicircles at base-pair resolution. We observe that negative superhelical stress induces local variation in the canonical B-form DNA structure by introducing kinks and defects that affect global minicircle structure and flexibility. We probe how these local and global conformational changes affect DNA interactions through the binding of triplex-forming oligonucleotides to DNA minicircles. We show that the energetics of triplex formation is governed by a delicate balance between electrostatics and bonding interactions. Our results provide mechanistic insight into how DNA supercoiling can affect molecular recognition, that may have broader implications for DNA interactions with other molecular species.

TopoStats - A program for automated tracing of biomolecules from AFM images.

We present TopoStats, a Python toolkit for automated editing and analysis of Atomic Force Microscopy images. The program automates identification and tracing of individual molecules in circular and linear conformations without user input. TopoStats was able to identify and trace a range of molecules within AFM images, finding, on average, ~90% of all individual molecules and molecular assemblies within a wide field of view, and without the need for prior processing. DNA minicircles of varying size, DNA origami rings and pore forming proteins were identified and accurately traced with contour lengths of traces typically within 10 nm of the predicted contour length. TopoStats was also able to reliably identify and trace linear and enclosed circular molecules within a mixed population. The program is freely available via GitHub (https://github.com/afm-spm/TopoStats) and is intended to be modified and adapted for use if required.

Exploring Step-by-Step Assembly of Nanoparticle:Cytochrome Biohybrid Photoanodes.

Coupling light-harvesting semiconducting nanoparticles (NPs) with redox enzymes has been shown to create artificial photosynthetic systems that hold promise for the synthesis of solar fuels. High quantum yields require efficient electron transfer from the nanoparticle to the redox protein, a property that can be difficult to control. Here, we have compared binding and electron transfer between dye-sensitized TiO2 nanocrystals or CdS quantum dots and two decaheme cytochromes on photoanodes. The effect of NP surface chemistry was assessed by preparing NPs capped with amine or carboxylic acid functionalities. For the TiO2 nanocrystals, binding to the cytochromes was optimal when capped with a carboxylic acid ligand, whereas for the CdS QDs, better adhesion was observed for amine capped ligand shells. When using TiO2 nanocrystals, dye-sensitized with a phosphonated bipyridine Ru(II) dye, photocurrents are observed that are dependent on the redox state of the decaheme, confirming that electrons are transferred from the TiO2 nanocrystals to the surface via the decaheme conduit. In contrast, when CdS NPs are used, photocurrents are not dependent on the redox state of the decaheme, consistent with a model in which electron transfer from CdS to the photoanode bypasses the decaheme protein. These results illustrate that although the organic shell of NPs nanoparticles crucially affects coupling with proteinaceous material, the coupling can be difficult to predict or engineer.