Dissociation of an antiviral compound from the internal pocket of human rhinovirus 14 capsid

Five of Five VHHs Neutralizing Poliovirus Bind the Receptor-Binding Site

Nanobodies, or VHHs, that recognize poliovirus type 1 have previously been selected and characterized as candidates for antiviral agents or reagents for standardization of vaccine quality control. In this study, we present high-resolution cryo-electron microscopy reconstructions of poliovirus with five neutralizing VHHs. All VHHs bind the capsid in the canyon at sites that extensively overlap the poliovirus receptor-binding site. In contrast, the interaction involves a unique (and surprisingly extensive) surface for each of the five VHHs. Five regions of the capsid were found to participate in binding with all five VHHs. Four of these five regions are known to alter during the expansion of the capsid associated with viral entry. Interestingly, binding of one of the VHHs, PVSS21E, resulted in significant changes of the capsid structure and thus seems to trap the virus in an early stage of expansion.

IMPORTANCE

We describe the cryo-electron microscopy structures of complexes of five neutralizing VHHs with the Mahoney strain of type 1 poliovirus at resolutions ranging from 3.8 to 6.3Å. All five VHHs bind deep in the virus canyon at similar sites that overlap extensively with the binding site for the receptor (CD155). The binding surfaces on the VHHs are surprisingly extensive, but despite the use of similar binding surfaces on the virus, the binding surface on the VHHs is unique for each VHH. In four of the five complexes, the virus remains essentially unchanged, but for the fifth there are significant changes reminiscent of but smaller in magnitude than the changes associated with cell entry, suggesting that this VHH traps the virus in a previously undescribed early intermediate state. The neutralizing mechanisms of the VHHs and their potential use as quality control agents for the end game of poliovirus eradication are discussed.

Recent Developments in Molecular Simulation Approaches to Study Spherical Virus Capsids

Viruses are a particularly challenging systems to study via molecular simulation methods. Virus capsids typically consist of over 100 subunit proteins and reach dimensions of over 100 nm; solvated viruses capsid systems can be over 1 million atoms in size. In this review, I will present recent developments which have attempted to overcome the significant computational expense to perform simulations which can inform experimental studies, make useful predictions about biological phenomena and calculate material properties relevant to nanotechnology design efforts.

Long-distance correlations of rhinovirus capsid dynamics contribute to uncoating and antiviral activity

Human rhinovirus (HRV) and other members of the enterovirus genus bind small-molecule antiviral compounds in a cavity buried within the viral capsid protein VP1. These compounds block the release of the viral protein VP4 and RNA from inside the capsid during the uncoating process. In addition, the antiviral compounds prevent "breathing" motions, the transient externalization of the N-terminal regions of VP1 and VP4 from the inside of intact viral capsid. The site for externalization of VP1/VP4 or release of RNA is likely between protomers, distant to the binding cavity for antiviral compounds. Molecular dynamics simulations were conducted to explore how the antiviral compound, WIN 52084, alters properties of the HRV 14 capsid through long-distance effect. We developed an approach to analyze capsid dynamics in terms of correlated radial motion and the shortest paths of correlated motions. In the absence of WIN, correlated radial motion is observed between residues separated by as much as 85 Å, a remarkably long distance. The most frequently populated path segments of the network were localized near the fivefold symmetry axis and included those connecting the N termini of VP1 and VP4 with other regions, in particular near twofold symmetry axes, of the capsid. The results provide evidence that the virus capsid exhibits concerted long-range dynamics, which have not been previously recognized. Moreover, the presence of WIN destroys this radial correlation network, suggesting that the underlying motions contribute to a mechanistic basis for the initial steps of VP1 and VP4 externalization and uncoating.

Microscopic Symmetry Imposed by Rotational Symmetry Boundary Conditions in Molecular Dynamics Simulation

A large number of viral capsids, as well as other macromolecular assemblies, have icosahedral structure or structures with other rotational symmetries. This symmetry can be exploited during molecular dynamics (MD) to model in effect the full viral capsid using only a subset of primary atoms plus copies of image atoms generated from rotational symmetry boundary conditions (RSBC). A pure rotational symmetry operation results in both primary and image atoms at short range, and within nonbonded interaction distance of each other, so that nonbonded interactions can not be specified by the minimum image convention and explicit treatment of image atoms is required. As such an unavoidable consequence of RSBC is that the enumeration of nonbonded interactions in regions surrounding certain rotational axes must include both a primary atom and its copied image atom, thereby imposing microscopic symmetry for some forces. We examined the possibility of artifacts arising from this imposed microscopic symmetry of RSBC using two simulation systems: a water shell and human rhinovirus 14 (HRV14) capsid with explicit water. The primary unit was a pentamer of the icosahedron, which has the advantage of direct comparison of icosahedrally equivalent spatial regions, for example regions near a 2-fold symmetry axis with imposed symmetry and a 2-fold axis without imposed symmetry. Analysis of structural and dynamic properties of water molecules and protein atoms found similar behavior near symmetry axes with imposed symmetry and where the minimum image convention fails compared with that in other regions in the simulation system, even though an excluded volume effect was detected for water molecules near the axes with imposed symmetry. These results validate the use of RSBC for icosahedral viral capsids or other rotationally symmetric systems.

The Q cycle of cytochrome bc complexes: a structure perspective

Aspects of the crystal structures of the hetero-oligomeric cytochrome bc(1) and b(6)f ("bc") complexes relevant to their electron/proton transfer function and the associated redox reactions of the lipophilic quinones are discussed. Differences between the b(6)f and bc(1) complexes are emphasized. The cytochrome bc(1) and b(6)f dimeric complexes diverge in structure from a core of subunits that coordinate redox groups consisting of two bis-histidine coordinated hemes, a heme b(n) and b(p) on the electrochemically negative (n) and positive (p) sides of the complex, the high potential [2Fe-2S] cluster and c-type heme at the p-side aqueous interface and aqueous phase, respectively, and quinone/quinol binding sites on the n- and p-sides of the complex. The bc(1) and b(6)f complexes diverge in subunit composition and structure away from this core. b(6)f Also contains additional prosthetic groups including a c-type heme c(n) on the n-side, and a chlorophyll a and β-carotene. Common structure aspects; functions of the symmetric dimer. (I) Quinone exchange with the bilayer. An inter-monomer protein-free cavity of approximately 30Å along the membrane normal×25Å (central inter-monomer distance)×15Å (depth in the center), is common to both bc(1) and b(6)f complexes, providing a niche in which the lipophilic quinone/quinol (Q/QH(2)) can be exchanged with the membrane bilayer. (II) Electron transfer. The dimeric structure and the proximity of the two hemes b(p) on the electrochemically positive side of the complex in the two monomer units allow the possibility of two alternate routes of electron transfer across the complex from heme b(p) to b(n): intra-monomer and inter-monomer involving electron cross-over between the two hemes b(p). A structure-based summary of inter-heme distances in seven bc complexes, representing mitochondrial, chromatophore, cyanobacterial, and algal sources, indicates that, based on the distance parameter, the intra-monomer pathway would be favored kinetically. (III) Separation of quinone binding sites. A consequence of the dimer structure and the position of the Q/QH(2) binding sites is that the p-side QH(2) oxidation and n-side Q reduction sites are each well separated. Therefore, in the event of an overlap in residence time by QH(2) or Q molecules at the two oxidation or reduction sites, their spatial separation would result in minimal steric interference between extended Q or QH(2) isoprenoid chains. (IV) Trans-membrane QH(2)/Q transfer. (i) n/p-side QH(2)/Q transfer may be hindered by lipid acyl chains; (ii) the shorter less hindered inter-monomer pathway across the complex would not pass through the center of the cavity, as inferred from the n-side antimycin site on one monomer and the p-side stigmatellin site on the other residing on the same surface of the complex. (V) Narrow p-side portal for QH(2)/Q passage. The [2Fe-2S] cluster that serves as oxidant, and whose histidine ligand serves as a H(+) acceptor in the oxidation of QH(2), is connected to the inter-monomer cavity by a narrow extended portal, which is also occupied in the b(6)f complex by the 20 carbon phytyl chain of the bound chlorophyll.

A Comparison of Three Perturbation Molecular Dynamics Methods for Modeling Conformational Transitions

Targeted, steered, and biased molecular dynamics (MD) are widely used methods for studying transition processes of biomolecules. They share the common feature of adding external perturbations along a conformational progress variable to guide the transition in a predefined direction in conformational space, yet differ in how these perturbations are applied. In the present paper, we report a comparison of these three methods on generating transition paths for two different processes: the unfolding of the B domain of protein A and a conformational transition of the catalytic domain of a Src kinase Lyn. Transition pathways were calculated with different simulation parameters including the choice of progress variable and the simulation length or biasing force constant. A comparison of the generated paths based on structural similarity finds that the three perturbation MD methods generate similar transition paths for a given progress variable in most cases. On the other hand, the path depends more strongly on the choice of progress variable used to move the system between the initial and final states. Potentials of mean force (PMF) were calculated starting from unfolding trajectories to estimate the relative probabilities of the paths. A lower PMF was found for the lowest biasing force constant with BMD.

Antiviral compounds discovered by virtual screening of small-molecule libraries against dengue virus E protein

Infection by the mosquito-borne dengue virus causes dengue fever and the sometimes fatal dengue hemorrhagic fever. The increasing number of dengue infections per year suggests that the virus is becoming more virulent and its transmission is expanding. Nevertheless, no effective treatment for dengue infection currently exists. In a search for antiviral agents effective against dengue virus, we investigated the potential of targeting a structural protein site rather than an enzymatic one. Using this approach, we now report the discovery of a small molecule ligand that inhibits viral growth. Our results also provide the first evidence that the binding site, a pocket located at the hinge between domains 1 and 2 of the envelope protein (E protein) on the virus surface, is a valid target for antiviral therapy. Ligand candidates were identified from libraries of approximately 142,000 compounds using a computational high-throughput screening protocol targeting this pocket of the E protein. Cell-based assays were conducted on 23 top-ranked compounds. Among four with good antiviral activity profiles, the compound P02 was found to inhibit viral reproduction at micromolar concentrations. Using saturation transfer difference NMR spectroscopy, we also show that the compound binds virus and competes for binding E protein with the known ligand N-octyl-beta-D-glucoside. Together, the results are consistent with an inhibition mechanism against maturation or host-cell entry mediated by ligand binding to the E-protein pocket. P02 is a promising lead compound for future development of an effective treatment against dengue virus and related flaviviruses.

Free-energy reconstruction from experiments performed under different biasing programs

Recently developed nonequilibrium statistical physics relationships, including Jarzynski's equality and the Crooks fluctuation theorem, have been used to calculate equilibrium thermodynamic properties using data from both laboratory and computational experiments. Although Jarzynski's derivation does not include an explicit time dependency, prior work utilizing the relationship to reconstruct free-energy surfaces has combined data from experiments performed under identical conditions. Here, a formalism is developed for combining data from a variety of biasing protocols, as in dynamic force spectroscopy experiments. The method is then demonstrated on data from simulations conducted under a wide range of pulling velocities and with a random biasing protocol.

The folding mechanism of collagen-like model peptides explored through detailed molecular simulations

Collagen has a unique folding mechanism that begins with the formation of a triple-helical structure near its C terminus followed by propagation of this structure to the N terminus. To elucidate factors that affect the folding of collagen, we explored the folding pathway of collagen-like model peptides using detailed molecular simulations with explicit solvent. Using biased molecular dynamics we examined the latter stages of folding of a peptide model of native collagen, (Pro-Hyp-Gly)10, and a peptide that models a Gly --> Ser mutation found in several forms of osteogenesis imperfecta, (Pro-Hyp-Gly)3-Pro-Hyp-Ser-(Pro-Hyp-Gly)6. Starting from an unfolded state that contains a C-terminal nucleated trimer, (Pro-Hyp-Gly)10 folds to a structure where two of the three chains associate through water-mediated hydrogen bonds and the third is relatively separated from this dimer. Calculated free-energy profiles for folding from this intermediate to the final triple-helical structure suggest that further folding occurs at a rate of approximately one Pro-Hyp-Gly triplet per msec. In contrast, after 6 nsec of biased dynamics, the region N-terminal to the Ser residue in (Pro-Hyp-Gly)3-Pro-Hyp-Ser-(Pro-Hyp-Gly)6 folds to a structure where the three chains form close contacts near the N terminus, away from the mutation site. Further folding to an ideal triple-helical structure at the site of the mutation is unfavorable as the free energy of a triple-helical conformation at this position is more than 20 kcal/mol higher than that of a structure with unassociated chains. These data provide insights into the folding pathway of native collagen and the events underlying the formation of misfolded structures.

VP4 protein from human rhinovirus 14 is released by pressure and locked in the capsid by the antiviral compound WIN

Rhinoviruses are the major causative agents of the common cold in humans. Here, we studied the stability of human rhinovirus type 14 (HRV14) under conditions of high hydrostatic pressure, low temperature, and urea in the absence and presence of an antiviral drug. Capsid dissociation and changes in the protein conformation were monitored by fluorescence spectroscopy, light scattering, circular dichroism, gel filtration chromatography, mass spectrometry and infectivity assays. The data show that high pressure induces the dissociation of HRV14 and that this process is inhibited by WIN 52084. MALDI-TOF mass spectrometry experiments demonstrate that VP4, the most internal viral protein, is released from the capsid by pressure treatment. This release of VP4 is concomitant with loss of infectivity. Our studies also show that at least one antiviral effect of the WIN drugs involves the locking of VP4 inside the capsid by blocking the dynamics associated with cell attachment.

Structure of a membrane-binding domain from a non-enveloped animal virus: insights into the mechanism of membrane permeability and cellular entry

The gamma(1)-peptide is a 21-residue lipid-binding domain from the non-enveloped Flock House virus (FHV). Unlike enveloped viruses, the entry of non-enveloped viruses into cells is believed to occur without membrane fusion. In this study, we performed NMR experiments to establish the solution structure of a membrane-binding peptide from a small non-enveloped icosahedral virus. The three-dimensional structure of the FHV gamma(1)-domain was determined at pH 6.5 and 4.0 in a hydrophobic environment. The secondary and tertiary structures were evaluated in the context of the capacity of the peptide for permeabilizing membrane vesicles of different lipid composition, as measured by fluorescence assays. At both pH values, the peptide has a kinked structure, similar to the fusion domain from the enveloped viruses. The secondary structure was similar in three different hydrophobic environments as follows: water/trifluoroethanol, SDS, and membrane vesicles of different compositions. The ability of the peptide to induce vesicle leakage was highly dependent on the membrane composition. Although the gamma-peptide shares some structural properties to fusion domains of enveloped viruses, it did not induce membrane fusion. Our results suggest that small protein components such as the gamma-peptide in nodaviruses (such as FHV) and VP4 in picornaviruses have a crucial role in conducting nucleic acids through cellular membranes and that their structures resemble the fusion domains of membrane proteins from enveloped viruses.

Src kinase activation: A switched electrostatic network

Src tyrosine kinases are essential in numerous cell signaling pathways, and improper functioning of these enzymes has been implicated in many diseases. The activity of Src kinases is regulated by conformational activation, which involves several structural changes within the catalytic domain (CD): the orientation of two lobes of CD; rearrangement of the activation loop (A-loop); and movement of an alpha-helix (alphaC), which is located at the interface between the two lobes, into or away from the catalytic cleft. Conformational activation was investigated using biased molecular dynamics to explore the transition pathway between the active and the down-regulated conformation of CD for the Src-kinase family member Lyn kinase, and to gain insight into the interdependence of these changes. Lobe opening is observed to be a facile motion, whereas movement of the A-loop motion is more complex requiring secondary structure changes as well as communication with alphaC. A key result is that the conformational transition involves a switch in an electrostatic network of six polar residues between the active and the down-regulated conformations. The exchange between interactions links the three main motions of the CD. Kinetic experiments that would demonstrate the contribution of the switched electrostatic network to the enzyme mechanism are proposed. Possible implications for regulation conferred by interdomain interactions are also discussed.