Using enveloping distribution sampling to compute the folding free enthalpy of a β-peptide with a very unstable folded conformation in solution: The advantage of focused sampling using EDS

Methods for Classical-Mechanical Molecular Simulation in Chemistry: Achievements, Limitations, Perspectives

More than a half century ago it became feasible to simulate, using classical-mechanical equations of motion, the dynamics of molecular systems on a computer. Since then classical-physical molecular simulation has become an integral part of chemical research. It is widely applied in a variety of branches of chemistry and has significantly contributed to the development of chemical knowledge. It offers understanding and interpretation of experimental results, semiquantitative predictions for measurable and nonmeasurable properties of substances, and allows the calculation of properties of molecular systems under conditions that are experimentally inaccessible. Yet, molecular simulation is built on a number of assumptions, approximations, and simplifications which limit its range of applicability and its accuracy. These concern the potential-energy function used, adequate sampling of the vast statistical-mechanical configurational space of a molecular system and the methods used to compute particular properties of chemical systems from statistical-mechanical ensembles. During the past half century various methodological ideas to improve the efficiency and accuracy of classical-physical molecular simulation have been proposed, investigated, evaluated, implemented in general simulation software or were abandoned. The latter because of fundamental flaws or, while being physically sound, computational inefficiency. Some of these methodological ideas are briefly reviewed and the most effective methods are highlighted. Limitations of classical-physical simulation are discussed and perspectives are sketched.

Effect of Oxidative Damage on the Stability and Dimerization of Superoxide Dismutase 1

During their life cycle, proteins are subject to different modifications involving reactive oxygen species. Such oxidative damage to proteins may lead to the formation of insoluble aggregates and cytotoxicity and is associated with age-related disorders including neurodegenerative diseases, cancer, and diabetes. Superoxide dismutase 1 (SOD1), a key antioxidant enzyme in human cells, is particularly susceptible to such modifications. Moreover, this homodimeric metalloenzyme has been directly linked to both familial and sporadic amyotrophic lateral sclerosis (ALS), a devastating, late-onset motor neuronal disease, with more than 150 ALS-related mutations in the SOD1 gene. Importantly, oxidatively damaged SOD1 aggregates have been observed in both familial and sporadic forms of the disease. However, the molecular mechanisms as well as potential implications of oxidative stress in SOD1-induced cytotoxicity remain elusive. In this study, we examine the effects of oxidative modification on SOD1 monomer and homodimer stability, the key molecular properties related to SOD1 aggregation. We use molecular dynamics simulations in combination with thermodynamic integration to study microscopic-level site-specific effects of oxidative "mutations" at the dimer interface, including lysine, arginine, proline and threonine carbonylation, and cysteine oxidation. Our results show that oxidative damage of even single residues at the interface may drastically destabilize the SOD1 homodimer, with several modifications exhibiting a comparable effect to that of the most drastic ALS-causing mutations known. Additionally, we show that the SOD1 monomer stability decreases upon oxidative stress, which may lead to partial local unfolding and consequently to increased aggregation propensity. Importantly, these results suggest that oxidative stress may play a key role in development of ALS, with the mutations in the SOD1 gene being an additional factor.

A comparison of pathway-independent and pathway-dependent methods in the calculation of conformational free enthalpy differences

The multistep umbrella sampling method, which belongs to pathway-dependent methods to calculate conformational free enthalpy differences, is used to calculate the free enthalpy difference between a right-handed 2.710/12 -helix and a left-handed 314 -helix of a hexa-β-peptide in methanol solution. The same conformational free enthalpy difference was previously investigated using pathway-independent methods such as direct counting and enveloping distribution sampling. Our results show that the pathway-dependent simulations are sensitive to the choice of the pathway and its parameter values. A pathway based on restraining distances of hydrogen-bonding atom pairs shows poor sampling for two different values of the restraining force constant. Another pathway based on restraining backbone dihedral angles did smoothly sample the transition between the two helical conformations, but only with a proper choice of the restraining force constant. The results illustrate that if, and only if, a proper pathway and proper parameters are chosen, the multistep umbrella sampling can be almost 50 times more efficient than the pathway-independent methods in this case. The analysis illustrates the advantages and pitfalls of the much used multistep umbrella sampling methodology.

Practical Aspects of Free-Energy Calculations: A Review

Free-energy calculations in the framework of classical molecular dynamics simulations are nowadays used in a wide range of research areas including solvation thermodynamics, molecular recognition, and protein folding. The basic components of a free-energy calculation, that is, a suitable model Hamiltonian, a sampling protocol, and an estimator for the free energy, are independent of the specific application. However, the attention that one has to pay to these components depends considerably on the specific application. Here, we review six different areas of application and discuss the relative importance of the three main components to provide the reader with an organigram and to make nonexperts aware of the many pitfalls present in free energy calculations.