Small metal ions play critical roles in numerous biological processes. Of particular interest is how metalloenzymes are allosterically regulated by the binding of specific ions. Understanding how ion binding affects these biological processes requires atomic models that accurately treat the microscopic interactions with the protein ligands. Theoretical approaches at different levels of sophistication can contribute to a deeper understanding of these systems, although computational models must strike a balance between accuracy and efficiency in order to enable long molecular dynamics simulations. In this study, we present a systematic effort to optimize the parameters of a polarizable force field based on classical Drude oscillators to accurately represent the interactions between ions (K⁺, Na⁺, Ca²⁺, and Cl^–) and coordinating amino-acid residues for a set of 30 biologically important proteins. By combining ab initio calculations and experimental thermodynamic data, we derive a polarizable force field that is consistent with a wide range of properties, including the geometries and interaction energies of gas-phase ion/protein-like model compound clusters, and the experimental solvation free-energies of the cations in liquids. The resulting models display significant improvements relative to the fixed-atomic-charge additive CHARMM C36 force field, particularly in their ability to reproduce the many-body electrostatic nonadditivity effects estimated from ab initio calculations. The analysis clarifies the fundamental limitations of the pairwise additivity assumption inherent in classical fixed-charge force fields, and shows its dramatic failures in the case of Ca²⁺ binding sites. These optimized polarizable models, amenable to computationally efficient large-scale MD simulations, set a firm foundation and offer a powerful avenue to study the roles of the ions in soluble and membrane transport proteins.

Despite decades of investigations, the principal mechanisms responsible for the high affinity and specificity of proteins for key physiological cations K⁺, Na⁺, and Ca²⁺ remain a hotly debated topic. At the core of the debate is an apparent need (or lack thereof) for an accurate description of the electrostatic response of the charge distribution in a protein to the binding of an ion. These effects range from partial electronic polarization of the directly ligating atoms to long-range effects related to partial charge transfer and electronic delocalization effects. While accurate modeling of cation recognition by metalloproteins warrants the use of quantum-mechanics (QM) calculations, the most popular approximations used in major biomolecular simulation packages rely on the implicit modeling of electronic polarization effects. That is, high-level QM computations for ion binding to proteins are desirable, but they are often unfeasible, because of the large size of the reactive-site models and the need to sample conformational space exhaustively at finite temperature. Several solutions to this challenge have been proposed in the field, ranging from the recently developed Drude polarizable force-field for simulations of metalloproteins to approximate tight-binding density functional theory (DFTB). To delineate the usefulness of different approximations, we examined the accuracy of three recent and commonly used theoretical models and numerical algorithms, namely, CHARMM C36, the latest developed Drude polarizable force fields, and DFTB3 with the latest 3OB parameters. We performed MD simulations for 30 cation-selective proteins with high-resolution X-ray structures to create ensembles of structures for analysis with different levels of theory, e.g., additive and polarizable force fields, DFTB3, and DFT. The results from DFT computations were used to benchmark CHARMM C36, Drude, and DFTB3 performance. The explicit modeling of quantum effects unveils the key electrostatic properties of the protein sites and the importance of specific ion-protein interactions. One of the most interesting findings is that secondary coordination shells of proteins are noticeably perturbed in a cation-dependent manner, showing significant delocalization and long-range effects of charge transfer and polarization upon binding Ca²⁺.

Recently, microbial prodigiosin (PG) has received much attention due to its numerous beneficial applications. The aim of this study was to establish the bioprocessing of marine chitinous wastes (MCWs) for the cost-effective preparation of PG. Of the MCWs, demineralized shrimp shell powders (de-SSP) were found to be a potential source of carbon/nitrogen (C/N) for PG production by bacterial fermentation using Serratia marcescens strains. Further, PG scale-up production was investigated in a 15 L bioreactor system, and the highest yield (6200 mg/L) was achieved during fermentation using 5 L of a novel-designed culture broth that included 1.60% C/N sources (a de-SSP/casein ratio of 7/3), 0.02% K₂SO_4, and 0.05% K₂HPO₄, with an initial pH of 6–7. Fermentation was conducted in the dark at 27.5 °C for 8.0 h. This study was the first to report on the utilization of shrimp wastes for cost-effective, large-scale (5 L/pilot) PG production with high productivity (6200 mg/L) in a short cultivation time. The combination of 0.02% K₂SO₄ and 0.05% K₂HPO₄ was also found to be a novel salt composition that significantly enhanced PG yield. The red compound was purified and confirmed as PG after analyzing its HPLC profile, mass, and UV/vis spectra. The purified PG was then tested for its bioactivities and showed effective anticancer activities, moderated antioxidant activities, and novel anti-NO effects.

Recent development of the Drude polarizable (Drude) force field (FF), based on the extension of an induced dipole model, has reached a milestone in the past few years providing a complete set of polarizable parameters for proteins, water, ions, and many lipid types. This FF enables stable simulations up to microseconds, surpassing the capability of other polarizable FFs. The quality of the Drude FF, however, has remained largely untested for modeling the secondary structures of small peptides in explicit solvents compared with classical non‐polarizable FFs. It is critical to benchmark the complex and mutually dependent dynamics of hydrogen‐bond (H‐bond) networks formed by water–water, protein–water, and protein–protein interactions that are expected to have a major impact on the stability of protein structures and their conformational space. Here, a direct comparison is presented between the current Drude FF and the CHARMM‐36 non‐polarizable classical FF for 1) the solvation free energy of mimetics for all amino acid side‐chain equivalents, 2) limited conformational space, 3) protein–water and protein–protein interactions, and 4) the comparative lifetimes of H‐bonds. The impact of counterions on the stabilization of secondary structure in model peptides is additionally discussed and compared between these FFs.