Interaction of Charged Amino-Acid Side Chains with Ions: An Optimization Strategy for Classical Force Fields

Peptides Used for Heavy Metal Remediation: A Promising Approach

In recent years, heavy metal pollution has become increasingly prominent, severely damaging ecosystems and biodiversity, and posing a serious threat to human health. However, the results of current methods for heavy metal restoration are not satisfactory, so it is urgent to find a new and effective method. Peptides are the units that make up proteins, with small molecular weights and strong biological activities. They can effectively repair proteins by forming complexes, reducing heavy metal ions, activating the plant's antioxidant defense system, and promoting the growth and metabolism of microorganisms. Peptides show great potential for the remediation of heavy metal contamination due to their special structure and properties. This paper reviews the research progress in recent years on the use of peptides to remediate heavy metal pollution, describes the mechanisms and applications of remediation, and provides references for the remediation of heavy metal pollution.

Recalibrating the calcium trap in amino acid carboxyl groups via classical molecular dynamics simulations

In order to use classical molecular dynamics to complement experiments accurately, it is important to use robust descriptions of the system. The interactions between biomolecules, like aspartic and glutamic acid, and dissolved ions are often studied using standard biomolecular force-fields, where the interactions between biomolecules and cations are often not parameterized explicitly. In this study, we have employed metadynamics simulations to investigate different interactions of Ca with aspartic and glutamic acid and constructed the free energy profiles of Ca²⁺-carboxylate association. Starting from a generally accepted, AMBER-based force field, the association was substantially over and under-estimated, depending on the choice of water model (TIP3P and SPC/fw, respectively). To rectify this discrepancy, we have replaced the default calcium parameters. Additionally, we modified the σ_ij value in the hetero-atomic Lennard-Jones interaction by 0.5% to further improve the interaction between Ca and carboxylate, based on comparison with the experimentally determined association constant for Ca with the carboxylate group of L-aspartic acid. The corrected description retrieved the structural properties of the ion pair in agreement with the original biomolecule - Ca²⁺ interaction in AMBER, whilst also producing an association constant comparable to experimental observations. This refined force field was then used to investigate the interactions between amino acids, calcium and carbonate ions during biogenic and biomimetic calcium carbonate mineralisation.

How Accurate Is the Egg-Box Model in Describing the Binding of Calcium to Polygalacturonate? A Molecular Dynamics Simulation Study

We performed molecular dynamics (MD) simulations of octameric galacturonate, GalA₈, chains in the presence of Ca²⁺ in a ratio of R = [Ca²⁺]/[GalA] = 0.25 in order to determine to which extent the popular "egg-box model" (EBM) is able to describe the association between Ca²⁺ cations and polygalacturonate (polyGalA) chains. To this aim, we slightly revised the empirical parameters for the interaction between Ca²⁺ and the carboxylate oxygen atoms of GalA units so as to reproduce the experimental Ca²⁺-GalA association constant. We also defined an ad hoc order parameter, referred to as the egg-box score (EBS), that quantifies any deviation of the local coordination geometry of calcium cations with respect to an "ideal" EBM coordination geometry. The results reveal that the local coordination geometry of Ca²⁺ cations bound to polyGalA chains differs from that of the EBM. Moreover, polyGalA chains exhibit significant conformational disorder, and the cross-link angles formed between polyGalA chain axes are broadly distributed. Overall, the present study suggests that the EBM fails to describe accurately the association modes between calcium and polyGalA chains at a molar ratio R of 0.25.

Impacts of targeting different hydration free energy references on the development of ion potentials

Hydration free energy (HFE) as the most important solvation parameter is often targeted in ion model development, even though the reported values differ by dozens of kcal mol^-1 mainly due to the experimentally undetermined HFE of the proton ΔG°(H⁺). The choice of ΔG°(H⁺) obviously affects the hydration of single ions and the relative HFE between the ions with different (magnitude or sign) charges, and the impacts of targeted HFEs on the ion solvation and ion-ion interactions are largely unrevealed. Here we designed point charge models of K⁺, Mg²⁺, Al³⁺, and Cl^- ions targeting a variety of HFE references and then investigated the HFE influences on the simulations of dilute and concentrated ion solutions and of the salt ion pairs in gas, liquid, and solid phases. Targeting one more property of ion-water oxygen distances (IOD) leaves the ion-water binding distance invariant, while the binding strength increases with the decreasing (more negative) HFE of ions as a result of a decrease in ΔG°(H⁺) for the cation and an increase in ΔG°(H⁺) for the anion. The increase in ΔG°(H⁺) leads to strengthened cation-anion interactions and thus to close ion-ion contacts, low osmotic pressures, and small activity derivatives in concentrated ion solutions as well as too stable ion pairs of the salts in different phases. The ion diffusivity and water exchange rates around the ions are simply not HFE dependent but rather more complex. Targeting both the aqueous IOD and salt crystal properties of KCl was also attempted and the comparison between different models indicates the complexity and challenge in obtaining a balanced performance between different phases using classical force fields. Our results also support that a real ΔG°(H⁺) value of -259.8 kcal mol^-1 recommended by Hünenberger and Reif guides ion models to reproduce ion-water and ion-ion interactions reasonably at relatively low salt concentrations. Simulations of a metalloprotein show that a relatively more positive ΔG°(H⁺) for Mg²⁺ model is better for a reasonable description of the metal binding network.

Role of water-bridged interactions in metal ion coupled protein allostery

Allosteric communication between distant parts of proteins controls many cellular functions, in which metal ions are widely utilized as effectors to trigger the allosteric cascade. Due to the involvement of strong coordination interactions, the energy landscape dictating the metal ion binding is intrinsically rugged. How metal ions achieve fast binding by overcoming the landscape ruggedness and thereby efficiently mediate protein allostery is elusive. By performing molecular dynamics simulations for the Ca²⁺ binding mediated allostery of the calmodulin (CaM) domains, each containing two Ca²⁺ binding helix-loop-helix motifs (EF-hands), we revealed the key role of water-bridged interactions in Ca²⁺ binding and protein allostery. The bridging water molecules between Ca²⁺ and binding residue reduces the ruggedness of ligand exchange landscape by acting as a lubricant, facilitating the Ca²⁺ coupled protein allostery. Calcium-induced rotation of the helices in the EF-hands, with the hydrophobic core serving as the pivot, leads to exposure of hydrophobic sites for target binding. Intriguingly, despite being structurally similar, the response of the two symmetrically arranged EF-hands upon Ca²⁺ binding is asymmetric. Breakage of symmetry is needed for efficient allosteric communication between the EF-hands. The key roles that water molecules play in driving allosteric transitions are likely to be general in other metal ion mediated protein allostery.

Natural proteins often utilize allostery in executing a variety of functions. Metal ions are typical cofactors to trigger the allosteric cascade. In this work, using the Ca²⁺ sensor protein calmodulin as the model system, we revealed crucial roles of water-bridged interactions in the metal ion coupled protein allostery. The coordination of the Ca²⁺ to the binding site involves an intermediate in which the water molecule bridges the Ca²⁺ and the liganding residue. The bridging water reduces the free energy barrier height of ligand exchange, therefore facilitating the ligand exchange and allosteric coupling by acting as a lubricant. We also showed that the response of the two symmetrically arranged EF-hand motifs of CaM domains upon Ca²⁺ binding is asymmetric, which is directly attributed to the differing dehydration process of the Ca²⁺ ions and is needed for efficient allosteric communication.

Simulating the binding of key organic functional groups to aqueous calcium carbonate species

The interaction of organic molecules with mineral systems is relevant to a wide variety of scientific problems both in the environment and minerals processing. In this study, the coordination of small organics that contain the two most relevant functional groups for biomineralisation of calcium carbonate, namely carboxylate and ammonium, with the corresponding mineral ions are examined in aqueous solution. Specifically, two force fields have been examined based on rigid-ion or polarisable models, with the latter being within the AMOEBA formalism. Here the parameters for the rigid-ion model are determined to target the accurate reproduction of the hydration structure and solvation thermodynamics, while both force fields are designed to be compatible with the corresponding recently published models for aqueous calcium carbonate. The application of these force fields to ion pairing in aqueous solution is studied in order to quantitatively determine the extent of association.

Ion Pair Supramolecular Structure Identified by ATR-FTIR Spectroscopy and Simulations in Explicit Solvent*

The present work uses ATR-FTIR spectroscopy assisted by simulations in explicit solvent and frequency calculations to investigate the supramolecular structure of carboxylate alkali-metal ion pairs in aqueous solutions. ATR-FTIR spectra in the 0.25-4.0 M concentration range displayed cation-specific behaviors, which enabled the measurement of the appearance concentration thresholds of contact ion pairs between 1.9 and 2.6 M depending on the cation. Conformational explorations performed using a non-local optimization method associated to a polarizable force-field (AMOEBA), followed by high quantum chemistry level (RI-B97-D3/dhf-TZVPP) optimizations, mode-dependent scaled harmonic frequency calculations and electron density analyses, were used to identify the main supramolecular structures contributing to the experimental spectra. A thorough analysis enables us to reveal the mechanisms responsible for the spectroscopic sensitivity of the carboxylate group and the respective role played by the cation and the water molecules, highlighting the necessity of combining advanced experimental and theoretical techniques to provide a fair and accurate description of ion pairing.

Three Reasons Why Aspartic Acid and Glutamic Acid Sequences Have a Surprisingly Different Influence on Mineralization

Understanding the role of polymers rich in aspartic acid (Asp) and glutamic acid (Glu) is the key to gaining precise control over mineralization processes. Despite their chemical similarity, experiments revealed a surprisingly different influence of Asp and Glu sequences. We conducted molecular dynamics simulations of Asp and Glu peptides in the presence of calcium and chloride ions to elucidate the underlying phenomena. In line with experimental differences, in our simulations, we indeed find strong differences in the way the peptides interact with ions in solution. The investigated Asp pentapeptide tends to pull a lot of ions into its vicinity, and many structures with clusters of calcium and chloride ions on the surface of the peptide can be observed. Under the same conditions, comparatively fewer ions can be found in proximity of the investigated Glu pentapeptide, and the structures are characterized by single calcium ions bound to multiple carboxylate groups. Based on our simulation data, we identified three reasons contributing to these differences, leading to a new level of understanding additive-ion interactions.

Sulfation and Calcium Favor Compact Conformations of Chondroitin in Aqueous Solutions

The effects of sulfation and calcium cations (Ca²⁺) on the atomic-resolution conformational properties of chondroitin carbohydrate polymers in aqueous solutions are not well studied owing to experimental challenges. Here, we compare all-atom explicit-solvent molecular dynamics simulations results for pairs of O-type (nonsulfated) and A-type (GlcNAc 4-O-sulfated) chondroitin 20-mers in 140 mM NaCl with and without Ca²⁺ and find that both sulfation and Ca²⁺ favor more compact polymer conformations. We also show that subtle differences in force-field parametrization can have dramatic effects on Ca²⁺ binding to chondroitin carboxylate and sulfate functional groups and thereby determine Ca²⁺-mediated intra- and interstrand association. In addition to providing an atomic-resolution picture of the interaction of Ca²⁺ with sulfated and nonsulfated chondroitin polymers, the molecular dynamics data emphasize the importance of careful force-field parametrization to balance ion-water and ion-chondroitin interactions and suggest additional parametrization efforts to tune interactions involving sulfate.

Investigating the calcium binding characteristics of black bean protein hydrolysate

The black bean protein has been widely utilized to prepare hydrolysates with different bioactive properties. Herein, we hydrolyzed the black bean protein to prepare hydrolysate with calcium binding activity and characterized its behavior. Our results showed that ficin was superior in obtaining hydrolysate with calcium binding capacity in comparison with trypsin, alcalase and bromelain. In particular, the optimal capacity of ficin hydrolysate reached 77.54 ± 1.61 μg mg-1, where the optimal hydrolysis conditions of ficin were a temperature of 70 °C, a pH value of 6.2, an enzyme concentration of 1.61% and a time of 3 h. This might be due to high proportions of aspartic acid and glutamic acid (35.59%). Further spectral analysis evidenced the formation of hydrolysate-calcium complexes, demonstrating that the interaction between hydrolysate and calcium ions primarily occur on carboxyl oxygen atoms and amino nitrogen atoms. These findings provide a possible utilization of black bean hydrolysate to serve as a calcium supplement nutraceutical to enhance the absorption and bioavailability.

Evidence for an egg-box-like structure in iron(ii)–polygalacturonate hydrogels: a combined EXAFS and molecular dynamics simulation study

The coordination of Fe(ii) with polygalacturonic acid (polyGalA) in Fe(ii)–polyGalA hydrogels exhibits an octahedral geometry that follows the “egg-box model”.

Challenges and advances in atomistic simulations of potassium and sodium ion channel gating and permeation

Ion channels are implicated in many essential physiological events such as electrical signal propagation and cellular communication. The advent of K⁺ and Na⁺ ion channel structure determination has facilitated numerous investigations of molecular determinants of their behaviour. At the same time, rapid development of computer hardware and molecular simulation methodologies has made computational studies of large biological molecules in all-atom representation tractable. The concurrent evolution of experimental structural biology with biomolecular computer modelling has yielded mechanistic details of fundamental processes unavailable through experiments alone, such as ion conduction and ion channel gating. This review is a short survey of the atomistic computational investigations of K⁺ and Na⁺ ion channels, focusing on KcsA and several voltage-gated channels from the K_V and Na_V families, which have garnered many successes and engendered several long-standing controversies regarding the nature of their structure-function relationship. We review the latest advancements and challenges facing the field of molecular modelling and simulation regarding the structural and energetic determinants of ion channel function and their agreement with experimental observations.

Binding of Divalent Cations to Insulin: Capillary Electrophoresis and Molecular Simulations

In the present study, we characterize the binding of divalent cations to insulin in aqueous salt solutions by means of capillary electrophoresis and molecular dynamics simulations. The results show a strong pH dependence. At low pH, at which all the carboxylate groups are protonated and the protein has an overall positive charge, all the cations exhibit only weak and rather unspecific interactions with insulin. In contrast, at close to neutral pH, when all the carboxylate groups are deprotonated and negatively charged, the charge-neutralizing effect of magnesium, calcium, and zinc, in particular, on the electrophoretic mobility of insulin is significant. This is also reflected in the results of molecular dynamics simulations showing accumulation of cations at the protein surface, which becomes smaller in magnitude upon effective inclusion of electronic polarization via charge rescaling.

Complexation of carboxylate on smectite surfaces

We report a first principles molecular dynamics (FPMD) study of carboxylate complexation on clay surfaces. By taking acetate as a model carboxylate, we investigate its inner-sphere complexes adsorbed on clay edges (including (010) and (110) surfaces) and in interlayer space. Simulations show that acetate forms stable monodentate complexes on edge surfaces and a bidentate complex with Ca²⁺ in the interlayer region. The free energy calculations indicate that the complexation on edge surfaces is slightly more stable than in interlayer space. By integrating pK_as and desorption free energies of Al coordinated water calculated previously (X. Liu, X. Lu, E. J. Meijer, R. Wang and H. Zhou, Geochim. Cosmochim. Acta, 2012, 81, 56-68; X. Liu, J. Cheng, M. Sprik, X. Lu and R. Wang, Geochim. Cosmochim. Acta, 2014, 140, 410-417), the pH dependence of acetate complexation has been revealed. It shows that acetate forms inner-sphere complexes on (110) in a very limited mildly acidic pH range while it can complex on (010) in the whole common pH range. The results presented in this study form a physical basis for understanding the geochemical processes involving clay-organics interactions.

Relating side chain organization of PNIPAm with its conformation in aqueous methanol

Combining nuclear magnetic resonance (NMR), dynamic light scattering (DLS), and μs long all-atom simulations with two million particles, we establish a delicate correlation between increased side chain organization of PNIPAm and its collapse in aqueous methanol mixtures. We find that the preferential binding of methanol with PNIPAm side chains, bridging distal monomers along the polymer backbone, results in increased organization. Furthermore, methanol-PNIPAm preferential binding is dominated by hydrogen bonding. Our findings reveal that the collapse of PNIPAm is dominated by enthalpic interactions and that the standard poor solvent (entropic) effects play no major role.

DrugMiner: comparative analysis of machine learning algorithms for prediction of potential druggable proteins

Application of computational methods in drug discovery has received increased attention in recent years as a way to accelerate drug target prediction. Based on 443 sequence-derived protein features, we applied the most commonly used machine learning methods to predict whether a protein is druggable as well as to opt for superior algorithm in this task. In addition, feature selection procedures were used to provide the best performance of each classifier according to the optimum number of features. When run on all features, Neural Network was the best classifier, with 89.98% accuracy, based on a k-fold cross-validation test. Among all the algorithms applied, the optimum number of most-relevant features was 130, according to the Support Vector Machine-Feature Selection (SVM-FS) algorithm. This study resulted in the discovery of new drug target which potentially can be employed in cell signaling pathways, gene expression, and signal transduction. The DrugMiner web tool was developed based on the findings of this study to provide researchers with the ability to predict druggable proteins. DrugMiner is freely available at www.DrugMiner.org.

Why does high pressure destroy co-non-solvency of PNIPAm in aqueous methanol?

It is well known that poly(N-isopropylacrylamide) (PNIPAm) exhibits an interesting, yet puzzling, phenomenon of co-non-solvency. Co-non-solvency occurs when two competing good solvents for PNIPAm, such as water and alcohol, are mixed together. As a result, the same PNIPAm collapses within intermediate mixing ratios. This complex conformational transition is driven by preferential binding of methanol with PNIPAm. Interestingly, co-non-solvency can be destroyed when applying high hydrostatic pressures. In this work, using a large scale molecular dynamics simulation employing high pressures, we propose a microscopic picture behind the suppression of the co-non-solvency phenomenon. Based on thermodynamic and structural analysis, our results suggest that the preferential binding of methanol with PNIPAm gets partially lost at high pressures, making the background fluid reasonably homogeneous for the polymer. This is consistent with the hypothesis that the co-non-solvency phenomenon is driven by preferential binding and is not based on depletion effects.

Improving the description of interactions between Ca2+ and protein carboxylate groups, including γ-carboxyglutamic acid: revised CHARMM22* parameters

We show that the CHARMM22* force-field over-binds the interaction between aqueous carboxylates and Ca²⁺, and introduce a modification that can recover experimentally-determined binding free energies for these systems.

Molecular simulation of oligo-glutamates in a calcium-rich aqueous solution: insights into peptide-induced polymorph selection

Advanced simulation methods show how glutamate oligomers prestructure Ca ions and induce structural motifs in correspondence with calciumoxalate pseudopolymorphs.

Molecular dynamics simulations of peptides at the air-water interface: influencing factors on peptide-templated mineralization

Biomineralization is the intricate, biomedically highly relevant process by which living organisms deposit minerals on biological matrices to stiffen tissues and build skeletal structures and shells. Rapaport and coworkers ( J. Am. Chem. Soc. 2000 , 122 , 12523 ; Adv. Funct. Mater. 2008 , 18 , 2889 ; Acta Biomater. 2012 , 8 , 2466 ) have designed a class of self-assembling amphiphilic peptides that are capable of forming hydrogels and attracting ions from the environment, generating structures akin to the extracellular matrix and promoting bone regeneration. The air-water interface serves both in experiment and in simulations as a model hydrophobic surface to mimic the cell's organic-aqueous interface and to investigate the organization of the peptide matrix into ordered β-pleated monolayers and the subsequent onset of biomineral formation. To obtain insight into the underlying molecular mechanism, we have used molecular dynamics simulations to study the effect of peptide sequence on aggregate stability and ion-peptide interactions. We find-in excellent agreement with experimental observations-that the nature of the peptide termini (proline vs phenylalanine) affect the aggregate order, while the nature of the acidic side chains (aspartic vs glutamic acid) affect the aggregate's stability in the presence of ions. These simulations provide valuable microscopic insight into the way ions and peptide templates mutually affect each other during the early stages of biomineralization preceding nucleation.

Role of electrostatic interactions in binding of peptides and intrinsically disordered proteins to their folded targets. 1. NMR and MD characterization of the complex between the c-Crk N-SH3 domain and the peptide Sos

Intrinsically disordered proteins (IDPs) often rely on electrostatic interactions to bind their structured targets. To obtain insight into the mechanism of formation of the electrostatic encounter complex, we investigated the binding of the peptide Sos (PPPVPPRRRR), which serves as a minimal model for an IDP, to the c-Crk N-terminal SH3 domain. Initially, we measured ¹⁵N relaxation rates at two magnetic field strengths and determined the binding shifts for the complex of Sos with wild-type SH3. We have also recorded a 3 μs molecular dynamics (MD) trajectory of this complex using the Amber ff99SB*-ILDN force field. The comparison of the experimental and simulated data shows that MD simulation consistently overestimates the strength of salt bridge interactions at the binding interface. The series of simulations using other advanced force fields also failed to produce any satisfactory results. To address this issue, we have devised an empirical correction to the Amber ff99SB*-ILDN force field whereby the Lennard-Jones equilibrium distance for the nitrogen-oxygen pair across the Arg-to-Asp and Arg-to-Glu salt bridges has been increased by 3%. Implementing this correction resulted in a good agreement between the simulations and the experiment. Adjusting the strength of salt bridge interactions removed a certain amount of strain contained in the original MD model, thus improving the binding of the hydrophobic N-terminal portion of the peptide. The arginine-rich C-terminal portion of the peptide, freed from the effect of the overstabilized salt bridges, was found to interconvert more rapidly between its multiple conformational states. The modified MD protocol has also been successfully used to simulate the entire binding process. In doing so, the peptide was initially placed high above the protein surface. It then arrived at the correct bound pose within ∼2 Å of the crystallographic coordinates. This simulation allowed us to analyze the details of the dynamic binding intermediate, i.e., the electrostatic encounter complex. However, an experimental characterization of this transient, weakly populated state remains out of reach. To overcome this problem, we designed the double mutant of c-Crk N-SH3 in which mutations Y186L and W169F abrogate tight Sos binding and shift the equilibrium toward the intermediate state resembling the electrostatic encounter complex. The results of the combined NMR and MD study of this engineered system will be reported in the next part of this paper.