Search

In the preceding paper by Stork and Tavan, [J. Chem. Phys. 126, 165105 (2007)], the authors have reformulated an electrostatic theory which treats proteins surrounded by dielectric solvent continua and approximately solves the associated Poisson equation [B. Egwolf and P. Tavan, J. Chem. Phys. 118, 2039 (2003)]. The resulting solution comprises analytical expressions for the electrostatic reaction field (RF) and potential, which are generated within the protein by the polarization of the surrounding continuum. Here the field and potential are represented in terms of Gaussian RF dipole densities localized at the protein atoms. Quite like in a polarizable force field, also the RF dipole at a given protein atom is induced by the partial charges and RF dipoles at the other atoms. Based on the reformulated theory, the authors have suggested expressions for the RF forces, which obey Newton’s third law. Previous continuum approaches, which were also built on solutions of the Poisson equation, used to violate the reactio principle required by this law, and thus were inapplicable to molecular dynamics (MD) simulations. In this paper, the authors suggest a set of techniques by which one can surmount the few remaining hurdles still hampering the application of the theory to MD simulations of soluble proteins and peptides. These techniques comprise the treatment of the RF dipoles within an extended Lagrangian approach and the optimization of the atomic RF polarizabilities. Using the well-studied conformational dynamics of alanine dipeptide as the simplest example, the authors demonstrate the remarkable accuracy and efficiency of the resulting RF-MD approach. [ABSTRACT FROM AUTHOR]

The authors reformulate and revise an electrostatic theory treating proteins surrounded by dielectric solvent continua [B. Egwolf and P. Tavan, J. Chem. Phys. 118, 2039 (2003)] to make the resulting reaction field (RF) forces compatible with Newton’s third law. Such a compatibility is required for their use in molecular dynamics (MD) simulations, in which the proteins are modeled by all-atom molecular mechanics force fields. According to the original theory the RF forces, which are due to the electric field generated by the solvent polarization and act on the partial charges of a protein, i.e., the so-called qE forces, can be quite accurately computed from Gaussian RF dipoles localized at the protein atoms. Using a slightly different approximation scheme also the RF energies of given protein configurations are obtained. However, because the qE forces do not account for the dielectric boundary pressure exerted by the solvent continuum on the protein, they do not obey the principle that actio equals reactio as required by Newton’s third law. Therefore, their use in MD simulations is severely hampered. An analysis of the original theory has led the authors now to a reformulation removing the main difficulties. By considering the RF energy, which represents the dominant electrostatic contribution to the free energy of solvation for a given protein configuration, they show that its negative configurational gradient yields mean RF forces obeying the reactio principle. Because the evaluation of these mean forces is computationally much more demanding than that of the qE forces, they derive a suggestion how the qE forces can be modified to obey Newton’s third law. Various properties of the thus established theory, particularly issues of accuracy and of computational efficiency, are discussed. A sample application to a MD simulation of a peptide in solution is described in the following paper [M. Stork and P. Tavan, J. Chem. Phys., 126, 165106 (2007). [ABSTRACT FROM AUTHOR]

The molecular structures of amyloid fibers characterizing neurodegenerative diseases such as Huntington's or transmissible spongiform encephalopathies are unknown. Recently, x-ray diffraction patterns of poly-Gln fibers and electron microscopy images of two-dimensional crystals formed from building blocks of prion rods have suggested that the corresponding amyloid fibers are generated by the aggregation of parallel beta-helices. To explore this intriguing concept, we study the stability of small beta-helices in aqueous solution by molecular dynamics simulations. In particular, for the Huntington aggregation nucleus, which is thought to be formed of poly-Gln polymers, we show that three-coiled beta-helices are unstable at the suggested circular geometries and stable at a triangular shape with 18 residues per coil. Moreover, we demonstrate that individually unstable two-coiled triangular poly-Gln beta-helices become stabilized upon dimerization, suggesting that seeded aggregation of Huntington amyloids requires dimers of at least 36 Gln repeats (or monomers of approximately 54 Gln) for the formation of sufficiently stable aggregation nuclei. An analysis of our results and of sequences occurring in native beta-helices leads us to the proposal of a revised model for the PrP(Sc) aggregation nucleus.