Browse
Submit Data
Databases
API
Help

Dataset Information

28 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Microsecond-timescale MD simulation of EGFR minor mutation predicts the structural flexibility of EGFR kinase core that reflects EGFR inhibitor sensitivity.

ABSTRACT: Approximately 15-30% of patients with lung cancer harbor mutations in the EGFR gene. Major EGFR mutations (>90% of EGFR-mutated lung cancer) are highly sensitive to EGFR tyrosine kinase inhibitors (TKIs). Many uncommon EGFR mutations have been identified, but little is known regarding their characteristics, activation, and sensitivity to various EGFR-TKIs, including allosteric inhibitors. We encountered a case harboring an EGFR-L747P mutation, originally misdiagnosed with EGFR-del19 mutation using a routine diagnostic EGFR mutation test, which was resistant to EGFR-TKI gefitinib. Using this minor mutation and common EGFR-activating mutations, we performed the binding free energy calculations and microsecond-timescale molecular dynamic (MD) simulations, revealing that the L747P mutation considerably stabilizes the active conformation through a salt-bridge formation between K745 and E762. We further revealed why several EGFR inhibitors, including the allosteric inhibitor, were ineffective. Our computational structural analysis strategy would be beneficial for future drug development targeting the EGFR minor mutations.

SUBMITTER: Yoshizawa T

PROVIDER: S-EPMC8052404 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Flexibility of the bacterial chaperone trigger factor in microsecond-timescale molecular dynamics simulations.

Project description:The bacterial chaperone trigger factor (TF) is the first chaperone to be encountered by a nascent protein chain as it emerges from the ribosome exit tunnel. Experimental results suggest that TF possesses considerable conformational flexibility, and in an attempt to provide an atomic-level view of this flexibility, we have performed independent 1.5-μs molecular dynamics simulations of TF in explicit solvent using two different simulation force fields (OPLS-AA/L and AMBER ff99SB-ILDN). Both simulations indicate that TF possesses tremendous flexibility, with huge excursions from the crystallographic conformation caused by reorientations of the protein's constituent domains; both simulations also predict the formation of extensive contacts between TF's PPIase domain and the Arm 1 domain that is involved in nascent-chain binding. In the OPLS simulation, however, TF rapidly settles into a very compact conformation that persists for at least 1 μs, whereas in the AMBER simulation, it remains highly dynamic; additional simulations in which the two force fields were swapped suggest that these differences are at least partly attributable to sampling issues. The simulation results provide potential rationalizations of a number of experimental observations regarding TF's conformational behavior and have implications for using simulations to model TF's function on translating ribosomes.

| S-EPMC3736688 | biostudies-literature

Multiscale Simulation Reveals Passive Proton Transport Through SERCA on the Microsecond Timescale.

Project description:The sarcoplasmic reticulum Ca2+-ATPase (SERCA) transports two Ca2+ ions from the cytoplasm to the reticulum lumen at the expense of ATP hydrolysis. In addition to transporting Ca2+, SERCA facilitates bidirectional proton transport across the sarcoplasmic reticulum to maintain the charge balance of the transport sites and to balance the charge deficit generated by the exchange of Ca2+. Previous studies have shown the existence of a transient water-filled pore in SERCA that connects the Ca2+ binding sites with the lumen, but the capacity of this pathway to sustain passive proton transport has remained unknown. In this study, we used the multiscale reactive molecular dynamics method and free energy sampling to quantify the free energy profile and timescale of the proton transport across this pathway while also explicitly accounting for the dynamically coupled hydration changes of the pore. We find that proton transport from the central binding site to the lumen has a microsecond timescale, revealing a novel passive cytoplasm-to-lumen proton flow beside the well-known inverse proton countertransport occurring in active Ca2+ transport. We propose that this proton transport mechanism is operational and serves as a functional conduit for passive proton transport across the sarcoplasmic reticulum.

| S-EPMC7474205 | biostudies-literature

Microsecond timescale MD simulations at the transition state of <i>Pm</i>HMGR predict remote allosteric residues.

Project description:Understanding the mechanisms of enzymatic catalysis requires a detailed understanding of the complex interplay of structure and dynamics of large systems that is a challenge for both experimental and computational approaches. More importantly, the computational demands of QM/MM simulations mean that the dynamics of the reaction can only be considered on a timescale of nanoseconds even though the conformational changes needed to reach the catalytically active state happen on a much slower timescale. Here we demonstrate an alternative approach that uses transition state force fields (TSFFs) derived by the quantum-guided molecular mechanics (Q2MM) method that provides a consistent treatment of the entire system at the classical molecular mechanics level and allows simulations at the microsecond timescale. Application of this approach to the second hydride transfer transition state of HMG-CoA reductase from Pseudomonas mevalonii (PmHMGR) identified three remote residues, R396, E399 and L407, (15-27 Å away from the active site) that have a remote dynamic effect on enzyme activity. The predictions were subsequently validated experimentally via site-directed mutagenesis. These results show that microsecond timescale MD simulations of transition states are possible and can predict rather than just rationalize remote allosteric residues.

| S-EPMC8115266 | biostudies-literature

Sequence-specific backbone resonance assignments and microsecond timescale molecular dynamics simulation of human eosinophil-derived neurotoxin.

Project description:Eight active canonical members of the pancreatic-like ribonuclease A (RNase A) superfamily have been identified in human. All structural homologs share similar RNA-degrading functions, while also cumulating other various biological activities in different tissues. The functional homologs eosinophil-derived neurotoxin (EDN, or RNase 2) and eosinophil cationic protein (ECP, or RNase 3) are known to be expressed and secreted by eosinophils in response to infection, and have thus been postulated to play an important role in host defense and inflammatory response. We recently initiated the biophysical and dynamical investigation of several vertebrate RNase homologs and observed that clustering residue dynamics appear to be linked with the phylogeny and biological specificity of several members. Here we report the 1H, 13C and 15N backbone resonance assignments of human EDN (RNase 2) and its molecular dynamics simulation on the microsecond timescale, providing means to pursue this comparative atomic-scale functional and dynamical analysis by NMR and computation over multiple time frames.

| S-EPMC5589483 | biostudies-literature

Recovery of the poisoned topoisomerase II for DNA religation: coordinated motion of the cleavage core revealed with the microsecond atomistic simulation.

Project description:Type II topoisomerases resolve topological problems of DNA double helices by passing one duplex through the reversible double-stranded break they generated on another duplex. Despite the wealth of information in the cleaving operation, molecular understanding of the enzymatic DNA ligation remains elusive. Topoisomerase poisons are widely used in anti-cancer and anti-bacterial therapy and have been employed to entrap the intermediates of topoisomerase IIβ with religatable DNA substrate. We removed drug molecules from the structure and conducted molecular dynamics simulations to investigate the enzyme-mediated DNA religation. The drug-unbound intermediate displayed transitions toward the resealing-compliant configuration: closing distance between the cleaved DNA termini, B-to-A transformation of the double helix, and restoration of the metal-binding motif. By mapping the contact configurations and the correlated motions between enzyme and DNA, we identified the indispensable role of the linker preceding winged helix domain (WHD) in coordinating the movements of TOPRIM, the nucleotide-binding motifs, and the bound DNA substrate during gate closure. We observed a nearly vectorial transition in the recovery of the enzyme and identified the previously uncharacterized roles of Asn508 and Arg677 in DNA rejoining. Our findings delineate the dynamic mechanism of the DNA religation conducted by type II topoisomerases.

| S-EPMC4538842 | biostudies-literature

Microsecond-timescale simulations suggest 5-HT-mediated preactivation of the 5-HT<sub>3A</sub> serotonin receptor.

Project description:Aided by efforts to improve their speed and efficiency, molecular dynamics (MD) simulations provide an increasingly powerful tool to study the structure-function relationship of pentameric ligand-gated ion channels (pLGICs). However, accurate reporting of the channel state and observation of allosteric regulation by agonist binding with MD remains difficult due to the timescales necessary to equilibrate pLGICs from their artificial and crystalized conformation to a more native, membrane-bound conformation in silico. Here, we perform multiple all-atom MD simulations of the homomeric 5-hydroxytryptamine 3A (5-HT3A) serotonin receptor for 15 to 20 μs to demonstrate that such timescales are critical to observe the equilibration of a pLGIC from its crystalized conformation to a membrane-bound conformation. These timescales, which are an order of magnitude longer than any previous simulation of 5-HT3A, allow us to observe the dynamic binding and unbinding of 5-hydroxytryptamine (5-HT) (i.e., serotonin) to the binding pocket located on the extracellular domain (ECD) and allosteric regulation of the transmembrane domain (TMD) from synergistic 5-HT binding. While these timescales are not long enough to observe complete activation of 5-HT3A, the allosteric regulation of ion gating elements by 5-HT binding is indicative of a preactive state, which provides insight into molecular mechanisms that regulate channel activation from a resting state. This mechanistic insight, enabled by microsecond-timescale MD simulations, will allow a careful examination of the regulation of pLGICs at a molecular level, expanding our understanding of their function and elucidating key structural motifs that can be targeted for therapeutic regulation.

| S-EPMC6955379 | biostudies-literature

Starting-structure dependence of nanosecond timescale intersubstate transitions and reproducibility of MD-derived order parameters.

Project description:Factors affecting the accuracy of molecular dynamics (MD) simulations are investigated by comparing generalized order parameters for backbone NH vectors of the B3 immunoglobulin-binding domain of streptococcal protein G (GB3) derived from simulations with values obtained from NMR spin relaxation (Yao L, Grishaev A, Cornilescu G, Bax A, J Am Chem Soc 2010;132:4295-4309.). Choices for many parameters of the simulations, such as buffer volume, water model, or salt concentration, have only minor influences on the resulting order parameters. In contrast, seemingly minor conformational differences in starting structures, such as orientations of sidechain hydroxyl groups, resulting from applying different protonation algorithms to the same structure, have major effects on backbone dynamics. Some, but not all, of these effects are mitigated by increased sampling in simulations. Most discrepancies between simulated and experimental results occur for residues located at the ends of secondary structures and involve large amplitude nanosecond timescale transitions between distinct conformational substates. These transitions result in autocorrelation functions for bond vector reorientation that do not converge when calculated over individual simulation blocks, typically of length similar to the overall rotational diffusion time. A test for convergence before averaging the order parameters from different blocks results in better agreement between order parameters calculated from different sets of simulations and with NMR-derived order parameters. Thus, MD-derived order parameters are more strongly affected by transitions between conformational substates than by fluctuations within individual substates themselves, while conformational differences in the starting structures affect the frequency and scale of such transitions.

| S-EPMC3557728 | biostudies-literature

Reorientation and dimerization of the membrane-bound antimicrobial peptide PGLa from microsecond all-atom MD simulations.

Project description:The membrane-active antimicrobial peptide PGLa from Xenopus laevis is known from solid-state (2)H-, (15)N-, and (19)F-NMR spectroscopy to occupy two distinct ?-helical surface adsorbed states in membranes: a surface-bound S-state with a tilt angle of ~95° at low peptide/lipid molar ratio (P/L = 1:200), and an obliquely tilted T-state with a tilt angle of 127° at higher peptide concentration (P/L = 1:50). Using a rapid molecular-dynamics insertion protocol in combination with microsecond-scale simulation, we have characterized the structure of both states in detail. As expected, the amphiphilic peptide resides horizontally on the membrane surface in a monomeric form at a low P/L, whereas the T-state is seen in the simulations to be a symmetric antiparallel dimer, with close contacts between small glycine and alanine residues at the interface. The computed tilt angles and azimuthal rotations, as well as the quadrupolar splittings predicted from the simulations agree with the experimental NMR data. The simulations reveal many structural details previously inaccessible, such as the immersion depth of the peptide in the membrane and the packing of the dimerization interface. The study highlights the ability and limitations of current state-of-the-art multimicrosecond all-atom simulations of membrane-active peptides to complement experimental data from solid-state NMR.

| S-EPMC3414899 | biostudies-other

Concerted interconversion between ionic lock substates of the beta(2) adrenergic receptor revealed by microsecond timescale molecular dynamics.

Project description:The recently solved crystallographic structures for the A(2A) adenosine receptor and the beta(1) and beta(2) adrenergic receptors have shown important differences between members of the class-A G-protein-coupled receptors and their archetypal model, rhodopsin, such as the apparent breaking of the ionic lock that stabilizes the inactive structure. Here, we characterize a 1.02 mus all-atom simulation of an apo-beta(2) adrenergic receptor that is missing the third intracellular loop to better understand the inactive structure. Although we find that the structure is remarkably rigid, there is a rapid influx of water into the core of the protein, as well as a slight expansion of the molecule relative to the crystal structure. In contrast to the x-ray crystal structures, the ionic lock rapidly reforms, although we see an activation-precursor-like event wherein the ionic lock opens for approximately 200 ns, accompanied by movements in the transmembrane helices associated with activation. When the lock reforms, we see the structure return to its inactive conformation. We also find that the ionic lock exists in three states: closed (or locked), semi-open with a bridging water molecule, and open. The interconversion of these states involves the concerted motion of the entire protein. We characterize these states and the concerted motion underlying their interconversion. These findings may help elucidate the connection between key local events and the associated global structural changes during activation.

| S-EPMC2800975 | biostudies-literature

Coupling between Histone Conformations and DNA Geometry in Nucleosomes on a Microsecond Timescale: Atomistic Insights into Nucleosome Functions.

Project description:An octamer of histone proteins wraps about 200bp of DNA into two superhelical turns to form nucleosomes found in chromatin. Although the static structure of the nucleosomal core particle has been solved, details of the dynamic interactions between histones and DNA remain elusive. We performed extensively long unconstrained, all-atom microsecond molecular dynamics simulations of nucleosomes including linker DNA segments and full-length histones in explicit solvent. For the first time, we were able to identify and characterize the rearrangements in nucleosomes on a microsecond timescale including the coupling between the conformation of the histone tails and the DNA geometry. We found that certain histone tail conformations promoted DNA bulging near its entry/exit sites, resulting in the formation of twist defects within the DNA. This led to a reorganization of histone-DNA interactions, suggestive of the formation of initial nucleosome sliding intermediates. We characterized the dynamics of the histone tails upon their condensation on the core and linker DNA and showed that tails may adopt conformationally constrained positions due to the insertion of "anchoring" lysines and arginines into the DNA minor grooves. Potentially, these phenomena affect the accessibility of post-translationally modified histone residues that serve as important sites for epigenetic marks (e.g., at H3K9, H3K27, H4K16), suggesting that interactions of the histone tails with the core and linker DNA modulate the processes of histone tail modifications and binding of the effector proteins. We discuss the implications of the observed results on the nucleosome function and compare our results to different experimental studies.

| S-EPMC4738025 | 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