I. Quantal effects in biochemical cooperativity and a proposed mechanism for the differentiation of calcium signaling in synaptic plasticity. II. Evolutionary algorithms for the optimization of methods in computational chemistry

In Part 1 of this thesis, we propose that biochemical cooperativity is a fundamentally non-ideal process. We show quantal effects underlying biochemical cooperativity and highlight apparent ergodic breaking at small volumes. The apparent ergodic breaking manifests itself in a divergence of deterministic and stochastic models. We further predict that this divergence of deterministic and stochastic results is a failure of the deterministic methods rather than an issue of stochastic simulations.

Ergodic breaking at small volumes may allow these molecular complexes to function as switches to a greater degree than has previously been shown. We propose that this ergodic breaking is a phenomenon that the synapse might exploit to differentiate Ca$^{2+}$ signaling that would lead to either the strengthening or weakening of a synapse. Techniques such as lattice-based statistics and rule-based modeling are tools that allow us to directly confront this non-ideality. A natural next step to understanding the chemical physics that underlies these processes is to consider \textit{in silico} specifically atomistic simulation methods that might augment our modeling efforts.

In the second part of this thesis, we use evolutionary algorithms to optimize \textit{in silico} methods that might be used to describe biochemical processes at the subcellular and molecular levels. While we have applied evolutionary algorithms to several methods, this thesis will focus on the optimization of charge equilibration methods. Accurate charges are essential to understanding the electrostatic interactions that are involved in ligand binding, as frequently discussed in the first part of this thesis.

Biochemical studies of postsynaptic density signaling proteins with a focus on synaptic GTPase activating protein and PDZ domains

Memory storage in the brain involves adjustment of the strength of existing synapses and formation of new neural networks. A key process underlying memory formation is synaptic plasticity, the ability of excitatory synapses to strengthen or weaken their connections in response to patterns of activity between their connected neurons. Synaptic plasticity is governed by the precise pattern of Ca²⁺ influx through postsynaptic N-methyl-D-aspartate-type glutamate receptors (NMDARs), which can lead to the activation of the small GTPases Ras and Rap. Differential activation of Ras and Rap acts to modulate synaptic strength by promoting the insertion or removal of 2-amino-3-(3-hydroxy-5-methyl-isoxazol-4-yl)propanoic acid receptors (AMPARs) from the synapse. Synaptic GTPase activating protein (synGAP) regulates AMPAR levels by catalyzing the inactivation of GTP-bound (active) Ras or Rap. synGAP is positioned in close proximity to the cytoplasmic tail regions of the NMDAR through its association with the PDZ domains of PSD-95. SynGAP’s activity is regulated by the prominent postsynaptic protein kinase, Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) and cyclin-dependent kinase 5 (CDK5), a known binding partner of CaMKII. Modulation of synGAP’s activity by phosphorylation may alter the ratio of active Ras to Rap in spines, thus pushing the spine towards the insertion or removal of AMPARs, subsequently strengthening or weakening the synapse. To date, all biochemical studies of the regulation of synGAP activity by protein kinases have utilized impure preparations of membrane bound synGAP. Here we have clarified the effects of phosphorylation of synGAP on its Ras and Rap GAP activities by preparing and utilizing purified, soluble recombinant synGAP, Ras, Rap, CaMKII, CDK5, PLK2, and CaM. Using mass spectrometry, we have confirmed the presence of previously identified CaMKII and CDK5 sites in synGAP, and have identified novel sites of phosphorylation by CaMKII, CDK5, and PLK2. We have shown that the net effect of phosphorylation of synGAP by CaMKII, CDK5, and PLK2 is an increase in its GAP activity toward HRas and Rap1. In contrast, there is no effect on its GAP activity toward Rap2. Additionally, by assaying the GAP activity of phosphomimetic synGAP mutants, we have been able to hypothesize the effects of CDK5 phosphorylation at specific sites in synGAP. In the course of this work, we also found, unexpectedly, that synGAP is itself a Ca²⁺/CaM binding protein. While Ca²⁺/CaM binding does not directly affect synGAP activity, it causes a conformational change in synGAP that increases the rate of its phosphorylation and exposes additional phosphorylation sites that are inaccessible in the absence of Ca²⁺/CaM.

The postsynaptic density (PSD) is an electron-dense region in excitatory postsynaptic neurons that contains a high concentration of glutamate receptors, cytoskeletal proteins, and associated signaling enzymes. Within the PSD, three major classes of scaffolding molecules function to organize signaling enzymes and glutamate receptors. PDZ domains present in the Shank and PSD-95 scaffolds families serve to physically link AMPARs and NMDARs to signaling molecules in the PSD. Because of the specificity and high affinity of PDZ domains for their ligands, I reasoned that these interacting pairs could provide the core components of an affinity chromatography system, including affinity resins, affinity tags, and elution agents. I show that affinity columns containing the PDZ domains of PSD-95 can be used to purify active PDZ domain-binding proteins to very high purity in a single step. Five heterologously expressed neuronal proteins containing endogenous PDZ domain ligands (NMDAR GluN2B subunit Tail, synGAP, neuronal nitric oxide synthase PDZ domain, cysteine rich interactor of PDZ three and cypin) were purified using PDZ domain resin, with synthetic peptides having the sequences of cognate PDZ domain ligands used as elution agents. I also show that conjugation of PDZ domain-related affinity tags to Proteins Of Interest (POIs) that do not contain endogenous PDZ domains or ligands does not alter protein activity and enables purification of the POIs on PDZ domain-related affinity resins.

Atomic-level structure and deformation in metallic glasses

Metallic glasses (MGs) are a relatively new class of materials discovered in 1960 and lauded for its high strengths and superior elastic properties. Three major obstacles prevent their widespread use as engineering materials for nanotechnology and industry: 1) their lack of plasticity mechanisms for deformation beyond the elastic limit, 2) their disordered atomic structure, which prevents effective study of their structure-to-property relationships, and 3) their poor glass forming ability, which limits bulk metallic glasses to sizes on the order of centimeters. We focused on understanding the first two major challenges by observing the mechanical properties of nanoscale metallic glasses in order to gain insight into its atomic-level structure and deformation mechanisms. We found that anomalous stable plastic flow emerges in room-temperature MGs at the nanoscale in wires as little as ~100 nanometers wide regardless of fabrication route (ion-irradiated or not). To circumvent experimental challenges in characterizing the atomic-level structure, extensive molecular dynamics simulations were conducted using approximated (embedded atom method) potentials to probe the underlying processes that give rise to plasticity in nanowires. Simulated results showed that mechanisms of relaxation via the sample free surfaces contribute to tensile ductility in these nanowires. Continuing with characterizing nanoscale properties, we studied the fracture properties of nano-notched MGnanowires and the compressive response of MG nanolattices at cryogenic (~130 K) temperatures. We learned from these experiments that nanowires are sensitive to flaws when the (amorphous) microstructure does not contribute stress concentrations, and that nano-architected structures with MG nanoribbons are brittle at low temperatures except when elastic shell buckling mechanisms dominate at low ribbon thicknesses (~20 nm), which instead gives rise to fully recoverable nanostructures regardless of temperature. Finally, motivated by understanding structure-to-property relationships in MGs, we studied the disordered atomic structure using a combination of in-situ X-ray tomography and X-ray diffraction in a diamond anvil cell and molecular dynamics simulations. Synchrotron X-ray experiments showed the progression of the atomic-level structure (in momentum space) and macroscale volume under increasing hydrostatic pressures. Corresponding simulations provided information on the real space structure, and we found that the samples displayed fractal scaling (r^d ∝ V, d < 3) at short length scales (< ~8 Å), and exhibited a crossover to a homogeneous scaling (d = 3) at long length scales. We examined this underlying fractal structure of MGs with parallels to percolation clusters and discuss the implications of this structural analogy to MG properties and the glass transition phenomenon.