Engineering, characterization and phage display of hepatitis C virus NS3 protease and NS4A cofactor peptide as a single-chain protein

The polyprotein encoded by hepatitis C virus (HCV) genomic RNA is processed into functional polypeptides by both host- and virus-encoded proteases. The HCV-encoded NS3 protease and its cofactor peptide NS4A form a non-covalent complex, which participates in processing the viral polyprotein. This proteolytic activity is believed to be essential for virus proliferation and thus the NS3 protease is a prime target for developing anti-HCV pharmacological agents. Recent X-ray crystallography structural studies have revealed the nature of this non- covalent complex between NS3 protease and the 'active' central segment of NS4A, providing the opportunity to design a single-chain polypeptide. To this end, the DNA sequence encoding for the NS4A peptide (residues 21-34) was genetically fused via a short linker, capable of making a beta-turn, to the N-terminus of the NS3 protease domain. This engineered single-chain NS3-protease (scNS3) is fully active with kinetic parameters virtually identical with those of the NS3/ NS4A non-covalent complex. Moreover, the scNS3 protease can be displayed on filamentous phage and affinity selected using an immobilized specific inhibitor. The scNS3 expressed as a soluble protein and in a phage-display format facilitates enzyme engineering for further structural studies and in vitro selection of potential drug- resistant mutants. These are important steps towards developing effective anti-protease compounds.      

Affinity selection of a camelized V(H) domain antibody inhibitor of hepatitis C virus NS3 protease

The HCV genome encodes, within the NS3 gene, a serine protease whose activity specifically cleaves the viral polyprotein precursor. Proteolytic processing of HCV polyprotein precursor by the viral NS3 proteinase is essential for virion maturation and designing specific inhibitors of this protease as possible anti-viral agents is a desirable and practical objective. With a view to studying both the function of HCV NS3 protease and to designing inhibitors of this enzyme, we directed our interest towards engineering macromolecular inhibitors of the viral protease catalytic activity. We describe here the affinity-selection and biochemical characterization of one inhibitor, cV(H)E2, a 'camelized' variable domain antibody fragment, isolated from a phage displayed synthetic repertoire, which is a potent and selective inhibitor of proteolysis by the NS3 enzyme. In addition to being useful as a biological probe to study the function of HCV protease, this inhibitor can serve as a potential pharmacophore model to design antivirals. Moreover, the results suggest a way of engineering improved human-derived small recognition units tailored for enzyme inhibition.      

Structural Modeling and Biochemical Characterization of Recombinant KPN_02809, a Zinc-Dependent Metalloprotease from Klebsiella pneumoniae MGH 78578

Klebsiella pneumoniae is a Gram-negative, cylindrical rod shaped opportunistic pathogen that is found in the environment as well as existing as a normal flora in mammalian mucosal surfaces such as the mouth, skin, and intestines. Clinically it is the most important member of the family of Enterobacteriaceae that causes neonatal sepsis and nosocomial infections. In this work, a combination of protein sequence analysis, structural modeling and molecular docking simulation approaches were employed to provide an understanding of the possible functions and characteristics of a hypothetical protein (KPN_02809) from K. pneumoniae MGH 78578. The computational analyses showed that this protein was a metalloprotease with zinc binding motif, HEXXH. To verify this result, a ypfJ gene which encodes for this hypothetical protein was cloned from K. pneumoniae MGH 78578 and the protein was overexpressed in Escherichia coli BL21 (DE3). The purified protein was about 32 kDa and showed maximum protease activity at 30 °C and pH 8.0. The enzyme activity was inhibited by metalloprotease inhibitors such as EDTA, 1,10-phenanthroline and reducing agent, 1,4-dithiothreitol (DTT). Each molecule of KPN_02809 protein was also shown to bind one zinc ion. Hence, for the first time, we experimentally confirmed that KPN_02809 is an active enzyme with zinc metalloprotease activity.     

Influenza A Virus NS1 Protein Structural Flexibility Analysis According to Its Structural Polymorphism Using Computational Approaches

Influenza A viruses are highly contagious RNA viruses that cause respiratory tract infections in humans and animals. Their non-structural protein NS1, a homodimer of two 230-residue chains, is the main viral factor in counteracting the antiviral defenses of the host cell. Its RNA-binding domain is an obligate dimer that is connected to each of the two effector domains by a highly flexible unstructured linker region of ten amino acids. The flexibility of NS1 is a key property that allows its effector domains and its RNA binding domain to interact with several protein partners or RNAs. The three-dimensional structures of full-length NS1 dimers revealed that the effector domains could adopt three distinct conformations as regards their mutual interactions and their orientation relative to the RNA binding domain (closed, semi-open and open). The origin of this structural polymorphism is currently being investigated and several hypotheses are proposed, among which one posits that it is a strain-specific property. In the present study, we explored through computational molecular modeling the dynamic and flexibility properties of NS1 from three important influenza virus A strains belonging to three distinct subtypes (H1N1, H6N6, H5N1), for which at least one conformation is available in the Protein Data Bank. In order to verify whether NS1 is stable in three forms for the three strains, we constructed homology models if the corresponding forms were not available in the Protein Data Bank. Molecular dynamics simulations were performed in order to predict the stability over time of the three distinct sequence variants of NS1, in each of their three distinct conformations. Our results favor the co-existence of three stable structural forms, regardless of the strain, but also suggest that the length of the linker, along with the presence of specific amino acids, modulate the dynamic properties and the flexibility of NS1.     

Transient structural ordering of the RNA-binding domain of carnation mottle virus p7 movement protein modulates nucleic acid binding

Plant viral movement proteins bind to RNA and participate in the intra- and intercellular movement of the RNAs from plant viruses. However, the role and magnitude of the conformational changes associated with the formation of RNA-protein complexes are not yet defined. Here we describe studies on the relevance of a preexisting nascent alpha-helix at the C terminus of the RNA-binding domain of p7, a movement protein from carnation mottle virus, to RNA binding. Synthetic peptide analogues and single amino acid mutation at the RNA-binding domain of recombinant p7 protein were used to correlate the transient structural order in aqueous solution with RNA-binding potential.     

Molecular Dynamics Studies on the Structural Stability Prediction of SARS-CoV-2 Variants Including Multiple Mutants

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has caused the Coronavirus Disease (COVID-19) pandemic worldwide. The spike protein in SARS-CoV-2 fuses with and invades cells in the host respiratory system by binding to angiotensin-converting enzyme 2 (ACE2). The spike protein, however, undergoes continuous mutation from a D614G single mutant to an omicron variant, including multiple mutants. In this study, variants, including multiple mutants (double, triple mutants, B.1.620, delta, alpha, delta_E484Q, mu, and omicron) were investigated in patients. The 3D structure of the full-length spike protein was used in conformational analysis depending on the SARS-CoV-2 variants. The structural stability of the variant types was analyzed based on the distance between the receptor-binding domain (RBD) of each chain in the spike protein and the binding free energy between the spike protein and bound ACE2 in the one-, two-, and three-open-complex forms using molecular dynamics (MD) simulation. Omicron variants, the most prevalent in the recent history of the global pandemic, which consist of 32 mutations, showed higher stability in all open-complex forms compared with that of the wild type and other variants. We suggest that the conformational stability of the spike protein is the one of the important determinants for the differences in viral infectivity among variants, including multiple mutants.     

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.     

Probing the Structure and Function of the Cytosolic Domain of the Human Zinc Transporter ZnT8 with Nickel(II) Ions

The human zinc transporter ZnT8 provides the granules of pancreatic β-cells with zinc (II) ions for assembly of insulin hexamers for storage. Until recently, the structure and function of human ZnTs have been modelled on the basis of the 3D structures of bacterial zinc exporters, which form homodimers with each monomer having six transmembrane α-helices harbouring the zinc transport site and a cytosolic domain with an α,β structure and additional zinc-binding sites. However, there are important differences in function as the bacterial proteins export an excess of zinc ions from the bacterial cytoplasm, whereas ZnT8 exports zinc ions into subcellular vesicles when there is no apparent excess of cytosolic zinc ions. Indeed, recent structural investigations of human ZnT8 show differences in metal binding in the cytosolic domain when compared to the bacterial proteins. Two common variants, one with tryptophan (W) and the other with arginine (R) at position 325, have generated considerable interest as the R-variant is associated with a higher risk of developing type 2 diabetes. Since the mutation is at the apex of the cytosolic domain facing towards the cytosol, it is not clear how it can affect zinc transport through the transmembrane domain. We expressed the cytosolic domain of both variants of human ZnT8 and have begun structural and functional studies. We found that (i) the metal binding of the human protein is different from that of the bacterial proteins, (ii) the human protein has a C-terminal extension with three cysteine residues that bind a zinc(II) ion, and (iii) there are small differences in stability between the two variants. In this investigation, we employed nickel(II) ions as a probe for the spectroscopically silent Zn(II) ions and utilised colorimetric and fluorimetric indicators for Ni(II) ions to investigate metal binding. We established Ni(II) coordination to the C-terminal cysteines and found differences in metal affinity and coordination in the two ZnT8 variants. These structural differences are thought to be critical for the functional differences regarding the diabetes risk. Further insight into the assembly of the metal centres in the cytosolic domain was gained from potentiometric investigations of zinc binding to synthetic peptides corresponding to N-terminal and C-terminal sequences of ZnT8 bearing the metal-coordinating ligands. Our work suggests the involvement of the C-terminal cysteines, which are part of the cytosolic domain, in a metal chelation and/or acquisition mechanism and, as now supported by the high-resolution structural work, provides the first example of metal-thiolate coordination chemistry in zinc transporters.     

Analysis of the Structural Mechanism of ATP Inhibition at the AAA1 Subunit of Cytoplasmic Dynein-1 Using a Chemical “Toolkit”

Dynein is a ~1.2 MDa cytoskeletal motor protein that carries organelles via retrograde transport in eukaryotic cells. The motor protein belongs to the ATPase family of proteins associated with diverse cellular activities and plays a critical role in transporting cargoes to the minus end of the microtubules. The motor domain of dynein possesses a hexameric head, where ATP hydrolysis occurs. The presented work analyzes the structure–activity relationship (SAR) of dynapyrazole A and B, as well as ciliobrevin A and D, in their various protonated states and their 46 analogues for their binding in the AAA1 subunit, the leading ATP hydrolytic site of the motor domain. This study exploits in silico methods to look at the analogues’ effects on the functionally essential subsites of the motor domain of dynein 1, since no similar experimental structural data are available. Ciliobrevin and its analogues bind to the ATP motifs of the AAA1, namely, the walker-A (W-A) or P-loop, the walker-B (W-B), and the sensor I and II. Ciliobrevin A shows a better binding affinity than its D analogue. Although the double bond in ciliobrevin A and D was expected to decrease the ligand potency, they show a better affinity to the AAA1 binding site than dynapyrazole A and B, lacking the bond. In addition, protonation of the nitrogen atom in ciliobrevin A and D, as well as dynapyrazole A and B, at the N9 site of ciliobrevin and the N7 of the latter increased their binding affinity. Exploring ciliobrevin A geometrical configuration suggests the E isomer has a superior binding profile over the Z due to binding at the critical ATP motifs. Utilizing the refined structure of the motor domain obtained through protein conformational search in this study exhibits that Arg1852 of the yeast cytoplasmic dynein could involve in the “glutamate switch” mechanism in cytoplasmic dynein 1 in lieu of the conserved Asn in AAA+ protein family.     

F?rster Resonance Energy Transfer (FRET) as a Tool for Dissecting the Molecular Mechanisms for Maturation of the Shigella Type III Secretion Needle Tip Complex

Förster resonance energy transfer (FRET) provides a powerful tool for monitoring intermolecular interactions and a sensitive technique for studying Å-level protein conformational changes. One system that has particularly benefited from the sensitivity and diversity of FRET measurements is the maturation of the Shigella type III secretion apparatus (T3SA) needle tip complex. The Shigella T3SA delivers effector proteins into intestinal cells to promote bacterial invasion and spread. The T3SA is comprised of a basal body that spans the bacterial envelope and a needle with an exposed tip complex that matures in response to environmental stimuli. FRET measurements demonstrated bile salt binding by the nascent needle tip protein IpaD and also mapped resulting structural changes which led to the recruitment of the translocator IpaB. At the needle tip IpaB acts as a sensor for host cell contact but prior to secretion, it is stored as a heterodimeric complex with the chaperone IpgC. FRET analyses showed that chaperone binding to IpaB’s N-terminal domain causes a conformational change in the latter. These FRET analyses, with other biophysical methods, have been central to understanding T3SA maturation and will be highlighted, focusing on the details of the FRET measurements and the relevance to this particular system.     

Synthesis and in vitro Evaluation of West Nile Virus Protease Inhibitors Based on the 2‐{6‐[2‐(5‐Phenyl‐4H‐{1,2,4]triazol‐3‐ylsulfanyl)acetylamino]benzothiazol‐2‐ylsulfanyl}acetamide Scaffold

In recent years, clinical symptoms resulting from West Nile virus (WNV) infection have worsened in severity, with an increased frequency in neuroinvasive diseases among the elderly. As there are presently no successful therapies against WNV for use in humans, continual efforts to develop new chemotherapeutics against this virus are highly desired. The viral NS2B‐NS3 protease is a promising target for viral inhibition due to its importance in viral replication and its unique substrate preference. In this study, a WNV NS2B‐NS3 protease inhibitor with a 2‐{6‐[2‐(5‐phenyl‐4H‐[1,2,4]triazol‐3‐ylsulfanyl)acetylamino]benzothiazol‐2‐ylsulfanyl}acetamide scaffold was identified during screening. Optimization of this initial hit by synthesis and screening of a focused compound library with this scaffold led to the identification of a novel uncompetitive inhibitor ( 1 a24 , IC₅₀=3.4±0.2 μM ) of the WNV NS2B‐NS3 protease. Molecular docking of 1 a24 into the WNV protease showed that the compound interferes with productive interactions of the NS2B cofactor with the NS3 protease and is an allosteric inhibitor of the WNV NS3 protease.     

An Engineered Calmodulin‐Based Allosteric Switch for Peptide Biosensing

This work describes the development of a new platform for allosteric protein engineering that takes advantage of the ability of calmodulin to change conformation upon binding to peptide and protein ligands. The switch we have developed consists of a fusion protein in which calmodulin is genetically inserted into the sequence of TEM1 β‐lactamase. In this approach, calmodulin acts as the input domain, whose ligand‐dependent conformational changes control the activity of the β‐lactamase output domain. The new allosteric enzyme exhibits up to 120 times higher catalytic activity in the activated (peptide bound) state compared to the inactive (no peptide bound) state in vitro. Activation of the enzyme is ligand‐dependent—peptides with higher affinities for wild‐type calmodulin exhibit increased switch activity. Calmodulin's ability to “turn on” the activity of β‐lactamase makes this a potentially valuable scaffold for the directed evolution of highly specific biosensors for detecting toxins and other clinically relevant biomarkers.     

Imidazolidinones and Imidazolidine-2,4-diones as Antiviral Agents

Imidazolidinones and imidazolidine-2,4-diones are important classes of heterocyclic compounds that possess potent activities against several viruses such as dengue virus, enterovirus, hepatitis C virus (HCV), and human immunodeficiency virus (HIV). The first imidazolidinone derivative as an anti-HIV agent was reported in 1996. Imidazolidinones inhibit HIV aspartic protease activity, and also act as CCR5 co-receptor antagonists. Significant effort has been devoted to the design of various imidazolidinone analogues that are active against drug-resistant HIV strains, with fewer side effects. Different scaffolds have been designed through both rational drug design strategies and computer-aided drug design. Imidazolidinones have been found to be potent against HIV, and preclinical studies are currently in progress. There are some reports of imidazolidinones as having both anti-HCV and anti-dengue virus activity, and more research has yet to be done along these lines. These compounds inhibit NS3 serine protease of HCV, and NS2B-NS3 protease of dengue virus. Pyridyl-imidazolidinones possess very specific and potent activity against human enterovirus 71 (EV71) by targeting the EV71 capsid protein VP1, and inhibiting viral adsorption and/or viral RNA uncoating.     

NMR chemical shift perturbation study of the N-terminal domain of Hsp90 upon binding of ADP,AMP-PNP,geldanamycin, and radicicol

Hsp90 is one of the most abundant chaperone proteins in the cytosol. In an ATP-dependent manner it plays an essential role in the folding and activation of a range of client proteins involved in signal transduction and cell cycle regulation. We used NMR shift perturbation experiments to obtain information on the structural implications of the binding of AMP-PNP (adenylyl-imidodiphosphate-a non-hydrolysable ATP analogue), ADP and the inhibitors radicicol and geldanamycin. Analysis of (1)H,(15)N correlation spectra showed a specific pattern of chemical shift perturbations at N210 (ATP binding domain of Hsp90, residues 1-210) upon ligand binding. This can be interpreted qualitatively either as a consequence of direct ligand interactions or of ligand-induced conformational changes within the protein. All ligands show specific interactions in the binding site, which is known from the crystal structure of the N-terminal domain of Hsp90. For AMP-PNP and ADP, additional shift perturbations of residues outside the binding pocket were observed and can be regarded as a result of conformational rearrangement upon binding. According to the crystal structures, these regions are the first alpha-helix and the "ATP-lid" ranging from amino acids 85 to 110. The N-terminal domain is therefore not a passive nucleotide-binding site, as suggested by X-ray crystallography, but responds to the binding of ATP in a dynamic way with specific structural changes required for the progression of the ATPase cycle.     

Production of an Active,Human Membrane Protein in Saccharomyces cerevisiae: Full-Length FICD

The human Fic domain-containing protein (FICD) is a type II endoplasmic reticulum (ER) membrane protein that is important for the maintenance of ER proteostasis. Structural and in vitro biochemical characterisation of FICD AMPylase and deAMPylase activity have been restricted to the soluble ER-luminal domain produced in Escherichia coli. Information about potentially important features, such as structural motifs, modulator binding sites or other regulatory elements, is therefore missing for the approximately 100 N-terminal residues including the transmembrane region of FICD. Expressing and purifying the required quantity and quality of membrane proteins is demanding because of the low yields and poor stability often observed. Here, we produce full-length FICD by combining a Saccharomyces cerevisiae-based platform with green fluorescent protein (GFP) tagging to optimise the conditions for expression, solubilisation and purification. We subsequently employ these conditions to purify milligram quantities of His-tagged FICD per litre of culture, and show that the purified, detergent-solubilised membrane protein is an active deAMPylating enzyme. Our work provides a straightforward methodology for producing not only full-length FICD, but also other membrane proteins in S. cerevisiae for structural and biochemical characterisation.     

NMR backbone assignment of a protein kinase catalytic domain by a combination of several approaches: application to the catalytic subunit of cAMP-dependent protein kinase

Protein phosphorylation is one of the most important mechanisms used for intracellular regulation in eukaryotic cells. Currently, one of the best-characterized protein kinases is the catalytic subunit of cAMP-dependent protein kinase or protein kinase A (PKA). PKA has the typical bilobular structure of kinases, with the active site consisting of a cleft between the two structural lobes. For full kinase activity, the catalytic subunit has to be phosphorylated. The catalytic subunit of PKA has two main phosphorylation sites: Thr197 and Ser338. Binding of ATP or inhibitors to the ATP site induces large structural changes. Here we describe the partial backbone assignment of the PKA catalytic domain by NMR spectroscopy, which represents the first NMR assignment of any protein kinase catalytic domain. Backbone resonance assignment for the 42 kDa protein was accomplished by an approach employing 1) triply ((2)H,(13)C,(15)N) labeled protein and classical NMR assignment experiments, 2) back-calculation of chemical shifts from known X-ray structures, 3) use of paramagnetic adenosine derivatives as spin-labels, and 4) selective amino acid labeling. Interpretation of chemical-shift perturbations allowed mapping of the interaction surface with the protein kinase inhibitor H7. Furthermore, structural conformational changes were observed by comparison of backbone amide shifts obtained by 2D (1)H,(15)N TROSY of an inactive Thr197Ala mutant with the wild-type enzyme.     

Integrated Self-Assembly of the Mms6 Magnetosome Protein to Form an Iron-Responsive Structure

A common feature of biomineralization proteins is their self-assembly to produce a surface consistent in size with the inorganic crystals that they produce. Mms6, a small protein of 60 amino acids from Magnetospirillum magneticum strain AMB-1 that promotes the in vitro growth of superparamagnetic magnetite nanocrystals, assembles in aqueous solution to form spherical micelles that could be visualized by TEM and AFM. The results reported here are consistent with the view that the N and C-terminal domains interact with each other within one polypeptide chain and across protein units in the assembly. From studies to determine the amino acid residues important for self-assembly, we identified the unique GL repeat in the N-terminal domain with additional contributions from amino acids in other positions, throughout the molecule. Analysis by CD spectroscopy identified a structural change in the iron-binding C-terminal domain in the presence of Fe³⁺. A change in the intrinsic fluorescence of tryptophan in the N-terminal domain showed that this structural change is transmitted through the protein. Thus, self-assembly of Mms6 involves an interlaced structure of intra- and inter-molecular interactions that results in a coordinated structural change in the protein assembly with iron binding.