Search

Aggregation of amyloid-β (Aβ) peptides from soluble monomers to insoluble amyloid fibrils has been hypothesized to be one of the crucial steps in the progression of Alzheimer’s disease (AD). Due to the disordered nature of Aβ peptides, identifying aggregation inhibitors as potential drug candidates against AD has been a great challenge. In this communication, we report an atomistic simulation study of the inhibition mechanism of two small molecules, homotaurine and scyllo-inositol, which are AD drug candidates currently under investigation. Using a replica exchange molecular dynamics method to extensively explore the conformational space of Aβ42 monomer with and without the presence of the small molecule agents, we found that both drug candidates reduce the β-strand propensity of the C-terminus region (I31-A42) and promote a conformational change of Aβ42 monomer toward a more collapsed phase through a non-specific binding mechanism. These findings provide atomistic-level insights into understanding of the inhibitory mechanisms of the two potential small-molecule drug candidates for AD treatment in the future.

Molecular dynamics simulations using physics-based atomistic force fields have been increasingly used to characterize the heterogeneous structural ensembles of intrinsically disordered proteins (IDPs). To evaluate the accuracy of the latest atomistic explicit-solvent force fields in modeling larger IDPs with nontrivial structural features, we focus on the 61-residue N-terminal transactivation domain (TAD) of tumor suppressor p53, an important protein in cancer biology that has been extensively studied, and abundant experimental data is available for evaluation of simulated ensembles. We performed extensive replica exchange with solute tempering simulations, in excess of 1.0 μs/replica, to generate disordered structural ensembles of p53-TAD using six latest explicit solvent protein force fields. Multiple local and long-range structural properties, including chain dimension, residual secondary structures, and transient long-range contacts, were analyzed and compared against available experimental data. The results show that IDPs such as p53-TAD remain highly challenging for atomistic simulations due to conformational complexity and difficulty in achieving adequate convergence. Structural ensembles of p53-TAD generated using various force fields differ significantly from each other. The a99SB-disp force field demonstrates the best agreement with experimental data at all levels and proves to be suitable for simulating unbound p53-TAD and how its conformational properties may be modulated by phosphorylation and other cellular signals or cancer-associated mutations. Feasibility of such detailed structural characterization is a key step toward establishing the sequence-disordered ensemble-function-disease relationship of p53 and other biologically important IDPs.

The generalized Born (GB) theory is a prime choice for implicit treatment of solvent that provides a favorable balance between efficiency and accuracy for reliable simulation of protein conformational equilibria. In GB, the dielectric boundary is a key physical property that needs to be properly described. While it is widely accepted that the molecular surface (MS) should provide the most physical description, most existing GB models are based on van der Waals (vdW)-like surfaces for computational simplicity and efficiency. A simple and effective approximation to molecular volume is explored here using atom-centered dielectric functions within the context of a generalized Born model with simple switching (GBSW). The new model, termed GBSW/MS2, is as efficient as the original vdW-like-surface-based GBSW model, but is able to reproduce the Born radii calculated from the "exact" Poisson-Boltzmann theory with a correlation of 0.95. More importantly, examination of the potentials of mean force of hydrogen-bonding and charge-charge interactions demonstrates that GBSW/MS2 correctly captures the first desolvation peaks, a key signature of true MS. Physical parameters including atomic input radii and peptide backbone torsion were subsequently optimized on the basis of solvation free energies of model compounds, potentials of mean force of their interactions, and conformational equilibria of a set of helical and β-hairpin model peptides. The resulting GBSW/MS2 protein force field reasonably recapitulates the structures and stabilities of these model peptides. Several remaining limitations and possible future developments are also discussed.

Incorporating receptor flexibility into molecular docking should improve results for flexible proteins. However, the incorporation of explicit all-atom flexibility with molecular dynamics for the entire protein chain may also introduce significant error and “noise” that could decrease docking accuracy and deteriorate the ability of a scoring function to rank native-like poses. We address this apparent paradox by comparing the success of several flexible receptor models in cross-docking and multiple receptor ensemble docking for p38α mitogen-activated protein (MAP) kinase. Explicit all-atom receptor flexibility has been incorporated into a CHARMM-based molecular docking method (CDOCKER) using both molecular dynamics (MD) and torsion angle molecular dynamics (TAMD) for the refinement of predicted protein-ligand binding geometries. These flexible receptor models have been evaluated, and the accuracy and efficiency of TAMD sampling is directly compared to MD sampling. Several flexible receptor models are compared, encompassing flexible side chains, flexible loops, multiple flexible backbone segments, and treatment of the entire chain as flexible. We find that although including side chain and some backbone flexibility is required for improved docking accuracy as expected, docking accuracy also diminishes as additional and unnecessary receptor flexibility is included into the conformational search space. Ensemble docking results demonstrate that including protein flexibility leads to to improved agreement with binding data for 227 active compounds. This comparison also demonstrates that a flexible receptor model enriches high affinity compound identification without significantly increasing the number of false positives from low affinity compounds.