Effect of magnesium addition on the vacancy-Mg binding free enthalpy in Al-Cr alloys

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Enthalpy/entropy compensation: Influence of DNA flanking sequence on the binding of 7-amino actinomycin D to its primary binding site in short DNA duplexes

The effect of the context of the flanking sequence on ligand binding to DNA oligonucleotides that contain consensus binding sites was investigated for the binding of the intercalator 7-amino actinomycin D. Seven self-complementary DNA oligomers each containing a centrally located primary binding site, 5'-A-G-C-T-3', flanked on either side by the sequences (AT)(n) or (AA)(n) (with n = 2, 3, 4) and AA(AT)(2), were studied. For different flanking sequences, (AA)(n)-series or (AT)(n)-series, differential fluorescence enhancements of the ligand due to binding were observed. Thermodynamic studies indicated that the flanking sequences not only affected DNA stability and secondary structure but also modulated ligand binding to the primary binding site. The magnitude of the ligand binding affinity to the primary site was inversely related to the sequence dependent stability. The enthalpy of ligand binding was directly measured by isothermal titration calorimetry, and this made it possible to parse the binding free energy into its energetic and entropic terms.

Computational study of ligand binding in lipid transfer proteins: Structures, interfaces, and free energies of protein-lipid complexes

Plant nonspecific lipid transfer proteins (nsLTPs) bind a wide variety of lipids, which allows them to perform disparate functions. Recent reports on their multifunctionality in plant growth processes have posed new questions on the versatile binding abilities of these proteins. The lack of binding specificity has been customarily explained in qualitative terms on the basis of a supposed structural flexibility and nonspecificity of hydrophobic protein-ligand interactions. We present here a computational study of protein-ligand complexes formed between five nsLTPs and seven lipids bound in two different ways in every receptor protein. After optimizing geometries inmolecular dynamics calculations, we computed Poisson- Boltzmann electrostatic potentials, solvation energies, properties of the protein-ligand interfaces, and estimates of binding free energies of the resulting complexes. Our results provide the first quantitative information on the ligand abilities of nsLTPs, shed new light into protein-lipid interactions, and reveal new features which supplement commonly held assumptions on their lack of binding specificity.

Metal-free graphitic carbon nitride as mechano-catalyst for hydrogen evolution reaction

Graphitic carbon nitride (g-C3N4), as a promising metal-free catalyst for photo-catalytic and electrochemical water splitting, has recently attracted tremendous research interest. However, the underlying catalytic mechanism for the hydrogen evolution reaction (HER) is not fully understood. By using density functional theory calculations, here we have established that the binding free energy of hydrogen atom (ΔGH∗0) on g-C3N4 is very sensitive to mechanical strain, leading to substantial tuning of the HER performance of g-C3N4 at different coverages. The experimentally-observed high HER activity in N-doped graphene supported g-C3N4 (Zheng et al., 2014) is actually attributed to electron-transfer induced strain. A more practical strategy to induce mechanical strain in g-C3N4 is also proposed by doping a bridge carbon atom in g-C3N4 with an isoelectronic silicon atom. The calculated ΔGH∗0 on the Si-doped g-C3N4 is ideal for HER. Our results indicate that g-C3N4 would be an excellent metal-free mechano-catalyst for HER and this finding is expected to guide future experiments to efficiently split water into hydrogen based on the g-C3N4 materials.

Clustering and solute-vacancy binding energies in Al-4.4% Zn alloy with 0.01% ternary additions

Isochronal and isothermal ageing experiments have been carried out to determine the influence of 0.01 at. % addition of a second solute on the clustering rate in the quenched Al-4,4 a/o Zn alloy. The influence of quenching and ageing temperatures has been interpreted to obtain the apparent vacancy formation and vacancy migration energies in the various ternary alloys. Using a vacancy-aided clustering model the following values of binding free energy have been evaluated: Ce-0.18; Dy-0.24; Fe-0.18; Li-0.25; Mn-0.27; Nb-0.18; Pt-0.23; Sb-0.21; Si-0.30; Y-0.25; and Yb-0.23 (± 0.02 eV). These binding energy values refer to that between a solute atom and a single vacancy. The values of vacancy migration energy (c. 0.4 eV) and the experimental activation energy for solute diffusion (c. 1.1 eV) are unaffected by the presence of the ternary atoms in the Al-Zn alloy.

Docking and free energy simulations to predict conformational domains involved in hCG-LH receptor interactions using recombinant antibodies

Single chain fragment variables (ScFvs) have been extensively employed in studying the protein-protein interactions. ScFvs derived from phage display libraries have an additional advantage of being generated against a native antigen, circumventing loss of information on conformational epitopes. In the present study, an attempt has been made to elucidate human chorionic gonadotropin (hCG)-luteinizing hormone (LH) receptor interactions by using a neutral and two inhibitory ScFvs against hCG. The objective was to dock a computationally derived model of these ScFvs onto the crystal structure of hCG and understand the differential roles of the mapped epitopes in hCG-LH receptor interactions. An anti-hCG ScFv, whose epitope was mapped previously using biochemical tools, served as the positive control for assessing the quality of docking analysis. To evaluate the role of specific side chains at the hCG-ScFv interface, binding free energy as well as residue interaction energies of complexes in solution were calculated using molecular mechanics Poisson-Boltzmann/surface area method after performing the molecular dynamic simulations on the selected hCG-ScFv models and validated using biochemical and SPR analysis. The robustness of these calculations was demonstrated by comparing the theoretically determined binding energies with the experimentally obtained kinetic parameters for hCG-ScFv complexes. Superimposition of hCG-ScFv model onto a model of hCG complexed with the 51-266 residues of LH receptor revealed importance of the residues previously thought to be unimportant for hormone binding and response. This analysis provides an alternate tool for understanding the structure-function analysis of ligand-receptor interactions. Proteins 2011;79:3108-3122. (C) 2011 Wiley-Liss, Inc.

Binding of BIS like and other ligands with the GSK-3 beta kinase: a combined docking and MM-PBSA study

Binding of several bisindolylmaleimide (BIS) like (BIS-3, BIS-8 and UCN1) and other ligands (H89, SB203580 and Y27632) with the glycogen synthase kinase-3 (GSK-3 beta) has been studied using combined docking, molecular dynamics and Poisson-Boltzmann surface area analysis approaches. The study generated novel binding modes of these ligands that can rationalize why some ligands inhibit GSK-3 beta while others do not. The relative binding free energies associated with these binding modes are in agreement with the corresponding measured specificities. This study further provides useful insight regarding possible existence of multiple conformations of some ligands like H89 and BIS-8. It is also found that binding modes of BIS-3, BIS-8 and UCN1 with GSK-3 beta and PDK1 kinases are similar. These new insights are expected to be useful for future rational design of novel, more potent GSK-3 beta-specific inhibitors as promising therapeutics.

Energetic basis of molecular recognition in a DNA aptamer

The thermal stability and ligand binding properties of the L-argininamide-binding DNA aptamer (5'-GATCGAAACGTAGCGCCTTCGATC3') were studied by spectroscopic and calorimetric methods. Differential calorimetric studies showed that the uncomplexed aptamer melted in a two-state reaction with a melting temperature T-m = 50.2 +/- 0.2 degrees C and a folding enthalpy Delta H degrees(fold) = -49.0 +/- 2.1 kcal mol(-1). These values agree with values of T-m = 49.6 degrees C and Delta H degrees(fold) = -51.2 kcal mol(-1) predicted for a simple hairpin structure. Melting of the uncomplexed aptamer was dependent upon salt concentration, but independent of strand concentration. The T of aptamer melting was found to increase as L-argininamide concentrations increased. Analysis of circular dichroism titration data using a single-site binding model resulted in the determination of a binding free energy Delta G degrees(bind) = -5.1 kcal mol(-1). Isothermal titration calorimetry studies revealed an exothermic binding reaction with Delta H degrees(bind) = -8.7 kcal mol(-1). Combination of enthalpy and free energy produce ail unfavorable entropy of -T Delta S degrees = +3.6 kcal mol(-1). A molar heat capacity change of -116 cal mol(-1) K-1 was determined from calorimetric measurements at four temperatures over the range of 15-40 degrees C. Molecular dynamics simulations were used to explore the structures of the unligated and ligated aptamer structures.

Docking and molecular dynamics studies toward the binding of new natural phenolic marine inhibitors and aldose reductase

Phenolic marine natural product is a kind of new potential aldose reductase inhibitors (ARIs). In order to investigate the binding mode and inhibition mechanism, molecular docking and dynamics studies were performed to explore the interactions of six phenolic inhibitors with human aldose reductase (hALR2). Considering physiological environment, all the neutral and other two ionized states of each phenolic inhibitor were adopted in the simulation. The calculations indicate that all the inhibitors are able to form stable hydrogen bonds with the hALR2 active pocket which is mainly constructed by residues TYR48, HIS110 and TRP111, and they impose the inhibition effect by occupying the active space. In all inhibitors, only La and its two ionized derivatives La_ion1 and La_ion2, in which neither of the ortho-hydrogens of 3-hydroxyl is substituted by Br, bind with hALR2 active residues using the terminal 3-hydroxyl. While, all the other inhibitors, at least one of whose ortho-sites of 3- and 6-hydroxyls are substituted by Br substituent which take much electron-withdrawing effect and steric hindrance, bind with hALR2 through the lactone group. This means that the Br substituent can effectively regulate the binding modes of phenolic inhibitors. Although the lactone bound inhibitors have relatively high RMSD values, our dynamics study shows that both binding modes are of high stability. For each inhibitor molecule, the ionization does not change its original binding mode, but it does gradually increase the binding free energy, which reveals that besides hydrogen bonds, the electrostatic effect is also important to the inhibitor–hALR2 interaction.

Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.

Binding studies of a protonated dioxatetraazamacrocycle with carboxylate substrates

The tetraprotonated form of the dioxatetraazamacrocycle, 6,19-dioxa-3,9,16,22-tetraaza[22.2.2.2(11,14)]-triaconta-1(26),11,13,24, 27,29-hexaene, (H4L1)(4+), was used as the receptor for binding studies with carboxylate anionic substrates of different shapes, sizes, and charges [succinate (suc(2-)), cyclo- hexanetricarboxylate (cta(3-)), phthalate (ph(2-)), isophthalate (iph(2-)), terephthalate (tph(2-)), and benezenetricarboxylate (btc(3-))]. Association constants were determined by potentiometry in aqueous solution at 298.2 K and 0.10 M KCl and by H-1 NMR titration in D2O. The strongest association was found for the btc3- anion at 5-7 pH region. From both techniques it was possible to establish the binding preference trend of the receptor for the different substrates, and the H-1 NMR spectroscopy gave important suggestions about the type of interactions between partners and the location of the substrates in the supramolecular entities formed. The effective binding constants at pH 6 follow the order: btc(3-)>iph(2-)>cta(3-) =ph(2-)>tph(2-)>suc(2-). All the studies suggest that the anionic substrates bind to the receptor via N-H center dot center dot center dot O = C hydrogen bonds and electrostatic interactions, and the aromatic substrates can also establish pi-pi stacking interactions. The crystal structures of (H4L1)(4+) and its supramolecular assemblies with ph(2-) and tph(2-) were determined by X-ray diffraction. The last two structures showed that the association process in solid state occurs via multiple N-H center dot center dot center dot O = C hydrogen bonds with the anionic substrate located outside the macrocyclic cavity of the receptor. Molecular dynamics simulations carried out for the association of (H4L1)(4+) with tph(2-) and btC(3-) in water solution established at atomic level the existence of all interactions suggested by the experimental studies, which act cooperatively in the binding process. Furthermore, the binding free energies were estimated and the values are in agreement with the experimental ones, indicating that the binding of these two anionic substrates occurs into the receptor cavity. However, the tph(2-) has also propensity to leave the macrocyclic cavity and its molecular recognition can also happen at the top of the receptor. (C) 2008 Elsevier Ltd. All rights reserved.

Binding of an oligomeric ellagitannin series to bovine serum albumin (BSA): analysis by isothermal titration calorimetry (ITC)

A unique series of oligomeric ellagitannins was used to study their interactions with bovine serum albumin (BSA) by isothermal titration calorimetry. Oligomeric ellagitannins, ranging from monomer to heptamer and a mixture of octamer–undecamers, were isolated as individual pure compounds. This series allowed studying the effects of oligomer size and other structural features. The monomeric to trimeric ellagitannins deviated most from the overall trends. The interactions of ellagitannin oligomers from tetramers to octa–undecamers with BSA revealed strong similarities. In contrast to the equilibrium binding constant, enthalpy showed an increasing trend from the dimer to larger oligomers. It is likely that first the macrocyclic part of the ellagitannin binds to the defined binding sites on the protein surface and then the “flexible tail” of the ellagitannin coats the protein surface. The results highlight the importance of molecular flexibility to maximize binding between the ellagitannin and protein surfaces.

Flexibility and inhibitor binding in Cdc25 phosphatases

Cdc25 phosphatases involved in cell cycle checkpoints are now active targets for the development of anti-cancer therapies. Rational drug design would certainly benefit from detailed structural information for Cdc25s. However, only apo- or sulfate-bound crystal structures of the Cdc25 catalytic domain have been described so far. Together with previously available crystalographic data, results from molecular dynamics simulations, bioinformatic analysis, and computer-generated conformational ensembles shown here indicate that the last 30-40 residues in the C-terminus of Cdc25B are partially unfolded or disordered in solution. The effect of C-terminal flexibility upon binding of two potent small molecule inhibitors to Cdc25B is then analyzed by using three structural models with variable levels of flexibility, including an equilibrium distributed ensemble of Cdc25B backbone conformations. The three Cdc25B structural models are used in combination with flexible docking, clustering, and calculation of binding free energies by the linear interaction energy approximation to construct and validate Cdc25B-inhibitor complexes. Two binding sites are identified on top and beside the Cdc25B active site. The diversity of interaction modes found increases with receptor flexibility. Backbone flexibility allows the formation of transient cavities or compact hydrophobic units on the surface of the stable, folded protein core that are unexposed or unavailable for ligand binding in rigid and densely packed crystal structures. The present results may help to speculate on the mechanisms of small molecule complexation to partially unfolded or locally disordered proteins.

Path Optimization in Free Energy Calculations

Free energy calculations are a computational method for determining thermodynamic quantities, such as free energies of binding, via simulation.

Currently, due to computational and algorithmic limitations, free energy calculations are limited in scope.

In this work, we propose two methods for improving the efficiency of free energy calculations.

First, we expand the state space of alchemical intermediates, and show that this expansion enables us to calculate free energies along lower variance paths.

We use Q-learning, a reinforcement learning technique, to discover and optimize paths at low computational cost.

Second, we reduce the cost of sampling along a given path by using sequential Monte Carlo samplers.

We develop a new free energy estimator, pCrooks (pairwise Crooks), a variant on the Crooks fluctuation theorem (CFT), which enables decomposition of the variance of the free energy estimate for discrete paths, while retaining beneficial characteristics of CFT.

Combining these two advancements, we show that for some test models, optimal expanded-space paths have a nearly 80% reduction in variance relative to the standard path.

Additionally, our free energy estimator converges at a more consistent rate and on average 1.8 times faster when we enable path searching, even when the cost of path discovery and refinement is considered.

Understanding Co-translational Protein Targeting and Lithium Dendrite Formation through Free Energy Simulations and Coarse-Grained Models

We describe the application of alchemical free energy methods and coarse-grained models to study two key problems: (i) co-translational protein targeting and insertion to direct membrane proteins to the endoplasmic reticulum for proper localization and folding, (ii) lithium dendrite formation during recharging of lithium metal batteries. We show that conformational changes in the signal recognition particle, a central component of the protein targeting machinery, confer additional specificity during the the recognition of signal sequences. We then develop a three-dimensional coarse-grained model to study the long-timescale dynamics of membrane protein integration at the translocon and a framework for the calculation of binding free energies between the ribosome and translocon. Finally, we develop a coarse-grained model to capture the dynamics of lithium deposition and dissolution at the electrode interface with time-dependent voltages to show that pulse plating and reverse pulse plating methods can mitigate dendrite growth.