Search

Single-particle cryogenic electron microscopy (cryo-EM) has become a powerful mainstay tool in high resolution structural biology thanks to advances in hardware, software and sample preparation technology. In my thesis, I utilized this technique to unravel the function of various challenging biological macromolecules. My first focus was bacterial ribosomal biogenesis: understanding how bacteria assemble their ribosomes. Ribosomes are the factories of the cell, responsible for manufacturing all proteins. Ribosomes themselves are huge, with the bacterial version made of 52 proteins and 4566 RNA nucleotides. How these components assemble has long been a mystery. Early groundbreaking work sketched out a biogenesis pathway using purified components in vitro – but under non-physiological conditions. We sought to understand how the bacterial ribosome – specifically the large subunit 50S – is built inside the cell. To achieve this, we engineered a conditional knock-out bacterial strain that lacked one specific ribosomal protein (L17). This caused the cells to accumulate incomplete intermediates along the 50S biogenesis pathway. These intermediates were purified and examined with mass spectrometry and single-particle cryo-EM. Two major hurdles arose in this project: firstly, the biogenesis intermediates exhibited a preferred orientation when vitrified for cryo-EM analysis. This means that instead of showing many different views required for reconstruction of the 3D structure, the intermediates only adopted one view on the cryo-EM grid. To overcome this problem, we engineered a method to induce additional views on the microscope by tilting the stage. Using another test protein that also exhibited preferred orientation (hemagglutinin), we optimized and characterized this new tilt methodology and showed it was generally applicable to overcoming preferred orientation, regardless of type of specimen. We also created a software tool, called 3DFSC (3dfsc.salk.edu), for other microscopists to calculate the degree of directional anisotropy in their structures due to preferred orientation. Using this tilt strategy finally enabled the structural elucidation of our 50S intermediates. The second challenge in the project was the large amount of heterogeneity present in the sample. Through hierarchical 3D classification schemes using the latest software tools, we obtained 14 different 50S intermediate structures, all from imaging a single cryo-EM grid. By analyzing the missing components of each intermediate, and corroborating these observations with mass spectrometry data, we outlined the first in vivo 50S assembly pathway, and showed that ribosome assembly occurs step-wise and in parallel pathways. My second focus was on pushing the resolution limits of single-particle cryo-EM using adeno-associated virus (AAV) serotype 2 homogeneous virus-like particles (VLPs) that lack DNA. Exploiting several technical advances to improve resolution, including use of gold grids, per-particle CTF refinement, and correction for Ewald sphere curvature, we managed to obtain a 1.86 Å resolution reconstruction of the AAV2L336C variant VLP, the highest resolution icosahedral virus reconstruction solved by single-particle cryo-EM to date. Using our structure, we were able to show improvements using Ewald sphere curvature correction and shed light on the mechanistic basis as to why the L336C mutation resulted in defects in genome packaging and infectivity compared to the WT viral particles. My third focus was the understanding of small membrane proteins involved in infectious diseases. Membrane proteins are a challenge to work with due to the need for them to be extracted from the lipid bilayer for studies as compared to soluble proteins. Infectious diseases have a huge burden on society, with the top three infectious agents accounting for 2.7 million deaths in 2016. The third most deadly infectious disease is malaria, a mosquito-borne parasite which kills 450,000 people annually. One drug used early on for treating malaria was chloroquine but its usefulness waned due to development of resistance. Chloroquine resistance is mediated by the chloroquine resistance transporter (PfCRT). Although small (49 kDa) for single-particle cryo-EM, we solved its structure by using fragment antibody technology to add mass and help with image alignment and 3D reconstruction. The 3.2 Å structure resembles other drug metabolite transporters, and the chloroquine resistance mutations map to a ring around the central cavity, suggesting this central pore as the drug binding site. Tuberculosis (TB) is the top killer, above malaria and HIV/AIDS, being responsible for 1.3 million deaths. In TB, a common antibiotic target is the bacterium’s cell wall synthesis machinery. One family of such enzymes is the arabinosyltransferases, which synthesize the critical arabinose sugars. Using single-particle cryo-EM, we solved two high resolution structures of one such essential enzyme, AftD. Due to the low yield of the protein, a picoliter automated sample dispensing robot was crucial to allow for initial cryo-EM analysis. We then performed mutagenesis studies in M. smegmatis, a TB model organism, which uncovered the critical amino acid residues in the active site and determined that a bound acyl-carrier-protein was likely involved in allosteric inhibition of AftD’s active site. Another member of the family, EmbB, is the target of a widely used frontline TB drug called ethambutol. We have solved the high resolution structures of the apo and putative drug-bound states of EmbB, allowing us to map out, for the first time, both the active site and drug-resistance mutations of this crucial enzyme. The atomic structures of the functional pockets of Mycobacterial AftD and malarial PfCRT will hopefully enable structure-based drug design to improve existing drugs or potentially even develop new treatments against these infectious maladies. In conclusion, the continual and breathtaking improvements in single-particle cryo-EM methodology has been instrumental in allowing the elucidation of the aforementioned biological macromolecules from ribosome biogenesis intermediates, to AAV2 vehicle, Plasmodium drug resistance transporter to mycobacterial glycosyltransferases – structures of which help explain biological function.

Proteins and DNA are essential to life as we know it and understanding their function is understanding their structure and dynamics. The importance of the latter is being appreciated more in recent years and has led to the development of novel interdisciplinary techniques and approaches to studying protein function. Three techniques to study protein structure and dynamics have been used and combined in different ways in the context of this thesis and have led to a better understanding of the three systems described herein. X-ray crystallography is the oldest and still arguably most popular technique to study macromolecular structures. Nuclear magnetic resonance (NMR) spectroscopy is a not much younger technique that is a powerful tool not only to probe molecular structure but also dynamics. The last technique described herein are molecular dynamics (MD) simulations, which are only just growing out of their infancy. MD simulations are computer simulations of macromolecules based on structures solved by X-ray crystallography or NMR spectroscopy, that can give mechanistic insight into dynamic processes of macromolecules whose amplitudes can be estimated by the former two techniques. MD simulations of the model protein GB3 (B3 immunoglobulin-binding domain of streptococcal protein G) were conducted to identify origins of discrepancies between order parameters derived from different sets of MD simulations and NMR relaxation experiments.The results highlight the importance of time scales as well as sampling when comparing MD simulations to NMR experiments. Discrepancies are seen for unstructured regions like loops and termini and often correspond to nanosecond time scale transitions between conformational substates that are either over- or undersampled in simulation. Sampling biases can be somewhat remedied by running longer (microsecond time scale) simulations. However, some discrepancies persist over even very long trajectories. We show that these discrepancies can be due to the choice of the starting structure and more specifically even differences in protonation procedures. A test for convergence on the nanosecond time scale is shown to be able to correct for many of the observed discrepancies. Next, MD simulations were used to predict in vitro thermostability of members of the bacterial Ribonuclease HI (RNase H) family of endonucleases. Thermodynamic stability is a central requirement for protein function and a goal of protein engineering is improvement of stability, particularly for applications in biotechnology. The temperature dependence of the generalized order parameter, S, for four RNase H homologs, from psychrotrophic, mesophilic and thermophilic organisms, is highly correlated with experimentally determined melting temperatures and with calculated free energies of folding at the midpoint temperature of the simulations. This study provides an approach for in silico mutational screens to improve thermostability of biologically and industrially relevant enzymes. Lastly, we used a combination of X-ray crystallography, NMR spectroscopy and MD simulations to study specificity of the interaction between Drosophila Hox proteins and their DNA target sites. Hox proteins are transcription factors specifying segment identity during embryogenesis of bilaterian animals. The DNA binding homeodomains have been shown to confer specificity to the different Hox paralogs, while being very similar in sequence and structure. Our results underline earlier findings about the importance of the N-terminal arm and linker region of Hox homeodomains, the cofactor Exd, as well as DNA shape, for specificity. A comparison of predicted DNA shapes based on sequence alone with the shapes observed for different DNA target sequences in four crystal structures when in complex with the Drosophila Hox protein AbdB and the cofactor Exd, shows that a combined ”induced fit”/”conformational selection” mechanism is the most likely mechanism by which Hox homeodomains recognize DNA shape and achieve specificity. The minor groove widths for all sequences is close to identical for all ternary complexes found in the different crystal structures, whereas predicted shapes vary between the different DNA sequences. The sequences that have shown higher affinity to AbdB in vitro have a predicted DNA shape that matches the observed DNA shape in the ternary complexes more closely than the sequences that show low in vitro affinity to AbdB. This strongly suggests that the AbdB-Exd complex selects DNA sequences with a higher propensity to adopt the final shape in their unbound form, leading to higher affinity. An additional AbdB monomer binding site with a strongly preformed binding competent shape is observed for one of the oligomers in the reverse complement strand of one of the canonical (weak) Hox-Exd complex binding site. The shape preference seems strong enough for AbdB monomer binding to compete with AbdB-Exd dimer binding to that same oligomer, suggested by the presence of both binding modes in the same crystal. The monomer binding site is essentially able to compete with the dimer binding site, even though binding with the cofactor is not possible, because its shape is very close to the ideal shape. A comparison of different crystal structures solved herein and in the literature as well as a set of molecular dynamics simulations was performed and led to insights about the importance of residues in the Hox N-terminal arm for the preference of certain Hox paralogs to certain DNA shapes. Taken together all these insights contribute to our understanding of Hox specificity in particular as well as protein-DNA interactions in general.