Elizabeth A. Proctor

Structural and dynamic determinants of protein-peptide recognition.

Protein-peptide interactions play important roles in many cellular processes, including signal transduction, trafficking, and immune recognition. Protein conformational changes upon binding, an ill-defined peptide binding surface, and the large number of peptide degrees of freedom make the prediction of protein-peptide interactions particularly challenging. To address these challenges, we perform rapid molecular dynamics simulations in order to examine the energetic and dynamic aspects of protein-peptide binding. We find that, in most cases, we recapitulate the native binding sites and native-like poses of protein-peptide complexes. Inclusion of electrostatic interactions in simulations significantly improves the prediction accuracy. Our results also highlight the importance of protein conformational flexibility, especially side-chain movement, which allows the peptide to optimize its conformation. Our findings not only demonstrate the importance of sufficient sampling of the protein and peptide conformations, but also reveal the possible effects of electrostatics and conformational flexibility on peptide recognition.

Glutathionylation at Cys-111 Induces Dissociation of Wild Type and FALS Mutant SOD1 Dimers

Mutation of the ubiquitous cytosolic enzyme Cu/Zn superoxide dismutase (SOD1) is hypothesized to cause familial amyotrophic lateral sclerosis (FALS) through structural destabilization leading to misfolding and aggregation. Considering the late onset of symptoms as well as the phenotypic variability among patients with identical SOD1 mutations, it is clear that nongenetic factor(s) impact ALS etiology and disease progression. Here we examine the effect of Cys-111 glutathionylation, a physiologically prevalent post-translational oxidative modification, on the stabilities of wild type SOD1 and two phenotypically diverse FALS mutants, A4V and I112T. Glutathionylation results in profound destabilization of SOD1(WT) dimers, increasing the equilibrium dissociation constant K(d) to ~10-20 μM, comparable to that of the aggressive A4V mutant. SOD1(A4V) is further destabilized by glutathionylation, experiencing an ~30-fold increase in K(d). Dissociation kinetics of glutathionylated SOD1(WT) and SOD1(A4V) are unchanged, as measured by surface plasmon resonance, indicating that glutathionylation destabilizes these variants by decreasing association rate. In contrast, SOD1(I112T) has a modestly increased dissociation rate but no change in K(d) when glutathionylated. Using computational structural modeling, we show that the distinct effects of glutathionylation on different SOD1 variants correspond to changes in composition of the dimer interface. Our experimental and computational results show that Cys-111 glutathionylation induces structural rearrangements that modulate stability of both wild type and FALS mutant SOD1. The distinct sensitivities of SOD1 variants to glutathionylation, a modification that acts in part as a coping mechanism for oxidative stress, suggest a novel mode by which redox regulation and aggregation propensity interact in ALS.

Discrete molecular dynamics

Discrete molecular dynamics (DMD) has emerged as a simplification of traditional molecular dynamics (MD). DMD employs discrete step function potentials in place of the continuous potential used in traditional MD. As a result, the simulation engine solves the ballistic equations of motion for only those particles participating in a collision, instead of solving Newtons equations of motion for every particle in the system. Because fewer calculations are performed per time step, the DMD technique allows for longer time and length scales to become accessible in the simulation of large biomolecules. The use of coarse‐grained models extends the computational advantage of this method. Although some accuracy is sacrificed to speed, because of the usefulness of DMD to the simulation of many particles at longer timescales, the technique has seen application to diverse molecular systems.

Fast screening of protein surfaces using geometric invariant fingerprints

We develop a rapid and efficient method for the comparison of protein local surface similarities using geometric invariants (fingerprints). By combining fast fingerprint comparison with explicit alignment, we successfully screen the entire Protein Data Bank for proteins that possess local surface similarities. Our method is independent of sequence and fold similarities, and has potential application to protein structure annotation and protein-protein interface design.

Nonnative SOD1 trimer is toxic to motor neurons in a model of amyotrophic lateral sclerosis

Significance Protein aggregation is a hallmark of neurodegenerative disease and is hypothesized to cause neuron death. Despite extensive study of disease-associated aggregating proteins, mechanisms of neuron death remain a mystery, and no cures or effective treatments yet exist. Here, we demonstrate the toxicity of a small aggregate of the Cu,Zn superoxide dismutase (SOD1) protein, associated with amyotrophic lateral sclerosis (ALS). We present an experimentally verified structural model of this toxic species and show that SOD1 mutants designed to promote formation of this aggregate increase cell death, providing a direct link between aggregate presence and neuron death. Knowledge of toxic species and the ability to manipulate their formation provides a valuable direction for pursuit of therapeutic strategies in ALS. Since the linking of mutations in the Cu,Zn superoxide dismutase gene (sod1) to amyotrophic lateral sclerosis (ALS) in 1993, researchers have sought the connection between SOD1 and motor neuron death. Disease-linked mutations tend to destabilize the native dimeric structure of SOD1, and plaques containing misfolded and aggregated SOD1 have been found in the motor neurons of patients with ALS. Despite advances in understanding of ALS disease progression and SOD1 folding and stability, cytotoxic species and mechanisms remain unknown, greatly impeding the search for and design of therapeutic interventions. Here, we definitively link cytotoxicity associated with SOD1 aggregation in ALS to a nonnative trimeric SOD1 species. We develop methodology for the incorporation of low-resolution experimental data into simulations toward the structural modeling of metastable, multidomain aggregation intermediates. We apply this methodology to derive the structure of a SOD1 trimer, which we validate in vitro and in hybridized motor neurons. We show that SOD1 mutants designed to promote trimerization increase cell death. Further, we demonstrate that the cytotoxicity of the designed mutants correlates with trimer stability, providing a direct link between the presence of misfolded oligomers and neuron death. Identification of cytotoxic species is the first and critical step in elucidating the molecular etiology of ALS, and the ability to manipulate formation of these species will provide an avenue for the development of future therapeutic strategies.

Discrete Molecular Dynamics Distinguishes Nativelike Binding Poses from Decoys in Difficult Targets

Virtual screening is one of the major tools used in computer-aided drug discovery. In structure-based virtual screening, the scoring function is critical to identifying the correct docking pose and accurately predicting the binding affinities of compounds. However, the performance of existing scoring functions has been shown to be uneven for different targets, and some important drug targets have proven especially challenging. In these targets, scoring functions cannot accurately identify the native or near-native binding pose of the ligand from among decoy poses, which affects both the accuracy of the binding affinity prediction and the ability of virtual screening to identify true binders in chemical libraries. Here, we present an approach to discriminating native poses from decoys in difficult targets for which several scoring functions failed to correctly identify the native pose. Our approach employs Discrete Molecular Dynamics simulations to incorporate protein-ligand dynamics and the entropic effects of binding. We analyze a collection of poses generated by docking and find that the residence time of the ligand in the native and nativelike binding poses is distinctly longer than that in decoy poses. This finding suggests that molecular simulations offer a unique approach to distinguishing the native (or nativelike) binding pose from decoy poses that cannot be distinguished using scoring functions that evaluate static structures. The success of our method emphasizes the importance of protein-ligand dynamics in the accurate determination of the binding pose, an aspect that is not addressed in typical docking and scoring protocols.

Applications of Discrete Molecular Dynamics in biology and medicine

Discrete Molecular Dynamics (DMD) is a physics-based simulation method using discrete energetic potentials rather than traditional continuous potentials, allowing microsecond time scale simulations of biomolecular systems to be performed on personal computers rather than supercomputers or specialized hardware. With the ongoing explosion in processing power even in personal computers, applications of DMD have similarly multiplied. In the past two years, researchers have used DMD to model structures of disease-implicated protein folding intermediates, study assembly of protein complexes, predict protein-protein binding conformations, engineer rescue mutations in disease-causative protein mutants, design a protein conformational switch to control cell signaling, and describe the behavior of polymeric dispersants for environmental cleanup of oil spills, among other innovative applications.

Cu,Zn-Superoxide Dismutase without Zn Is Folded but Catalytically Inactive

Amyotrophic lateral sclerosis has been linked to the gain of aberrant function of superoxide dismutase, Cu,Zn-SOD1 upon protein misfolding. The mechanism of SOD1 misfolding is thought to involve mutations leading to the loss of Zn, followed by protein unfolding and aggregation. We show that the removal of Zn from SOD1 may not lead to an immediate unfolding but immediately deactivates the enzyme through a combination of subtle structural and electronic effects. Using quantum mechanics/discrete molecular dynamics, we showed that both Zn-less wild-type (WT)-SOD1 and its D124N mutant that does not bind Zn have at least metastable folded states. In those states, the reduction potential of Cu increases, leading to the presence of detectable amounts of Cu(I) instead of Cu(II) in the active site, as confirmed experimentally. The Cu(I) protein cannot participate in the catalytic Cu(I)-Cu(II) cycle. However, even without the full reduction to Cu(I), the Cu site in the Zn-less variants of SOD1 is shown to be catalytically incompetent: unable to bind superoxide in a way comparable to the WT-SOD1. The changes are more radical and different in the D124N Zn-less mutant than in the Zn-less WT-SOD1, suggesting D124N being perhaps not the most adequate model for Zn-less SOD1. Overall, Zn in SOD1 appears to be influencing the Cu site directly by adjusting its reduction potential and geometry. Thus, the role of Zn in SOD1 is not just structural, as was previously thought; it is a vital part of the catalytic machinery.

Identification of neurotoxic cytokines by profiling Alzheimer’s disease tissues and neuron culture viability screening

Alzheimer’s disease (AD) therapeutics based on the amyloid hypothesis have shown minimal efficacy in patients, suggesting that the activity of amyloid beta (Aβ) represents only one aspect of AD pathogenesis. Since neuroinflammation is thought to play an important role in AD, we hypothesized that cytokines may play a direct role in promoting neuronal death. Here, we profiled cytokine expression in a small cohort of human AD and control brain tissues. We identified AD-associated cytokines using partial least squares regression to correlate cytokine expression with quantified pathologic disease state and then used neuron cultures to test whether cytokines up-regulated in AD tissues could affect neuronal viability. This analysis identified cytokines that were associated with the pathological severity. Of the top correlates, only TNF-α reduced viability in neuron culture when applied alone. VEGF also reduced viability when applied together with Aβ, which was surprising because VEGF has been viewed as a neuro-protective protein. We found that this synthetic pro-death effect of VEGF in the context of Aβ was commensurate with VEGFR-dependent changes in multiple signaling pathways that govern cell fate. Our findings suggest that profiling of tissues combined with a culture-based screening approach can successfully identify new mechanisms driving neuronal death.

Traffic within the Cytochrome b6f Lipoprotein Complex: Gating of the Quinone Portal

The cytochrome bc complexes b6f and bc1 catalyze proton-coupled quinol/quinone redox reactions to generate a transmembrane proton electrochemical gradient. Quinol oxidation on the electrochemically positive (p) interface of the complex occurs at the end of a narrow quinol/quinone entry/exit Qp portal, 11 Å long in bc complexes. Superoxide, which has multiple signaling functions, is a by-product of the p-side quinol oxidation. Although the transmembrane core and the chemistry of quinone redox reactions are conserved in bc complexes, the rate of superoxide generation is an order of magnitude greater in the b6f complex, implying that functionally significant differences in structure exist between the b6f and bc1 complexes on the p-side. A unique structure feature of the b6f p-side quinol oxidation site is the presence of a single chlorophyll-a molecule whose function is unrelated to light harvesting. This study describes a cocrystal structure of the cytochrome b6f complex with the quinol analog stigmatellin, which partitions in the Qp portal of the bc1 complex, but not effectively in b6f. It is inferred that the Qp portal is partially occluded in the b6f complex relative to bc1. Based on a discrete molecular-dynamics analysis, occlusion of the Qp portal is attributed to the presence of the chlorophyll phytyl tail, which increases the quinone residence time within the Qp portal and is inferred to be a cause of enhanced superoxide production. This study attributes a novel (to our knowledge), structure-linked function to the otherwise enigmatic chlorophyll-a in the b6f complex, which may also be relevant to intracellular redox signaling.