KEAP1 Cancer Mutants: A Large-Scale Molecular Dynamics Study of Protein Stability.
Ontology highlight
ABSTRACT:
SUBMITTER: Wilson CJ
PROVIDER: S-EPMC8161161 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Surfactant protein D (SP-D), a mammalian C-type lectin, is the primary innate inhibitor of influenza A virus (IAV) in the lung. Interactions of SP-D with highly branched viral N-linked glycans on hemagglutinin (HA), an abundant IAV envelope protein and critical virulence factor, promote viral aggregation and neutralization through as yet unknown molecular mechanisms. Two truncated human SP-D forms, wild-type (WT) and double mutant D325A+R343V, representing neck and carbohydrate recognition domains are compared in this study. Whereas both WT and D325A+R343V bind to isolated glycosylated HA, WT does not inhibit IAV in neutralization assays; in contrast, D325A+R343V neutralization compares well with that of full-length native SP-D. To elucidate the mechanism for these biochemical observations, we have determined crystal structures of D325A+R343V in the presence and absence of a viral nonamannoside (Man9). On the basis of the D325A+R343V-Man9 structure and other crystallographic data, models of complexes between HA and WT or D325A+R343V were produced and subjected to molecular dynamics. Simulations reveal that whereas WT and D325A+R343V both block the sialic acid receptor site of HA, the D325A+R343V complex is more stable, with stronger binding caused by additional hydrogen bonds and hydrophobic interactions with HA residues. Furthermore, the blocking mechanism of HA differs for WT and D325A+R343V because of alternate glycan binding modes. The combined results suggest a mechanism through which the mode of SP-D-HA interaction could significantly influence viral aggregation and neutralization. These studies provide the first atomic-level molecular view of an innate host defense lectin inhibiting its viral glycoprotein target.
| S-EPMC3927399 | biostudies-literature
Project description:We show how to apply the method of temperature-accelerated molecular dynamics (TAMD) in collective variables [Maragliano L, Vanden-Eijnden E (2006) Chem Phys Lett 426:168-175] to sample the conformational space of multidomain proteins in all-atom, explicitly solvated molecular dynamics simulations. The method allows the system to hyperthermally explore the free-energy surface in a set of collective variables computed at the physical temperature. As collective variables, we pick Cartesian coordinates of centers of contiguous subdomains. The method is applied to the GroEL subunit, a 55-kDa, three-domain protein, and HIV-1 gp120. For GroEL, the method induces in about 40 ns conformational changes that recapitulate the t --> r('') transition and are not observed in unaccelerated molecular dynamics: The apical domain is displaced by 30 A, with a twist of 90 degrees relative to the equatorial domain, and the root-mean-squared deviation relative to the r('') conformer is reduced from 13 to 5 A, representing fairly high predictive capability. For gp120, the method predicts both counterrotation of inner and outer domains and disruption of the so-called bridging sheet. In particular, TAMD on gp120 initially in the CD4-bound conformation visits conformations that deviate by 3.6 A from the gp120 conformer in complex with antibody F105, again reflecting good predictive capability. TAMD generates plausible all-atom models of the so-far structurally uncharacterized unliganded conformation of HIV-1 gp120, which may prove useful in the development of inhibitors and immunogens. The fictitious temperature employed also gives a rough estimate of 10 kcal/mol for the free-energy barrier between conformers in both cases.
| S-EPMC2841907 | biostudies-other
Project description:Cryo-electron microscopy (cryoEM) has become a well established technique with the potential to produce structures of large and dynamic supramolecular complexes that are not amenable to traditional approaches for studying structure and dynamics. The size and low resolution of such molecular systems often make structural modelling and molecular dynamics simulations challenging and computationally expensive. This, together with the growing wealth of structural data arising from cryoEM and other structural biology methods, has driven a trend in the computational biophysics community towards the development of new pipelines for analysing global dynamics using coarse-grained models and methods. At the centre of this trend has been a return to elastic network models, normal mode analysis (NMA) and ensemble analyses such as principal component analysis, and the growth of hybrid simulation methodologies that make use of them. Here, this field is reviewed with a focus on ProDy, the Python application programming interface for protein dynamics, which has been developed over the last decade. Two key developments in this area are highlighted: (i) ensemble NMA towards extracting and comparing the signature dynamics of homologous structures, aided by the recent SignDy pipeline, and (ii) pseudoatom fitting for more efficient global dynamics analyses of large and low-resolution supramolecular assemblies from cryoEM, revisited in the CryoDy pipeline. It is believed that such a renewal and extension of old models and methods in new pipelines will be critical for driving the field forward into the next cryoEM revolution.
| S-EPMC8972803 | biostudies-literature
Project description:HIV-1 integrase is one of the three essential enzymes required for viral replication and has great potential as a novel target for anti-HIV drugs. Although tremendous efforts have been devoted to understanding this protein, the conformation of the catalytic core domain around the active site, particularly the catalytic loop overhanging the active site, is still not well characterized by experimental methods due to its high degree of flexibility. Recent studies have suggested that this conformational dynamics is directly correlated with enzymatic activity, but the details of this dynamics is not known. In this study, we conducted a series of extended-time molecular dynamics simulations and locally enhanced sampling simulations of the wild-type and three loop hinge mutants to investigate the conformational dynamics of the core domain. A combined total of >480 ns of simulation data was collected which allowed us to study the conformational changes that were not possible to observe in the previously reported short-time molecular dynamics simulations. Among the main findings are a major conformational change (>20 A) in the catalytic loop, which revealed a gatinglike dynamics, and a transient intraloop structure, which provided a rationale for the mutational effects of several residues on the loop including Q(148), P(145), and Y(143). Further, clustering analyses have identified seven major conformational states of the wild-type catalytic loop. Their implications for catalytic function and ligand interaction are discussed. The findings reported here provide a detailed view of the active site conformational dynamics and should be useful for structure-based inhibitor design for integrase.
| S-EPMC1305464 | biostudies-literature
Project description:As the length of molecular dynamics (MD) trajectories grows with increasing computational power, so does the importance of clustering methods for partitioning trajectories into conformational bins. Of the methods available, the vast majority require users to either have some a priori knowledge about the system to be clustered or to tune clustering parameters through trial and error. Here we present non-parametric uses of two modern clustering techniques suitable for first-pass investigation of an MD trajectory. Being non-parametric, these methods require neither prior knowledge nor parameter tuning. The first method, HDBSCAN, is fast-relative to other popular clustering methods-and is able to group unstructured or intrinsically disordered systems (such as intrinsically disordered proteins, or IDPs) into bins that represent global conformational shifts. HDBSCAN is also useful for determining the overall stability of a system-as it tends to group stable systems into one or two bins-and identifying transition events between metastable states. The second method, iMWK-Means, with explicit rescaling followed by K-Means, while slower than HDBSCAN, performs well with stable, structured systems such as folded proteins and is able to identify higher resolution details such as changes in relative position of secondary structural elements. Used in conjunction, these clustering methods allow a user to discern quickly and without prior knowledge the stability of a simulated system and identify both local and global conformational changes.
| S-EPMC5719493 | biostudies-literature
Project description:Intrinsically disordered linkers provide multi-domain proteins with degrees of conformational freedom that are often essential for function. These highly dynamic assemblies represent a significant fraction of all proteomes, and deciphering the physical basis of their interactions represents a considerable challenge. Here we describe the difficulties associated with mapping the large-scale domain dynamics and describe two recent examples where solution state methods, in particular NMR spectroscopy, are used to investigate conformational exchange on very different timescales.
| S-EPMC5020063 | biostudies-literature
Project description:Manassantin A is a natural product that has been shown to have anticancer activity in cell-based assays, but has a largely unknown mode-of-action. Described here is the use of two different energetics-based approaches to identify protein targets of manassantin A. Using the stability of proteins from rates of oxidation technique with an isobaric mass tagging strategy (iTRAQ-SPROX) and the pulse proteolysis technique with a stable isotope labeling with amino acids in cell culture strategy (SILAC-PP), over 1000 proteins in a MDA-MB-231 cell lysate grown under hypoxic conditions were assayed for manassantin A interactions (both direct and indirect). A total of 28 protein hits were identified with manassantin A-induced thermodynamic stability changes. Two of the protein hits (filamin A and elongation factor 1?) were identified using both experimental approaches. The remaining 26 hit proteins were only assayed in either the iTRAQ-SPROX or the SILAC-PP experiment. The 28 potential protein targets of manassantin A identified here provide new experimental avenues along which to explore the molecular basis of manassantin A's mode of action. The current work also represents the first application iTRAQ-SPROX and SILAC-PP to the large-scale analysis of protein-ligand binding interactions involving a potential anticancer drug with an unknown mode-of-action.
| S-EPMC5270646 | biostudies-literature
Project description: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.
| S-EPMC9106056 | biostudies-literature
Project description:Advances in biomolecular simulation methods and access to large scale computer resources have led to a massive increase in the amount of data generated. The key enablers have been optimization and parallelization of the simulation codes. However, much of the software used to analyze trajectory data from these simulations is still run in serial, or in some cases many threads via shared memory. Here, we describe the addition of multiple levels of parallel trajectory processing to the molecular dynamics simulation analysis software CPPTRAJ. In addition to the existing OpenMP shared-memory parallelism, CPPTRAJ now has two additional levels of message passing (MPI) parallelism involving both across-trajectory processing and across-ensemble processing. All three levels of parallelism can be simultaneously active, leading to significant speed ups in data analysis of large datasets on the NCSA Blue Waters supercomputer by better leveraging the many available nodes and its parallel file system. © 2018 Wiley Periodicals, Inc.
| S-EPMC7313716 | biostudies-literature
Project description:Large-scale conformational changes in proteins that happen often on biological time scales may be relatively rare events on the molecular dynamics time scale. We have implemented an approach to targeted molecular dynamics called end-point targeted molecular dynamics that transforms proteins between two specified conformational states through the use of nonharmonic "soft" restraints. A key feature of the method is that the protein is free to discover its own conformational pathway through the plethora of possible intermediate states. The method is applied to the Shaker K(v)1.2 potassium channel in implicit solvent. The rate of cycling between the open and closed states was varied to explore how slow the cycling rate needed to be to ensure that microscopic reversibility along the transition pathways was well approximated. Results specific to the K(+) channel include: 1), a variation in backbone torsion angles of residues near the Pro-Val-Pro motif in the inner helix during both opening and closing; 2), the identification of possible occlusion sites in the closed channel located among Pro-Val-Pro residues and downstream; 3), a difference in the opening and closing pathways of the channel; and 4), evidence of a transient intermediate structural substate. The results also show that likely intermediate conformations during the opening-closing process can be generated in computationally tractable simulation times.
| S-EPMC2480670 | biostudies-literature