Mapping the early steps in the pH-induced conformational conversion of the prion protein

Under certain conditions, the prion protein (PrP) undergoes a conformational change from the normal cellular isoform, PrPC, to PrPSc, an infectious isoform capable of causing neurodegenerative diseases in many mammals. Conversion can be triggered by low pH, and in vivo this appears to take place in an endocytic pathway and/or caveolae-like domains. It has thus far been impossible to characterize the conformational change at high resolution by experimental methods. Therefore, to investigate the effect of acidic pH on PrP conformation, we have performed 10-ns molecular dynamics simulations of PrPC in water at neutral and low pH. The core of the protein is well maintained at neutral pH. At low pH, however, the protein is more dynamic, and the sheet-like structure increases both by lengthening of the native β-sheet and by addition of a portion of the N terminus to widen the sheet by another two strands. The side chain of Met-129, a polymorphic codon in humans associated with variant Creutzfeldt–Jakob disease, pulls the N terminus into the sheet. Neutralization of Asp-178 at low pH removes interactions that inhibit conversion, which is consistent with the Asp-178–Asn mutation causing human prion diseases.

Protein folding from a highly disordered denatured state: The folding pathway of chymotrypsin inhibitor 2 at atomic resolution

Previous experimental and theoretical studies have produced high-resolution descriptions of the native and folding transition states of chymotrypsin inhibitor 2 (CI2). In similar fashion, here we use a combination of NMR experiments and molecular dynamics simulations to examine the conformations populated by CI2 in the denatured state. The denatured state is highly unfolded, but there is some residual native helical structure along with hydrophobic clustering in the center of the chain. The lack of persistent nonnative structure in the denatured state reduces barriers that must be overcome, leading to fast folding through a nucleation–condensation mechanism. With the characterization of the denatured state, we have now completed our description of the folding/unfolding pathway of CI2 at atomic resolution.

Comparison of the early stages of forced unfolding for fibronectin type III modules

The structural changes accompanying stretch-induced early unfolding events were investigated for the four type III fibronectin (FN-III) modules, FN-III7, FN-III8, FN-III9, and FN-III10 by using steered molecular dynamics. Simulations revealed that two main energy barriers, I and II, have to be overcome to initiate unraveling of FN-III's tertiary structure. In crossing the first barrier, the two opposing β-sheets of FN-III are rotated against each other such that the β-strands of both β-sheets align parallel to the force vector (aligned state). All further events in the unfolding pathway proceed from this intermediate state. A second energy barrier has to be overcome to break the first major cluster of hydrogen bonds between adjacent β-strands. Simulations revealed that the height of barrier I varied significantly among the four modules studied, being largest for FN-III7 and lowest for FN-III10, whereas the height of barrier II showed little variation. Key residues affecting the mechanical stability of FN-III modules were identified. These results suggest that FN-III modules can be prestretched into an intermediate state with only minor changes to their tertiary structures. FN-III10, for example, extends 12 Å from the native “twisted” to the intermediate aligned state, and an additional 10 Å from the aligned state to further unfolding where the first β-strand is peeled away. The implications of the existence of intermediate states regarding the elasticity of fibrillar fibers and the stretch-induced exposure of cryptic sites are discussed.

The arginine finger of RasGAP helps Gln-61 align the nucleophilic water in GAP-stimulated hydrolysis of GTP

The Ras family of GTPases is a collection of molecular switches that link receptors on the plasma membrane to signaling pathways that regulate cell proliferation and differentiation. The accessory GTPase-activating proteins (GAPs) negatively regulate the cell signaling by increasing the slow intrinsic GTP to GDP hydrolysis rate of Ras. Mutants of Ras are found in 25–30% of human tumors. The most dramatic property of these mutants is their insensitivity to the negative regulatory action of GAPs. All known oncogenic mutants of Ras map to a small subset of amino acids. Gln-61 is particularly important because virtually all mutations of this residue eliminate sensitivity to GAPs. Despite its obvious importance for carcinogenesis, the role of Gln-61 in the GAP-stimulated GTPase activity of Ras has remained a mystery. Our molecular dynamics simulations of the p21ras–p120GAP–GTP complex suggest that the local structure around the catalytic region can be different from that revealed by the x-ray crystal structure. We find that the carbonyl oxygen on the backbone of the arginine finger supplied in trans by p120GAP (Arg-789) interacts with a water molecule in the active site that is forming a bridge between the NH2 group of the Gln-61 and the γ-phosphate of GTP. Thus, Arg-789 may play a dual role in generating the nucleophile as well as stabilizing the transition state for P—O bond cleavage.

Direct evidence for modified solvent structure within the hydration shell of a hydrophobic amino acid.

Neutron scattering experiments are used to determine scattering profiles for aqueous solutions of hydrophobic and hydrophilic amino acid analogs. Solutions of hydrophobic solutes show a shift in the main diffraction peak to smaller angle as compared with pure water, whereas solutions of hydrophilic solutes do not. The same difference for solutions of hydrophobic and hydrophilic side chains is also predicted by molecular dynamics simulations. The neutron scattering curves of aqueous solutions of hydrophobic amino acids at room temperature are qualitatively similar to differences between the liquid molecular structure functions measured for ambient and supercooled water. The nonpolar solute-induced expansion of water structure reported here is also complementary to recent neutron experiments where compression of aqueous solvent structure has been observed at high salt concentration.

Shared fine specificity between T-cell receptors and an antibody recognizing a peptide/major histocompatibility class I complex.

Cytotoxic T cells recognize mosaic structures consisting of target peptides embedded within self-major histocompatibility complex (MHC) class I molecules. This structure has been described in great detail for several peptide-MHC complexes. In contrast, how T-cell receptors recognize peptide-MHC complexes have been less well characterized. We have used a complete set of singly substituted analogs of a mouse MHC class I, Kk-restricted peptide, influenza hemagglutinin (Ha)255-262, to address the binding specificity of this MHC molecule. Using the same peptide-MHC complexes we determined the fine specificity of two Ha255-262-specific, Kk-restricted T cells, and of a unique antibody, pSAN, specific for the same peptide-MHC complex. Independently, a model of the Ha255-262-Kk complex was generated through homology modeling and molecular mechanics refinement. The functional data and the model corroborated each other showing that peptide residues 1, 3, 4, 6, and 7 were exposed on the MHC surface and recognized by the T cells. Thus, the majority, and perhaps all, of the side chains of the non-primary anchor residues may be available for T-cell recognition, and contribute to the stringent specificity of T cells. A striking similarity between the specificity of the T cells and that of the pSAN antibody was found and most of the peptide residues, which could be recognized by the T cells, could also be recognized by the antibody.

Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

Interresidue hydrogen bonding in a peptide nucleic acid.RNA heteroduplex.

Previous molecular mechanics calculations suggest that strands of peptide nucleic acids (PNAs) and complementary oligonucleotides form antiparallel duplexes stabilized by interresidue hydrogen bonds. In the computed structures, the amide carbonyl oxygen nearest the nucleobase (O7') forms an interresidue hydrogen bond with the backbone amide proton of the following residue, (n + 1)H1'. Of the 10 published two dimensional 1H NMR structures of a hexameric PNA.RNA heteroduplex. PNA(GAACTC).r(GAGUUC), 9 exhibit two to five potential interresidue hydrogen bonds. In our minimized average structure, created from the coordinates of these 10 NMR structures, three of the five possible interresidue hydrogen bond sites within the PNA backbone display the carbonyl oxygen (O7') and the amide proton (n + 1)H1' distances and N1'-H1'-(n - 1)O7' angles optimal for hydrogen bond formation. The finding of these interresidue hydrogen bonds supports the results of our previous molecular mechanics calculations.

Mechanical annealing of metallic electrodes at the atomic scale

The process of creating an atomically defined and robust metallic tip is described and quantified using measurements of contact conductance between gold electrodes and numerical simulations. Our experiments show how the same conductance behavior can be obtained for hundreds of cycles of formation and rupture of the nanocontact by limiting the indentation depth between the two electrodes up to a conductance value of approximately 5G0 in the case of gold. This phenomenon is rationalized using molecular dynamics simulations together with density functional theory transport calculations which show how, after repeated indentations (mechanical annealing), the two metallic electrodes are shaped into tips of reproducible structure. These results provide a crucial insight into fundamental aspects relevant to nanotribology or scanning probe microscopies.

Mechanical, electrical, and magnetic properties of Ni nanocontacts

The dynamic deformation upon stretching of Ni nanowires as those formed with mechanically controllable break junctions or with a scanning tunneling microscope is studied both experimentally and theoretically. Molecular dynamics simulations of the breaking process are performed. In addition, and in order to compare with experiments, we also compute the transport properties in the last stages before failure using the first-principles implementation of Landauer's formalism included in our transport package ALACANT.

Understanding the structure of the first atomic contact in gold

We have studied experimentally jump-to-contact (JC) and jump-out-of-contact (JOC) phenomena in gold electrodes. JC can be observed at first contact when two metals approach each other, while JOC occurs in the last contact before breaking. When the indentation depth between the electrodes is limited to a certain value of conductance, a highly reproducible behaviour in the evolution of the conductance can be obtained for hundreds of cycles of formation and rupture. Molecular dynamics simulations of this process show how the two metallic electrodes are shaped into tips of a well-defined crystallographic structure formed through a mechanical annealing mechanism. We report a detailed analysis of the atomic configurations obtained before contact and rupture of these stable structures and obtained their conductance using first-principles quantum transport calculations. These results help us understand the values of conductance obtained experimentally in the JC and JOC phenomena and improve our understanding of atomic-sized contacts and the evolution of their structural characteristics.

Mécanismes de transfert de proton d’une réaction acido-basique en phase aqueuse : une étude ab-initio

Les réactions de transfert de proton se retrouvent abondamment dans la nature et sont des processus cruciaux dans plusieurs réactions chimiques et biologiques, qui se produisent souvent en milieu aqueux. Les mécanismes régissant ces échanges de protons sont complexes et encore mal compris, suscitant un intérêt des chercheurs en vue d’une meilleure compréhension fondamentale du processus de transfert. Le présent manuscrit présente une étude mécanistique portant sur une réaction de transfert de proton entre un acide (phénol fonctionnalisé) et une base (ion carboxylate) en phase aqueuse. Les résultats obtenus sont basés sur un grand nombre de simulations de dynamique moléculaire ab-initio réalisées pour des systèmes de type « donneur-pont-accepteur », où le pont se trouve à être une unique molécule d’eau, permettant ainsi l’élaboration d’un modèle cinétique détaillé pour le système étudié. La voie de transfert principalement observée est un processus ultra-rapide (moins d’une picoseconde) passant par la formation d’une structure de type « Eigen » (H9O4+) pour la molécule d’eau pontante, menant directement à la formation des produits. Une seconde structure de la molécule d’eau pontante est également observée, soit une configuration de type « Zündel » (H5O2+) impliquant l’accepteur de proton (l’ion carboxylate) qui semble agir comme un cul-de-sac pour la réaction de transfert de proton.

Mécanismes de transfert de proton d’une réaction acido-basique en phase aqueuse : une étude ab-initio

Les réactions de transfert de proton se retrouvent abondamment dans la nature et sont des processus cruciaux dans plusieurs réactions chimiques et biologiques, qui se produisent souvent en milieu aqueux. Les mécanismes régissant ces échanges de protons sont complexes et encore mal compris, suscitant un intérêt des chercheurs en vue d’une meilleure compréhension fondamentale du processus de transfert. Le présent manuscrit présente une étude mécanistique portant sur une réaction de transfert de proton entre un acide (phénol fonctionnalisé) et une base (ion carboxylate) en phase aqueuse. Les résultats obtenus sont basés sur un grand nombre de simulations de dynamique moléculaire ab-initio réalisées pour des systèmes de type « donneur-pont-accepteur », où le pont se trouve à être une unique molécule d’eau, permettant ainsi l’élaboration d’un modèle cinétique détaillé pour le système étudié. La voie de transfert principalement observée est un processus ultra-rapide (moins d’une picoseconde) passant par la formation d’une structure de type « Eigen » (H9O4+) pour la molécule d’eau pontante, menant directement à la formation des produits. Une seconde structure de la molécule d’eau pontante est également observée, soit une configuration de type « Zündel » (H5O2+) impliquant l’accepteur de proton (l’ion carboxylate) qui semble agir comme un cul-de-sac pour la réaction de transfert de proton.

Developments in the Theory of X-ray Absorption and Compton Scattering.

Thesis (Ph.D.)--University of Washington, 2016-06