Use of MCSS to design small targeted libraries: application to picornavirus ligands.

Computational methods were used to design structure-based combinatorial libraries of antipicornaviral capsid-binding ligands. The multiple copy simultaneous search (MCSS) program was employed to calculate functionality maps for many diverse functional groups for both the poliovirus and rhinovirus capsid structures in the region of the known drug binding pocket. Based on the results of the MCSS calculations, small combinatorial libraries consisting of 10s or 100s of three-monomer compounds were designed and synthesized. Ligand binding was demonstrated by a noncell-based mass spectrometric assay, a functional immuno-precipitation assay, and crystallographic analysis of the complexes of the virus with two of the candidate ligands. The P1/Mahoney poliovirus strain was used in the experimental studies. A comparison showed that the MCSS calculations had correctly identified the observed binding site for all three monomer units in one ligand and for two out of three in the other ligand. The correct central monomer position in the second ligand was reproduced in calculations in which the several key residues lining the pocket were allowed to move. This study validates the computational methodology. It also illustrates that subtle changes in protein structure can lead to differences in docking results and points to the importance of including target flexibility, as well as ligand flexibility, in the design process.

Ab initio phasing of high-symmetry macromolecular complexes: successful phasing of authentic poliovirus data to 3.0 A resolution.

A genetic algorithm-based computational method for the ab initio phasing of diffraction data from crystals of symmetric macromolecular structures, such as icosahedral viruses, has been implemented and applied to authentic data from the P1/Mahoney strain of poliovirus. Using only single-wavelength native diffraction data, the method is shown to be able to generate correct phases, and thus electron density, to 3.0 A resolution. Beginning with no advance knowledge of the shape of the virus and only approximate knowledge of its size, the method uses a genetic algorithm to determine coarse, low-resolution (here, 20.5 A) models of the virus that obey the known non-crystallographic symmetry (NCS) constraints. The best scoring of these models are subjected to refinement and NCS-averaging, with subsequent phase extension to high resolution (3.0 A). Initial difficulties in phase extension were overcome by measuring and including all low-resolution terms in the transform. With the low-resolution data included, the method was successful in generating essentially correct phases and electron density to 6.0 A in every one of ten trials from different models identified by the genetic algorithm. Retrospective analysis revealed that these correct high-resolution solutions converged from a range of significantly different low-resolution phase sets (average differences of 59.7 degrees below 24 A). This method represents an efficient way to determine phases for icosahedral viruses, and has the advantage of producing phases free from model bias. It is expected that the method can be extended to other protein systems with high NCS.

The structure and oligomerization of the yeast arginine methyltransferase, Hmt1.

Protein methylation at arginines is ubiquitous in eukaryotes and affects signal transduction, gene expression and protein sorting. Hmt1/Rmt1, the major arginine methyltransferase in yeast, catalyzes methylation of arginine residues in several mRNA-binding proteins and facilitates their export from the nucleus. We now report the crystal structure of Hmt1 at 2.9 A resolution. Hmt1 forms a hexamer with approximate 32 symmetry. The surface of the oligomer is dominated by large acidic cavities at the dimer interfaces. Mutation of dimer contact sites eliminates activity of Hmt1 both in vivo and in vitro. Mutating residues in the acidic cavity significantly reduces binding and methylation of the substrate Npl3.

The crystal structure of an unusual processivity factor, herpes simplex virus UL42, bound to the C terminus of its cognate polymerase.

Herpes simplex virus DNA polymerase is a heterodimer composed of a catalytic subunit, Pol, and an unusual processivity subunit, UL42, which, unlike processivity factors such as PCNA, directly binds DNA. The crystal structure of a complex of the C-terminal 36 residues of Pol bound to residues 1-319 of UL42 reveals remarkable similarities between UL42 and PCNA despite contrasting biochemical properties and lack of sequence homology. Moreover, the Pol-UL42 interaction resembles the interaction between the cell cycle regulator p21 and PCNA. The structure and previous data suggest that the UL42 monomer interacts with DNA quite differently than does multimeric toroidal PCNA. The details of the structure lead to a model for the mechanism of UL42, provide the basis for drug design, and allow modeling of other proteins that lack sequence homology with UL42 or PCNA.

Three-dimensional structure of poliovirus receptor bound to poliovirus.

Poliovirus initiates infection by binding to its cellular receptor (Pvr). We have studied this interaction by using cryoelectron microscopy to determine the structure, at 21-A resolution, of poliovirus complexed with a soluble form of its receptor (sPvr). This density map aided construction of a homology-based model of sPvr and, in conjunction with the known crystal structure of the virus, allowed delineation of the binding site. The virion does not change significantly in structure on binding sPvr in short incubations at 4 degrees C. We infer that the binding configuration visualized represents the initial interaction that is followed by structural changes in the virion as infection proceeds. sPvr is segmented into three well-defined Ig-like domains. The two domains closest to the virion (domains 1 and 2) are aligned and rigidly connected, whereas domain 3 diverges at an angle of approximately 60 degrees. Two nodules of density on domain 2 are identified as glycosylation sites. Domain 1 penetrates the "canyon" that surrounds the 5-fold protrusion on the capsid surface, and its binding site involves all three major capsid proteins. The inferred pattern of virus-sPvr interactions accounts for most mutations that affect the binding of Pvr to poliovirus.

Molecular tectonic model of virus structural transitions: the putative cell entry states of poliovirus.

Upon interacting with its receptor, poliovirus undergoes conformational changes that are implicated in cell entry, including the externalization of the viral protein VP4 and the N terminus of VP1. We have determined the structures of native virions and of two putative cell entry intermediates, the 135S and 80S particles, at approximately 22-A resolution by cryo-electron microscopy. The 135S and 80S particles are both approximately 4% larger than the virion. Pseudoatomic models were constructed by adjusting the beta-barrel domains of the three capsid proteins VP1, VP2, and VP3 from their known positions in the virion to fit the 135S and 80S reconstructions. Domain movements of up to 9 A were detected, analogous to the shifting of tectonic plates. These movements create gaps between adjacent subunits. The gaps at the sites where VP1, VP2, and VP3 subunits meet are plausible candidates for the emergence of VP4 and the N terminus of VP1. The implications of these observations are discussed for models in which the externalized components form a transmembrane pore through which viral RNA enters the infected cell.

Crystal structure of a brominated RNA helix with four mismatched base pairs: An investigation into RNA conformational variability.

The X-ray crystal structure of a brominated RNA helix with four mismatched base pairs and sequence r(UG(Br)C(Br)CAGUUCGCUGGC)(2) was determined to 2.1 A using the methods of multiwavelength anomalous diffraction (MAD) applied to the bromine K-absorption edge. There are three molecules in the asymmetric unit with unique crystal-packing environments, revealing true conformational variability at high resolution for this sequence. The structure shows that the sequence itself does not define a consistent pattern of solvent molecules, with the exception of the mismatched base pairs, implying that specific RNA-protein interactions would occur only with the nucleotides. There are a number of significant tertiary interactions, some of which are a result of the brominated base pairs and others that are directly mediated by the RNA 2' hydroxyl groups. The mismatched base pairs exhibit a solvent network as well as a stacking pattern with their nearest neighbors that validate previous thermodynamic analysis.

Crystal structure of a 14 bp RNA duplex with non-symmetrical tandem GxU wobble base pairs.

Adjacent GxU wobble base pairs are frequently found in rRNA. Atomic structures of small RNA motifs help to provide a better understanding of the effects of various tandem mismatches on duplex structure and stability, thereby providing better rules for RNA structure prediction and validation. The crystal structure of an RNA duplex containing the sequence r(GGUAUUGC-GGUACC)2 has been solved at 2.1 A resolution using experimental phases. Novel refinement strategies were needed for building the correct solvent model. At present, this is the only short RNA duplex structure containing 5'-U-U-3'/3'-G-G-5' non-symmetric tandem GxU wobble base pairs. In the 14mer duplex, the six central base pairs are all displaced away from the helix axis, yielding significant changes in local backbone conformation, helix parameters and charge distribution that may provide specific recognition sites for biologically relevant ligand binding. The greatest deviations from A-form helix occur where the guanine of a wobble base pair stacks over a purine from the opposite strand. In this vicinity, the intra-strand phosphate distances increase significantly, and the major groove width increases up to 3 A. Structural comparisons with other short duplexes containing symmetrical tandem GxU or GxT wobble base pairs show that nearest-neighbor sequence dependencies govern helical twist and the occurrence of cross-strand purine stacks.

The crystal structure of Dps, a ferritin homolog that binds and protects DNA.

The crystal structure of Dps, a DNA-binding protein from starved E. coli that protects DNA from oxidative damage, has been solved at 1.6 A resolution. The Dps monomer has essentially the same fold as ferritin, which forms a 24-mer with 432 symmetry, a hollow core and pores at the three-fold axes. Dps forms a dodecamer with 23 (tetrahedral) point group symmetry which also has a hollow core and pores at the three-folds. The structure suggests a novel DNA-binding motif and a mechanism for DNA protection based on the sequestration of Fe ions.

Structure determination of echovirus 1.

The atomic structure of echovirus 1 (a member of the enterovirus genus of the picornavirus family) has been determined using cryo-crystallography and refined to 3.55 A resolution. Echovirus 1 crystallizes in space group P22121 with a = 352.45, b = 472.15 and c = 483.20 A. The crystals contain one full virus particle in the asymmetric unit allowing for 60-fold noncrystallographic symmetry averaging. The diffraction pattern shows strong pseudo-B-centering with reflections with h + l = 2n + 1 being systematically weak or absent below about 6 A resolution. The size of the unit cell and presence of pseudo-B-centering placed strong constraints on the allowed packing of the icosahedral particle in the crystal lattice. These constraints greatly facilitated the determination of the orientation and position of the virus by reducing the dimensionality of the search, but interactions between the crystallographic and noncrystallographic symmetries rendered the choice of space group ambiguous until very late in the structure determination. This structure determination provides a striking example of the power of packing analysis in molecular replacement and illustrates how subtle interactions between crystallographic and noncrystallographic symmetries can be resolved.

Structural studies of poliovirus mutants that overcome receptor defects.

In order to better understand the process of cell entry for non-enveloped viruses, we have solved the crystal structures of five poliovirus mutants which can infect cells expressing mutant poliovirus receptors. Four of these structures have been solved from frozen crystals using cryocrystallographic data collection methods. The mutations have a range of structural consequences, from small local perturbations to significant loop rearrangements. All of the mutant viruses are more labile to conversion to an apparent cell entry intermediate, suggesting that these mutant viruses could compensate for the suboptimal receptors by lowering the thermal energy required to undergo the receptor-mediated conformational change.

Ligand-induced conformational changes in poliovirus-antiviral drug complexes.

Crystal structures of the Mahoney strain of type 1 poliovirus complexed with the antiviral compounds R80633 and R77975 were determined at 2.9 A resolution. These compounds block infection by preventing conformational changes required for viral uncoating. In various drug-poliovirus complexes reported earlier, no significant conformational changes were found in the structures of the capsid proteins. In the structures reported here, the strain of virus is relatively insensitive to these antivirals. Correspondingly, significant conformational changes are necessary to accommodate the drug. These conformational changes affect both the immediate vicinity of the drug binding site, and more distant loops located near the fivefold axis. In addition, small but concerted shifts of the centers of mass of the major capsid proteins consistently have been detected whose magnitudes are correlated inversely with the effectiveness of the drugs. Collectively, the drug complexes appear to sample the conformational repertoire of poliovirus near equilibrium, and thus provide a possible model for the earliest stages of viral uncoating during infection.

A pseudo-cell based approach to efficient crystallographic refinement of viruses.

Strategies have been developed for the inexpensive refinement of atomic models of viruses and of other highly symmetric structures. These methods, which have been used in the refinement of several strains of poliovirus, focus on an arbitrary-sized parallelepiped (termed the 'protomer' box) containing a single complete averaged copy of the structural motif which forms the protein capsid, together with the fragments of other symmetry-related copies of the motif which are located in its immediate neighborhood. The Fourier transform of the protomer box provides reference structure factors for stereochemically restrained crystallographic refinement of the atomic model parameters. The phases of the reference structure factors are based on the averaged map, and are not permitted to change during the refinement. It is demonstrated that models refined using the protomer box methods do not differ significantly from models refined by more expensive full-cell calculations.

A genetic algorithm for the ab initio phasing of icosahedral viruses.

Genetic algorithms have been investigated as computational tools for the de novo phasing of low-resolution X-ray diffraction data from crystals of icosahedral viruses. Without advance knowledge of the shape of the virus and only approximate knowledge of its size, the virus can be modeled as the symmetry expansion of a short list of nearly tetrahedrally arranged lattice points which coarsely, but uniformly, sample the icosahedrally unique volume. The number of lattice points depends on an estimate of the non-redundant information content at the working resolution limit. This parameterization permits a simple matrix formulation of the model evaluation calculation, resulting in a highly efficient survey of the space of possible models. Initially, one bit per parameter is sufficient, since the assignment of ones and zeros to the lattice points yields a physically reasonable low-resolution image of the virus. The best candidate solutions identified by the survey are refined to relax the constraints imposed by the coarseness of the modeling, and then trials whose intensity-based statistics are comparatively good in all resolution ranges are chosen. This yields an acceptable starting point for symmetry-based direct phase extension about half the time. Improving efficiency by incorporating the selection criterion directly into the genetic algorithm's fitness function is discussed.

An enzyme-substrate complex involved in bacterial cell wall biosynthesis.

The crystal structure of UDP-N-acetylenolpyruvylglucosamine reductase in the presence of its substrate, enolpyruvyl-UDP-N-acetylglucosamine, has been solved to 2.7 A resolution. This enzyme is responsible for the synthesis of UDP-N-acetylmuramic acid in bacterial cell wall biosynthesis and consequently provides an attractive target for the design of antibacterial agents. The structure reveals a novel flavin binding motif, shows a striking alignment of the flavin with the substrate, and suggests a catalytic mechanism for the reduction of this unusual enol ether.

Binding of the antiviral drug WIN51711 to the sabin strain of type 3 poliovirus: structural comparison with drug binding in rhinovirus 14.

The crystal structure of the Sabin strain of type 3 poliovirus (P3/Sabin) complexed with the antiviral drug WIN51711 has been determined at 2.9 A resolution. Drugs of this kind are known to inhibit the uncoating of the virus during infection, by stabilizing the capsid against receptor-induced conformational changes. The electron density for the bound drug is very well defined so that its position and orientation are unambiguous. The drug binds in a nearly extended conformation, slightly bent in the middle, in a blind pocket formed predominantly by hydrophobic residues in the core of the beta-barrel of capsid protein VP1. Comparisons between this structure, the corresponding drug complex in human rhinovirus 14 (HRV 14), and the native structures of both viruses demonstrate that the binding of WIN51711 has markedly different effects on the structures of these two viruses. Unlike HRV14, wherein large conformational changes are observed in the coat protein after drug binding, the binding of this drug in poliovirus does not induce any significant conformational changes in the structure of the capsid protein, though the drug has a greater inhibitory effect in P3/Sabin than in HRV14. The implications of this result for the mechanism of capsid stabilization are discussed.

Structure of the complex between the Fab fragment of a neutralizing antibody for type 1 poliovirus and its viral epitope.

The crystal structure of the complex between the Fab fragment of C3, a neutralizing antibody for poliovirus, and a peptide corresponding to the viral epitope has been determined at 3.0 A resolution. Although this antibody was originally raised to heat inactivated (noninfectious) virus particles, it strongly neutralizes the Mahoney strain of type 1 poliovirus. Eleven peptide residues are well-defined in the electron-density map and form two type I beta-turns in series. At the carboxyl end, the peptide is bound snugly in the antibody-combining site and adopts a conformation that differs significantly from the structure of the corresponding residues in the virus. Structural comparisons between the peptide in the complex and the viral epitope suggests that on binding to infectious virions, this antibody may induce structural changes important for neutralization.

Role and mechanism of the maturation cleavage of VP0 in poliovirus assembly: structure of the empty capsid assembly intermediate at 2.9 A resolution.

The crystal structure of the P1/Mahoney poliovirus empty capsid has been determined at 2.9 A resolution. The empty capsids differ from mature virions in that they lack the viral RNA and have yet to undergo a stabilizing maturation cleavage of VP0 to yield the mature capsid proteins VP4 and VP2. The outer surface and the bulk of the protein shell are very similar to those of the mature virion. The major differences between the 2 structures are focused in a network formed by the N-terminal extensions of the capsid proteins on the inner surface of the shell. In the empty capsids, the entire N-terminal extension of VP1, as well as portions corresponding to VP4 and the N-terminal extension of VP2, are disordered, and many stabilizing interactions that are present in the mature virion are missing. In the empty capsid, the VP0 scissile bond is located some 20 A away from the positions in the mature virion of the termini generated by VP0 cleavage. The scissile bond is located on the rim of a trefoil-shaped depression in the inner surface of the shell that is highly reminiscent of an RNA binding site in bean pod mottle virus. The structure suggests plausible (and ultimately testable) models for the initiation of encapsidation, for the RNA-dependent autocatalytic cleavage of VP0, and for the role of the cleavage in establishing the ordered N-terminal network and in generating stable virions.

Structures of poliovirus complexes with anti-viral drugs: implications for viral stability and drug design.

BACKGROUND: Picornaviruses, such as the structurally related polioviruses and rhinoviruses, are important human pathogens which have been the target of major drug development efforts. Receptor-mediated uncoating and thermal inactivation of poliovirus and rhinovirus are inhibited by agents that bind to each virus by inserting into a pocket in the beta barrel of the viral capsid protein, VP1. This pocket, which is normally empty in human rhinovirus-14 (HRV14), is occupied by an unknown natural ligand in poliovirus. Structural studies of HRV14-drug complexes have shown that drug binding causes large, localized changes in the conformation of VP1. RESULTS: We report the crystal structures of six complexes between poliovirus and capsid-binding, antiviral drugs, including complexes of four different drugs with the Sabin vaccine strain of type 3 poliovirus, and complexes of one of these drugs with two other poliovirus strains that contain sequence differences in the drug-binding site. In each complex, the changes in capsid structure associated with drug binding are limited to minor adjustments in the conformations of a few side chains lining the binding site. CONCLUSIONS: The minor structural changes caused by drug binding suggest a model of drug action in which it is the conformational changes prevented by the bound drug, rather than obvious conformational changes induced by drug binding, which exert the biological effect. Our results, along with additional structures of rhinovirus-drug complexes, suggest possible improvements in drug design, and provide important clues about the nature of the conformational changes that are involved in the uncoating process.

Three-dimensional structure of Theiler virus.

Theiler murine encephalomyelitis virus strains are categorized into two groups, a neurovirulent group that rapidly kills the host, and a demyelinating group that causes a generally nonlethal infection of motor neurons followed by a persistent infection of the white matter with demyelinating lesions similar to those found in multiple sclerosis. The three-dimensional structure of the DA strain, a member of the demyelinating group, has been determined at 2.8 A resolution. As in other picornaviruses, the icosahedral capsid is formed by the packing of wedge-shaped eight-stranded antiparallel beta barrels. The surface of Theiler virus has large star-shaped plateaus at the fivefold axes and broad depressions spanning the twofold axes. Several unusual structural features are clustered near one edge of the depression. These include two finger-like loops projecting from the surface (one formed by residues 78-85 of VP1, and the other formed by residues 56-65 of VP3) and a third loop containing three cysteines (residues 87, 89, and 91 of VP3), which appear to be covalently modified. Most of the sequence differences between the demyelinating and neurovirulent groups that could play a role in determining pathogenesis map to the surface of the star-shaped plateau. The distribution of these sequence differences on the surface of the virion is consistent with models in which the differences in the pathogenesis of the two groups of Theiler viruses are the result of differences in immunological or receptor-mediated recognition processes.