Structural Investigation of DHICA Eumelanin Using Density Functional Theory and Classical Molecular Dynamics
Sepideh Soltani1, Shahin Sowlati-Hashjin2, Conrard Giresse Tetsassi Feugmo3
1Department of Physics and Astronomy, The University of Western Ontario, 1151 Richmond Street, London, ON N6A 3K7, Canada.
Molecules (Basel, Switzerland)
|December 11, 2022
Summary
This study used molecular dynamics simulations to investigate the self-assembly of eumelanin building blocks. Charged eumelanin forms soluble complexes with K+ ions, while uncharged eumelanin aggregates in aqueous solutions.
Area of Science:
- Biomaterials Science
- Polymer Chemistry
- Computational Biophysics
Background:
- Eumelanin, a key pigment in various tissues, is a complex polymer primarily composed of 5,6-dihydroxyindole-2-carboxylic acid (DHICA) and 5,6-dihydroxyindole (DHI) units.
- DHICA is the predominant building block in natural eumelanin, forming protomolecules with indole-2-carboxylic acid-5,6-quinone (ICAQ) and pyrrole-2,3,5-tricarboxylic acid (PTCA).
Purpose of the Study:
- To investigate the self-assembly behavior of DHICA-eumelanin protomolecules in aqueous solutions using molecular dynamics (MD) simulations.
- To explore the influence of molecular charge and counterions (K+) on DHICA-eumelanin aggregation and solubility at various concentrations.
Main Methods:
- Multi-microsecond molecular dynamics (MD) simulations were employed to model DHICA-eumelanin self-assembly.
- Density functional theory (DFT) calculations were used for initial parameterization of the DHICA-eumelanin protomolecule.
- Simulations were conducted on uncharged DHICA-eumelanin, charged DHICA-eumelanin (physiological pH), and binary mixtures.
Main Results:
- Uncharged DHICA-eumelanin spontaneously aggregated in aqueous solution, trapping water molecules within the aggregates.
- Charged DHICA-eumelanin (net charge of -4) exhibited increased solubility in the presence of K+ counterions, which preferentially bound to carboxylate groups.
- Binary mixtures showed aggregation of uncharged protomolecules, with charged protomolecules associating with the surface of these aggregates via K+ ions.
Conclusions:
- The charge state of DHICA-eumelanin significantly impacts its self-assembly and solubility in aqueous environments.
- Potassium ions play a crucial role in solubilizing charged DHICA-eumelanin by interacting with its carboxylate moieties, preventing aggregation.
- These findings provide insights into the structural organization and solution behavior of eumelanin, relevant for its biological functions and potential applications.
Related Concept Videos
Crystal Field Theory - Octahedral Complexes
27.1K
Crystal Field Theory
To explain the observed behavior of transition metal complexes (such as colors), a model involving electrostatic interactions between the electrons from the ligands and the electrons in the unhybridized d orbitals of the central metal atom has been developed. This electrostatic model is crystal field theory (CFT). It helps to understand, interpret, and predict the colors, magnetic behavior, and some structures of coordination compounds of transition metals.
CFT focuses on...
To explain the observed behavior of transition metal complexes (such as colors), a model involving electrostatic interactions between the electrons from the ligands and the electrons in the unhybridized d orbitals of the central metal atom has been developed. This electrostatic model is crystal field theory (CFT). It helps to understand, interpret, and predict the colors, magnetic behavior, and some structures of coordination compounds of transition metals.
CFT focuses on...
27.1K
UV–Vis Spectroscopy: Woodward–Fieser Rules
25.0K
UV–Visible absorption spectra of conjugated dienes arise from the lowest energy π → π* transitions. The light-absorbing part of the molecule is called the chromophore, and the substituents directly attached to the chromophore are called auxochromes. A strong correlation exists between the absorption maxima, λmax, and the structure of a conjugated π system. The Woodward–Fieser rules predict the value of λmax for a given...
25.0K


