Biological processes occur on space and time scales that are often unreachable for fully atomistic simulations. Therefore, simplified or coarse grain (CG) models for the theoretical study of these systems are frequently used. In this context, the accurate description of solvation properties remains an important and challenging field. In the present work, we report a new CG model based on the transient tetrahedral structures observed in pure water. Our representation lumps approximately 11 WATer molecules into FOUR tetrahedrally interconnected beads, hence the name WAT FOUR (WT4). Each bead carries a partial charge allowing the model to explicitly consider long-range electrostatics, generating its own dielectric permittivity and obviating the shortcomings of a uniform dielectric constant. We obtained a good representation of the aqueous environment for most biologically relevant temperature conditions in the range from 278 to 328 K. The model is applied to solvate simple CG electrolytes developed in this work (Na + , K + , and Cl -) and a recently published simplified representation of nucleic acids. In both cases, we obtained a good resemblance of experimental data and atomistic simulations. In particular, the solvation structure around DNA, partial charge neutralization by counterions, preference for sodium over potassium, and ion mediated minor groove narrowing as reported from X-ray crystallography are well reproduced by the present scheme. The set of parameters presented here opens the possibility of reaching the multimicroseconds time scale, including explicit solvation, ionic specificity, and long-range electrostatics, keeping nearly atomistic resolution with significantly reduced computational cost.

The interplay between dopamine and α-synuclein (AS) plays a central role in Parkinson's disease (PD). PD results primarily from a severe and selective devastation of dopaminergic neurons in substantia nigra pars compacta. The neuropathological hallmark of the disease is the presence of intraneuronal proteinaceous inclusions known as Lewy bodies within the surviving neurons, enriched in filamentous AS. In vitro, dopamine inhibits AS fibril formation, but the molecular determinants of this inhibition remain obscure. Here we use molecular dynamic (MD) simulations to investigate the binding of dopamine and several of its derivatives onto conformers representative of an NMR ensemble of AS structures in aqueous solution. Within the limitations inherent to MD simulations of unstructured proteins, our calculations suggest that the ligands bind to the 125YEMPS129 region, consistent with experimental findings. The ligands are further stabilized by long-range electrostatic interactions with glutamate 83 (E83) in the NAC region. These results suggest that by forming these interactions with AS, dopamine may affect AS aggregation and fibrillization properties. To test this hypothesis, we investigated in vitro the effects of dopamine on the aggregation of mutants designed to alter or abolish these interactions. We found that point mutations in the 125YEMPS129 region do not affect AS aggregation, which is consistent with the fact that dopamine interacts non-specifically with this region. In contrast, and consistent with our modeling studies, the replacement of glutamate by alanine at position 83 (E83A) abolishes the ability of dopamine to inhibit AS fibrillization.

The action of dopamine on the aggregation of the unstructured alpha-synuclein (α-syn) protein may be linked to the pathogenesis of Parkinson's disease. Dopamine and its oxidation derivatives may inhibit α-syn aggregation by non-covalent binding. Exploiting this fact, we applied an integrated computational and experimental approach to find alternative ligands that might modulate the fibrillization of α-syn. Ligands structurally and electrostatically similar to dopamine were screened from an established library. Five analogs were selected for in vitro experimentation from the similarity ranked list of analogs. Molecular dynamics simulations showed they were, like dopamine, binding non-covalently to α-syn and, although much weaker than dopamine, they shared some of its binding properties. In vitro fibrillization assays were performed on these five dopamine analogs. Consistent with our predictions, analyses by atomic force and transmission electron microscopy revealed that all of the selected ligands affected the aggregation process, albeit to a varying and lesser extent than dopamine, used as the control ligand. The in silico/in vitro approach presented here emerges as a possible strategy for identifying ligands interfering with such a complex process as the fibrillization of an unstructured protein.