Demet Akten - Academia.edu (original) (raw)
Papers by Demet Akten
Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) catalyzing the sixth step of glycolysis has been... more Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) catalyzing the sixth step of glycolysis has been investigated for allosteric features that might be used as potential target for specific inhibition of Staphylococcus aureus (S.aureus). X-ray structure of bacterial enzyme for which a tunnel-like opening passing through the center previously proposed as an allosteric site has been subjected to six independent 500 ns long Molecular Dynamics simulations. Harmonic bond restraints were employed at key residues to underline the allosteric feature of this region. A noticeable reduction was observed in the mobility of NAD+ binding domains when restrictions were applied. Also, a substantial decrease in cross-correlations between distant C$\alpha$ fluctuations was detected throughout the structure. Mutual information (MI) analysis revealed a similar decrease in the degree of correspondence in positional fluctuations in all directions everywhere in the receptor. MI between backbone and side-chain torsional variations changed its distribution profile and decreased considerably around the catalytic sites when restraints were employed. Principal component analysis clearly showed that the restrained state sampled a narrower range of conformations than apo state, especially in the first principal mode due to restriction in the conformational flexibility of NAD+ binding domain. Clustering the trajectory based on catalytic site residues displayed a smaller repertoire of conformations for restrained state compared to apo. Representative snapshots subjected to k-shortest pathway analysis revealed the impact of bond restraints on the allosteric communication which displayed distinct optimal and suboptimal pathways for two states, where observed frequencies of critical residues Gln51 and Val283 at the proposed site changed considerably.
Applied Sciences
Two independent 1.5 μs long MD simulations were conducted for the fully atomistic model of the hu... more Two independent 1.5 μs long MD simulations were conducted for the fully atomistic model of the human beta2-adrenergic receptor (β2AR) in a complex with a G protein to investigate the signal transmission in a fully active state via mutual information and transfer entropy based on α-carbon displacements and rotameric states of backbone and side-chain torsion angles. Significant correlations between fluctuations in α-Carbon displacements were mostly detected between transmembrane (TM) helices, especially TM5 and TM6 located at each end of ICL3 and TM7. Signal transmission across β2-AR was quantified by shared mutual information; a high amount of correspondence was distinguished in almost all loop regions when rotameric states were employed. Moreover, polar residues, especially Arg, made the most contribution to signal transmission via correlated side-chain rotameric fluctuations as they were more frequently observed in loop regions than hydrophobic residues. Furthermore, transfer entro...
The Journal of Physical Chemistry B
Chemical Biology & Drug Design
This study investigates the structural distinctiveness of orthosteric ligand‐binding sites of sev... more This study investigates the structural distinctiveness of orthosteric ligand‐binding sites of several human β2 adrenergic receptor (β2‐AR) conformations that have been obtained from a set of independent molecular dynamics (MD) simulations in the presence of intracellular loop 3 (ICL3). A docking protocol was established in order to classify each receptor conformation via its binding affinity to selected ligands with known efficacy. This work's main goal was to reveal many subtle features of the ligand‐binding site, presenting alternative conformations, which might be considered as either active‐ or inactive‐like but mostly specific for that ligand. Agonists, inverse agonists, and antagonists were docked to each MD conformer with distinct binding pockets, using different docking tools and scoring functions. Mostly favored receptor conformation persistently observed in all docking/scoring evaluations was classified as active or inactive based on the type of ligand's biological effect. Classified MD conformers were further tested for their ability to discriminate agonists from inverse agonists/antagonists, and several conformers were proposed as important targets to be used in virtual screening experiments that were often limited to a single X‐ray structure.
Journal of Computer-Aided Molecular Design
Molecular Informatics, 2011
We aim to uncover the binding modes of benzothiazoles, which have been reported as specific inhib... more We aim to uncover the binding modes of benzothiazoles, which have been reported as specific inhibitors of triosephosphate isomerase from the parasite Trypanosoma cruzi (TcTIM), by performing blind dockings on both TcTIM and human TIM (hTIM). Detailed analysis of binding sites and specific interactions are carried out based on ensemble dockings to multiple receptor conformers obtained from molecular dynamics simulations. In TcTIM dimer dockings, the inhibitors preferentially bind to the tunnel-shaped cavity formed at the interface of the subunits, whereas non-inhibitors mostly choose other sites. In contrast, TcTIM monomer binding interface and hTIM dimer interface do not present a specific binding site for the inhibitors. These findings point to the importance of the tunnel and of the dimeric form for inhibition of TcTIM. Specific interactions of the inhibitors and their sulfonate-free derivatives with the receptor residues indicate the significance of sulfonate group for binding affinity and positioning on the TcTIM dimer interface. One of the inhibitors also binds to the active site, which may explain its relatively higher inhibition effect on hTIM.
Journal of Biomolecular Structure and Dynamics, 2009
BMC Structural Biology, 2016
Frontiers in Molecular Biosciences
Journal of Biomolecular Structure and Dynamics
Drug repositioning has recently become one of the widely used drug design approaches in proposing... more Drug repositioning has recently become one of the widely used drug design approaches in proposing alternative compounds with potentially fewer side effects. In this study, structure-based pharmacophore modelling and docking was used to screen existing drug molecules to bring forward potential modulators for ligand-binding domain of human glucocorticoid receptor (hGR). There exist several drug molecules targeting hGR, yet their apparent side effects still persist. Our goal was to disclose new compounds via screening existing drug compounds to bring forward fast and explicit solutions. The so-called shared pharmacophore model was created using the most persistent pharmacophore features shared by several crystal structures of the receptor. The shared model was first used to screen a small database of 75 agonists and 300 antagonists/decoys, and exhibited a successful outcome in its ability to distinguish agonists from antagonists/decoys. Then, it was used to screen a database of over 5000 molecules composed of FDA-approved, worldwide used and investigational drug compounds. A total of 110 compounds satisfying the pharmacophore requirements were subjected to different docking experiments for further assessment of their binding ability. In the final hit list of 54 compounds which fulfilled all scoring criteria, 19 of them were nonsteroidal and when further investigated, each presented a unique scaffold with little structural resemblance to any known nonsteroidal GR modulators. Independent 100 ns long MD simulations conducted on three selected drug candidates in complex with hGR displayed stable conformations incorporating several hydrogen bonds common to all three compounds and the reference molecule dexamethasone.Communicated by Ramaswamy H. Sarma.
Frontiers in Molecular Biosciences
Proteins: Structure, Function, and Bioinformatics
Journal of Molecular Graphics and Modelling, 2014
Novel high affinity compounds for human β2-adrenergic receptor (β2-AR) were searched among the cl... more Novel high affinity compounds for human β2-adrenergic receptor (β2-AR) were searched among the clean drug-like subset of ZINC database consisting of 9,928,465 molecules that satisfy the Lipinski's rule of five. The screening protocol consisted of a high-throughput pharmacophore screening followed by an extensive amount of docking and rescoring. The pharmacophore model was composed of key features shared by all five inactive states of β2-AR in complex with inverse agonists and antagonists. To test the discriminatory power of the pharmacophore model, a small-scale screening was initially performed on a database consisting of 117 compounds of which 53 antagonists were taken as active inhibitors and 64 agonists as inactive inhibitors. Accordingly, 7.3% of the ZINC database subset (729,413 compounds) satisfied the pharmacophore requirements, along with 44 antagonists and 17 agonists. Afterwards, all these hit compounds were docked to the inactive apo form of the receptor using various docking and scoring protocols. Following each docking experiment, the best pose was further evaluated based on the existence of key residues for antagonist binding in its vicinity. After final evaluations based on the human intestinal absorption (HIA) and the blood brain barrier (BBB) penetration properties, 62 hit compounds have been clustered based on their structural similarity and as a result four scaffolds were revealed. Two of these scaffolds were also observed in three high affinity compounds with experimentally known Ki values. Moreover, novel chemical compounds with distinct structures have been determined as potential β2-AR drug candidates.
Simulation Methods for Polymers, 2004
Energ Fuel, 2003
We present a molecular model for the adsorption of CO 2 , N 2 , H 2 , and their mixtures in dehyd... more We present a molecular model for the adsorption of CO 2 , N 2 , H 2 , and their mixtures in dehydrated zeolite Na-4A. The interatomic potentials for this model were developed by comparing the results of grand canonical Monte Carlo (GCMC) simulations of single-component adsorption at room temperature with experimental measurements. GCMC simulation is also used to assess the adsorption selectivity of CO 2 /N 2 and CO 2 /H 2 mixtures, as a function of temperature and gasphase composition. At room temperature, Na-4A is strongly selective for CO 2 over both N 2 and H 2 , although this selectivity decreases slightly as the gas-phase pressure increases. Ideal adsorbed solution theory is shown to give accurate predictions of the adsorption selectivity at low CO 2 partial pressures, provided that a functional form that accurately describes the CO 2 singlecomponent isotherm is used. The adsorption properties of CO 2 /N 2 mixtures in Na-4A are compared to the same mixtures in silicalite.
Journal of biomolecular structure & dynamics, Jan 6, 2015
Peptide bond hydrolysis of bovine serum albumin (BSA) by chymotrypsin and trypsin was investigate... more Peptide bond hydrolysis of bovine serum albumin (BSA) by chymotrypsin and trypsin was investigated by employing time-resolved fluorescence spectroscopy. As a fluorescent cross-linking reagent, N-(1-pyrenyl) maleimide (PM) was attached to BSA, through all free amine groups of arginine, lysine, and/or single free thiol (Cys34). Time-resolved fluorescence spectroscopy was used to monitor fluorescence decays analyzed by exponential series method to obtain the changes in lifetime distributions. After the exposure of synthesized protein substrate PM-BSA to chymotrypsin and trypsin, it is observed that each protease produced a distinct change in the lifetime distribution profile, which was attributed to distinct chemical environments created by short peptide fragments in each hydrolysate. The persistence of excimer emission at longer lifetime regions for chymotrypsin, as opposed to trypsin, suggested the presence of small-scale hydrophobic clusters that might prevent some excimers from bei...
Biophysical Journal, 2015
The tunnel region at triosephosphate isomerase (TIM)'s di... more The tunnel region at triosephosphate isomerase (TIM)'s dimer interface, distant from its catalytic site, is a target site for certain benzothiazole derivatives that inhibit TIM's catalytic activity in Trypanosoma cruzi, the parasite that causes Chagas disease. We performed multiple 100-ns molecular-dynamics (MD) simulations and elastic network modeling (ENM) on both apo and complex structures to shed light on the still unclear inhibitory mechanism of one such inhibitor, named bt10. Within the time frame of our MD simulations, we observed stabilization of aromatic clusters at the dimer interface and enhancement of intersubunit hydrogen bonds in the presence of bt10, which point to an allosteric effect rather than destabilization of the dimeric structure. The collective dynamics dictated by the topology of TIM is known to facilitate the closure of its catalytic loop over the active site that is critical for substrate entrance and product release. We incorporated the ligand's effect on vibrational dynamics by applying mixed coarse-grained ENM to each one of 54,000 MD snapshots. Using this computationally efficient technique, we observed altered collective modes and positive shifts in eigenvalues due to the constraining effect of bt10 binding. Accordingly, we observed allosteric changes in the catalytic loop's dynamics, flexibility, and correlations, as well as the solvent exposure of catalytic residues. A newly (to our knowledge) introduced technique that performs residue-based ENM scanning of TIM revealed the tunnel region as a key binding site that can alter global dynamics of the enzyme.
Biophysical Journal, 2014
Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) catalyzing the sixth step of glycolysis has been... more Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) catalyzing the sixth step of glycolysis has been investigated for allosteric features that might be used as potential target for specific inhibition of Staphylococcus aureus (S.aureus). X-ray structure of bacterial enzyme for which a tunnel-like opening passing through the center previously proposed as an allosteric site has been subjected to six independent 500 ns long Molecular Dynamics simulations. Harmonic bond restraints were employed at key residues to underline the allosteric feature of this region. A noticeable reduction was observed in the mobility of NAD+ binding domains when restrictions were applied. Also, a substantial decrease in cross-correlations between distant C$\alpha$ fluctuations was detected throughout the structure. Mutual information (MI) analysis revealed a similar decrease in the degree of correspondence in positional fluctuations in all directions everywhere in the receptor. MI between backbone and side-chain torsional variations changed its distribution profile and decreased considerably around the catalytic sites when restraints were employed. Principal component analysis clearly showed that the restrained state sampled a narrower range of conformations than apo state, especially in the first principal mode due to restriction in the conformational flexibility of NAD+ binding domain. Clustering the trajectory based on catalytic site residues displayed a smaller repertoire of conformations for restrained state compared to apo. Representative snapshots subjected to k-shortest pathway analysis revealed the impact of bond restraints on the allosteric communication which displayed distinct optimal and suboptimal pathways for two states, where observed frequencies of critical residues Gln51 and Val283 at the proposed site changed considerably.
Applied Sciences
Two independent 1.5 μs long MD simulations were conducted for the fully atomistic model of the hu... more Two independent 1.5 μs long MD simulations were conducted for the fully atomistic model of the human beta2-adrenergic receptor (β2AR) in a complex with a G protein to investigate the signal transmission in a fully active state via mutual information and transfer entropy based on α-carbon displacements and rotameric states of backbone and side-chain torsion angles. Significant correlations between fluctuations in α-Carbon displacements were mostly detected between transmembrane (TM) helices, especially TM5 and TM6 located at each end of ICL3 and TM7. Signal transmission across β2-AR was quantified by shared mutual information; a high amount of correspondence was distinguished in almost all loop regions when rotameric states were employed. Moreover, polar residues, especially Arg, made the most contribution to signal transmission via correlated side-chain rotameric fluctuations as they were more frequently observed in loop regions than hydrophobic residues. Furthermore, transfer entro...
The Journal of Physical Chemistry B
Chemical Biology & Drug Design
This study investigates the structural distinctiveness of orthosteric ligand‐binding sites of sev... more This study investigates the structural distinctiveness of orthosteric ligand‐binding sites of several human β2 adrenergic receptor (β2‐AR) conformations that have been obtained from a set of independent molecular dynamics (MD) simulations in the presence of intracellular loop 3 (ICL3). A docking protocol was established in order to classify each receptor conformation via its binding affinity to selected ligands with known efficacy. This work's main goal was to reveal many subtle features of the ligand‐binding site, presenting alternative conformations, which might be considered as either active‐ or inactive‐like but mostly specific for that ligand. Agonists, inverse agonists, and antagonists were docked to each MD conformer with distinct binding pockets, using different docking tools and scoring functions. Mostly favored receptor conformation persistently observed in all docking/scoring evaluations was classified as active or inactive based on the type of ligand's biological effect. Classified MD conformers were further tested for their ability to discriminate agonists from inverse agonists/antagonists, and several conformers were proposed as important targets to be used in virtual screening experiments that were often limited to a single X‐ray structure.
Journal of Computer-Aided Molecular Design
Molecular Informatics, 2011
We aim to uncover the binding modes of benzothiazoles, which have been reported as specific inhib... more We aim to uncover the binding modes of benzothiazoles, which have been reported as specific inhibitors of triosephosphate isomerase from the parasite Trypanosoma cruzi (TcTIM), by performing blind dockings on both TcTIM and human TIM (hTIM). Detailed analysis of binding sites and specific interactions are carried out based on ensemble dockings to multiple receptor conformers obtained from molecular dynamics simulations. In TcTIM dimer dockings, the inhibitors preferentially bind to the tunnel-shaped cavity formed at the interface of the subunits, whereas non-inhibitors mostly choose other sites. In contrast, TcTIM monomer binding interface and hTIM dimer interface do not present a specific binding site for the inhibitors. These findings point to the importance of the tunnel and of the dimeric form for inhibition of TcTIM. Specific interactions of the inhibitors and their sulfonate-free derivatives with the receptor residues indicate the significance of sulfonate group for binding affinity and positioning on the TcTIM dimer interface. One of the inhibitors also binds to the active site, which may explain its relatively higher inhibition effect on hTIM.
Journal of Biomolecular Structure and Dynamics, 2009
BMC Structural Biology, 2016
Frontiers in Molecular Biosciences
Journal of Biomolecular Structure and Dynamics
Drug repositioning has recently become one of the widely used drug design approaches in proposing... more Drug repositioning has recently become one of the widely used drug design approaches in proposing alternative compounds with potentially fewer side effects. In this study, structure-based pharmacophore modelling and docking was used to screen existing drug molecules to bring forward potential modulators for ligand-binding domain of human glucocorticoid receptor (hGR). There exist several drug molecules targeting hGR, yet their apparent side effects still persist. Our goal was to disclose new compounds via screening existing drug compounds to bring forward fast and explicit solutions. The so-called shared pharmacophore model was created using the most persistent pharmacophore features shared by several crystal structures of the receptor. The shared model was first used to screen a small database of 75 agonists and 300 antagonists/decoys, and exhibited a successful outcome in its ability to distinguish agonists from antagonists/decoys. Then, it was used to screen a database of over 5000 molecules composed of FDA-approved, worldwide used and investigational drug compounds. A total of 110 compounds satisfying the pharmacophore requirements were subjected to different docking experiments for further assessment of their binding ability. In the final hit list of 54 compounds which fulfilled all scoring criteria, 19 of them were nonsteroidal and when further investigated, each presented a unique scaffold with little structural resemblance to any known nonsteroidal GR modulators. Independent 100 ns long MD simulations conducted on three selected drug candidates in complex with hGR displayed stable conformations incorporating several hydrogen bonds common to all three compounds and the reference molecule dexamethasone.Communicated by Ramaswamy H. Sarma.
Frontiers in Molecular Biosciences
Proteins: Structure, Function, and Bioinformatics
Journal of Molecular Graphics and Modelling, 2014
Novel high affinity compounds for human β2-adrenergic receptor (β2-AR) were searched among the cl... more Novel high affinity compounds for human β2-adrenergic receptor (β2-AR) were searched among the clean drug-like subset of ZINC database consisting of 9,928,465 molecules that satisfy the Lipinski's rule of five. The screening protocol consisted of a high-throughput pharmacophore screening followed by an extensive amount of docking and rescoring. The pharmacophore model was composed of key features shared by all five inactive states of β2-AR in complex with inverse agonists and antagonists. To test the discriminatory power of the pharmacophore model, a small-scale screening was initially performed on a database consisting of 117 compounds of which 53 antagonists were taken as active inhibitors and 64 agonists as inactive inhibitors. Accordingly, 7.3% of the ZINC database subset (729,413 compounds) satisfied the pharmacophore requirements, along with 44 antagonists and 17 agonists. Afterwards, all these hit compounds were docked to the inactive apo form of the receptor using various docking and scoring protocols. Following each docking experiment, the best pose was further evaluated based on the existence of key residues for antagonist binding in its vicinity. After final evaluations based on the human intestinal absorption (HIA) and the blood brain barrier (BBB) penetration properties, 62 hit compounds have been clustered based on their structural similarity and as a result four scaffolds were revealed. Two of these scaffolds were also observed in three high affinity compounds with experimentally known Ki values. Moreover, novel chemical compounds with distinct structures have been determined as potential β2-AR drug candidates.
Simulation Methods for Polymers, 2004
Energ Fuel, 2003
We present a molecular model for the adsorption of CO 2 , N 2 , H 2 , and their mixtures in dehyd... more We present a molecular model for the adsorption of CO 2 , N 2 , H 2 , and their mixtures in dehydrated zeolite Na-4A. The interatomic potentials for this model were developed by comparing the results of grand canonical Monte Carlo (GCMC) simulations of single-component adsorption at room temperature with experimental measurements. GCMC simulation is also used to assess the adsorption selectivity of CO 2 /N 2 and CO 2 /H 2 mixtures, as a function of temperature and gasphase composition. At room temperature, Na-4A is strongly selective for CO 2 over both N 2 and H 2 , although this selectivity decreases slightly as the gas-phase pressure increases. Ideal adsorbed solution theory is shown to give accurate predictions of the adsorption selectivity at low CO 2 partial pressures, provided that a functional form that accurately describes the CO 2 singlecomponent isotherm is used. The adsorption properties of CO 2 /N 2 mixtures in Na-4A are compared to the same mixtures in silicalite.
Journal of biomolecular structure & dynamics, Jan 6, 2015
Peptide bond hydrolysis of bovine serum albumin (BSA) by chymotrypsin and trypsin was investigate... more Peptide bond hydrolysis of bovine serum albumin (BSA) by chymotrypsin and trypsin was investigated by employing time-resolved fluorescence spectroscopy. As a fluorescent cross-linking reagent, N-(1-pyrenyl) maleimide (PM) was attached to BSA, through all free amine groups of arginine, lysine, and/or single free thiol (Cys34). Time-resolved fluorescence spectroscopy was used to monitor fluorescence decays analyzed by exponential series method to obtain the changes in lifetime distributions. After the exposure of synthesized protein substrate PM-BSA to chymotrypsin and trypsin, it is observed that each protease produced a distinct change in the lifetime distribution profile, which was attributed to distinct chemical environments created by short peptide fragments in each hydrolysate. The persistence of excimer emission at longer lifetime regions for chymotrypsin, as opposed to trypsin, suggested the presence of small-scale hydrophobic clusters that might prevent some excimers from bei...
Biophysical Journal, 2015
The tunnel region at triosephosphate isomerase (TIM)'s di... more The tunnel region at triosephosphate isomerase (TIM)'s dimer interface, distant from its catalytic site, is a target site for certain benzothiazole derivatives that inhibit TIM's catalytic activity in Trypanosoma cruzi, the parasite that causes Chagas disease. We performed multiple 100-ns molecular-dynamics (MD) simulations and elastic network modeling (ENM) on both apo and complex structures to shed light on the still unclear inhibitory mechanism of one such inhibitor, named bt10. Within the time frame of our MD simulations, we observed stabilization of aromatic clusters at the dimer interface and enhancement of intersubunit hydrogen bonds in the presence of bt10, which point to an allosteric effect rather than destabilization of the dimeric structure. The collective dynamics dictated by the topology of TIM is known to facilitate the closure of its catalytic loop over the active site that is critical for substrate entrance and product release. We incorporated the ligand's effect on vibrational dynamics by applying mixed coarse-grained ENM to each one of 54,000 MD snapshots. Using this computationally efficient technique, we observed altered collective modes and positive shifts in eigenvalues due to the constraining effect of bt10 binding. Accordingly, we observed allosteric changes in the catalytic loop's dynamics, flexibility, and correlations, as well as the solvent exposure of catalytic residues. A newly (to our knowledge) introduced technique that performs residue-based ENM scanning of TIM revealed the tunnel region as a key binding site that can alter global dynamics of the enzyme.
Biophysical Journal, 2014