Browse
Submit Data
Databases
API
Help

Dataset Information

26 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Microsecond timescale MD simulations at the transition state of PmHMGR predict remote allosteric residues.

ABSTRACT: 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.

SUBMITTER: Quinn TR

PROVIDER: S-EPMC8115266 | 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

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

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

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

Project description: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.

| S-EPMC8052404 | biostudies-literature

Allosteric inhibition of TEM-1 β lactamase: Microsecond molecular dynamics simulations provide mechanistic insights.

Project description:β-lactam antibiotics target DD-transpeptidases, enzymes that perform the last step of bacterial cell-wall synthesis. To block the antimicrobial activity of these antibiotics, bacteria have evolved lactamases that render them inert. Among these, TEM-1, a class A lactamase, has been extensively studied. In 2004, Horn et al. described a novel allosteric TEM-1 inhibitor, FTA, that binds distant from the TEM-1 orthosteric (penicillin-binding) pocket. TEM-1 has subsequently become a model for the study of allostery. In the present work, we perform molecular dynamics simulations of FTA-bound and FTA-absent TEM-1, totaling ~3 μS, that provide new insight into TEM-1 inhibition. In one of the simulations, bound FTA assumed a conformation different than that observed crystallographically. We provide evidence that the alternate pose is physiologically plausible and describe how it impacts our understanding of TEM-1 allostery.

| S-EPMC10044111 | biostudies-literature

Microsecond MD simulations of human CYP2D6 wild-type and five allelic variants reveal mechanistic insights on the function.

Project description:Characterization of cytochrome P450 2D6 (CYP2D6) and the impact of the major identified allelic variants on the activity of one of the most dominating drug-metabolising enzymes is essential to increase drug safety and avoid adverse reactions. Microsecond molecular dynamics simulations have been performed to capture the dynamic signatures of this complex enzyme and five allelic variants with diverse enzymatic activity. In addition to the apo simulations, three substrates (bufuralol, veliparib and tamoxifen) and two inhibitors (prinomastat and quinidine) were included to explore their influence on the structure and dynamical features of the enzyme. Our results indicate that the altered enzyme activity can be attributed to changes in the hydrogen bonding network within the active site, and local structural differences in flexibility, position and shape of the binding pocket. In particular, the increased (CYP2D6*53) or the decreased (CYP2D6*17) activity seems to be related to a change in dynamics of mainly the BC loop due to a modified hydrogen bonding network around this region. In addition, the smallest active site volume was found for CYP2D6*4 (no activity). CYP2D6*2 (normal activity) showed no major differences in dynamic behaviour compared to the wild-type.

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

Small molecule interactions with the SARS-CoV-2 main protease: In silico all-atom microsecond MD simulations, PELE Monte Carlo simulations, and determination of in vitro activity inhibition.

Project description:The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the ongoing COVID-19 pandemic. With some notable exceptions, safe and effective vaccines, which are now being widely distributed globally, have largely begun to stabilise the situation. However, emerging variants of concern and vaccine hesitancy are apparent obstacles to eradication. Therefore, the need for the development of potent antivirals is still of importance. In this context, the SARS-CoV-2 main protease (Mpro) is a critical target and numerous clinical trials, predominantly in the private domain, are currently in progress. Here, our aim was to extend our previous studies, with hypericin and cyanidin-3-O-glucoside, as potential inhibitors of the SARS-CoV-2 Mpro. Firstly, we performed all-atom microsecond molecular dynamics simulations, which highlight the stability of the ligands in the Mpro active site over the duration of the trajectories. We also invoked PELE Monte Carlo simulations which indicate that both hypericin and cyanidin-3-O-glucoside preferentially interact with the Mpro active site and known allosteric sites. For further validation, we performed an in vitro enzymatic activity assay that demonstrated that hypericin and cyanidin-3-O-glucoside inhibit Mpro activity in a dose-dependent manner at biologically relevant (μM) concentrations. However, both ligands are much less potent than the well-known covalent antiviral GC376, which was used as a positive control in our experiments. Nevertheless, the biologically relevant activity of hypericin and cyanidin-3-O-glucoside is encouraging. In particular, a synthetic version of hypericin has FDA orphan drug designation, which could simplify potential clinical evaluation in the context of COVID-19.

| S-EPMC8504156 | biostudies-literature

Microsecond molecular dynamics simulations of lipid mixing.

Project description:Molecular dynamics (MD) simulations of membranes are often hindered by the slow lateral diffusion of lipids and the limited time scale of MD. In order to study the dynamics of mixing and characterize the lateral distribution of lipids in converged mixtures, we report microsecond-long all-atom MD simulations performed on the special-purpose machine Anton. Two types of mixed bilayers, POPE:POPG (3:1) and POPC:cholesterol (2:1), as well as a pure POPC bilayer, were each simulated for up to 2 μs. These simulations show that POPE:POPG and POPC:cholesterol are each fully miscible at the simulated conditions, with the final states of the mixed bilayers similar to a random mixture. By simulating three POPE:POPG bilayers at different NaCl concentrations (0, 0.15, and 1 M), we also examined the effect of salt concentration on lipid mixing. While an increase in NaCl concentration is shown to affect the area per lipid, tail order, and lipid lateral diffusion, the final states of mixing remain unaltered, which is explained by the largely uniform increase in Na(+) ions around POPE and POPG. Direct measurement of water permeation reveals that the POPE:POPG bilayer with 1 M NaCl has reduced water permeability compared with those at zero or low salt concentration. Our calculations provide a benchmark to estimate the convergence time scale of all-atom MD simulations of lipid mixing. Additionally, equilibrated structures of POPE:POPG and POPC:cholesterol, which are frequently used to mimic bacterial and mammalian membranes, respectively, can be used as starting points of simulations involving these membranes.

| S-EPMC4196744 | biostudies-literature

Novel covalent and non-covalent complex-based pharmacophore models of SARS-CoV-2 main protease (M<sup>pro</sup>) elucidated by microsecond MD simulations.

Project description:As the world enters its second year of the pandemic caused by SARS-CoV-2, intense efforts have been directed to develop an effective diagnosis, prevention, and treatment strategies. One promising drug target to design COVID-19 treatments is the SARS-CoV-2 Mpro. To date, a comparative understanding of Mpro dynamic stereoelectronic interactions with either covalent or non-covalent inhibitors (depending on their interaction with a pocket called S1' or oxyanion hole) has not been still achieved. In this study, we seek to fill this knowledge gap using a cascade in silico protocol of docking, molecular dynamics simulations, and MM/PBSA in order to elucidate pharmacophore models for both types of inhibitors. After docking and MD analysis, a set of complex-based pharmacophore models was elucidated for covalent and non-covalent categories making use of the residue bonding point feature. The highest ranked models exhibited ROC-AUC values of 0.93 and 0.73, respectively for each category. Interestingly, we observed that the active site region of Mpro protein-ligand complex undergoes large conformational changes, especially within the S2 and S4 subsites. The results reported in this article may be helpful in virtual screening (VS) campaigns to guide the design and discovery of novel small-molecule therapeutic agents against SARS-CoV-2 Mpro protein.

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