Molecular dynamics of ß-hairpin models of epigenetic recognition motifs.

The conformations and stabilities of the ß-hairpin model peptides of Waters (Riemen, A. J.; Waters, M. L. Biochemistry 2009, 48, 1525; Hughes, R. M.; Benshoff, M. L.; Waters, M. L. Chemistry 2007, 13, 5753) have been experimentally characterized as a function of lysine ε-methylation. These models were developed to explore molecular recognition of known epigenetic recognition motifs. This system offered an opportunity to computationally examine the role of cation-π interactions, desolvation of the ε-methylated ammonium groups, and aromatic/aromatic interactions on the observed differences in NMR spectra. AMOEBA, a second-generation force field (Ponder, J. W.; Wu, C.; Ren, P.; Pande, V. S.; Chodera, J. D.; Schnieders, M. J.; Haque, I.; Mobley, D. L.; Lambrecht, D. S.; DiStasio, R. A., Jr.; Head-Gordon, M.; Clark, G. N.; Johnson, M. E.; Head-Gordon, T. J. Phys. Chem. B 2010, 114, 2549), was chosen as it includes both multipole electrostatics and polarizability thought to be essential to accurately characterize such interactions. Independent parametrization of ε-methylated amines was required from which aqueous solvation free energies were estimated and shown to agree with literature values. Molecular dynamics simulations (100 ns) using the derived parameters with model peptides, such as Ac-R-W-V-W-V-N-G-Orn-K(Me)(n)-I-L-Q-NH(2), where n = 0, 1, 2, or 3, were conducted in explicit solvent. Distances between the centers of the indole rings of the two-tryptophan residues, 2 and 4, and the ε-methylated ammonium group on Lys-9 as well as the distance between the N- and C-termini were monitored to estimate the strength and orientation of the cation-π and aromatic/aromatic interactions. In agreement with the experimental data, the stability of the ß-hairpin increased significantly with lysine ε-methylation. The ability of MD simulations to reproduce the observed NOEs for the four peptides was further estimated for the monopole-based force fields, AMBER, CHARMM, and OPLSAA. AMOEBA correctly predicted over 80% of the observed NOEs for all 4 peptides, while the three-monopole force fields were 40-50% predictive in only 2 cases and approximately 10% in the other 10 examples. Preliminary analysis suggests that the decreased cost of desolvation of the substituted ammonium group significantly compensated for the reduced cation-π interaction resulting from the increased separation due to steric bulk of the ε-methylated amines.

Drug-like leads for steric discrimination between substrate and inhibitors of human acetylcholinesterase.

Protection of the enzyme acetylcholinesterase (AChE) from the toxic effects of organophosphate insecticides and chemical warfare agents (OPs) may be provided by inhibitors that bind at the peripheral binding site (P-site) near the mouth of the active-site gorge. Compounds that bind to this site may selectively block access to the acylation site (A-site) catalytic serine for OPs, but not acetylcholine itself. To identify such compounds, we employed a virtual screening approach using AutoDock 4.2 and AutoDock Vina, confirmed using compounds experimentally known to bind specifically to either the A-site or P-site. Both programs demonstrated the ability to correctly predict the binding site. Virtual screening of the NCI Diversity Set II was conducted using the apo form of the enzyme, and with acetylcholine bound at the crystallographic locations in the A-site only and in both and A- and P-sites. The docking identified 32 compounds that were obtained for testing, and one was demonstrated to bind specifically to the P-site in an inhibitor competition assay.

Molecular modeling of the cytoskeleton.

Molecular modeling techniques have truly come of age in recent decades, and here we cover several of the most commonly used techniques, namely molecular dynamics, Brownian dynamics, and molecular docking. In each case, we explain the physical basis and limitations of the various techniques and then illustrate their application to various problems related to the cytoskeleton. This set of studies covers a relatively wide range of examples and is comprehensive enough to clearly see how these techniques could be applied to other systems. Finally, we cover several related methodologies that expand on these basic techniques to allow for more detailed and specific simulation and analysis.

Nucleotide effects on the structure and dynamics of actin.

Adenosine 5'-triphosphate or ATP is the primary energy source within the cell, releasing its energy via hydrolysis into adenosine 5'-diphosphate or ADP. Actin is an important ATPase involved in many aspects of cellular function, and the binding and hydrolysis of ATP regulates its polymerization into actin filaments as well as its interaction with a host of actin-associated proteins. Here we study the dynamics of monomeric actin in ATP, ADP-Pi, and ADP states via molecular dynamics simulations. As observed in some crystal structures we see that the DNase-I loop is an alpha-helix in the ADP state but forms an unstructured coil domain in the ADP-Pi and ATP states. We also find that this secondary structure change is reversible, and by mimicking nucleotide exchange we can observe the transition between the helical and coil states. Apart from the DNase-I loop, we also see several key structural differences in the nucleotide binding cleft as well as in the hydrophobic cleft between subdomains 1 and 3 where WH2-containing proteins have been shown to interact. These differences provide a structural basis for understanding the observed differences between the various nucleotide states of actin and provide some insight into how ATP regulates the interaction of actin with itself and other proteins.

Origin of co-conformational selectivity in a [3]rotaxane.

Co-conformational selectivity and structure-energy relationships in a [3]rotaxane are investigated with a recently developed multiple-sampling and statistical analysis procedure for modeling interlocked molecules and mechanical molecular devices. The results presented confirm the experimentally observed co-conformational selectivity. The theoretical calculations reveal that ring-ring interactions are very small and ring-shaft inter-component interactions decide the co-conformational preference. In particular, it is found that stronger ring binding at the central binding station on the shaft than at either of the two terminal binding stations gives rise to the observed co-conformational preference. Analysis of radius of gyration data shows that co-conformational isomerism is not strongly correlated to coiling of the shaft.