Browse
Submit Data
Databases
API
Help

Dataset Information

30 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Molecular dynamics simulations reveal the HIV-1 Vpu transmembrane protein to form stable pentamers.

ABSTRACT: The human immunodeficiency virus type I (HIV-1) Vpu protein is 81 residues long and has two cytoplasmic and one transmembrane (TM) helical domains. The TM domain oligomerizes to form a monovalent cation selective ion channel and facilitates viral release from host cells. Exactly how many TM domains oligomerize to form the pore is still not understood, with experimental studies indicating the existence of a variety of oligomerization states. In this study, molecular dynamics (MD) simulations were performed to investigate the propensity of the Vpu TM domain to exist in tetrameric, pentameric, and hexameric forms. Starting with an idealized ?-helical representation of the TM domain, a thorough search for the possible orientations of the monomer units within each oligomeric form was carried out using replica-exchange MD simulations in an implicit membrane environment. Extensive simulations in a fully hydrated lipid bilayer environment on representative structures obtained from the above approach showed the pentamer to be the most stable oligomeric state, with interhelical van der Waals interactions being critical for stability of the pentamer. Atomic details of the factors responsible for stable pentamer structures are presented. The structural features of the pentamer models are consistent with existing experimental information on the ion channel activity, existence of a kink around the Ile17, and the location of tetherin binding residues. Ser23 is proposed to play an important role in ion channel activity of Vpu and possibly in virus propagation.

SUBMITTER: Padhi S

PROVIDER: S-EPMC3819244 | biostudies-literature | 2013

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

HIV-1 Vpu protein forms stable oligomers in aqueous solution via its transmembrane domain self-association.

Project description:We report our findings on the assembly of the HIV-1 protein Vpu into soluble oligomers. Vpu is a key to HIV-1 protein. It has been considered exclusively a single-pass membrane protein. However, we revealed that this protein forms stable oligomers in aqueous solution, which is an interesting and rather unique observation, as the number of proteins transitioning between soluble and membrane embedded states is limited. Therefore, we undertook a study to characterize these oligomers by utilizing protein engineering, size exclusion chromatography, cryoEM and electron paramagnetic resonance (EPR) spectroscopy. We found that Vpu oligomerizes via its N-terminal transmembrane domain (TM). CryoEM analyses suggest that the oligomeric state most likely is a hexamer or hexamer-to-heptamer equilibrium. Both cryoEM and EPR suggest that, within the oligomer, the distant C-terminal region of Vpu is highly flexible. To the best of our knowledge, this is the first comprehensive study on soluble Vpu. We propose that these oligomers are stabilized via possibly hydrophobic interactions between Vpu TMs. Our findings contribute valuable information about this protein properties and about protein supramolecular complexes formation. The acquired knowledge could be further used in protein engineering, and could also help to uncover possible physiological function of these Vpu oligomers.

| S-EPMC10197565 | biostudies-literature

Sidekick for Membrane Simulations: Automated Ensemble Molecular Dynamics Simulations of Transmembrane Helices.

Project description:The interactions of transmembrane (TM) α-helices with the phospholipid membrane and with one another are central to understanding the structure and stability of integral membrane proteins. These interactions may be analyzed via coarse grained molecular dynamics (CGMD) simulations. To obtain statistically meaningful analysis of TM helix interactions, large (N ca. 100) ensembles of CGMD simulations are needed. To facilitate the running and analysis of such ensembles of simulations, we have developed Sidekick, an automated pipeline software for performing high throughput CGMD simulations of α-helical peptides in lipid bilayer membranes. Through an end-to-end approach, which takes as input a helix sequence and outputs analytical metrics derived from CGMD simulations, we are able to predict the orientation and likelihood of insertion into a lipid bilayer of a given helix of a family of helix sequences. We illustrate this software via analyses of insertion into a membrane of short hydrophobic TM helices containing a single cationic arginine residue positioned at different positions along the length of the helix. From analyses of these ensembles of simulations, we estimate apparent energy barriers to insertion which are comparable to experimentally determined values. In a second application, we use CGMD simulations to examine the self-assembly of dimers of TM helices from the ErbB1 receptor tyrosine kinase and analyze the numbers of simulation repeats necessary to obtain convergence of simple descriptors of the mode of packing of the two helices within a dimer. Our approach offers a proof-of-principle platform for the further employment of automation in large ensemble CGMD simulations of membrane proteins.

| S-EPMC4871227 | biostudies-literature

Continuous Constant pH Molecular Dynamics Simulations of Transmembrane Proteins.

Project description:Many membrane channels, transporters, and receptors utilize a pH gradient or proton coupling to drive functionally relevant conformational transitions. Conventional molecular dynamics simulations employ fixed protonation states, thus neglecting the coupling between protonation and conformational equilibria. Here we describe the membrane-enabled hybrid-solvent continuous constant pH molecular dynamics method for capturing atomic details of proton-coupled conformational dynamics of transmembrane proteins. Example protocols from our recent application studies of proton channels and ion/substrate transporters are discussed.

| S-EPMC8062021 | biostudies-literature

Dynamics of lipids, cholesterol, and transmembrane ?-helices from microsecond molecular dynamics simulations.

Project description:Extensive all-atom molecular dynamics (?24 ?s total) allowed exploration of configurational space and calculation of lateral diffusion coefficients of the components of a protein-embedded, cholesterol-containing model bilayer. The three model membranes are composed of an ?50/50 (by mole) dipalmitoylphosphatidylcholine (DPPC)/cholesterol bilayer and contained an ?-helical transmembrane protein (HIV-1 gp41 TM). Despite the high concentration of cholesterol, normal Brownian motion was observed and the calculated diffusion coefficients (on the order of 10(-9) cm(2)/s) are consistent with experiments. Diffusion is sensitive to a variety of parameters, and a temperature difference of ?4 K from thermostat artifacts resulted in 2-10-fold differences in diffusion coefficients and significant differences in lipid order, membrane thickness, and unit cell area. Also, the specific peptide sequence likely underlies the consistently observed faster diffusion in one leaflet. Although the simulations here present molecular dynamics (MD) an order of magnitude longer than those from previous studies, the three systems did not approach ergodicity. The distributions of cholesterol and DPPC around the peptides changed on the microsecond time scale, but not significantly enough to thoroughly explore configurational space. These simulations support conclusions of other recent microsecond MD in that even longer time scales are needed for equilibration of model membranes and simulations of more realistic cellular or viral bilayers.

| S-EPMC4254001 | biostudies-literature

Coarse grained molecular dynamics simulations of transmembrane protein-lipid systems.

Project description:Many biological cellular processes occur at the micro- or millisecond time scale. With traditional all-atom molecular modeling techniques it is difficult to investigate the dynamics of long time scales or large systems, such as protein aggregation or activation. Coarse graining (CG) can be used to reduce the number of degrees of freedom in such a system, and reduce the computational complexity. In this paper the first version of a coarse grained model for transmembrane proteins is presented. This model differs from other coarse grained protein models due to the introduction of a novel angle potential as well as a hydrogen bonding potential. These new potentials are used to stabilize the backbone. The model has been validated by investigating the adaptation of the hydrophobic mismatch induced by the insertion of WALP-peptides into a lipid membrane, showing that the first step in the adaptation is an increase in the membrane thickness, followed by a tilting of the peptide.

| S-EPMC2904924 | biostudies-literature

Transmembrane domain determinants of CD4 Downregulation by HIV-1 Vpu.

Project description:The transmembrane domains (TMDs) of integral membrane proteins do not merely function as membrane anchors but play active roles in many important biological processes. The downregulation of the CD4 coreceptor by the Vpu protein of HIV-1 is a prime example of a process that is dependent on specific properties of TMDs. Here we report the identification of Trp22 in the Vpu TMD and Gly415 in the CD4 TMD as critical determinants of Vpu-induced targeting of CD4 to endoplasmic reticulum (ER)-associated degradation (ERAD). The two residues participate in different aspects of ERAD targeting. Vpu Trp22 is required to prevent assembly of Vpu into an inactive, oligomeric form and to promote CD4 polyubiquitination and subsequent recruitment of the VCP-UFD1L-NPL4 dislocase complex. In the presence of a Vpu Trp22 mutant, CD4 remains integrally associated with the ER membrane, suggesting that dislocation from the ER into the cytosol is impaired. CD4 Gly415, on the other hand, contributes to CD4-Vpu interactions. We also identify two residues, Val20 and Ser23, in the Vpu TMD that mediate retention of Vpu and, by extension, CD4 in the ER. These findings highlight the exploitation of several TMD-mediated mechanisms by HIV-1 Vpu in order to downregulate CD4 and thus promote viral pathogenesis.

| S-EPMC3255848 | biostudies-literature

Heat-Stable Hazelnut Profilin: Molecular Dynamics Simulations and Immunoinformatics Analysis.

Project description:Heat treatment can modify the allergenic potential, reducing allergenicity in specific proteins. Profilins are one of the important hazelnut allergens; these proteins are considered panallergens due to their high capacity for cross-reactivity with other allergens. In the present work, we evaluated the thermostability of hazelnut profilin, combining molecular dynamics simulation and immunoinformatic techniques. This approach helped us to have reliable results in immunogenicity studies. We modeled Cor a 2 profilin and applied annealing simulation, equilibrium, and production simulation at constant temperatures ranging from 300 to 500 K using Gromacs software. Despite the hazelnut profilins being able to withstand temperatures of up to 400 K, this does not seem to reduce its allergenicity. We have found that profilin subjected to temperatures of 450 and 500 K could generate cross-reactivity with other food allergens. In conclusion, we note a remarkable thermostability of Cor a 2 at 400 K which avoids its structural unfolding.

| S-EPMC7464029 | biostudies-literature

Molecular Dynamics Simulations of Influenza A Virus NS1 Reveal a Remarkably Stable RNA-Binding Domain Harboring Promising Druggable Pockets.

Project description:The non-structural protein NS1 of influenza A viruses is considered to be the major antagonist of the interferon system and antiviral defenses of the cell. It could therefore represent a suitable target for novel antiviral strategies. As a first step towards the identification of small compounds targeting NS1, we here investigated the druggable potential of its RNA-binding domain since this domain is essential to the biological activities of NS1. We explored the flexibility of the full-length protein by running molecular dynamics simulations on one of its published crystal structures. While the RNA-binding domain structure was remarkably stable along the simulations, we identified a flexible site at the two extremities of the "groove" that is delimited by the antiparallel ?-helices that make up its RNA-binding interface. This groove region is able to form potential binding pockets, which, in 60% of the conformations, meet the druggability criteria. We characterized these pockets and identified the residues that contribute to their druggability. All the residues involved in the druggable pockets are essential at the same time to the stability of the RNA-binding domain and to the biological activities of NS1. They are also strictly conserved across the large sequence diversity of NS1, emphasizing the robustness of this search towards the identification of broadly active NS1-targeting compounds.

| S-EPMC7290946 | biostudies-literature

Molecular Dynamics Simulations Reveal the Conformational Transition of GH33 Sialidases.

Project description:Sialidases are increasingly used in the production of sialyloligosaccharides, a significant component of human milk oligosaccharides. Elucidating the catalytic mechanism of sialidases is critical for the rational design of better biocatalysts, thereby facilitating the industrial production of sialyloligosaccharides. Through comparative all-atom molecular dynamics simulations, we investigated the structural dynamics of sialidases in Glycoside Hydrolase family 33 (GH33). Interestingly, several sialidases displayed significant conformational transition and formed a new cleft in the simulations. The new cleft was adjacent to the innate active site of the enzyme, which serves to accommodate the glycosyl acceptor. Furthermore, the residues involved in the specific interactions with the substrate were evolutionarily conserved in the whole GH33 family, highlighting their key roles in the catalysis of GH33 sialidases. Our results enriched the catalytic mechanism of GH33 sialidases, with potential implications in the rational design of sialidases.

| S-EPMC10095477 | biostudies-literature

Elucidating activation and deactivation dynamics of VEGFR-2 transmembrane domain with coarse-grained molecular dynamics simulations.

Project description:The vascular endothelial growth factor receptor 2 (VEGFR-2) is a member of receptor tyrosine kinases (RTKs) and is a dimeric membrane protein that functions as a primary regulator of angiogenesis. As is usual with RTKs, spatial alignment of its transmembrane domain (TMD) is essential toward VEGFR-2 activation. Experimentally, the helix rotations within TMD around their own helical axes are known to participate importantly toward the activation process in VEGFR-2, but the detailed dynamics of the interconversion between the active and inactive TMD forms have not been clearly elucidated at the molecular level. Here, we attempt to elucidate the process by using coarse grained (CG) molecular dynamics (MD) simulations. We observe that inactive dimeric TMD in separation is structurally stable over tens of microseconds, suggesting that TMD itself is passive and does not allow spontaneous signaling of VEGFR-2. By starting from the active conformation, we reveal the mechanism of TMD inactivation through analyzing the CG MD trajectories. We observe that interconversions between a left-handed overlay and a right-handed one are essential for the process of going from an active TMD structure to the inactive form. In addition, our simulations find that the helices can rotate properly when the overlaying structure of the helices interconverts and when the crossing angle of the two helices changes by larger than ~40 degrees. As the activation right after the ligand attachment on VEGFR-2 will take place in the reverse manner of this inactivation process, these structural aspects will also appear importantly for the activation process. The rather large change in helix configuration for activation also explains why VEGFR-2 rarely self-activate and how the activating ligand structurally drive the whole VEGFR-2. This mechanism of TMD activation / inactivation within VEGFR-2 may help in further understanding the overall activation processes of other RTKs.

| S-EPMC9934429 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data