No abstract

No abstract

The G protein-coupled receptor (GPCR) family of proteins comprises signaling proteins that mediate cellular responses to various hormones and neurotransmitters, and serves as a prime target for drug discovery. Towards our goal of discovering secondary metabolites from natural sources that can function as neuronal drugs, we evaluated the modulatory effect of eckol on various GPCRs via cell-based functional assays. In addition, we conducted in silico predictions to obtain molecular insights into the functional effects of eckol. Functional assays revealed that eckol had a concentration-dependent agonist effect on dopamine D3 and D4 receptors. The half maximal effective concentration (EC50) of eckol for the dopamine D3 and D4 receptors was 48.62 ± 3.21 and 42.55 ± 2.54 µM, respectively, while the EC50 values of dopamine as a reference agonist for these two receptors were 2.9 and 3.3 nM, respectively. In silico studies revealed that a low binding energy in addition to hydrophilic, hydrophobic, π–alkyl, and π–π T-shaped interactions are potential mechanisms by which eckol binds to the dopamine receptors to exert its agonist effects. Molecular dynamics (MD) simulation revealed that Phe346 of the dopamine receptors is important for binding of eckol, similar to eticlopride and dopamine. Our results collectively suggest that eckol is a potential D3/D4 agonist for the management of neurodegenerative diseases, such as Parkinson’s disease.

We report a fully general technique addressing a long standing challenge of calculating conformational free energy differences between various states of a polymer chain from simulations using explicit solvent force fields. The main feature of our method is a special mapping variable, a path coordinate, which continuously connects two conformations. The path variable has been designed to preserve locality in the phase space near the path endpoints. We avoid the problem of sampling the unfolded states by creating an artificial confinement “tube” in the phase space that prevents the molecule from unfolding without affecting the calculation of the desired free energy difference. We applied our technique to compute the free energy difference between two native-like conformations of the small protein Trp-cage using the CHARMM force field with explicit solvent. We verified this result by comparing it with an independent, significantly more expensive calculation. Overall, the present study suggests that the new method of computing free energy differences between polymer chain conformations is accurate and highly computationally efficient.

Many biologically interesting functions such as allosteric switching or protein-ligand binding are determined by the kinetics and mechanisms of transitions between various conformational substates of the native basin of globular proteins. To advance our understanding of these processes, we constructed a two-dimensional free energy surface (FES) of the native basin of a small globular protein, Trp-cage. The corresponding order parameters were defined using two native substructures of Trp-cage. These calculations were based on extensive explicit water all-atom molecular dynamics simulations. Using the obtained two-dimensional FES, we studied the transition kinetics between two Trp-cage conformations, finding that switching process shows a borderline behavior between diffusive and weakly-activated dynamics. The transition is well-characterized kinetically as a biexponential process. We also introduced a new one-dimensional reaction coordinate for the conformational transition, finding reasonable qualitative agreement with the two-dimensional kinetics results. We investigated the distribution of all the 38 native nuclear magnetic resonance structures on the obtained FES, analyzing interactions that stabilize specific low-energy conformations. Finally, we constructed a FES for the same system but with simple dielectric model of water instead of explicit water, finding that the results were surprisingly similar in a small region centered on the native conformations. The dissimilarities between the explicit and implicit model on the larger-scale point to the important role of water in mediating interactions between amino acid residues.