Dynamics Free Energy Calculations Research Articles

Structural Studies of the 3',3'-cGAMP Riboswitch Induced by Cognate and Noncognate Ligands Using Molecular Dynamics Simulation.

Riboswtich RNAs can control gene expression through the structural change induced by the corresponding small-molecule ligands. Molecular dynamics simulations and free energy calculations on the aptamer domain of the 3′,3′-cGAMP riboswitch in the ligand-free, cognate-bound and noncognate-bound states were performed to investigate the structural features of the 3′,3′-cGAMP riboswitch induced by the 3′,3′-cGAMP ligand and the specificity of ligand recognition. The results revealed that the aptamer of the 3′,3′-cGAMP riboswitch in the ligand-free state has a smaller binding pocket and a relatively compact structure versus that in the 3′,3′-cGAMP-bound state. The binding of the 3′,3′-cGAMP molecule to the 3′,3′-cGAMP riboswitch induces the rotation of P1 helix through the allosteric communication from the binding sites pocket containing the J1/2, J1/3 and J2/3 junction to the P1 helix. Simultaneously, these simulations also revealed that the preferential binding of the 3′,3′-cGAMP riboswitch to its cognate ligand, 3′,3′-cGAMP, over its noncognate ligand, c-di-GMP and c-di-AMP. The J1/2 junction in the 3′,3′-cGAMP riboswitch contributing to the specificity of ligand recognition have also been found.

3D-QSAR-aided design of potent c-Met inhibitors using molecular dynamics simulation and binding free energy calculation

Mesenchymal-epithelial transition factor (c-Met) is a member of receptor tyrosine kinase. It involves in various cellular signaling pathways which includes proliferation, motility, migration, and invasion. Over-expression of c-Met has been reported in various cancers. Hence, it is an ideal therapeutic target for cancer. The main objective of the study is to identify crucial residues involved in the inhibition of c-Met kinase and to design a series of potent imidazo [4,5-b] pyrazine derivatives as c-Met inhibitors. Docking was used to identify important active site residues involved in the inhibition of c-Met kinase which was further validated by 100 ns of molecular dynamics simulation and free energy calculation using molecular mechanics generalized born surface area. Furthermore, binding energy decomposition identified that residues Tyr1230, Met1211, Asp1222, Tyr1159, Met1160, Val1092, Ala1108, and Leu1157 contributed favorably to the binding stability of compound 32. Receptor-guided Comparative Molecular Field Analysis (CoMFA) (q2 = 0.751, NOC = 6, r2 = 0.933) and Comparative Molecular Similarity Indices Analysis (COMSIA) (q2 = 0.744, NOC = 6, r2 = 0.950) models were generated based on the docked conformation of the most active compound 32. The robustness of these models was tested using various validation techniques and found to be predictive. The results of CoMFA and CoMSIA contour maps exposed the regions favorable to enhance the activity. Based on this information, 27 novel c-Met inhibitors were designed. These designed compounds exhibited potent activity than the most active compound of the existing dataset.Communicated by Ramaswamy H. Sarma

The role of LasR active site amino acids in the interaction with the Acyl Homoserine Lactones (AHLs) analogues: A computational study

The present work combines molecular docking calculations, 3D-QSAR, molecular dynamics simulations and free binding energy calculations (MM/PBSA and MM/GBSA) in a set of 28 structural analogues of acyl homoserine lactones with Quorum Sensing antagonist activity. The aim of this work is to understand how ligand binds and is affected by the molecular microenvironment in the active site of the LasR receptor for pseudomonas aeruginosa. We also study the stability of the interaction to find key structural characteristics that explain the antagonist activities of this set of ligands. This information is relevant for the rational modification or design of molecules and their identification as powerful LasR modulators.The analysis of molecular docking simulations shows that the 28 analogues have a similar binding mode compared to the native ligand. The carbonyl groups belonging to the lactone ring and the amide group of the acyl chain are oriented towards the amino acids forming hydrogen bond like interactions. The difference in antagonist activity is due to location and orientation of the LasR side chains within the hydrophobic pocket in its binding site. Additionally, we carried out molecular dynamics simulations to understand the conformational changes in the ligand-receptor interaction and the stability of each complex. Results show a direct relationship among the interaction energies of the ligands and the activities as an antagonist of the LasR receptor.

Hydrophobic alkyl chains substituted to the 8-position of cyclic nucleotides enhance activation of CNG and HCN channels by an intricate enthalpy - entropy compensation

Cyclic nucleotide-gated (CNG) and hyperpolarization-activated cyclic nucleotide-gated (HCN) channels are tetrameric non-specific cation channels in the plasma membrane that are activated by either cAMP or cGMP binding to specific binding domains incorporated in each subunit. Typical apparent affinities of these channels for these cyclic nucleotides range from several hundred nanomolar to tens of micromolar. Here we synthesized and characterized novel cAMP and cGMP derivatives by substituting either hydrophobic alkyl chains or similar-sized more hydrophilic heteroalkyl chains to the 8-position of the purine ring with the aim to obtain full agonists of higher potency. The compounds were tested in homotetrameric CNGA2, heterotetrameric CNGA2:CNGA4:CNGB1b and homotetrameric HCN2 channels. We show that nearly all compounds are full agonists and that longer alkyl chains systematically increase the apparent affinity, at the best more than 30 times. The effects are stronger in CNG than HCN2 channels which, however, are constitutively more sensitive to cAMP. Kinetic analyses reveal that the off-rate is significantly slowed by the hydrophobic alkyl chains. Molecular dynamics simulations and free energy calculations suggest that an intricate enthalpy - entropy compensation underlies the higher apparent affinity of the derivatives with the longer alkyl chains, which is shown to result from a reduced loss of configurational entropy upon binding.

Identification of (4-(9H-fluoren-9-yl) piperazin-1-yl) methanone derivatives as falcipain 2 inhibitors active against Plasmodium falciparum cultures

BackgroundFalcipain 2 (FP-2) is the hemoglobin-degrading cysteine protease of Plasmodium falciparum most extensively targeted to develop novel antimalarials. However, no commercial antimalarial drugs based on FP-2 inhibition are available yet due to the low selectivity of most FP-2 inhibitors against the human cysteine proteases. MethodsA structure-based virtual screening (SVBS) using Maybridge HitFinder™ compound database was conducted to identify potential FP-2 inhibitors. In vitro enzymatic and cell-growth inhibition assays were performed for the top-scoring compounds. Docking, molecular dynamics (MD) simulations and free energy calculations were employed to study the interaction of the best hits with FP-2 and other related enzymes. Results and conclusionsTwo hits based on 4-(9H-fluoren-9-yl) piperazin-1-yl) methanone scaffold, HTS07940 and HTS08262, were identified as inhibitors of FP-2 (half-maximal inhibitory concentration (IC50) = 64 μM and 14.7 μM, respectively) without a detectable inhibition against the human off-target cathepsin K (hCatK). HTS07940 and HTS08262 inhibited the growth of the multidrug-resistant P. falciparum strain FCR3 in culture (half-maximal inhibitory concentrations (IC50) = 2.91 μM and 34 μM, respectively) and exhibited only moderate cytotoxicity against HeLa cells (Half-maximal cytotoxic concentration (CC50) = 133 μM and 350 μM, respectively). Free energy calculations reproduced the experimental affinities of the hits for FP-2 and explained the selectivity with respect to hCatK. General significanceTo the best of our knowledge, HTS07940 stands among the most selective FP-2 inhibitors identified by SBVS reported so far, displaying moderate antiplasmodial activity and low cytotoxicity against human cells. Hence, this compound constitutes a promising lead for the design of more potent and selective FP-2 inhibitors.

The influence of residue in the position of 116 on the inhibitory potency of TH588 for MTH1

As one of the first-in-class inhibitor, TH588 was found to be efficient in the suppression of MutT homolog1 (MTH1). A recent work shows that the inhibitory potency of TH588 against human MTH1 (hsMTH1) is approximately 20-fold over that of mouse MTH1 (mmMTH1) and identifies residue in position 116 in MTH1 has an important contribution to TH588 affinity. But the effect of residue Leu or Met in position 116 on the binding affinity remains unclear. In this study, molecular dynamics (MD) simulations and free energy calculations were used to elucidate the mechanism about the effect of residue 116 to the different inhibitory potency of TH588 against MTH1. The binding free energy of TH588 in M116 complexes predicated by the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) is much lower than that in L116 complexes, which is consistent with the experiment results. The analysis of the individual energy terms suggests that the non-polar interactions are important for distinguishing the binding of TH588. The MD results show that the Leu116 disrupts the interactions between Asn33 and TH588, thus induces the conformational changes of Asn33 as well as TH588. The altered interactions between TH588 and mmMTH1 change the flexibility of TH588, which could induce the remarkable conformational fluctuation of mmMTH1. The conformations of the two loops covering the binding pocket have obvious influence on the opening or closure of the active site. The more open binding site may explain the lower inhibitor potency of TH588 against mmMTH1 than hsMTH1. Our results provide mechanistic insight into the effect of different residue Leu or Met in position 116 on the binding affinity of TH588 for MTH1, which is expected to contribute to the further rational design of more potent inhibitors.

Structure and permeability of ionomers studied by atomistic molecular simulation combined with the theory of solutions in the energy representation.

Ionomers play a key role in forming the catalyst layer of polymer electrolyte fuel cells. In the present work, we performed atomistic molecular dynamics simulations and free-energy calculations with the energy-representation method for sulfonated polyethersulfone (SPES) and its derivatives toward the rational design of ionomers for carbon alloy catalysts. It was observed that H2O aggregates strongly in the branched SPES systems with fluorocarbons and is located homogeneously in the systems without fluorocarbons. The O2 permeability was then examined within the framework of the solubility-diffusion mechanism. The permeability was seen to be large for the branched SPES with fluorocarbons, indicating that the performance of ionomers as a permeation medium for O2 may be tuned by the flexibility and branching of the polymer chain.

Negatively Charged Gangliosides Promote Membrane Association of Amphipathic Neurotransmitters

Lipophilic neurotransmitters (NTs) such as dopamine are chemical messengers enabling neurotransmission by adhering onto the extracellular surface of the post-synaptic membrane in a synapse, followed by binding to their receptors. Previous studies have shown that the strength of the NT–membrane association is dependent on the lipid composition of the membrane. Negatively charged lipids such as phosphatidylserine, phosphatidylglycerol, and phosphatidic acid have been indicated to promote NT–membrane binding, however these anionic lipids reside almost exclusively in the intracellular leaflet of the post-synaptic membrane instead of the extracellular leaflet facing the synaptic cleft. Meanwhile, the extracellular leaflet is relatively rich in biologically relevant anionic gangliosides such as monosialotetrahexosylganglioside (GM1), yet the role of gangliosides in NT–membrane association is not clear. Here, we explored the role of GM1 in modulating the binding of dopamine and histamine (as amphipathic/cationic NTs) as well as acetylcholine (as a hydrophilic/cationic NT) with the post-synaptic membrane surface. Atomistic molecular dynamics simulations and free energy calculations indicated that GM1 fosters membrane association of histamine and dopamine. For acetylcholine, this effect was not observed. The in silico results suggest that gangliosides form a charge-based vestibule in front of the post-synaptic membrane, attracting amphipathic NTs to the vicinity of the membrane. The results also stress the importance to understand the significance of the structural details of NTs, as exemplified by the GM1–acetylcholine interaction. In a larger context, the NT–membrane adherence, coupled to lateral diffusion in the membrane plane, is proposed to improve neurotransmission efficiency by advancing NT entry into the membrane-embedded ligand-binding sites.

Assigning crystallographic electron densities with free energy calculations-The case of the fluoride channel Fluc.

Approximately 90% of the structures in the Protein Data Bank (PDB) were obtained by X-ray crystallography or electron microscopy. Whereas the overall quality of structure is considered high, thanks to a wide range of tools for structure validation, uncertainties may arise from density maps of small molecules, such as organic ligands, ions or water, which are non-covalently bound to the biomolecules. Even with some experience and chemical intuition, the assignment of such disconnected electron densities is often far from obvious. In this study, we suggest the use of molecular dynamics (MD) simulations and free energy calculations, which are well-established computational methods, to aid in the assignment of ambiguous disconnected electron densities. Specifically, estimates of (i) relative binding affinities, for instance between an ion and water, (ii) absolute binding free energies, i.e., free energies for transferring a solute from bulk solvent to a binding site, and (iii) stability assessments during equilibrium simulations may reveal the most plausible assignments. We illustrate this strategy using the crystal structure of the fluoride specific channel (Fluc), which contains five disconnected electron densities previously interpreted as four fluoride and one sodium ion. The simulations support the assignment of the sodium ion. In contrast, calculations of relative and absolute binding free energies as well as stability assessments during free MD simulations suggest that four of the densities represent water molecules instead of fluoride. The assignment of water is compatible with the loss of these densities in the non-conductive F82I/F85I mutant of Fluc. We critically discuss the role of the ion force fields for the calculations presented here. Overall, these findings indicate that MD simulations and free energy calculations are helpful tools for modeling water and ions into crystallographic density maps.

Molecular dynamics simulation and binding free energy studies of novel leads belonging to the benzofuran class inhibitors of Mycobacterium tuberculosis Polyketide Synthase 13

In this work, the binding mechanism of new Polyketide Synthase 13 (Pks13) inhibitors has been studied through molecular dynamics simulation and free energy calculations. The drug Tam1 and its analogs, belonging to the benzofuran class, were submitted to 100 ns simulations, and according to the results obtained for root mean square deviation, all the simulations converged from approximately 30 ns. For the analysis of backbone flotation, the root mean square fluctuations were plotted for the Cα atoms; analysis revealed that the greatest fluctuation occurred in the residues that are part of the protein lid domain. The binding free energy value (ΔGbind) obtained for the Tam16 lead molecule was of −51.43 kcal/mol. When comparing this result with the ΔGbind values for the remaining analogs, the drug Tam16 was found to be the highest ranked: this result is in agreement with the experimental results obtained by Aggarwal and collaborators, where it was verified that the IC50 for Tam16 is the smallest necessary to inhibit the Pks13 (IC50 = 0.19 μM). The energy decomposition analysis suggested that the residues which most interact with inhibitors are: Ser1636, Tyr1637, Asn1640, Ala1667, Phe1670, and Tyr1674, from which the greatest energy contribution to Phe1670 was particularly notable. For the lead molecule Tam16, a hydrogen bond with the hydroxyl of the phenol not observed in the other analogs induced a more stable molecular structure. Aggarwal and colleagues reported this hydrogen bonding as being responsible for the stability of the molecule, optimizing its physic-chemical, toxicological, and pharmacokinetic properties.

Role of protein structure and the role of individual fingers in zinc finger protein–DNA recognition: a molecular dynamics simulation study and free energy calculations

Molecular dynamics and MM_GBSA energy calculations on various zinc finger proteins containing three and four fingers bound to their target DNA gave insights into the role of each finger in the DNA binding process as part of the protein structure. The wild type Zif 268 (PDB code: 1AAY) gave a ΔG value of - 76.1 (14)kcal/mol. Zinc fingers ZF1, ZF2 and ZF3 were mutated in one experiment and in another experiment one finger was cut and the rest of the protein was studied for binding. The ΔΔG values for the Zinc Finger protein with both ZF1 and ZF2 mutated was + 80kcal/mol, while mutating only ZF1 the ΔΔG value was + 52kcal/mol (relative to the wild type). Cutting ZF3 and studying the protein consisting only of ZF1 linked to ZF2 gave a ΔΔG value of + 68kcal/mol. Upon cutting ZF1, the resulting ZF2 linked to ZF3 protein gave a ΔΔG value of + 41kcal/mol. The above results shed light on the importance of each finger in the binding process, especially the role of ZF1 as the anchoring finger followed in importance by ZF2 and ZF3. The energy difference between the binding of the wild type protein Zif268 (1AAY) and that for individual finger binding to DNA according to the formula: ΔΔGlinkers, otherstructuralfactors = ΔGzif268 - (ΔGF1+F2+F3) gave a value = - 44.5kcal/mol. This stabilization can be attributed to the contribution of linkers and other structural factors in the intact protein in the DNA binding process. DNA binding energies of variant proteins of the wild type Zif268 which differ in their ZF1 amino acid sequence gave evidence of a good relationship between binding energy and recognition and specificity, this finding confirms the reported vital role of ZF1 in the ZF protein scanning and anchoring to the target DNA sequence. The role of hydrogen bonds in both specific and nonspecific amino acid-DNA contacts is discussed in relation to mutations. The binding energies of variant Zinc Finger proteins confirmed the role of ZF1 in the recognition, specificity and anchoring of the zinc finger protein to DNA.

Spontaneous insertion of GPI anchors into cholesterol-rich membrane domains

GPI-Anchored proteins (GPI-APs) can be exogenously transferred onto bilayer membranes both in vivo and in vitro, while the mechanism by which this transfer process occurs is unknown. In this work, we used atomistic molecular dynamics simulations and free energy calculations to characterize the essential influence of cholesterol on insertion of the GPI anchors into plasma membranes. We demonstrate, both dynamically and energetically, that in the presence of cholesterol, the tails of GPI anchors are able to penetrate inside the core of the lipid membrane spontaneously with a three-step mechanism, while in the absence of cholesterol no spontaneous insertion was observed. We ascribe the failure of insertion to the strong thermal fluctuation of lipid molecules in cholesterol-free bilayer, which generates a repulsive force in entropic origin. In the presence of cholesterol, however, the fluctuation of lipids is strongly reduced, thus decreasing the barrier for the anchor insertion. Based on this observation, we propose a hypothesis that addition of cholesterol creates vertical creases in membranes for the insertion of acyl chains. Moreover, we find that the GPI anchor could also spontaneously inserted into the boundary between cholesterol-rich and cholesterol-depleted domains. Our results shed light on the mechanism of cholesterol-mediated interaction between membrane proteins with acyl chain and plasma membranes in living cells.

Energetics of Preferential Binding of Retinoic Acid-Inducible Gene-I to Double-Stranded Viral RNAs with 5′ Tri-/Diphosphate over 5′ Monophosphate

Retinoic acid-inducible gene-I (RIG-I) is a cytosolic sensor protein that recognizes viral RNAs and triggers an innate immune response in cells. Panhandle-like base-paired blunt-ended 5′ ppp/pp-dsRNA is a characteristic feature of viral RNAs. Structural studies of RIG-I C-terminal domain bound 5′ ppp/pp-dsRNA complexes show the direct interaction between all the 5′ terminal phosphates (α, β, and γ) and protein, suggesting γ phosphate might be a major recognition determinant for RIG-I binding. Biochemical studies, however, suggest that 5′ pp-dsRNA is the minimal determinant for RIG-I binding and antiviral response. Despite biochemical and structural studies, the origin of viral RNA recognition by RIG-I is an unsolved problem. X-ray structures of RIG-I bound dsRNA not only provide atomic insight into the interaction network but also provide sufficiently good models for computational studies. We report structure-based molecular dynamics (MD) free energy calculations to quantitatively estimate the energetics of RIG-I binding to dsRNA containing 5′ ppp, 5′ pp, and 5′ p. The results suggest that RIG-I weakly discriminates between 5′ ppp-dsRNA and 5′ pp-dsRNA (favoring former) and strongly disfavors 5′ p-dsRNA with respect to the rest. Interestingly, direct interaction between γ phosphate of 5′ ppp-dsRNA and RIG-I is a robust feature of the MD simulations. dsRNA binding to RIG-I is associated with Mg2+ dissociation from the 5′ phosphate/s of dsRNA. The higher Mg2+ dissociation penalty from 5′ ppp-dsRNA with respect to 5′ pp-dsRNA offsets most of the favorable interaction between RIG-I and γ phosphate of 5′ ppp-dsRNA. This leads to weak discrimination between 5′ ppp-dsRNA and 5′ pp-dsRNA. 5′ p-dsRNA is discriminated strongly because of the loss of interaction with RIG-I.

Role of the Bound Phospholipids in the Structural Stability of Cholesteryl Ester Transfer Protein.

Cholesteryl ester transfer protein (CETP) facilitates the transfer of cholesteryl esters (CEs) from antiatherogenic high-density lipoproteins to proatherogenic low-density lipoproteins. Inhibition of CETP is therefore being pursued as a potential strategy to reduce cardiovascular risk. The crystal structure of CETP has revealed the existence of two neutral CEs and two charged phospholipids (PLs) in its hydrophobic tunnel. This is in direct contrast to the other lipid-binding proteins that contain only two bound lipids. Moreover, previous animal studies on mice showed no detectable PL-transfer activity of CETP. Thus, the role of bound PLs in CETP is completely unknown. Here, we employ molecular dynamics simulations and free-energy calculations to unravel the primary effects of bound PLs on CETP structure and dynamics and attempt to correlate the observed changes to its function. Our results suggest that the structure of CETP is elastic and can attain different conformations depending on the state of bound PLs. In solution, these PLs maintain CETP in a bent-untwisted conformation that can uphold neutral lipids in its core tunnel. Results also suggest that although both PLs complement each other in their action, the C-terminal PL (C-PL) imparts greater influence on CETP by virtue of its tighter binding. Our finding fits very well with the recent inhibitor-bound CETP crystal structure, where the inhibitor displaced the N-terminal PL for binding to CETP's central domain without disrupting the binding of C-PL. We speculate that the observed increased flexibility of CETP in the absence of PLs could play a crucial role in its binding with lipoproteins and subsequent lipid-transfer activity.

Ion Channel Sensing: Are Fluctuations the Crux of the Matter?

The nonselective cation channel TRPV1 is responsible for transducing noxious stimuli into action potentials propagating through peripheral nerves. It is activated by temperatures greater than 43 °C, while remaining completely nonconductive at temperatures lower than this threshold. The origin of this sharp response, which makes TRPV1 a biological temperature sensor, is not understood. Here we used molecular dynamics simulations and free energy calculations to characterize the molecular determinants of the transition between nonconductive and conductive states. We found that hydration of the pore and thus ion permeation depends critically on the polar character of its molecular surface: in this narrow hydrophobic enclosure, the motion of a polar side-chain is sufficient to stabilize either the dry or wet state. The conformation of this side-chain is in turn coupled to the hydration state of four peripheral cavities, which undergo a dewetting transition at the activation temperature.

Effect of Binding on Enantioselectivity of Epoxide Hydrolase.

Molecular dynamics simulations and free energy calculations have been used to investigate the effect of ligand binding on the enantioselectivity of an epoxide hydrolase (EH) from Aspergillus niger. Despite sharing a common mechanism, a wide range of alternative mechanisms have been proposed to explain the origin of enantiomeric selectivity in EHs. By comparing the interactions of ( R)- and ( S)-glycidyl phenyl ether (GPE) with both the wild type (WT, E = 3) and a mutant showing enhanced enantioselectivity to GPE (LW202, E = 193), we have examined whether enantioselectivity is due to differences in the binding pose, the affinity for the ( R)- or ( S)- enantiomers, or a kinetic effect. The two enantiomers were easily accommodated within the binding pockets of the WT enzyme and LW202. Free energy calculations suggested that neither enzyme had a preference for a given enantiomer. The two substrates sampled a wide variety of conformations in the simulations with the sterically hindered and unhindered carbon atoms of the GPE epoxide ring both coming in close proximity to the nucleophilic aspartic acid residue. This suggests that alternative pathways could lead to the formation of a ( S)- and ( R)-diol product. Together, the calculations suggest that the enantioselectivity is due to kinetic rather than thermodynamic effects and that the assumption that one substrate results in one product when interpreting the available experimental data and deriving E-values may be inappropriate in the case of EHs.

Nucleotide-dependent Conformational Dynamics and Energetics of Tubulin

Microtubules (MTs) are key components of the cytoskeleton and play a central role in cell division and development. They undergo stochastic switching between phases of growth and shrinkage driven by GTP hydrolysis in αβ-tubulin dimers. MT assembly requires that tubulin adopts a straight conformation in the lattice, while during disassembly it adopts a kinked conformation that is incompatible with the MT geometry. The mechanism by which GTP binding renders individual tubulin dimers polymerization-competent is still poorly understood. In particular, two models of MT assembly are being controversially discussed: in the allosteric model, GTP binding directly promotes tubulin straightening, whereas in the lattice model tubulin is kinked regardless of its nucleotide state, and it is the assembling MT itself that forces tubulin into a straight conformation. In order to distinguish between these two hypotheses, we have characterized the conformational dynamics and energetics of free tubulin in solution using atomistic molecular dynamics simulations and free energy calculations. Unexpectedly, and contradictory to the above models, we find that GTP-tubulin samples a broad range of intrinsic curvatures all of which are almost isoenergetic. In contrast, GDP-tubulin is confined to low-curvature regions of the conformational space. Furthermore, an additional energy barrier is seen in all simulated tubulin ensembles that is associated with collective rearrangements in the β-subunit of tubulin but not with the nucleotide state. This barrier separates different bending modi and may be related to the inhibitory mechanism of MT-associated proteins such as stathmin which suppresses tubulin polymerization in a nucleotide-independent manner. Overall, our results suggest a new combined conformational selection / induced fit model in which GTP binding renders tubulin more flexible at the intadimer interface and allows a specific type of tubulin bending that facilitates MT assembly.

Alanine mutation of the catalytic sites of Pantothenate Synthetase causes distinct conformational changes in the ATP binding region

The enzyme Pantothenate synthetase (PS) represents a potential drug target in Mycobacterium tuberculosis. Its X-ray crystallographic structure has demonstrated the significance and importance of conserved active site residues including His44, His47, Asn69, Gln72, Lys160 and Gln164 in substrate binding and formation of pantoyl adenylate intermediate. In the current study, molecular mechanism of decreased affinity of the enzyme for ATP caused by alanine mutations was investigated using molecular dynamics (MD) simulations and free energy calculations. A total of seven systems including wild-type + ATP, H44A + ATP, H47A + ATP, N69A + ATP, Q72A + ATP, K160A + ATP and Q164A + ATP were subjected to 50 ns MD simulations. Docking score, MM-GBSA and interaction profile analysis showed weak interactions between ATP (substrate) and PS (enzyme) in H47A and H160A mutants as compared to wild-type, leading to reduced protein catalytic activity. However, principal component analysis (PCA) and free energy landscape (FEL) analysis revealed that ATP was strongly bound to the catalytic core of the wild-type, limiting its movement to form a stable complex as compared to mutants. The study will give insight about ATP binding to the PS at the atomic level and will facilitate in designing of non-reactive analogue of pantoyl adenylate which will act as a specific inhibitor for PS.

The Interaction Modes of Haloimidazolium Salts in Solution.

We performed a comparative study on the interaction modes of 2-haloimidazolium salts with anions in solution, particularly with regard to halogen bonding, hydrogen bonding and anion-π interactions. The syntheses and solid-state analyses of a series of sterically and electronically modified 2-haloimidazolium structures are presented. Detailed isothermal titration calorimetry (ITC) measurements, quantum mechanics/molecular mechanics (QM/MM), classical molecular dynamics simulations (MD) and free-energy calculations together with NMR spectroscopy were used to elucidate the binding modes in solution. Our work reveals the absence of a potential anion-π interaction between the cationic imidazolium ring and the Lewis basic counteranion, and corroborates a formation of halogen bonding via the Lewis acidic iodine moiety and hydrogen bonding via the backbone hydrogen atoms, with repercussions in the field of organocatalysis.

Flavonoids as potent allosteric inhibitors of protein tyrosine phosphatase 1B: molecular dynamics simulation and free energy calculation

Protein tyrosine phosphatase 1B (PTP1B) is a member of the PTP superfamily which is considered to be a negative regulator of insulin receptor (IR) signaling pathway. PTP1B is a promising drug target for the treatment of type 2 diabetes, obesity, and cancer. The existence of allosteric site in PTP1B has turned the researcher’s attention to an alternate strategy for inhibition of this enzyme. Herein, the molecular interactions between the allosteric site of PTP1B with three non-competitive flavonoids, (MOR), (MOK), and (DPO) have been investigated. Three ligands were docked into allosteric site of the enzyme. The resulting protein–ligand complexes were used for molecular dynamics studies. Principal component and free-energy landscape (FEL) as well as cluster analyses were used to investigate the conformational and dynamical properties of the protein and identify representative enzyme substrates bounded to the inhibitors. Per residue energy decomposition analysis attributed dissimilar affinities of three inhibitors to the several hydrogen bonds and non-bonded interactions. In conclusion, our results exhibited an inhibitory pattern of the ligands against PTP1B.