Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Steered Molecular Dynamics Simulations Study on FABP4 Inhibitors.

ABSTRACT: Ordinary small molecule de novo drug design is time-consuming and expensive. Recently, computational tools were employed and proved their efficacy in accelerating the overall drug design process. Molecular dynamics (MD) simulations and a derivative of MD, steered molecular dynamics (SMD), turned out to be promising rational drug design tools. In this paper, we report the first application of SMD to evaluate the binding properties of small molecules toward FABP4, considering our recent interest in inhibiting fatty acid binding protein 4 (FABP4). FABP4 inhibitors (FABP4is) are small molecules of therapeutic interest, and ongoing clinical studies indicate that they are promising for treating cancer and other diseases such as metabolic syndrome and diabetes.

SUBMITTER: Tomarchio R

PROVIDER: S-EPMC10058326 | biostudies-literature | 2023 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Steered Molecular Dynamics Simulations Study on FABP4 Inhibitors.

Tomarchio Rosario R Patamia Vincenzo V Zagni Chiara C Crocetti Letizia L Cilibrizzi Agostino A Floresta Giuseppe G Rescifina Antonio A

Molecules (Basel, Switzerland) 20230317 6

Ordinary small molecule de novo drug design is time-consuming and expensive. Recently, computational tools were employed and proved their efficacy in accelerating the overall drug design process. Molecular dynamics (MD) simulations and a derivative of MD, steered molecular dynamics (SMD), turned out to be promising rational drug design tools. In this paper, we report the first application of SMD to evaluate the binding properties of small molecules toward FABP4, considering our recent interest i ...[more]

PMID: 36985701

Similar Datasets

Targeting N-Terminal Human Maltase-Glucoamylase to Unravel Possible Inhibitors Using Molecular Docking, Molecular Dynamics Simulations, and Adaptive Steered Molecular Dynamics Simulations.

Project description:There are multiple drugs for the treatment of type 2 diabetes, including traditional sulfonylureas biguanides, glinides, thiazolidinediones, α-glucosidase inhibitors, glucagon-like peptide-1 (GLP-1) receptor agonists, dipeptidyl peptidase IV (DPP-4) inhibitors, and sodium-glucose cotransporter 2 (SGLT2) inhibitors. α-Glucosidase inhibitors have been used to control postprandial glucose levels caused by type 2 diabetes since 1990. α-Glucosidases are rather crucial in the human metabolic system and are principally found in families 13 and 31. Maltase-glucoamylase (MGAM) belongs to glycoside hydrolase family 31. The main function of MGAM is to digest terminal starch products left after the enzymatic action of α-amylase; hence, MGAM becomes an efficient drug target for insulin resistance. In order to explore the conformational changes in the active pocket and unbinding pathway for NtMGAM, molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed for two NtMGAM-inhibitor [de-O-sulfonated kotalanol (DSK) and acarbose] complexes. MD simulations indicated that DSK bound to NtMGAM may influence two domains (inserted loop 1 and inserted loop 2) by interfering with the spiralization of residue 497-499. The flexibility of inserted loop 1 and inserted loop 2 can influence the volume of the active pocket of NtMGAM, which can affect the binding progress for DSK to NtMGAM. ASMD simulations showed that compared to acarbose, DSK escaped from NtMGAM easily with lower energy. Asp542 is an important residue on the bottleneck of the active pocket of NtMGAM and could generate hydrogen bonds with DSK continuously. Our theoretical results may provide some useful clues for designing new α-glucosidase inhibitors to treat type 2 diabetes.

| S-EPMC8435576 | biostudies-literature

Steered molecular dynamics simulations reveal critical residues for (un)binding of substrates, inhibitors and a product to the malarial M1 aminopeptidase.

Project description:Malaria is a life-threatening disease spread by mosquitoes. Plasmodium falciparum M1 alanyl aminopeptidase (PfM1-AAP) is a promising target for the treatment of malaria. The recently solved crystal structures of PfM1-AAP revealed that the buried active site can be accessed through two channel openings: a short N-terminal channel with the length of 8 Å and a long C-terminal channel with the length of 30 Å. It is unclear, however, how substrates and inhibitors migrate to the active site and a product of cleavage leaves. Here, we study the molecular mechanism of substrate and inhibitor migration to the active site and the product release using steered molecular dynamics simulations. We identified a stepwise passage of substrates and inhibitors in the C-terminal channel of PfM1-AAP, involving (I) ligand recognition at the opening of the channel, (II) ionic translation to the 'water reservoir', (III) ligand reorientation in the 'water reservoir' and (IV) passage in a suitable conformation into the active site. Endorsed by enzymatic analysis of functional recombinant PfM1-AAP and mutagenesis studies, our novel ligand-residue binding network analysis has identified the functional residues controlling ligand migration within the C-terminal channel of PfM1-AAP. Furthermore, from unbinding simulations of the Arg product we propose a charge repulsion as the driving force to expel the product out from the N-terminal channel of PfM1-AAP. Our work paves the way towards the design of a novel class of PfM1-AAP inhibitors based on preventing substrate entry to the active site.

| S-EPMC6239339 | biostudies-literature

Elucidation of the glucose transport pathway in glucose transporter 4 via steered molecular dynamics simulations.

Project description:BackgroundGLUT4 is a predominant insulin regulated glucose transporter expressed in major glucose disposal tissues such as adipocytes and muscles. Under the unstimulated state, GLUT4 resides within intracellular vesicles. Various stimuli such as insulin translocate this protein to the plasma membrane for glucose transport. In the absence of a crystal structure for GLUT4, very little is known about the mechanism of glucose transport by this protein. Earlier we proposed a homology model for GLUT4 and performed a conventional molecular dynamics study revealing the conformational rearrangements during glucose and ATP binding. However, this study could not explain the transport of glucose through the permeation tunnel.Methodology/principal findingsTo elucidate the molecular mechanism of glucose transport and its energetic, a steered molecular dynamics study (SMD) was used. Glucose was pulled from the extracellular end of GLUT4 to the cytoplasm along the pathway using constant velocity pulling method. We identified several key residues within the tunnel that interact directly with either the backbone ring or the hydroxyl groups of glucose. A rotation of glucose molecule was seen near the sugar binding site facilitating the sugar recognition process at the QLS binding site.Conclusions/significanceThis study proposes a possible glucose transport pathway and aids the identification of several residues that make direct interactions with glucose during glucose transport. Mutational studies are required to further validate the observation made in this study.

| S-EPMC3192114 | biostudies-literature

Computational Study for the Unbinding Routes of β-N-Acetyl-d-Hexosaminidase Inhibitor: Insight from Steered Molecular Dynamics Simulations.

Project description:β-N-Acetyl-d-hexosaminidase from Ostrinia furnacalis (OfHex1) is a new target for the design of insecticides. Although some of its inhibitors have been found, there is still no commercial drug available at present. The residence time of the ligand may be important for its pharmacodynamic effect. However, the unbinding routes of ligands from OfHex1 still remain largely unexplored. In the present study, we first simulated the six dissociation routes of N,N,N-trimethyl-d-glucosamine-chitotriomycin (TMG-chitotriomycin, a highly selective inhibitor of OfHex1) from the active pocket of OfHex1 by steered molecular dynamics simulations. By comparing the potential of mean forces (PMFs) of six routes, Route 1 was considered as the most possible route with the lowest energy barrier. Furthermore, the structures of six different states for Route 1 were snapshotted, and the key amino acid residues affecting the dissociated time were analyzed in the unbinding pathway. Moreover, we also analyzed the "open⁻close" mechanism of Glu368 and Trp448 and found that their conformational changes directly affected the dissociation of TMG-chitotriomycin. Our findings would be helpful to understanding and identifying novel inhibitors against OfHex1 from virtual screening or lead-optimization.

| S-EPMC6471479 | biostudies-literature

Steered molecular dynamics simulations on the "tail helix latch" hypothesis in the gelsolin activation process.

Project description:The molecular basis of the "tail helix latch" hypothesis in the gelsolin activation process has been studied by using the steered molecular dynamics simulations. In the present nanosecond scale simulations, the tail helix of gelsolin was pulled away from the S2 binding surface, and the required forces were calculated, from which the properties of binding between the tail helix and S2 domain and their dynamic unbinding processes were obtained. The force profile provides a detailed rupture mechanism that includes six major unbinding steps. In particular, the hydrogen bonds formed between Arg-207 and Asp-744 and between Arg-221 and Leu-753 are of the most important interaction pairs. The two hydrogen bond "clamps" stabilize the complex. The subsequent simulation on Arg-207-Ala (R207A) mutation of gelsolin indicated that this mutation facilitates the unbinding of the tail helix and that the contribution of the hydrogen bond between Arg-207 and Asp-744 to the binding is more than 50%, which offers a new clue for further mutagenesis study on the activation mechanism of gelsolin. Surrounding water molecules enhance the stability of the tail helix and facilitate the rupture process. Additionally, temperature also has a significant effect on the conformation of the arginine and arginine-related interactions, which revealed the molecular basis of the temperature dependence in gelsolin activation.

| S-EPMC1302184 | biostudies-other

Mechanically untying a protein slipknot: multiple pathways revealed by force spectroscopy and steered molecular dynamics simulations.

Project description:Protein structure is highly diverse when considering a wide range of protein types, helping to give rise to the multitude of functions that proteins perform. In particular, certain proteins are known to adopt a knotted or slipknotted fold. How such proteins undergo mechanical unfolding was investigated utilizing a combination of single molecule atomic force microscopy (AFM), protein engineering, and steered molecular dynamics (SMD) simulations to show the mechanical unfolding mechanism of the slipknotted protein AFV3-109. Our results reveal that the mechanical unfolding of AFV3-109 can proceed via multiple parallel unfolding pathways that all cause the protein slipknot to untie and the polypeptide chain to completely extend. These distinct unfolding pathways proceed via either a two- or three-state unfolding process involving the formation of a well-defined, stable intermediate state. SMD simulations predict the same contour length increments for different unfolding pathways as single molecule AFM results, thus providing a plausible molecular mechanism for the mechanical unfolding of AFV3-109. These SMD simulations also reveal that two-state unfolding is initiated from both the N- and C-termini, while three-state unfolding is initiated only from the C-terminus. In both pathways, the protein slipknot was untied during unfolding, and no tightened slipknot conformation was observed. Detailed analysis revealed that interactions between key structural elements lock the knotting loop in place, preventing it from shrinking and the formation of a tightened slipknot conformation. Our results demonstrate the bifurcation of the mechanical unfolding pathway of AFV3-109 and point to the generality of a kinetic partitioning mechanism for protein folding/unfolding.

| S-EPMC3663486 | biostudies-literature

Concerted dynamic motions of an FABP4 model and its ligands revealed by microsecond molecular dynamics simulations.

Project description:In this work, we investigate the dynamic motions of fatty acid binding protein 4 (FABP4) in the absence and presence of a ligand by explicitly solvated all-atom molecular dynamics simulations. The dynamics of one ligand-free FABP4 and four ligand-bound FABP4s is compared via multiple 1.2 μs simulations. In our simulations, the protein interconverts between the open and closed states. Ligand-free FABP4 prefers the closed state, whereas ligand binding induces a conformational transition to the open state. Coupled with opening and closing of FABP4, the ligand adopts distinct binding modes, which are identified and compared with crystal structures. The concerted dynamics of protein and ligand suggests that there may exist multiple FABP4-ligand binding conformations. Thus, this work provides details about how ligand binding affects the conformational preference of FABP4 and how ligand binding is coupled with a conformational change of FABP4 at an atomic level.

| S-EPMC4196735 | biostudies-literature

On the improvement of free-energy calculation from steered molecular dynamics simulations using adaptive stochastic perturbation protocols.

Project description:The potential of mean force (PMF) calculation in single molecule manipulation experiments performed via the steered molecular dynamics (SMD) technique is a computationally very demanding task because the analyzed system has to be perturbed very slowly to be kept close to equilibrium. Faster perturbations, far from equilibrium, increase dissipation and move the average work away from the underlying free energy profile, and thus introduce a bias into the PMF estimate. The Jarzynski equality offers a way to overcome the bias problem by being able to produce an exact estimate of the free energy difference, regardless of the perturbation regime. However, with a limited number of samples and high dissipation the Jarzynski equality also introduces a bias. In our previous work, based on the Brownian motion formalism, we introduced three stochastic perturbation protocols aimed at improving the PMF calculation with the Jarzynski equality in single molecule manipulation experiments and analogous computer simulations. This paper describes the PMF reconstruction results based on full-atom molecular dynamics simulations, obtained with those three protocols. We also want to show that the protocols are applicable with the second-order cumulant expansion formula. Our protocols offer a very noticeable improvement over the simple constant velocity pulling. They are able to produce an acceptable estimate of PMF with a significantly reduced bias, even with very fast perturbation regimes. Therefore, the protocols can be adopted as practical and efficient tools for the analysis of mechanical properties of biological molecules.

| S-EPMC4169427 | biostudies-literature

Mechanism of Cardiac Tropomyosin Transitions on Filamentous Actin As Revealed by All-Atom Steered Molecular Dynamics Simulations.

Project description:The three-state model of tropomyosin (Tm) positioning along filamentous actin allows for Tm to act as a gate for myosin head binding with actin. The blocked state of Tm prevents myosin binding, while the open state allows for strong binding. Intermediate to this transition is the closed state. The details of the transition from the blocked to the closed state and then finally to the open state by Tm have not been fully accessible to experiment. Utilizing steered molecular dynamics, we investigate the work required to move the Tm strand through the extant set of proposed transitions. We find that an azimuthal motion around the actin filament by Tm is most probable in spite of increased initial energy barrier from the topographical landscape of actin.

| S-EPMC6019281 | biostudies-literature

Nanomechanical properties of MscL α helices: A steered molecular dynamics study.

Project description:Gating of mechanosensitive (MS) channels is driven by a hierarchical cascade of movements and deformations of transmembrane helices in response to bilayer tension. Determining the intrinsic mechanical properties of the individual transmembrane helices is therefore central to understanding the intricacies of the gating mechanism of MS channels. We used a constant-force steered molecular dynamics (SMD) approach to perform unidirectional pulling tests on all the helices of MscL in M. tuberculosis and E. coli homologs. Using this method, we could overcome the issues encountered with the commonly used constant-velocity SMD simulations, such as low mechanical stability of the helix during stretching and high dependency of the elastic properties on the pulling rate. We estimated Young's moduli of the α-helices of MscL to vary between 0.2 and 12.5 GPa with TM2 helix being the stiffest. We also studied the effect of water on the properties of the pore-lining TM1 helix. In the absence of water, this helix exhibited a much stiffer response. By monitoring the number of hydrogen bonds, it appears that water acts like a 'lubricant' (softener) during TM1 helix elongation. These data shed light on another physical aspect underlying hydrophobic gating of MS channels, in particular MscL.

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