Inhibition of mugineic acid-ferric complex uptake in barley by copper, zinc and cobalt

The interfering effects of copper, zinc, and cobalt on the uptake of mugineic acid-ferric complex were studied in barley ( Hordeum vulgare , cv. Minorimugi) grown in nutrient solution. Short-term uptake experiments of 3 h were performed utilizing both ionic and mugineic acid-complex forms of each metal at two different concentrations. Copper was most effective in decreasing iron uptake when added in an ionic form at either concentration. The inhibition order at higher concentrations followed Cu(II) > Zn(II) ≥ Co(II), Co(III), which is consistent with the stability constants of these metal complexes with mugineic acid. The displacement of iron from its mugineic acid complex by these metals is suggested as a probable explanation for the decreased iron uptake. The inhibitory effect of metal complexes with mugineic acid on iron uptake was only found in cases with higher concentrations of Cu(II) and Zn(II) complexes. Deformation of the specific iron transport system in the plasma membrane due to their adsorption may be responsible for this effect.     

DFT study of the V(IV)/V(V) oxidation mechanism in the presence of N-hydroxyacetamide

The oxidation mechanism of V(IV)/V(V) in the presence of N-hydroxyacetamide (acetohydroxamic acid, HL) in aqueous solution has been investigated using density functional theory (DFT) calculations aiming to contribute to the understanding of this process at a molecular level. The mechanism proposed involves formation of the *OH, *OOH, H2O2 radicals and complexes formed from the interaction of these species with VOL2 complex. The Gibbs free energy of each step of the mechanism has been evaluated. The solvation energy has been estimated by means of united atoms Hartree-Fock/polarizable continuum method (UAHF/PCM). The Gibbs free energy of the global reaction of the V(IV)/V(V) oxidation has been estimated and compared with the available experimental equilibrium constant. The difference between the calculated and experimental estimates for the reaction energy of the global equation is about 1.5 kcal mol(-1). The thermodynamic profile of the reaction mechanism has been provided and discussed in terms of the possible intermediates. The influence of the ligand and the reaction rate in terms of the steady-state approximation has been briefly discussed.     

Theoretical investigations of the experimentally observed selectivity of a cobalt imprinted polymer

A cobalt imprinted polymer synthesised, for reducing the volume of radioactive waste generated during nuclear reactor decontaminations, using vinylbenzyl iminodiacetate (VbIDA) as the functional ligand, has been found to be selective for cobaltous ions over excess ferrous ions. The selectivity of the polymer has been investigated through theoretical calculation of the formation energies of complexes involved by using the ab-initio density functional theory (DFT) code SIESTA (Spanish Initiative for Electronic Simulations with Thousands of Atoms). The formation energies of complexes of Fe²⁺, Co²⁺, Cu²⁺ and Ni²⁺ with the free functional ligands as well as with ligands attached to the crosslinkers have been calculated. The calculations revealed that the ferrous forms an unstable complex with the ligands attached to the crosslinkers. The formation energy calculation results were found to corroborate the experimentally observed selectivity order.     

Hydroxylated phytosiderophore species possess an enhanced chelate stability and affinity for iron(III)

Graminaceous plant species acquire soil iron by the release of phytosiderophores and subsequent uptake of iron(III)-phytosiderophore complexes. As plant species differ in their ability for phytosiderophore hydroxylation prior to release, an electrophoretic method was set up to determine whether hydroxylation affects the net charge of iron(III)-phytosiderophore complexes, and thus chelate stability. At pH 7.0, non-hydroxylated (deoxymugineic acid) and hydroxylated (mugineic acid; epi-hydroxymugineic acid) phytosiderophores form single negatively charged iron(III) complexes, in contrast to iron(III)-nicotianamine. As the degree of phytosiderophore hydroxylation increases, the corresponding iron(III) complex was found to be less readily protonated. Measured pKa values of the amino groups and calculated free iron(III) concentrations in presence of a 10-fold chelator excess were also found to decrease with increasing degree of hydroxylation, confirming that phytosiderophore hydroxylation protects against acid-induced protonation of the iron(III)-phytosiderophore complex. These effects are almost certainly associated with intramolecular hydrogen bonding between the hydroxyl and amino functions. We conclude that introduction of hydroxyl groups into the phytosiderophore skeleton increases iron(III)-chelate stability in acid environments such as those found in the rhizosphere or the root apoplasm and may contribute to an enhanced iron acquisition.     

Density functional theory study on the interaction between keto-9H guanine and aspartic acid

A theoretical study was performed using density functional theory (DFT) to investigate hydrogen bonding interactions in signature complexes formed between keto-9H guanine (Gua) and aspartic acid (Asp) at neutral pH. Optimized geometries, binding energies and the theoretical IR spectra of guanine, aspartic acid and their corresponding complexes (Gua-Asp) were calculated using the B3LYP method and the 6-31+G(d) basis set. Stationary points found to be at local minima on the potential energy surface were verified by second derivative harmonic vibrational frequency calculations at the same level of theory. AIM theory was used to analyze the hydrogen bonding characteristics of these DNA base complex systems. Our results show that the binding motif for the most stable complex is strikingly similar to a Watson-Crick motif observed in the guanine-cytosine base pair. We have found a range of hydrogen bonding interactions between guanine and aspartic acid in the six complexes. This was further verified by theoretical IR spectra of ω(C-H---O-H) cm⁻¹ stretches for the Gua-Asp complexes. The electron density plot indicates strong hydrogen bonding as shown by the 2p _z dominant HOMO orbital character.     

Empirical calculation of the relative free energies of peptide binding to the molecular chaperone DnaK

We describe a methodology to calculate the relative free energies of protein-peptide complex formation. The interaction energy was decomposed into nonpolar, electrostatic and entropic contributions. A free energy-surface area relationship served to calculate the nonpolar free energy term. The electrostatic free energy was calculated with the finite difference Poisson-Boltzmann method and the entropic contribution was estimated from the loss in the conformational entropy of the peptide side chains. We applied this methodology to a series of DnaK*peptide complexes. On the basis of the single known crystal structure of the peptide-binding domain of DnaK with a bound heptapeptide, we modeled ten other DnaK*heptapeptide complexes with experimentally measured K(d) values from 0.06 microM to 11 microM, using molecular dynamics to refine the structures of the complexes. Molecular dynamic trajectories, after equilibration, were used for calculating the energies with greater accuracy. The calculated relative binding free energies were compared with the experimentally determined free energies. Linear scaling of the calculated terms was applied to fit them to the experimental values. The calculated binding free energies were between -7.1 kcal/mol and - 9.4 kcal/mol with a correlation coefficient of 0.86. The calculated nonpolar contributions are mainly due to the central hydrophobic binding pocket of DnaK for three amino acid residues. Negative electrostatic fields generated by the protein increase the binding affinity for basic residues flanking the hydrophobic core of the peptide ligand. Analysis of the individual energy contributions indicated that the nonpolar contributions are predominant compared to the other energy terms even for peptides with low affinity and that inclusion of the change in conformational entropy of the peptide side chains does not improve the discriminative power of the calculation. The method seems to be useful for predicting relative binding energies of peptide ligands of DnaK and might be applicable to other protein-peptide systems, particularly if only the structure of one protein-ligand complex is available.     

The binding of aluminum to mugineic acid and related compounds as studied by potentiometric titration

The phytosiderophores, mugineic acid (MA) and epi-hydroxymugineic acid (HMA), together with a related compound, nicotianamine (NA), were investigated for their ability to bind Al(III). Potentiometric titration analysis demonstrated that MA and HMA bind Al(III), in contrast to NA which does not under normal physiological conditions. With MA and HMA, in addition to the Al complex (AlL), the protonated (AlLH) and deprotonated (AlLH₋₁) complexes were identified from an analysis of titration curves, where L denotes the phytosiderophore form in which all the carboxylate functions are ionized. The equilibrium formation constants of the Al(III) phytosiderophore complexes are much smaller than those of the corresponding Fe(III) complexes. The higher selectivity of phytosiderophores for Fe(III) over Al(III) facilitates Fe(III) acquisition in alkaline conditions where free Al(III) levels are higher than free Fe(III) levels.     

DNA-lipid complexes: stability of honeycomb-like and spaghetti-like structures.

A molecular level theory is presented for the thermodynamic stability of two (similar) types of structural complexes formed by (either single strand or supercoiled) DNA and cationic liposomes, both involving a monolayer-coated DNA as the central structural unit. In the "spaghetti" complex the central unit is surrounded by another, oppositely curved, monolayer, thus forming a bilayer mantle. The "honeycomb" complex is a bundle of hexagonally packed DNA-monolayer units. The formation free energy of these complexes, starting from a planar cationic/neutral lipid bilayer and bare DNA, is expressed as a sum of electrostatic, bending, mixing, and (for the honeycomb) chain frustration contributions. The electrostatic free energy is calculated using the Poisson-Boltzmann equation. The bending energy of the mixed lipid layers is treated in the quadratic curvature approximation with composition-dependent bending rigidity and spontaneous curvature. Ideal lipid mixing is assumed within each lipid monolayer. We found that the most stable monolayer-coated DNA units are formed when the charged/neutral lipid composition corresponds (nearly) to charge neutralization; the optimal monolayer radius corresponds to close DNA-monolayer contact. These conclusions are also valid for the honeycomb complex, as the chain frustration energy is found to be negligible. Typically, the stabilization energies for these structures are on the order of 1 k(B)T/A of DNA length, reflecting mainly the balance between the electrostatic and bending energies. The spaghetti complexes are less stable due to the additional bending energy of the external monolayer. A thermodynamic analysis is presented for calculating the equilibrium lipid compositions when the complexes coexist with excess bilayer.     

Role of DNA conformation & energetic insights in Msx-1-DNA recognition as revealed by molecular dynamics studies on specific and nonspecific complexes

In most of homeodomain–DNA complexes, glutamine or lysine is present at 50th position and interacts with 5th and 6th nucleotide of core recognition region. Molecular dynamics simulations of Msx-1–DNA complex (Q50-TG) and its variant complexes, that is specific (Q50K-CC), nonspecific (Q50-CC) having mutation in DNA and (Q50K-TG) in protein, have been carried out. Analysis of protein–DNA interactions and structure of DNA in specific and nonspecific complexes show that amino acid residues use sequence-dependent shape of DNA to interact. The binding free energies of all four complexes were analysed to define role of amino acid residue at 50th position in terms of binding strength considering the variation in DNA on stability of protein–DNA complexes. The order of stability of protein–DNA complexes shows that specific complexes are more stable than nonspecific ones. Decomposition analysis shows that N-terminal amino acid residues have been found to contribute maximally in binding free energy of protein–DNA complexes. Among specific protein–DNA complexes, K50 contributes more as compared to Q50 towards binding free energy in respective complexes. The sequence dependence of local conformation of DNA enables Q50/Q50K to make hydrogen bond with nucleotide(s) of DNA. The changes in amino acid sequence of protein are accommodated and stabilized around TAAT core region of DNA having variation in nucleotides.     

Energetics of galactose- and glucose-aromatic amino acid interactions: implications for binding in galactose-specific proteins

An aromatic amino acid is present in the binding site of a number of sugar binding proteins. The interaction of the saccharide with the aromatic residue is determined by their relative position as well as orientation. The position-orientation of the saccharide relative to the aromatic residue was found to vary in different sugar-binding proteins. In the present study, interaction energies of the complexes of galactose (Gal) and of glucose (Glc) with aromatic residue analogs have been calculated by ab initio density functional (U-B3LYP/ 6-31G**) theory. The position-orientations of the saccharide with respect to the aromatic residue observed in various Gal-, Glc-, and mannose-protein complexes were chosen for the interaction energy calculations. The results of these calculations show that galactose can interact with the aromatic residue with similar interaction energies in a number of position-orientations. The interaction energy of Gal-aromatic residue analog complex in position-orientations observed for the bound saccharide in Glc/Man-protein complexes is comparable to the Glc-aromatic residue analog complex in the same position-orientation. In contrast, there is a large variation in interaction energies of complexes of Glc- and of Gal- with the aromatic residue analog in position-orientations observed in Gal-protein complexes. Furthermore, the conformation wherein the O6 atom is away from the aromatic residue is preferred for the exocyclic -CH2OH group in Gal-aromatic residue analog complexes. The implications of these results for saccharide binding in Gal-specific proteins and the possible role of the aromatic amino acid to ensure proper positioning and orientation of galactose in the binding site have been discussed.     

A structural and free energy analysis of Ag complexes to five small peptides

The complexes of Ag⁺ with the peptides MetGly, ProGly, GlyPro, GlyHis and GlyProAla were investigated using hybrid density functional theory at the B3LYP/DZVP level. The silver ion binding free energies at 298 K to each of these peptides was calculated to be 60.8, 52.0, 54.3, 71.2 and 63.3 kcal mol⁻¹, respectively. Structural information and relative free energies are presented for several isomers for each of the five complexes. Each of the global minima found for the five complexes is a charge-solvated ion. An important finding is that the Ag⁺-ProGly is the only complex where a salt bridge structure is energetically favored occurring at 4.0 kcal mol⁻¹ higher in free energy than the global minimum. The Ag⁺ ion in this salt bridge structure is attached to the carboxylate anion of zwitterionic ProGly in which the terminal amino nitrogen is protonated. For all the other complexes studied, the salt bridge structure occurs at much higher energies. All the dipeptide complexes with Ag⁺, but one, exhibit a di- or tri-coordinate metal where the sites of attachment are amino and carbonyl groups. However, the highest coordination numbers are not always the global minima due to steric costs. The global minimum of the Ag⁺-GlyProAla complex is the only structure found in this study where the metal is tetra-coordinated, binding to the terminal amino nitrogen and all three carbonyl oxygen atoms. Silver binding to sulphur and imidazole nitrogen atoms of MetGly and GlyHis, respectively, are present in the three most energetically favored species in each of these cases.     

Binding free energy calculations of galectin-3-ligand interactions

Galectins show remarkable binding specificity towards beta-galactosides. A recently developed method for calculating binding free energies between a protein and its substrates has been used to evaluate the binding specificity of galectin-3. Five disaccharides and a tetrasaccharide were used as the substrates. The calculated binding free energies agree quite well with the experimental data and the ranking of binding affinities is well reproduced. For all the six protein-ligand complexes it was observed that electrostatic interactions oppose binding whereas the non-polar contributions drive complex formation. The observed binding specificity of galectin-3 for galactosides rather than glucosides is discussed in light of our results.     

Effective regulation of iron acquisition in graminaceous plants. The role of mugineic acids as phytosiderophores

Discovery of mugineic acids as phytosiderophores has shown that some graminaceous monocotyledonous plants have a different iron acquisition strategy (strategy II) from dicotyledonous and nongraminaceous monocotyledonous plants (strategy I). The process of iron acquisition by strategy II plants can be divided into four main steps: biosynthesis, secretion, solubilization, and uptake, all of which are effectively regulated by different systems. The biosynthesis of mugineic acids is controlled by an on-off system which is operated under the control of iron demand in the plant. All mugineic acids share the same biosynthetic pathway from L-methionine to 2'-deoxymugineic acid, but the subsequent steps differ among plant species and even cultivars. The biosynthesis of mugineic acids is associated with the methionine recycling pathway. The secretion of mugineic acids shows a distinct diumal rhythm. Mugineic acids solubilize sparingly soluble inorganic iron by chelation and possess a high chelation affinity for iron, but not for other polyvalent ions such as Ca²⁺, Mg²⁺ and Al³⁺. The iron uptake process is regulated by a specific uptake system that transports the mugineic acid-Fe(III) complex as an intact molecule. This system specifically recognizes the mugineic acid-Fe(III) complexes, but not other mugineic acid-metal or synthetic chelator-Fe(III) complexes, suggesting that binding sites with strict recognition for stereostructure of the complex are located on the plasma membrane. All these regulatory systems are considered to represent an efficient strategy to acquire adequate amounts of iron and to avoid factors unfavorable for iron acquisition such as high pH, high concentrations of bicarbonate, Ca^2- and Mg²⁺, microbial degradation, and uptake of other metals that are common in calcareous soils.     

“Russian doll” complexes of [n]cycloparaphenylenes: a theoretical study

The formation of "Russian doll" complexes consisting of [n]cycloparaphenylenes was predicted using quantum chemistry tools. The electronic structures of multiple inclusion complexes containing up to four macrocycles were explored at the M06-2X/6-31G* level of theory. The binding energy between the macrocycles increases from the center to the periphery of the complex and can be >60?kcal?mol(-1) for macrocycles containing 14 and 19 repeating units. It has been demonstrated that additional electrostatic interactions originating from the asymmetric electron density distribution observed when comparing the concave and convex macrocycle sides are responsible for the high binding energies in these Russian doll complexes. Oxidation or reduction of the Russian doll complexes creates polarons that are delocalized across the complexes. In the case of polaron cations, most of the polarons are localized at the macrocycle with the smallest ionization potential; for polaron anions, the negative charge is localized across the outer rings of the complex. Because anion polarons are more delocalized than cation polarons, the relaxation energies of the polaron anions were found to be smaller than those of the polaron cations.     

A study of the coordination shell of aluminum(III) and magnesium(II) in model protein environments: thermodynamics of the complex formation and metal exchange reactions

Al(III) toxicity in living organisms is based on competition with other metal cations. Mg(II) is one of the most affected cations, since the size similarity dominates over the charge identity. The slow ligand exchange rates for Al(III) render the ion useless as a metal ion at the active sites of enzymes and provide a mechanism by which Al(III) inhibits Mg(II) dependent biochemical processes. Al(III) cation interactions with relevant bioligands have been studied in a protein-model environment in gas and aqueous phases using density functional theory methods. The protein model consists of the metal cation bound to two chosen bioligands (functional groups of the amino acid side chains, one of them being always an acetate) and water molecules interacting with the cation to complete its first coordination shell. Analogous Mg(II) complexes are calculated and compared with the Al(III) ones. Formation energies of the complexes are calculated in both phases and magnesium/aluminum exchange energies evaluated. The effect of different dielectric media is also analyzed. The presence of an acetate ligand in the binding site is found to promote both, complex formation and metal exchange reactions. In addition, buried binding sites (with low dielectric constant) of the protein favor metal exchange, whereas fully solvated environments of high dielectric constant require the presence of two anionic ligands for metal exchange to occur.     

Calculating proton uptake/release and binding free energy taking into account ionization and conformation changes induced by protein-inhibitor association: application to plasmepsin, cathepsin D and endothiapepsin-pepstatin complexes

The protein-inhibitor binding energies of enzymes are often pH dependent, and binding induces either proton uptake or proton release. The proton uptake/release and the binding energy for three complexes with available experimental data were numerically studied: pepstatin-cathepsin D, pepstatin-plasmepsin II and pepstatin-endothiapepsin. Very good agreement with the experimental data was achieved when conformational changes were taken into account. The role of the desolvation energy and the conformational changes was revealed by modeling the complex, the separated molecules in the complex conformation and the free molecules. It was shown that the conformational changes induced by the complex formation are as important for the proton transfer as the loss of solvation energy caused by the burial of interface residues. The residues responsible for the proton transfer were identified and their contribution to the proton uptake/release calculated. These residues were found to be scattered along the whole protein rather than being localized only at the active site. In the case of cathepsin D, these residues were found to be highly conserved among the cathepsin D sequences of other species. It was shown that conformation and ionization changes induced by the complex formation are critical for the correct calculation of the binding energy. Taking into account the electrostatics and the van der Waals (vdW) energies within the Boltzmann distribution of energies and allowing ionization and conformation changes to occur makes the calculated binding energy more realistic and closer to the experimental value. The interplay between electrostatic and vdW forces makes the pH dependence of the binding energy smoother, because the vdW force acts in reaction to the changes of the electrostatic energy. It was found that a small fraction of the ionizable groups remain uncharged in both the free and complexed molecules. The sequence and structural position of these groups aligns well within the three proteases, suggesting that these may have specific role.     

Potential of mean force for protein-protein interaction studies.

Calculating protein-protein interaction energies is crucial for understanding protein-protein associations. On the basis of the methodology of mean-field potential, we have developed an empirical approach to estimate binding free energy for protein-protein interactions. This knowledge-based approach has been used to derive distance-dependent free energies of protein complexes from a nonredundant training set in the Protein Data Bank (PDB), with a careful treatment of homology. We calculate atom pair potentials for 16 pair interactions, which can reflect the importance of hydrophobic interactions and specific hydrogen-bonding interactions. The derived potentials for hydrogen-bonding interactions show a valley of favorable interactions at a distance of approximately 3 A, corresponding to that of an established hydrogen bond. For the test set of 28 protein complexes, the calculated energies have a correlation coefficient of 0.75 compared with experimental binding free energies. The performance of the method in ranking the binding energies of different protein-protein complexes shows that the energy estimation can be applied to value binding free energies for protein-protein associations.     

Theoretical study of the interaction between 5-methylcytosine and acrylamide

The hydrogen-bonded complexes between 5-methylcytosine and acrylamide have been investigated using the density function theory (DFT) method. Five stable complexes have been found with no imaginary frequencies. Complex C3 is the most stable one with interaction energies of -69.01?kJ?mol(-1) corrected for basis set superposition error (BSSE). The charge change in the process of these complexes formation has also been examined. The atoms in molecules (AIM) theory and natural bond orbital (NBO) method have been performed to investigate the hydrogen bonds involved in all the complexes. The electron density and its corresponding Laplacian at the bond and ring critical points have been analyzed. In C3 complex, there is the largest stabilization energy (18.17?kJ?mol(-1)) between N11-H12 antibonding orbital and lone electron pair of O17. It can be seen that the hydrogen bonds play a crucial role in the stability of all the complexes between 5-methylcytosine and acrylamide. The theoretical results could provide helpful information for other researchers in further work.     

Computational analysis of binding of P1 variants to trypsin

The binding of P1 variants of bovine pancreatic trypsin inhibitor (BPTI) to trypsin has been investigated by means of molecular dynamics simulations. The specific interaction formed between the amino acid at the primary binding (P1) position of the binding loop of BPTI and the specificity pocket of trypsin was estimated by use of the linear interaction energy (LIE) method. Calculations for 13 of the naturally occurring amino acids at the P1 position were carried out, and the results obtained were found to correlate well with the experimental binding free energies. The LIE calculations rank the majority of the 13 variants correctly according to the experimental association energies and the mean error between calculated and experimental binding free energies is only 0.38 kcal/mole, excluding the Glu and Asp variants, which are associated with some uncertainties regarding protonation and the possible presence of counter-ions. The three-dimensional structures of the complex with three of the P1 variants (Asn, Tyr, and Ser) included in this study have not at present been solved by any experimental techniques and, therefore, were modeled on the basis of experimental data from P1 variants of similar size. Average structures were calculated from the MD simulations, from which specific interactions explaining the broad variation in association energies were identified. The present study also shows that explicit treatment of the complex water-mediated hydrogen bonding network at the protein-protein interface is of crucial importance for obtaining reliable binding free energies. The successful reproduction of relative binding energies shows that this type of methodology can be very useful as an aid in rational design and redesign of biologically active macromolecules.     

Prediction of protein complexes using empirical free energy functions.

A long sought goal in the physical chemistry of macromolecular structure, and one directly relevant to understanding the molecular basis of biological recognition, is predicting the geometry of bimolecular complexes from the geometries of their free monomers. Even when the monomers remain relatively unchanged by complex formation, prediction has been difficult because the free energies of alternative conformations of the complex have been difficult to evaluate quickly and accurately. This has forced the use of incomplete target functions, which typically do no better than to provide tens of possible complexes with no way of choosing between them. Here we present a general framework for empirical free energy evaluation and report calculations, based on a relatively complete and easily executable free energy function, that indicate that the structures of complexes can be predicted accurately from the structures of monomers, including close sequence homologues. The calculations also suggest that the binding free energies themselves may be predicted with reasonable accuracy. The method is compared to an alternative formulation that has also been applied recently to the same data set. Both approaches promise to open new opportunities in macromolecular design and specificity modification.