Electron density sharpening as a general technique in crystallographic studies

Discovery of novel druggable pockets on polyomavirus VP1 through crystallographic fragment-based screening to develop capsid assembly inhibitors

Polyomaviruses are a family of ubiquitous double-stranded DNA viruses many of which are human pathogens. These include BK polyomavirus which causes severe urinary tract infection in immunocompromised patients and Merkel cell polyomavirus associated with aggressive cancers. The small genome of polyomaviruses lacks conventional drug targets, and no specific drugs are available at present. Here we focus on the main structural protein VP1 of BK polyomavirus which is responsible for icosahedral capsid formation. To provide a foundation towards rational drug design, we crystallized truncated VP1 pentamers and subjected them to a high-throughput screening for binding drug-like fragments through a direct X-ray analysis. To enable a highly performant screening, rigorous optimization of the crystallographic pipeline and processing with the latest generation PanDDA2 software were necessary. As a result, a total of 144 binding hits were established. Importantly, the hits are well clustered in six surface pockets. Three pockets are located on the outside of the pentamer and map on the regions where the 'invading' C-terminal arm of another pentamer is attached upon capsid assembly. Another set of three pockets is situated within the wide pore along the five-fold axis of the VP1 pentamer. These pockets are situated at the interaction interface with the minor capsid protein VP2 which is indispensable for normal functioning of the virus. Here we systematically analyse the three outside pockets which are highly conserved across various polyomaviruses, while point mutations in these pockets are detrimental for viral replication. We show that one of the pockets can accommodate antipsychotic drug trifluoperazine. For each pocket, we derive pharmacophore features which enable the design of small molecules preventing the interaction between VP1 pentamers and therefore inhibiting capsid assembly. Our data lay a foundation towards a rational development of first-in-class drugs targeting polyomavirus capsid.

Crystal structure of an RNA/DNA strand exchange junction

Short segments of RNA displace one strand of a DNA duplex during diverse processes including transcription and CRISPR-mediated immunity and genome editing. These strand exchange events involve the intersection of two geometrically distinct helix types-an RNA:DNA hybrid (A-form) and a DNA:DNA homoduplex (B-form). Although previous evidence suggests that these two helices can stack on each other, it is unknown what local geometric adjustments could enable A-on-B stacking. Here we report the X-ray crystal structure of an RNA-5'/DNA-3' strand exchange junction at an anisotropic resolution of 1.6 to 2.2 Å. The structure reveals that the A-to-B helical transition involves a combination of helical axis misalignment, helical axis tilting and compression of the DNA strand within the RNA:DNA helix, where nucleotides exhibit a mixture of A- and B-form geometry. These structural principles explain previous observations of conformational stability in RNA/DNA exchange junctions, enabling a nucleic acid architecture that is repeatedly populated during biological strand exchange events.

Catalytic and structural properties of ATP-dependent caprolactamase from Pseudomonas jessenii

Caprolactamase is the first enzyme in the caprolactam degradation pathway of Pseudomonas jessenii. It is composed of two subunits (CapA and CapB) and sequence-related to other ATP-dependent enzymes involved in lactam hydrolysis, like 5-oxoprolinases and hydantoinases. Low sequence similarity also exists with ATP-dependent acetone- and acetophenone carboxylases. The caprolactamase was produced in Escherichia coli, isolated by His-tag affinity chromatography, and subjected to functional and structural studies. Activity toward caprolactam required ATP and was dependent on the presence of bicarbonate in the assay buffer. The hydrolysis product was identified as 6-aminocaproic acid. Quantum mechanical modeling indicated that the hydrolysis of caprolactam was highly disfavored (ΔG0 '= 23 kJ/mol), which explained the ATP dependence. A crystal structure showed that the enzyme exists as an (αβ)2 tetramer and revealed an ATP-binding site in CapA and a Zn-coordinating site in CapB. Mutations in the ATP-binding site of CapA (D11A and D295A) significantly reduced product formation. Mutants with substitutions in the metal binding site of CapB (D41A, H99A, D101A, and H124A) were inactive and less thermostable than the wild-type enzyme. These residues proved to be essential for activity and on basis of the experimental findings we propose possible mechanisms for ATP-dependent lactam hydrolysis.

Guest-protein incorporation into solvent channels of a protein host crystal (hostal)

Soaking small molecules into the solvent channels of protein crystals is the most common method of obtaining crystalline complexes with ligands such as substrates or inhibitors. The solvent channels of some protein crystals are large enough to allow the incorporation of macromolecules, but soaking of protein guests into protein crystals has not been reported. Such protein host crystals (here given the name hostals) incorporating guest proteins may be useful for a wide range of applications in biotechnology, for example as cargo systems or for diffraction studies analogous to the crystal sponge method. The present study takes advantage of crystals of the Escherichia coli tryptophan repressor protein (ds-TrpR) that are extensively domain-swapped and suitable for incorporating guest proteins by diffusion, as they are robust and have large solvent channels. Confocal fluorescence microscopy is used to follow the migration of cytochrome c and fluorophore-labeled calmodulin into the solvent channels of ds-TrpR crystals. The guest proteins become uniformly distributed in the crystal within weeks and enriched within the solvent channels. X-ray diffraction studies on host crystals with high concentrations of incorporated guests demonstrate that diffraction limits of ∼2.5 Å can still be achieved. Weak electron density is observed in the solvent channels, but the guest-protein structures could not be determined by conventional crystallographic methods. Additional approaches that increase the ordering of guests in the host crystal are discussed that may support protein structure determination using the hostal system in the future. This host system may also be useful for biotechnological applications where crystallographic order of the guest is not required.

Structural Basis for tRNA Mimicry by a Bacterial Y RNA

Noncoding Y RNAs are present in both animal cells and many bacteria. In all species examined, Y RNAs tether the Ro60 protein to an effector protein to perform various cellular functions. Recently, a new Y RNA subfamily was identified in bacteria. Bioinformatic analyses of these YrlA (Y RNA-like A) RNAs predict that the effector-binding domain resembles tRNA. We present the structure of this domain, the overall folding of which is strikingly similar to canonical tRNAs. The tertiary interactions that are responsible for stabilizing tRNA are present in YrlA, making it a close tRNA mimic. However, YrlA lacks a free CCA end and contains a kink in the stem corresponding to the anticodon stem. Since nucleotides in the D and T stems are conserved among YrlAs, they may be an interaction site for an unknown factor. Our experiments identify YrlA RNAs as a new class of tRNA mimics.

Automated map sharpening by maximization of detail and connectivity

An algorithm for automatic map sharpening is presented that is based on optimization of the detail and connectivity of the sharpened map. The detail in the map is reflected in the surface area of an iso-contour surface that contains a fixed fraction of the volume of the map, where a map with high level of detail has a high surface area. The connectivity of the sharpened map is reflected in the number of connected regions defined by the same iso-contour surfaces, where a map with high connectivity has a small number of connected regions. By combining these two measures in a metric termed the `adjusted surface area', map quality can be evaluated in an automated fashion. This metric was used to choose optimal map-sharpening parameters without reference to a model or other interpretations of the map. Map sharpening by optimization of the adjusted surface area can be carried out for a map as a whole or it can be carried out locally, yielding a locally sharpened map. To evaluate the performance of various approaches, a simple metric based on map-model correlation that can reproduce visual choices of optimally sharpened maps was used. The map-model correlation is calculated using a model with B factors (atomic displacement factors; ADPs) set to zero. This model-based metric was used to evaluate map sharpening and to evaluate map-sharpening approaches, and it was found that optimization of the adjusted surface area can be an effective tool for map sharpening.

Structural characterization of a highly-potent V3-glycan broadly neutralizing antibody bound to natively-glycosylated HIV-1 envelope

Broadly neutralizing antibodies (bNAbs) isolated from HIV-1-infected individuals inform HIV-1 vaccine design efforts. Developing bNAbs with increased efficacy requires understanding how antibodies interact with the native oligomannose and complex-type N-glycan shield that hides most protein epitopes on HIV-1 envelope (Env). Here we present crystal structures, including a 3.8-Å X-ray free electron laser dataset, of natively glycosylated Env trimers complexed with BG18, the most potent V3/N332gp120 glycan-targeting bNAb reported to date. Our structures show conserved contacts mediated by common D gene-encoded residues with the N332gp120 glycan and the gp120 GDIR peptide motif, but a distinct Env-binding orientation relative to PGT121/10-1074 bNAbs. BG18's binding orientation provides additional contacts with N392gp120 and N386gp120 glycans near the V3-loop base and engages protein components of the V1-loop. The BG18-natively-glycosylated Env structures facilitate understanding of bNAb-glycan interactions critical for using V3/N332gp120 bNAbs therapeutically and targeting their epitope for immunogen design.

Covalently linked dengue virus envelope glycoprotein dimers reduce exposure of the immunodominant fusion loop epitope

A problem in the search for an efficient vaccine against dengue virus is the immunodominance of the fusion loop epitope (FLE), a segment of the envelope protein E that is buried at the interface of the E dimers coating mature viral particles. Anti-FLE antibodies are broadly cross-reactive but poorly neutralizing, displaying a strong infection enhancing potential. FLE exposure takes place via dynamic ‘breathing' of E dimers at the virion surface. In contrast, antibodies targeting the E dimer epitope (EDE), readily exposed at the E dimer interface over the region of the conserved fusion loop, are very potent and broadly neutralizing. We here engineer E dimers locked by inter-subunit disulfide bonds, and show by X-ray crystallography and by binding to a panel of human antibodies that these engineered dimers do not expose the FLE, while retaining the EDE exposure. These locked dimers are strong immunogen candidates for a next-generation vaccine.

The immunodominant epitope of dengue virus envelope protein (E) induces poorly neutralizing antibodies, which poses a problem for vaccine development. Here, the authors engineer covalently locked E dimers exposing an epitope that has been shown to induce potent and broadly neutralizing antibodies.

Photocyclic behavior of rhodopsin induced by an atypical isomerization mechanism

Vertebrate rhodopsin (Rh) contains 11-cis-retinal as a chromophore to convert light energy into visual signals. On absorption of light, 11-cis-retinal is isomerized to all-trans-retinal, constituting a one-way reaction that activates transducin (G_t) followed by chromophore release. Here we report that bovine Rh, regenerated instead with a six-carbon-ring retinal chromophore featuring a C¹¹=C¹² double bond locked in its cis conformation (Rh6mr), employs an atypical isomerization mechanism by converting 11-cis to an 11,13-dicis configuration for prolonged G_t activation. Time-dependent UV-vis spectroscopy, HPLC, and molecular mechanics analyses revealed an atypical thermal reisomerization of the 11,13-dicis to the 11-cis configuration on a slow timescale, which enables Rh6mr to function in a photocyclic manner similar to that of microbial Rhs. With this photocyclic behavior, Rh6mr repeatedly recruits and activates G_t in response to light stimuli, making it an excellent candidate for optogenetic tools based on retinal analog-bound vertebrate Rhs. Overall, these comprehensive structure-function studies unveil a unique photocyclic mechanism of Rh activation by an 11-cis-to-11,13-dicis isomerization.

Molecular dynamics-based refinement and validation for sub-5 Å cryo-electron microscopy maps

Two structure determination methods, based on the molecular dynamics flexible fitting (MDFF) paradigm, are presented that resolve sub-5 Å cryo-electron microscopy (EM) maps with either single structures or ensembles of such structures. The methods, denoted cascade MDFF and resolution exchange MDFF, sequentially re-refine a search model against a series of maps of progressively higher resolutions, which ends with the original experimental resolution. Application of sequential re-refinement enables MDFF to achieve a radius of convergence of ~25 Å demonstrated with the accurate modeling of β-galactosidase and TRPV1 proteins at 3.2 Å and 3.4 Å resolution, respectively. The MDFF refinements uniquely offer map-model validation and B-factor determination criteria based on the inherent dynamics of the macromolecules studied, captured by means of local root mean square fluctuations. The MDFF tools described are available to researchers through an easy-to-use and cost-effective cloud computing resource on Amazon Web Services.

DOI:http://dx.doi.org/10.7554/eLife.16105.001

To understand the roles that proteins and other large molecules play inside cells, it is important to determine their structures. One of the techniques that researchers can use to do this is called cryo-electron microscopy (cryo-EM), which rapidly freezes molecules to fix them in position before imaging them in fine detail.

The cryo-EM images are like maps that show the approximate position of atoms. These images must then be processed in order to build a three-dimensional model of the protein that shows how its atoms are arranged relative to each other. One computational approach called Molecular Dynamics Flexible Fitting (MDFF) works by flexibly fitting possible atomic structures into cryo-EM maps. Although this approach works well with relatively undetailed (or ‘low resolution’) cryo-EM images, it struggles to handle the high-resolution cryo-EM maps now being generated.

Singharoy, Teo, McGreevy et al. have now developed two MDFF methods – called cascade MDFF and resolution exchange MDFF – that help to resolve atomic models of biological molecules from cryo-EM images. Each method can refine poorly guessed models into ones that are consistent with the high-resolution experimental images. The refinement is achieved by interpreting a range of images that starts with a ‘fuzzy’ image. The contrast of the image is then progressively improved until an image is produced that has a resolution that is good enough to almost distinguish individual atoms.

The method works because each cryo-EM image shows not just one, but a collection of atomic structures that the molecule can take on, with the fuzzier parts of the image representing the more flexible parts of the molecule. By taking into account this flexibility, the large-scale features of the protein structure can be determined first from the fuzzier images, and increasing the contrast of the images allows smaller-scale refinements to be made to the structure.

The MDFF tools have been designed to be easy to use and are available to researchers at low cost through cloud computing platforms. They can now be used to unravel the structure of many different proteins and protein complexes including those involved in photosynthesis, respiration and protein synthesis.

DOI:http://dx.doi.org/10.7554/eLife.16105.002

Refinement of atomic models in high resolution EM reconstructions using Flex-EM and local assessment

As the resolutions of Three Dimensional Electron Microscopic reconstructions of biological macromolecules are being improved, there is a need for better fitting and refinement methods at high resolutions and robust approaches for model assessment. Flex-EM/MODELLER has been used for flexible fitting of atomic models in intermediate-to-low resolution density maps of different biological systems. Here, we demonstrate the suitability of the method to successfully refine structures at higher resolutions (2.5-4.5Å) using both simulated and experimental data, including a newly processed map of Apo-GroEL. A hierarchical refinement protocol was adopted where the rigid body definitions are relaxed and atom displacement steps are reduced progressively at successive stages of refinement. For the assessment of local fit, we used the SMOC (segment-based Manders' overlap coefficient) score, while the model quality was checked using the Qmean score. Comparison of SMOC profiles at different stages of refinement helped in detecting regions that are poorly fitted. We also show how initial model errors can have significant impact on the goodness-of-fit. Finally, we discuss the implementation of Flex-EM in the CCP-EM software suite.

Rational design of crystal contact-free space in protein crystals for analyzing spatial distribution of motions within protein molecules

Contacts with neighboring molecules in protein crystals inevitably restrict the internal motions of intrinsically flexible proteins. The resultant clear electron densities permit model building, as crystallographic snapshot structures. Although these still images are informative, they could provide biased pictures of the protein motions. If the mobile parts are located at a site lacking direct contacts in rationally designed crystals, then the amplitude of the movements can be experimentally analyzed. We propose a fusion protein method, to create crystal contact-free space (CCFS) in protein crystals and to place the mobile parts in the CCFS. Conventional model building fails when large amplitude motions exist. In this study, the mobile parts appear as smeared electron densities in the CCFS, by suitable processing of the X-ray diffraction data. We applied the CCFS method to a highly mobile presequence peptide bound to the mitochondrial import receptor, Tom20, and a catalytically relevant flexible segment in the oligosaccharyltransferase, AglB. These two examples demonstrated the general applicability of the CCFS method to the analysis of the spatial distribution of motions within protein molecules.

Insights into Cullin-RING E3 ubiquitin ligase recruitment: structure of the VHL-EloBC-Cul2 complex

The von Hippel-Lindau tumor suppressor protein (VHL) recruits a Cullin 2 (Cul2) E3 ubiquitin ligase to downregulate HIF-1α, an essential transcription factor for the hypoxia response. Mutations in VHL lead to VHL disease and renal cell carcinomas. Inhibition of this pathway to upregulate erythropoietin production is a promising new therapy to treat ischemia and chronic anemia. Here, we report the crystal structure of VHL bound to a Cul2 N-terminal domain, Elongin B, and Elongin C (EloC). Cul2 interacts with both the VHL BC box and cullin box and a novel EloC site. Comparison with other cullin E3 ligase structures shows that there is a conserved, yet flexible, cullin recognition module and that cullin selectivity is influenced by distinct electrostatic interactions. Our structure provides a structural basis for the study of the pathogenesis of VHL disease and rationale for the design of novel compounds that may modulate cullin-substrate receptor interactions.

Structural insight into HIV-1 restriction by MxB

The myxovirus resistance (Mx) proteins are interferon-induced dynamin GTPases that can inhibit a variety of viruses. Recently, MxB, but not MxA, was shown to restrict HIV-1 by an unknown mechanism that likely occurs in close proximity to the host cell nucleus and involves the viral capsid. Here, we present the crystal structure of MxB and reveal determinants involved in HIV-1 restriction. MxB adopts an extended antiparallel dimer and dimerization, but not higher-ordered oligomerization, is critical for restriction. Although MxB is structurally similar to MxA, the orientation of individual domains differs between MxA and MxB, and their antiviral functions rely on separate determinants, indicating distinct mechanisms for virus inhibition. Additionally, MxB directly binds the HIV-1 capsid, and this interaction depends on dimerization and the N terminus of MxB as well as the assembled capsid lattice. These insights establish a framework for understanding the mechanism by which MxB restricts HIV-1.

E pluribus unum, no more: from one crystal, many conformations

Several distinct computational approaches have recently been implemented to represent conformational heterogeneity from X-ray crystallography datasets that are averaged in time and space. As these modeling methods mature, newly discovered alternative conformations are being used to derive functional protein mechanisms. Room temperature X-ray data collection is emerging as a key variable for sampling functionally relevant conformations also observed in solution studies. Although concerns about radiation damage are warranted with higher temperature data collection, 'diffract and destroy' strategies on X-ray free electron lasers may permit radiation damage-free data collection. X-ray crystallography need not be confined to 'static unique snapshots'; these experimental and computational advances are revealing how the many conformations populated within a single crystal are used in biological mechanisms.

Structural basis of HIV-1 Vpu-mediated BST2 antagonism via hijacking of the clathrin adaptor protein complex 1

BST2/tetherin, an antiviral restriction factor, inhibits the release of enveloped viruses from the cell surface. Human immunodeficiency virus-1 (HIV-1) antagonizes BST2 through viral protein u (Vpu), which downregulates BST2 from the cell surface. We report the crystal structure of a protein complex containing Vpu and BST2 cytoplasmic domains and the core of the clathrin adaptor protein complex 1 (AP1). This, together with our biochemical and functional validations, reveals how Vpu hijacks the AP1-dependent membrane trafficking pathways to mistraffick BST2. Vpu mimics a canonical acidic dileucine-sorting motif to bind AP1 in the cytosol, while simultaneously interacting with BST2 in the membrane. These interactions enable Vpu to build on an intrinsic interaction between BST2 and AP1, presumably causing the observed retention of BST2 in juxtanuclear endosomes and stimulating its degradation in lysosomes. The ability of Vpu to hijack AP-dependent trafficking pathways suggests a potential common theme for Vpu-mediated downregulation of host proteins.

DOI:http://dx.doi.org/10.7554/eLife.02362.001

HIV is a retrovirus that attacks the immune system, making the body increasingly susceptible to opportunistic infections and disease and eventually leading to AIDS. While antiretroviral drugs have allowed people with AIDS to live longer, there is no cure or vaccine for HIV.

Two types of HIV exist, with HIV-1 being much more common and pathogenic than HIV-2. Like other ‘complex’ retroviruses, the HIV-1 genome contains genes that encode various proteins that allow the virus to disrupt the immune response of the host it is attacking. Viral protein u is a protein encoded by HIV-1 (but not HIV-2) that counteracts an antiviral protein called BST2 in the host. BST2, which is part of the host's innate immune response, prevents newly formed viruses from leaving the surface of infected cells. By counteracting BST2, viral protein u allows the virus to spread in the host more efficiently.

Like many proteins, newly produced BST2 is packaged inside structures called vesicles in a part of the cell called the trans-Golgi network, and then sent to its destination. Complexes formed by various proteins make sure that the vesicles take their cargo to their correct destinations within the cell. Two adaptor protein complexes—known as AP1 and AP2—are thought to be involved the transport of BST2. However, it is not known how viral protein u stops BST2 from reaching the cell surface, or how it decreases the amount of BST2 in the cell as a whole. Jia et al. show how viral protein u and BST2 jointly interact with AP1. This interaction leads to the mistrafficking and degradation of BST2 and the counteraction of its antiviral activity.

DOI:http://dx.doi.org/10.7554/eLife.02362.002