Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

What can we learn from molecular dynamics simulations for GPCR drug design?

ABSTRACT: Recent years have seen a tremendous progress in the elucidation of experimental structural information for G-protein coupled receptors (GPCRs). Although for the vast majority of pharmaceutically relevant GPCRs structural information is still accessible only by homology models the steadily increasing amount of structural information fosters the application of structure-based drug design tools for this important class of drug targets. In this article we focus on the application of molecular dynamics (MD) simulations in GPCR drug discovery programs. Typical application scenarios of MD simulations and their scope and limitations will be described on the basis of two selected case studies, namely the binding of small molecule antagonists to the human CC chemokine receptor 3 (CCR3) and a detailed investigation of the interplay between receptor dynamics and solvation for the binding of small molecules to the human muscarinic acetylcholine receptor 3 (hM3R).

SUBMITTER: Tautermann CS

PROVIDER: S-EPMC4334948 | biostudies-literature | 2015

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

What can we learn from molecular dynamics simulations for GPCR drug design?

Tautermann Christofer S CS Seeliger Daniel D Kriegl Jan M JM

Computational and structural biotechnology journal 20141210

Recent years have seen a tremendous progress in the elucidation of experimental structural information for G-protein coupled receptors (GPCRs). Although for the vast majority of pharmaceutically relevant GPCRs structural information is still accessible only by homology models the steadily increasing amount of structural information fosters the application of structure-based drug design tools for this important class of drug targets. In this article we focus on the application of molecular dynami ...[more]

PMID: 25709761

Similar Datasets

Molecular dynamics simulations and drug discovery.

Project description:This review discusses the many roles atomistic computer simulations of macromolecular (for example, protein) receptors and their associated small-molecule ligands can play in drug discovery, including the identification of cryptic or allosteric binding sites, the enhancement of traditional virtual-screening methodologies, and the direct prediction of small-molecule binding energies. The limitations of current simulation methodologies, including the high computational costs and approximations of molecular forces required, are also discussed. With constant improvements in both computer power and algorithm design, the future of computer-aided drug design is promising; molecular dynamics simulations are likely to play an increasingly important role.

| S-EPMC3203851 | biostudies-literature

Structure-based simulations reveal concerted dynamics of GPCR activation.

Project description:G protein-coupled receptors (GPCRs) are a vital class of proteins that transduce biological signals across the cell membrane. However, their allosteric activation mechanism is not fully understood; crystal structures of active and inactive receptors have been reported, but the functional pathway between these two states remains elusive. Here, we use structure-based (G?-like) models to simulate activation of two GPCRs, rhodopsin and the ?? adrenergic receptor (??AR). We used data-derived reaction coordinates that capture the activation mechanism for both proteins, showing that activation proceeds through quantitatively different paths in the two systems. Both reaction coordinates are determined from the dominant concerted motions in the simulations so the technique is broadly applicable. There were two surprising results. First, the main structural changes in the simulations were distributed throughout the transmembrane bundle, and not localized to the obvious areas of interest, such as the intracellular portion of Helix 6. Second, the activation (and deactivation) paths were distinctly nonmonotonic, populating states that were not simply interpolations between the inactive and active structures. These transitions also suggest a functional explanation for ??AR's basal activity: it can proceed through a more broadly defined path during the observed transitions.

| S-EPMC4177312 | biostudies-other

Neural relational inference to learn long-range allosteric interactions in proteins from molecular dynamics simulations.

Project description:Protein allostery is a biological process facilitated by spatially long-range intra-protein communication, whereby ligand binding or amino acid change at a distant site affects the active site remotely. Molecular dynamics (MD) simulation provides a powerful computational approach to probe the allosteric effect. However, current MD simulations cannot reach the time scales of whole allosteric processes. The advent of deep learning made it possible to evaluate both spatially short and long-range communications for understanding allostery. For this purpose, we applied a neural relational inference model based on a graph neural network, which adopts an encoder-decoder architecture to simultaneously infer latent interactions for probing protein allosteric processes as dynamic networks of interacting residues. From the MD trajectories, this model successfully learned the long-range interactions and pathways that can mediate the allosteric communications between distant sites in the Pin1, SOD1, and MEK1 systems. Furthermore, the model can discover allostery-related interactions earlier in the MD simulation trajectories and predict relative free energy changes upon mutations more accurately than other methods.

| S-EPMC8964751 | biostudies-literature

Can molecular dynamics simulations improve the structural accuracy and virtual screening performance of GPCR models?

Project description:The determination of G protein-coupled receptor (GPCR) structures at atomic resolution has improved understanding of cellular signaling and will accelerate the development of new drug candidates. However, experimental structures still remain unavailable for a majority of the GPCR family. GPCR structures and their interactions with ligands can also be modelled computationally, but such predictions have limited accuracy. In this work, we explored if molecular dynamics (MD) simulations could be used to refine the accuracy of in silico models of receptor-ligand complexes that were submitted to a community-wide assessment of GPCR structure prediction (GPCR Dock). Two simulation protocols were used to refine 30 models of the D3 dopamine receptor (D3R) in complex with an antagonist. Close to 60 μs of simulation time was generated and the resulting MD refined models were compared to a D3R crystal structure. In the MD simulations, the receptor models generally drifted further away from the crystal structure conformation. However, MD refinement was able to improve the accuracy of the ligand binding mode. The best refinement protocol improved agreement with the experimentally observed ligand binding mode for a majority of the models. Receptor structures with improved virtual screening performance, which was assessed by molecular docking of ligands and decoys, could also be identified among the MD refined models. Application of weak restraints to the transmembrane helixes in the MD simulations further improved predictions of the ligand binding mode and second extracellular loop. These results provide guidelines for application of MD refinement in prediction of GPCR-ligand complexes and directions for further method development.

| S-EPMC8186765 | biostudies-literature

Understanding Drug Skin Permeation Enhancers Using Molecular Dynamics Simulations.

Project description:Our skin constitutes an effective permeability barrier that protects the body from exogenous substances but concomitantly severely limits the number of pharmaceutical drugs that can be delivered transdermally. In topical formulation design, chemical permeation enhancers (PEs) are used to increase drug skin permeability. In vitro skin permeability experiments can measure net effects of PEs on transdermal drug transport, but they cannot explain the molecular mechanisms of interactions between drugs, permeation enhancers, and skin structure, which limits the possibility to rationally design better new drug formulations. Here we investigate the effect of the PEs water, lauric acid, geraniol, stearic acid, thymol, ethanol, oleic acid, and eucalyptol on the transdermal transport of metronidazole, caffeine, and naproxen. We use atomistic molecular dynamics (MD) simulations in combination with developed molecular models to calculate the free energy difference between 11 PE-containing formulations and the skin's barrier structure. We then utilize the results to calculate the final concentration of PEs in skin. We obtain an RMSE of 0.58 log units for calculated partition coefficients from water into the barrier structure. We then use the modified PE-containing barrier structure to calculate the PEs' permeability enhancement ratios (ERs) on transdermal metronidazole, caffeine, and naproxen transport and compare with the results obtained from in vitro experiments. We show that MD simulations are able to reproduce rankings based on ERs. However, strict quantitative correlation with experimental data needs further refinement, which is complicated by significant deviations between different measurements. Finally, we propose a model for how to use calculations of the potential of mean force of drugs across the skin's barrier structure in a topical formulation design.

| S-EPMC10428223 | biostudies-literature

Statistics for the analysis of molecular dynamics simulations: providing P values for agonist-dependent GPCR activation.

Project description:Molecular dynamics (MD) is the common computational technique for assessing efficacy of GPCR-bound ligands. Agonist efficacy measures the capability of the ligand-bound receptor of reaching the active state in comparison with the free receptor. In this respect, agonists, neutral antagonists and inverse agonists can be considered. A collection of MD simulations of both the ligand-bound and the free receptor are needed to provide reliable conclusions. Variability in the trajectories needs quantification and proper statistical tools for meaningful and non-subjective conclusions. Multiple-factor (time, ligand, lipid) ANOVA with repeated measurements on the time factor is proposed as a suitable statistical method for the analysis of agonist-dependent GPCR activation MD simulations. Inclusion of time factor in the ANOVA model is consistent with the time-dependent nature of MD. Ligand and lipid factors measure agonist and lipid influence on receptor activation. Previously reported MD simulations of adenosine A2a receptor (A2aR) are reanalyzed with this statistical method. TM6-TM3 and TM7-TM3 distances are selected as dependent variables in the ANOVA model. The ligand factor includes the presence or absence of adenosine whereas the lipid factor considers DOPC or DOPG lipids. Statistical analysis of MD simulations shows the efficacy of adenosine and the effect of the membrane lipid composition. Subsequent application of the statistical methodology to NECA A2aR agonist, with resulting P values in consistency with its pharmacological profile, suggests that the method is useful for ligand comparison and potentially for dynamic structure-based virtual screening.

| S-EPMC7672096 | biostudies-literature

GPCR Conformations: Implications for Rational Drug Design

Project description:G protein-coupled receptors (GPCRs) comprise a large class of transmembrane proteins that play critical roles in both normal physiology and pathophysiology. These critical roles offer targets for therapeutic intervention, as exemplified by the substantial fraction of current pharmaceutical agents that target members of this family. Tremendous contributions to our understanding of GPCR structure and dynamics have come from both indirect and direct structural characterization techniques. Key features of GPCR conformations derived from both types of characterization techniques are reviewed.

| S-EPMC4052540 | biostudies-literature

Improving the Efficiency of Ligand-Binding Protein Design with Molecular Dynamics Simulations.

Project description:Custom-designed ligand-binding proteins represent a promising class of macromolecules with exciting applications toward the design of new enzymes or the engineering of antibodies and small-molecule recruited proteins for therapeutic interventions. However, several challenges remain in designing a protein sequence such that the binding site organization results in high affinity interaction with a bound ligand. Here, we study the dynamics of explicitly solvated designed proteins through all-atom molecular dynamics (MD) simulations to gain insight into the causes that lead to the low affinity or instability of most of these designs, despite the prediction of their success by the computational design methodology. Simulations ranging from 500 to 1000 ns per replicate were conducted on 37 designed protein variants encompassing two distinct folds and a range of ligand affinities, resulting in more than 180 ?s of combined sampling. The simulations provide retrospective insights into the properties affecting ligand affinity that can prove useful in guiding further steps of design optimization. Features indicate that entropic components are particularly important for affinity, which are not easily incorporated in the empirical models often used in design protocols. Additionally, we demonstrate that the application of machine learning approaches built upon the output from the simulations can help discriminate between successful and failed binders, such that MD could act as a screening step in protein design, resulting in a more efficient process.

| S-EPMC7532806 | biostudies-literature

Fragment-based design of selective GPCR ligands guided by free energy simulations.

Project description:Fragment-based drug discovery relies on successful optimization of weakly binding ligands for affinity and selectivity. Herein, we explored strategies for structure-based evolution of fragments binding to a G protein-coupled receptor. Molecular dynamics simulations combined with rigorous free energy calculations guided synthesis of nanomolar ligands with up to >1000-fold improvements of binding affinity and close to 40-fold subtype selectivity.

| S-EPMC8603191 | biostudies-literature

Flexibility in a drug transport accessory protein: molecular dynamics simulations of MexA.

Project description:Drug resistance in gram-negative bacteria may be conferred via efflux through a tripartite complex of an inner membrane pump, an outer membrane pore, and a periplasmic adaptor protein. These are AcrB, TolC, and AcrA, respectively, in Escherichia coli. In Pseudomonas aerugonisa, their homologs are MexB, OprM, and MexA. Defining the interdomain dynamics of the adaptor protein is essential to understanding the mechanism of complex formation. Extended (25 ns) molecular dynamics simulations of MexA have been performed to determine such interdomain dynamics. Analysis of conformational drift demonstrates substantial motions of the three domains of MexA relative to one another. Principal components analysis reveals a hinge-bending motion and rotation of the alpha-helical hairpin relative to the other domains to be the two dominant motions. These two motions provide an element of considerable flexibility which is likely to be exploited in the adaptor function of MexA.

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