Activation of the West Nile virus NS3 protease: molecular dynamics evidence for a conformational selection mechanism

Quantum Biochemistry and MM-PBSA Description of the ZIKV NS2B-NS3 Protease: Insights into the Binding Interactions beyond the Catalytic Triad Pocket

The Zika virus protease NS2B-NS3 has a binding site formed with the participation of a H51-D75-S135 triad presenting two forms, active and inactive. Studies suggest that the inactive conformation is a good target for the design of inhibitors. In this paper, we evaluated the co-crystallized structures of the protease with the inhibitors benzoic acid (5YOD) and benzimidazole-1-ylmethanol (5H4I). We applied a protocol consisting of two steps: first, classical molecular mechanics energy minimization followed by classical molecular dynamics were performed, obtaining stabilized molecular geometries; second, the optimized/relaxed geometries were used in quantum biochemistry and molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) calculations to estimate the ligand interactions with each amino acid residue of the binding pocket. We show that the quantum-level results identified essential residues for the stabilization of the 5YOD and 5H4I complexes after classical energy minimization, matching previously published experimental data. The same success, however, was not observed for the MM-PBSA simulations. The application of quantum biochemistry methods seems to be more promising for the design of novel inhibitors acting on NS2B-NS3.

Ligand Accessibility Insights to the Dengue Virus NS3-NS2B Protease Assessed by Long-Timescale Molecular Dynamics Simulations

Dengue is a tropical disease caused by the dengue virus (DENV), with an estimate of 300 million new cases every year. Due to the limited vaccine efficiency and absence of effective antiviral treatment, new drug candidates are urgently needed. DENV NS3-NS2B protease complex is essential for viral post-translational processing and maturation, and this enzyme has been extensively studied as a relevant drug target. Crystal structures often underestimate NS3-NS2B flexibility, whereas they can adopt different conformational states depending on the bound substrate. We conducted molecular dynamics simulations (∼30 μs) with a non- and covalently bound inhibitor to understand the conformational changes in the DENV-3 NS3-NS2B complex. Our results show that the open-closing movement of the protease exposes multiple druggable subpockets that can be investigated in later drug discovery efforts.

MD simulations reveal alternate conformations of the oxyanion hole in the Zika virus NS2B/NS3 protease

Recent crystallography studies have shown that the binding site oxyanion hole plays an important role in inhibitor binding, but can exist in two conformations (active/inactive). We have undertaken molecular dynamics (MD) calculations to better understand oxyanion hole dynamics and thermodynamics. We find that the Zika virus (ZIKV) NS2B/NS3 protease maintains a stable closed conformation over multiple 100-ns conventional MD simulations in both the presence and absence of inhibitors. The S1, S2, and S3 pockets are stable as well. However, in two of eight simulations, the A132-G133 peptide bond in the binding pocket of S1' spontaneously flips to form a 3₁₀ -helix that corresponds to the inactive conformation of the oxyanion hole, and then maintains this conformation until the end of the 100-ns conventional MD simulations without inversion of the flip. This conformational change affects the S1' pocket in ZIKV NS2B/NS3 protease active site, which is important for small molecule binding. The simulation results provide evidence at the atomic level that the inactive conformation of the oxyanion hole is more favored energetically when no specific interactions are formed between substrate/inhibitor and oxyanion hole residues. Interestingly, however, transition between the active and inactive conformation of the oxyanion hole can be observed by boosting the valley potential in accelerated MD simulations. This supports a proposed induced-fit mechanism of ZIKV NS2B/NS3 protease from computational methods and provides useful direction to enhance inhibitor binding predictions in structure-based drug design.

All-atom virus simulations

The constant threat of viral disease can be combated by the development of novel vaccines and therapeutics designed to disrupt key features of virus structure or infection cycle processes. Such development relies on high-resolution characterization of viruses and their dynamical behaviors, which are often challenging to obtain solely by experiment. In response, all-atom molecular dynamics simulations are widely leveraged to study the structural components of viruses, leading to some of the largest simulation endeavors undertaken to date. The present work reviews exemplary all-atom simulation work on viruses, as well as progress toward simulating entire virions.

An Overview of Current Approaches Toward the Treatment and Prevention of West Nile Virus Infection

The persistence of West Nile virus (WNV) infections throughout the USA since its inception in 1999 and its continuous spread throughout the globe calls for an urgent need of effective treatments and prevention measures. Although the licensing of several WNV vaccines for veterinary use provides a proof of concept, similar efforts on the development of an effective vaccine for humans remain still unsuccessful. Increased understanding of biology and pathogenesis of WNV together with recent technological advancements have raised hope that an effective WNV vaccine may be available in the near future. In addition, rapid progress in the structural and functional characterization of WNV and other flaviviral proteins have provided a solid base for the design and development of several classes of inhibitors as potential WNV therapeutics. Moreover, the therapeutic monoclonal antibodies demonstrate an excellent efficacy against WNV in animal models and represent a promising class of WNV therapeutics. However, there are some challenges as to the design and development of a safe and efficient WNV vaccine or therapeutic. In this chapter, we discuss the current approaches, progress, and challenges toward the development of WNV vaccines, therapeutic antibodies, and antiviral drugs.

A conformational switch high-throughput screening assay and allosteric inhibition of the flavivirus NS2B-NS3 protease

The flavivirus genome encodes a single polyprotein precursor requiring multiple cleavages by host and viral proteases in order to produce the individual proteins that constitute an infectious virion. Previous studies have revealed that the NS2B cofactor of the viral NS2B-NS3 heterocomplex protease displays a conformational dynamic between active and inactive states. Here, we developed a conformational switch assay based on split luciferase complementation (SLC) to monitor the conformational change of NS2B and to characterize candidate allosteric inhibitors. Binding of an active-site inhibitor to the protease resulted in a conformational change of NS2B and led to significant SLC enhancement. Mutagenesis of key residues at an allosteric site abolished this induced conformational change and SLC enhancement. We also performed a virtual screen of NCI library compounds to identify allosteric inhibitors, followed by in vitro biochemical screening of the resultant candidates. Only three of these compounds, NSC135618, 260594, and 146771, significantly inhibited the protease of Dengue virus 2 (DENV2) in vitro, with IC₅₀ values of 1.8 μM, 11.4 μM, and 4.8 μM, respectively. Among the three compounds, only NSC135618 significantly suppressed the SLC enhancement triggered by binding of active-site inhibitor in a dose-dependent manner, indicating that it inhibits the conformational change of NS2B. Results from virus titer reduction assays revealed that NSC135618 is a broad spectrum flavivirus protease inhibitor, and can significantly reduce titers of DENV2, Zika virus (ZIKV), West Nile virus (WNV), and Yellow fever virus (YFV) on A549 cells in vivo, with EC₅₀ values in low micromolar range. In contrast, the cytotoxicity of NSC135618 is only moderate with CC₅₀ of 48.8 μM on A549 cells. Moreover, NSC135618 inhibited ZIKV in human placental and neural progenitor cells relevant to ZIKV pathogenesis. Results from binding, kinetics, Western blot, mass spectrometry and mutagenesis experiments unambiguously demonstrated an allosteric mechanism for inhibition of the viral protease by NSC135618.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Enzymatic detoxication, conformational selection, and the role of molten globule active sites

The role of conformational ensembles in enzymatic reactions remains unclear. Discussion concerning "induced fit" versus "conformational selection" has, however, ignored detoxication enzymes, which exhibit catalytic promiscuity. These enzymes dominate drug metabolism and determine drug-drug interactions. The detoxication enzyme glutathione transferase A1-1 (GSTA1-1), exploits a molten globule-like active site to achieve remarkable catalytic promiscuity wherein the substrate-free conformational ensemble is broad with barrierless transitions between states. A quantitative index of catalytic promiscuity is used to compare engineered variants of GSTA1-1 and the catalytic promiscuity correlates strongly with characteristics of the thermodynamic partition function, for the substrate-free enzymes. Access to chemically disparate transition states is encoded by the substrate-free conformational ensemble. Pre-steady state catalytic data confirm an extension of the conformational selection model, wherein different substrates select different starting conformations. The kinetic liability of the conformational breadth is minimized by a smooth landscape. We propose that "local" molten globule behavior optimizes detoxication enzymes.

Structured crowding and its effects on enzyme catalysis

Macromolecular crowding decreases the diffusion rate, shifts the equilibrium of protein-protein and protein-substrate interactions, and changes protein conformational dynamics. Collectively, these effects contribute to enzyme catalysis. Here we describe how crowding may bias the conformational change and dynamics of enzyme populations and in this way affect catalysis. Crowding effects have been studied using artificial crowding agents and in vivo-like environments. These studies revealed a correlation between protein dynamics and function in the crowded environment. We suggest that crowded environments be classified into uniform crowding and structured crowding. Uniform crowding represents random crowding conditions created by synthetic particles with a narrow size distribution. Structured crowding refers to the highly coordinated cellular environment, where proteins and other macromolecules are clustered and organized. In structured crowded environments the perturbation of protein thermal stability may be lower; however, it may still be able to modulate functions effectively and dynamically. Dynamic, allosteric enzymes could be more sensitive to cellular perturbations if their free energy landscape is flatter around the native state; on the other hand, if their free energy landscape is rougher, with high kinetic barriers separating deep minima, they could be more robust. Above all, cells are structured; and this holds both for the cytosol and for the membrane environment. The crowded environment is organized, which limits the search, and the crowders are not necessarily inert. More likely, they too transmit allosteric effects, and as such play important functional roles. Overall, structured cellular crowding may lead to higher enzyme efficiency and specificity.

Peptide binding to the PDZ3 domain by conformational selection

The PDZ domains, a large family of peptide recognition proteins, bind to the C-terminal segment of membrane ion channels and receptors thereby mediating their localization. The peptide binding process is not known in detail and seems to differ among different PDZ domains. For the third PDZ domain of the synaptic protein PSD-95 (PDZ3), a lock-and-key mechanism was postulated on the basis of the almost perfect overlap of the crystal structures in the presence and absence of its peptide ligand. Here, peptide binding to PDZ3 is investigated by explicit solvent molecular dynamics (MD) simulations (for a total of 1.3 μs) and the cut-based free energy profile method for determining free energy barriers and basins. The free energy landscape of apo PDZ3 indicates that there are multiple basins within the native state. These basins differ by the relative orientation of the α2 helix and β2 strand, the two secondary structure elements that make up the peptide binding site. Only the structure with the smallest aperture of the binding site is populated in the MD simulations of the complex whose analysis reveals that the peptide ligand binds to PDZ3 by selecting one of three conformations. Thus, the dynamical information obtained by the atomistic simulations increment the static, that is, partial, picture of the PDZ3 binding mechanism based on the X-ray crystallography data. Importantly, the simulation results show for the first time that conformational selection is a possible mechanism of peptide binding by PDZ domains in general.

Virtual ligand screening of the National Cancer Institute (NCI) compound library leads to the allosteric inhibitory scaffolds of the West Nile Virus NS3 proteinase

Viruses of the genus Flavivirus are responsible for significant human disease and mortality. The N-terminal domain of the flaviviral nonstructural (NS)3 protein codes for the serine, chymotrypsin-fold proteinase (NS3pro). The presence of the nonstructural (NS)2B cofactor, which is encoded by the upstream gene in the flaviviral genome, is necessary for NS3pro to exhibit its proteolytic activity. The two-component NS2B-NS3pro functional activity is essential for the viral polyprotein processing and replication. Both the structure and the function of NS2B-NS3pro are conserved in the Flavivirus family. Because of its essential function in the posttranslational processing of the viral polyprotein precursor, NS2B-NS3pro is a promising target for anti-flavivirus drugs. To identify selective inhibitors with the reduced cross-reactivity and off-target effects, we focused our strategy on the allosteric inhibitors capable of targeting the NS2B-NS3pro interface rather than the NS3pro active site. Using virtual ligand screening of the diverse, ∼275,000-compound library and the catalytic domain of the two-component West Nile virus (WNV) NS2B-NS3pro as a receptor, we identified a limited subset of the novel inhibitory scaffolds. Several of the discovered compounds performed as allosteric inhibitors and exhibited a nanomolar range potency in the in vitro cleavage assays. The inhibitors were also potent in cell-based assays employing the sub-genomic, luciferase-tagged WNV and Dengue viral replicons. The selectivity of the inhibitors was confirmed using the in vitro cleavage assays with furin, a human serine proteinase, the substrate preferences of which are similar to those of WNV NS2B-NS3pro. Conceptually, the similar in silico drug discovery strategy may be readily employed for the identification of inhibitors of other flaviviruses.

Enzyme dynamics point to stepwise conformational selection in catalysis

Recent data increasingly reveal that conformational dynamics are indispensable to enzyme function throughout the catalytic cycle, in substrate recruiting, chemical transformation, and product release. Conformational transitions may involve conformational selection and induced fit, which can be viewed as a special case in the catalytic network. NMR, X-ray crystallography, single-molecule FRET, and simulations clearly demonstrate that the free enzyme dynamics already encompass all the conformations necessary for substrate binding, preorganization, transition-state stabilization, and product release. Conformational selection and substate population shift at each step of the catalytic turnover can accommodate enzyme specificity and efficiency. Within such a framework, entropy can have a larger role in conformational dynamics than in direct energy transfer in dynamically promoted catalysis.

Flaviviral protease inhibitors identified by fragment-based library docking into a structure generated by molecular dynamics

Fragment-based docking was used to select a conformation for virtual screening from a molecular dynamics trajectory of the West Nile virus nonstructural 3 protease. This conformation was chosen from an ensemble of 100 molecular dynamics snapshots because it optimally accommodates benzene, the most common ring in known drugs, and two positively charged fragments (methylguanidinium and 2-phenylimidazoline). The latter fragments were used as probes because of the large number of hydrogen bond acceptors in the substrate binding site of the protease. Upon high-throughput docking of a diversity set of 18,694 molecules and pose filtering, only five compounds were chosen for experimental validation, and two of them are active in the low micromolar range in an enzymatic assay and a tryptophan fluorescence quenching assay. Evidence for specific binding to the protease active site is provided by nuclear magnetic resonance spectroscopy. The two inhibitors have different scaffolds (diphenylurea and diphenyl ester) and are promising lead candidates because they have a molecular weight of about 300 Da.

NMR study of complexes between low molecular mass inhibitors and the West Nile virus NS2B-NS3 protease

The two-component NS2B-NS3 protease of West Nile virus is essential for its replication and presents an attractive target for drug development. Here, we describe protocols for the high-yield expression of stable isotope-labelled samples in vivo and in vitro. We also describe the use of NMR spectroscopy to determine the binding mode of new low molecular mass inhibitors of the West Nile virus NS2B-NS3 protease which were discovered using high-throughput in vitro screening. Binding to the substrate-binding sites S1 and S3 is confirmed by intermolecular NOEs and comparison with the binding mode of a previously identified low molecular mass inhibitor. Our results show that all these inhibitors act by occupying the substrate-binding site of the protease rather than by an allosteric mechanism. In addition, the NS2B polypeptide chain was found to be positioned near the substrate-binding site, as observed previously in crystal structures of the protease in complex with peptide inhibitors or bovine pancreatic trypsin inhibitor. This indicates that the new low molecular mass compounds, although inhibiting the protease, also promote the proteolytically active conformation of NS2B, which is very different from the crystal structure of the protein without inhibitor.