Nonequilibrium molecular dynamics simulations of various mutants of hen egg white lysozyme have been performed at 300 K and 1 bar in the presence of both external static electric and low-frequency microwave (2.45 GHz) fields of varying intensity. Significant nonthermal field effects were noted, such as marked changes in the protein's secondary structure relative to the zero-field state, depending on the field conditions, mutation, and orientation with respect to the applied field. This occurred primarily as a consequence of alignment of the protein's total dipole moment with the external field, although the dipolar alignment of water molecules in both the solvation layer and the bulk was also found to be influential. Substantial differences in behavior were found for proteins with and without overall net charges, particularly with respect to translational motion. Localized motion and perturbation of hydrogen bonds were also found to be evident for charged residues.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease linked to the misfolding of Cu/Zn superoxide dismutase (SOD1).ALS-related defects in SOD1 result in a gain of toxic function that coincides with aberrant oligomerization. The structural events triggering oligomerization have remained enigmatic, however, as is the case in other protein-misfolding diseases. Here, we target the critical conformational change that defines the earliest step toward aggregation. Using nuclear spin relaxation dispersion experiments, we identified a short-lived (0.4 ms) and weakly populated (0.7%) conformation of metal-depleted SOD1 that triggers aberrant oligomerization. This excited state emanates from the folded ground state and is suppressed by metal binding, but is present in both the disulfide-oxidized and disulfide-reduced forms of the protein. Our results pinpoint a perturbed region of the excited-state structure that forms intermolecular contacts in the earliest nonnative dimer/ oligomer. The conformational transition that triggers oligomerization is a common feature of WT SOD1 and ALS-associated mutants that have widely different physicochemical properties. But compared with WT SOD1, the mutants have enhanced structural distortions in their excited states, and in some cases slightly higher excited-state populations and lower kinetic barriers, implying increased susceptibility to oligomerization. Our results provide a unified picture that highlights both (i) a common denominator among different SOD1 variants that may explain why diverse mutations cause the same disease, and (ii) a structural basis that may aid in understanding how different mutations affect disease propensity and progression.amyotrophic lateral sclerosis ͉ nuclear magnetic resonance relaxation ͉ protein misfolding

Nonequilibrium molecular dynamics simulations of a charge-neutral mutant of hen egg white lysozyme have been performed at 300 K and 1 bar in the presence of external microwave fields (2.45 to 100 GHz) of an rms electric field intensity of 0.05 V Å(-1). A systematic study was carried out of the distributions of persistence times and energies of each intraprotein hydrogen bond in between breakage and reformation, in addition to overall persistence over 20 ns simulations, vis-á-vis equilibrium, zero-field conditions. It was found that localized translational motion for formally charged residues led to greater disruption of associated hydrogen bonds, although induced rotational motion of strongly dipolar residues also led to a degree of hydrogen bond perturbation. These effects were most apparent in the solvent exposed exterior of hen egg white lysozyme, in which the intraprotein hydrogen bonds tend to be weaker.

FKBP12 (FK506 binding protein 12 kDa) is an important drug target. Nuclear magnetic resonance (NMR) order parameters, describing amplitudes of motion on the pico- to nanosecond time scale, can provide estimates of changes in conformational entropy upon ligand binding. Here we report backbone and methyl-axis order parameters of the apo and FK506-bound forms of FKBP12, based on N andH NMR relaxation. Binding of FK506 to FKBP12 results in localized changes in order parameters, notably for the backbone of residues E54 and I56 and the side chains of I56, I90, and I91, all positioned in the binding site. The order parameters increase slightly upon FK506 binding, indicating an unfavorable entropic contribution to binding of TΔ S = -18 ± 2 kJ/mol at 293 K. Molecular dynamics simulations indicate a change in conformational entropy, associated with all dihedral angles, of TΔ S = -26 ± 9 kJ/mol. Both these values are significant compared to the total entropy of binding determined by isothermal titration calorimetry and referenced to a reactant concentration of 1 mM ( TΔ S = -29 ± 1 kJ/mol). Our results reveal subtle differences in the response to ligand binding compared to that of the previously studied rapamycin-FKBP12 complex, despite the high degree of structural homology between the two complexes and their nearly identical ligand-FKBP12 interactions. These results highlight the delicate dependence of protein dynamics on drug interactions, which goes beyond the view provided by static structures, and reinforce the notion that protein conformational entropy can make important contributions to the free energy of ligand binding.

Non-equilibrium molecular dynamics simulations of a solvated 21-residue polyalanine (A21) peptide, featuring a high propensity for helix formation, have been performed at 300 K and 1 bar in the presence of external electromagnetic (e/m) fields in the microwave region (2.45 GHz) and an r.m.s. electric field intensity range of 0.01-0.05 V/Å. To investigate how the field presence affects transitions between the conformational states of a protein, we report 16 independent 40 ns-trajectories of A21 starting from both extended and fully folded states. We observe folding-behavior of the peptide consistent with prior simulation and experimental studies. The peptide displays a natural tendency to form stable elements of secondary structure which are stabilized by tertiary interactions with proximate regions of the peptide. Consistent with our earlier work, the presence of external e/m fields disrupts this behavior, involving a mechanism of localized dipolar alignment which serves to enhance intra-protein perturbations in hydrogen bonds (English, et al., J. Chem. Phys. 2010, 133, 091105), leading to more frequent transitions between shorter-lifetime states.