Molecular dynamics simulations on the first two helices of Vpu from HIV-1

Repurposing dye ligands as antivirals via a docking approach on viral membrane and globular proteins - SARS-CoV-2 and HPV-16

A series of dye ligands are docked to three different proteins, E and 3a of severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) and E6 of human papilloma virus type 16 (HPV-16) using three different software. A four-level selection algorithm is used based on nonparametric statistics of numerical key values such as the "rank" derived from (i) averaged estimated binding energies (EBEs) and (ii) absolute EBE value of each of the software, (iii) frequency of ranking and (iv) rank of the area-under-curve values (AUCs) from decoy docking. A series of repurposing drugs and known antivirals used in experimental studies are docked for comparison. One dye ligand is ranked best for all proteins using the selection algorithm levels i - iii. Another three dye ligands are ranked top for the proteins individually when using all four levels.

Rotational Dynamics of The Transmembrane Domains Play an Important Role in Peptide Dynamics of Viral Fusion and Ion Channel Forming Proteins—A Molecular Dynamics Simulation Study

Focusing on the transmembrane domains (TMDs) of viral fusion and channel-forming proteins (VCPs), experimentally available and newly generated peptides in an ideal conformation of the S and E proteins of severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) and SARS-CoV, gp41 and Vpu, both of human immunodeficiency virus type 1 (HIV-1), haemagglutinin and M2 of influenza A, as well as gB of herpes simplex virus (HSV), are embedded in a fully hydrated lipid bilayer and used in multi-nanosecond molecular dynamics simulations. It is aimed to identify differences in the dynamics of the individual TMDs of the two types of viral membrane proteins. The assumption is made that the dynamics of the individual TMDs are decoupled from their extra-membrane domains, and that the mechanics of the TMDs are distinct from each other due to the different mechanism of function of the two types of proteins. The diffusivity coefficient (DC) of the translational and rotational diffusion is decreased in the oligomeric state of the TMDs compared to those values when calculated from simulations in their monomeric state. When comparing the calculations for two different lengths of the TMD, a longer full peptide and a shorter purely TMD stretch, (i) the difference of the calculated DCs begins to level out when the difference exceeds approximately 15 amino acids per peptide chain, and (ii) the channel protein rotational DC is the most affected diffusion parameter. The rotational dynamics of the individual amino acids within the middle section of the TMDs of the fusion peptides remain high upon oligomerization, but decrease for the channel peptides, with an increasing number of monomers forming the oligomeric state, suggesting an entropic penalty on oligomerization for the latter.

Conserved residues within the HIV-1 Vpu transmembrane-proximal hinge region modulate BST2 binding and antagonism

Background

BST2 inhibits HIV-1 release by tethering nascent virions to the surface of infected cells. HIV-1 Vpu overcomes this restriction by removing BST2 from viral budding sites via BST2 intracellular trapping and sequestration, surface downregulation and/or displacement mechanisms. Vpu is composed of a short luminal tail, a transmembrane domain (TMD) and a cytoplasmic hinge region that is followed by two helices. BST2 counteraction relies on the ability of Vpu to physically bind BST2 through TMD interactions and recruit the clathrin-dependent trafficking machinery via a canonical acidic di-leucine signalling motif within the helix-2 of Vpu. The highly conserved Vpu transmembrane-proximal hinge region encompasses residues that resemble an acidic leucine-based trafficking motif, whose functional roles are currently ill-defined. In this study, we investigated the contribution of these residues towards Vpu-mediated BST2 antagonism.

Results

We show that while these conserved residues have no intrinsic activity on the cellular distribution of Vpu in the absence of BST2, they regulate the ability of Vpu to bind to BST2 and, consequently, govern both BST2-dependent trafficking properties of the protein as well as its co-localization with BST2. Moreover, these residues, particularly a glutamic acid residue positioned immediately following the TMD, are a determinant not only for efficient targeting of BST2, but also binding and degradation of CD4, another host membrane protein targeted by Vpu. Mechanistically, our data are consistent with a role of these residues in the maintenance of the Vpu TMD conformational configuration such that interactions with membrane-associated host targets are favoured.

Conclusions

Altogether, this work demonstrates an important regulatory role of the transmembrane-proximal Vpu hinge region residues towards enabling the protein to efficiently engage its target host proteins. Thus, this highly conserved, cytosolic Vpu hinge region may represent an attractive target for the development of anti-Vpu inhibitors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-017-0345-6) contains supplementary material, which is available to authorized users.

Membrane protein assembly: two cytoplasmic phosphorylated serine sites of Vpu from HIV-1 affect oligomerization

Viral protein U (Vpu) encoded by human immunodeficiency virus type 1 (HIV-1) is a short integral membrane protein which is known to self-assemble within the lipid membrane and associate with host factors during the HIV-1 infectivity cycle. In this study, full-length Vpu (M group) from clone NL4-3 was over-expressed in human cells and purified in an oligomeric state. Various single and double mutations were constructed on its phosphorylation sites to mimic different degrees of phosphorylation. Size exclusion chromatography of wild-type Vpu and mutants indicated that the smallest assembly unit of Vpu was a dimer and over time Vpu formed higher oligomers. The rate of oligomerization increased when (i) the degree of phosphorylation at serines 52 and 56 was decreased and (ii) when the ionic strength was increased indicating that the cytoplasmic domain of Vpu affects oligomerization. Coarse-grained molecular dynamic simulations with models of wild-type and mutant Vpu in a hydrated lipid bilayer supported the experimental data in demonstrating that, in addition to a previously known role in downregulation of host factors, the phosphorylation sites of Vpu also modulate oligomerization.

Viral channel forming proteins--How to assemble and depolarize lipid membranes in silico

Viral channel forming proteins (VCPs) have been discovered in the late 70s and are found in many viruses to date. Usually they are small and have to assemble to form channels which depolarize the lipid membrane of the host cells. Structural information is just about to emerge for just some of them. Thus, computational methods play a pivotal role in generating plausible structures which can be used in the drug development process. In this review the accumulation of structural data is introduced from a historical perspective. Computational performances and their predictive power are reported guided by biological questions such as the assembly, mechanism of function and drug–protein interaction of VCPs. An outlook of how coarse grained simulations can contribute to yet unexplored issues of these proteins is given. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

•

Early references about the discovery of viral channel forming proteins.

•

Latest structural information about the class of proteins.

•

Identification of structural motifs, assembly mechanism of function and drug action using computational methods.

•

Outlook for the use of coarse grained techniques to address assembly and integration into cellular processes.

Patch formation of a viral channel forming protein within a lipid membrane – Vpu of HIV-1

Dimer-first formation leads to larger assemblies with potentially relevant structures.

Folding Simulations of the Transmembrane Helix of Virus Protein U in an Implicit Membrane Model

Vpu is an 81-amino-acid auxiliary membrane protein encoded by human immunodeficiency virus type 1 (HIV-1). One of its roles is to amplify viral release by self-assembling in homo-oligomers to form functional water-filled pores enabling the flux of ions across the membrane. Various NMR and CD studies have shown that the transmembrane domain of Vpu has a helical conformation. With a recently developed implicit membrane model and an efficient Monte Carlo (MC) algorithm using concerted backbone rotations, we simulate the folding of the transmembrane domain of Vpu at atomic resolution. The implicit membrane environment is based on the generalized Born theory and enables very long time scale events, such as folding to be observed using detailed all-atom representation of the protein. Such studies are currently computationally unfeasible with fully explicit lipid bilayer molecular dynamics simulations. The correct helical transmembrane structure of Vpu is predicted from extended conformations and remains stably inserted. Tilt and kink angles agree well with experimental estimates from NMR measurements. The experimentally observed change in tilt angle in membranes of varying hydrophobic width is accurately reproduced. The extensive simulation of a pentamer of the Vpu transmembrane domain in the implicit membrane gives results similar to the ones reported previously for fully explicit bilayer simulations.

Structural determination of virus protein U from HIV-1 by NMR in membrane environments

Virus protein U (Vpu) from HIV-1, a small membrane protein composed of a transmembrane helical domain and two α-helices in an amphipathic cytoplasmic domain, down modulates several cellular proteins, including CD4, BST-2/CD317/tetherin, NTB-A, and CCR7. The interactions of Vpu with these proteins interfere with the immune system and enhance the release of newly synthesized virus particles. It is essential to characterize the structure and dynamics of Vpu in order to understand the mechanisms of the protein-protein interactions, and potentially to discover antiviral drugs. In this article, we describe investigations of the cytoplasmic domain of Vpu as well as full-length Vpu by NMR spectroscopy. These studies are complementary to earlier analysis of the transmembrane domain of Vpu. The results suggest that the two helices in the cytoplasmic domain form a U-shape. The length of the inter-helical loop in the cytoplasmic domain and the orientation of the third helix vary with the lipid composition, which demonstrate that the C-terminal helix is relatively flexible, providing accessibility for interaction partners.

In silico analysis reveals sequential interactions and protein conformational changes during the binding of chemokine CXCL-8 to its receptor CXCR1

Chemokine CXCL-8 plays a central role in human immune response by binding to and activate its cognate receptor CXCR1, a member of the G-protein coupled receptor (GPCR) family. The full-length structure of CXCR1 is modeled by combining the structures of previous NMR experiments with those from homology modeling. Molecular docking is performed to search favorable binding sites of monomeric and dimeric CXCL-8 with CXCR1 and a mutated form of it. The receptor-ligand complex is embedded into a lipid bilayer and used in multi ns molecular dynamics (MD) simulations. A multi-steps binding mode is proposed: (i) the N-loop of CXCL-8 initially binds to the N-terminal domain of receptor CXCR1 driven predominantly by electrostatic interactions; (ii) hydrophobic interactions allow the N-terminal Glu-Leu-Arg (ELR) motif of CXCL-8 to move closer to the extracellular loops of CXCR1; (iii) electrostatic interactions finally dominate the interaction between the N-terminal ELR motif of CXCL-8 and the EC-loops of CXCR1. Mutation of CXCR1 abrogates this mode of binding. The detailed binding process may help to facilitate the discovery of agonists and antagonists for rational drug design.

Mechanism of function of viral channel proteins and implications for drug development

Viral channel-forming proteins comprise a class of viral proteins which, similar to their host companions, are made to alter electrochemical or substrate gradients across lipid membranes. These proteins are active during all stages of the cellular life cycle of viruses. An increasing number of proteins are identified as channel proteins, but the precise role in the viral life cycle is yet unknown for the majority of them. This review presents an overview about these proteins with an emphasis on those with available structural information. A concept is introduced which aligns the transmembrane domains of viral channel proteins with those of host channels and toxins to give insights into the mechanism of function of the viral proteins from potential sequence identities. A summary of to date investigations on drugs targeting these proteins is given and discussed in respect of their mode of action in vivo.

Membrane raft association of the Vpu protein of human immunodeficiency virus type 1 correlates with enhanced virus release

The Vpu protein of human immunodeficiency virus type 1 (HIV-1) is known to enhance virion release from certain cell types. To accomplish this function, Vpu interacts with the restriction factor known as bone marrow stromal cell antigen 2 (BST-2)/tetherin. In this study, we analyzed whether the Vpu protein is associated with microdomains known as lipid or membrane rafts. Our results indicate that Vpu partially partitions into detergent-resistant membrane (DRM) fractions when expressed alone or in the context of simian-human immunodeficiency virus (SHIV) infection. The ability to be partitioned into rafts was observed with both subtype B and C Vpu proteins. The use of cholesterol lowering lovastatin/M-β-cyclodextrin and co-patching experiments confirmed that Vpu can be detected in cholesterol rich regions of membranes. Finally, we present data showing that raft association-defective transmembrane mutants of Vpu have impaired enhanced virus release function, but still maintain the ability to down-regulate CD4.

Structural implications of mutations assessed by molecular dynamics: Vpu1-32 from HIV-1

Structural pore models are generated for Vpu(1-32)WT from HIV-1 as well as for three mutants W23L, S24L and R31V. A computational methodology is employed which samples the whole conformational space of the pentameric assemblies of Vpu. The analysis of the related energy landscape reveals a small set of reasonable pore models, which are thoroughly investigated regarding their structural properties as well as their putative stability under native-like conditions. The models are also discussed in respect of earlier experimental findings about their channel activities. The study proposes functional pores reflecting the experimentally found conductance states of Vpu and its mutants.

Exploring the conformational space of Vpu from HIV-1: a versatile adaptable protein

The dynamic behavior of monomeric Vpu(1-32) from HIV-1 in different lipid environments has been studied. The peptide shows highly flexible behavior during the simulations and easily adapts to changing lipid environments as it experiences when travelling through the Golgi apparatus. Protein-lipid interactions do not show any significant correlation towards lipid type or thickness based on multiple 10 ns simulations. The averaged structure of a series of 16 independent simulations suggest kink around Ser-24, which compensates the polarity of its side chain by forming hydrogen bonds with the carbonyl backbone of adjacent amino acids towards the N-terminus.

Structural and functional characterization of human SGT and its interaction with Vpu of the human immunodeficiency virus type 1

The small glutamine-rich tetratricopeptide repeat protein (SGT) belongs to a family of cochaperones that interacts with both Hsp70 and Hsp90 via the so-called TPR domain. Here, we present the crystal structure of the TPR domain of human SGT (SGT-TPR), which shows that it contains typical features found in the structures of other TPR domains. Previous studies show that full-length SGT can bind to both Vpu and Gag of human immunodeficiency virus type 1 (HIV-1) and the overexpression of SGT in cells reduces the efficiency of HIV-1 particle release. We show that SGT-TPR can bind Vpu and reduce the amount of HIV-1 p24, which is the viral capsid, secreted from cells transfected with the HIV-1 proviral construct, albeit at a lower efficiency than full-length SGT. This indicates that the TPR domain of SGT is sufficient for the inhibition of HIV-1 particle release but the N- and/or C-terminus also have some contributions. The SGT binding site in Vpu was also identified by using peptide array and confirmed by GST pull-down assay.

Molecular characterization of the HIV type 1 subtype C accessory genes vif, vpr, and vpu

HIV-1 Vif, Vpr, and Vpu proteins have a profound effect on efficient viral replication and pathogenesis. This study describes the genotypic characterisation of vif , vpr and vpu from 20 South African HIV-1 subtype C primary isolates, and extensive analysis and comparison of known motifs. All HIV-1 subtype C Vif, Vpr and Vpu proteins revealed the presence of highly conserved structural and functional motifs similar to other sub-types, for example, the Vif-APOBEC3G interaction domains. However, several differences were noted when these sequences were compared to subtype B, such as the presence of the LRLL motif which has been implicated in targeting subtype C Vpu predominantly to the cell surface, instead of the Golgi apparatus. A better understanding of the structure/function relationship of these proteins may lead to the development of new classes of antiviral drugs. These results indicate that antiviral drugs that target the conserved functional domains within Vif, Vpr or Vpu could be active against all circulating subtypes, including HIV-1 subtype C.

Full Length Vpu from HIV-1: Combining Molecular Dynamics Simulations with NMR Spectroscopy

Based on structures made available by solution NMR, molecular models of the protein Vpu from HIV-1 were built and refined by 6 ns MD simulations in a fully hydrated lipid bilayer. Vpu is an 81 amino acid type I integral membrane protein encoded by the human immunodeficiency virus type-1 (HIV-1) and closely related simian immunodeficiency viruses (SIVs). Its role is to amplify viral release. Upon phosphorylation, the cytoplasmic domain adopts a more compact shape with helices 2 and 3 becoming almost parallel to each other. A loss of helicity for several residues belonging to the helices adjacent to both ends of the loop region containing serines 53 and 57 is observed. A fourth helix, present in one of the NMR-based structures of the cytoplasmic domain and located near the C-terminus, is lost upon phosphorylation.

Structure and Dynamics of the HIV-1 Vpu Transmembrane Domain Revealed by Solid-State NMR with Magic-Angle Spinning†

We report solid-state nuclear magnetic resonance (NMR) measurements on the peptide Vpu(1-40), comprising residues 1-40 of the 81-residue type 1 integral membrane protein Vpu encoded by the HIV-1 genome. On the basis of a combination of 13C and 15N NMR chemical shifts under magic-angle spinning (MAS), effects of local mobility on NMR signal intensities, site-specific MAS NMR line widths, and NMR-detected hydrogen-deuterium exchange, we develop a model for the structure and dynamics of the Vpu(1-40) monomer in phospholipid bilayer membranes. Our data are largely consistent with earlier structural studies of Vpu peptides by Opella and co-workers, in which solution NMR and solid-state NMR without MAS were used, but our data provide new information about local variations in the degree of mobility and structural order. In addition, our data indicate that the transmembrane alpha-helix of Vpu(1-40) extends beyond the hydrophobic core of the bilayer. We find no evidence for heterogeneity in the conformation and intermolecular contacts of the transmembrane alpha-helix, with the exception of two distinct chemical shifts observed for the C alpha and C beta atoms of A18 that may reflect distinct modes of helix-helix interaction. These results have possible implications for the supramolecular structure of Vpu oligomers that form cation-selective ion channels.

NMR studies for identifying phosphopeptide ligands of the HIV-1 protein Vpu binding to the F-box protein beta-TrCP

The human immunodeficiency virus type 1 (HIV-1) Vpu enhances viral particle release and, its interaction with the ubiquitin ligase SCF-beta-TrCP triggers the HIV-1 receptor CD4 degradation by the proteasome. The interaction between beta-TrCP protein and ligands containing the phosphorylated DpSGXXpS motif plays a key role for the development of severe disease states, such as HIV or cancer. This study examines the binding and conformation of phosphopeptides (P1, LIERAEDpSG and P2, EDpSGNEpSE) from HIV protein Vpu to beta-TrCP with the objective of defining the minimum length of peptide needed for effective binding. The screening step can be analyzed by NMR spectroscopy, in particular, saturation transfer NMR methods clearly identify the residues in the peptide that make direct contact with beta-TrCP protein when bound. An analysis of saturation transfer difference (STD) spectra provided clear evidence that the two peptides efficiently bound beta-TrCP receptor protein. To better characterize the ligand-protein interaction, the bound conformation of the phosphorylated peptides was determined using transferred NOESY methods, which gave rise to a well-defined structure. P1 and P2 can fold in a bend arrangement for the DpSG motif, showing the protons identified by STD-NMR as exposed in close proximity at the molecule surface. Ser phosphorylation allows electrostatic interaction and hydrogen bond with the amino acids of the beta-TrCP binding pocket. The upstream LIER hydrophobic region was also essential in binding to a hydrophobic pocket of the beta-TrCP WD domain. These findings are in good agreement with a recently published X-ray structure of a shorter beta-Catenin fragment with the beta-TrCP complex.

Computational analysis of mutations in the transmembrane region of Vpu from HIV-1

Vpu is an 81 amino acid integral membrane protein encoded by HIV-1. Its alpha-helical transmembrane (TM) domain (residues approximately 6-28) enhances virion release by oligomerizing into bundles and forming ion-conducting channels across the plasma membrane. Its cytoplasmic domain (residues approximately 29-81) is also alpha-helical and binds to the transmembrane protein CD4, inducing its degradation. Mutations within the TM domain have been found to abrogate enhanced particle release from the infected cell (Tiganos et al. Virology (1998) 251 96-107). A series of computational models of monomeric, pentameric and hexameric Vpu(1-31) mutants have been constructed, embedded in fully hydrated lipid bilayers and subjected to a 3 ns molecular dynamics (MD) simulation. None of the mutations has any destabilizing effect on the secondary and tertiary structure. One of the mutants, in which the position of a tryptophan residue within the TM domain is altered, is known not to induce CD4 degradation; an extended kinked model of this mutant has been generated (Vpu(1-52)IVW-k) and during subsequent MD simulations, the bend between the TM and a part of the cytoplasmic domain is found to unwind and a complex salt bridge involving Lys-37 is formed.

Interaction of amiloride and one of its derivatives with Vpu from HIV-1: a molecular dynamics simulation

Vpu is an 81-residue membrane protein, with a single transmembrane segment that is encoded by HIV-1 and is involved in the enhancement of virion release via formation of an ion channel. Cyclohexamethylene amiloride (Hma) has been shown to inhibit ion channel activity. In the present 12-ns simulation study a putative binding site of Hma blockers in a pentameric model bundle built of parallel aligned helices of the first 32 residues of Vpu was found near Ser-23. Hma orientates along the channel axis with its alkyl ring pointing inside the pore, which leads to a blockage of the pore.

Vpu from HIV-1 on an atomic scale: experiments and computer simulations

Vpu is an 81 amino acid protein encoded by HIV-1. Its role is to amplify viral release by two mechanisms: (i) docking to CD4 with the consequence of targeting CD4 for ubiquitine-mediated degradation, and (ii) formation of ion channels to enhance particle release. The intensive research on its in vivo function, combined with structural investigations, makes this viral membrane protein one of the better characterised membrane proteins. The wealth of structural information enables the use of computational methods to elucidate the mechanisms of function on an atomic scale. The discovery of Vpu and the development of structural models in a chronological order is summarised and first efforts on investigating the mechanics are outlined.