Fine structure of synapses on dendritic spines

Camillo Golgi's "Reazione Nera" led to the discovery of dendritic spines, small appendages originating from dendritic shafts. With the advent of electron microscopy (EM) they were identified as sites of synaptic contact. Later it was found that changes in synaptic strength were associated with changes in the shape of dendritic spines. While live-cell imaging was advantageous in monitoring the time course of such changes in spine structure, EM is still the best method for the simultaneous visualization of all cellular components, including actual synaptic contacts, at high resolution. Immunogold labeling for EM reveals the precise localization of molecules in relation to synaptic structures. Previous EM studies of spines and synapses were performed in tissue subjected to aldehyde fixation and dehydration in ethanol, which is associated with protein denaturation and tissue shrinkage. It has remained an issue to what extent fine structural details are preserved when subjecting the tissue to these procedures. In the present review, we report recent studies on the fine structure of spines and synapses using high-pressure freezing (HPF), which avoids protein denaturation by aldehydes and results in an excellent preservation of ultrastructural detail. In these studies, HPF was used to monitor subtle fine-structural changes in spine shape associated with chemically induced long-term potentiation (cLTP) at identified hippocampal mossy fiber synapses. Changes in spine shape result from reorganization of the actin cytoskeleton. We report that cLTP was associated with decreased immunogold labeling for phosphorylated cofilin (p-cofilin), an actin-depolymerizing protein. Phosphorylation of cofilin renders it unable to depolymerize F-actin, which stabilizes the actin cytoskeleton. Decreased levels of p-cofilin, in turn, suggest increased actin turnover, possibly underlying the changes in spine shape associated with cLTP. The findings reviewed here establish HPF as an appropriate method for studying the fine structure and molecular composition of synapses on dendritic spines.

Crystal structure of the sodium-proton antiporter NhaA dimer and new mechanistic insights

Sodium-proton antiporters rapidly exchange protons and sodium ions across the membrane to regulate intracellular pH, cell volume, and sodium concentration. How ion binding and release is coupled to the conformational changes associated with transport is not clear. Here, we report a crystal form of the prototypical sodium-proton antiporter NhaA from Escherichia coli in which the protein is seen as a dimer. In this new structure, we observe a salt bridge between an essential aspartic acid (Asp163) and a conserved lysine (Lys300). An equivalent salt bridge is present in the homologous transporter NapA, but not in the only other known crystal structure of NhaA, which provides the foundation of most existing structural models of electrogenic sodium-proton antiport. Molecular dynamics simulations show that the stability of the salt bridge is weakened by sodium ions binding to Asp164 and the neighboring Asp163. This suggests that the transport mechanism involves Asp163 switching between forming a salt bridge with Lys300 and interacting with the sodium ion. pKa calculations suggest that Asp163 is highly unlikely to be protonated when involved in the salt bridge. As it has been previously suggested that Asp163 is one of the two residues through which proton transport occurs, these results have clear implications to the current mechanistic models of sodium-proton antiport in NhaA.

Adsorption of biomedical coating molecules, amino acids, and short peptides on magnetite (110)

Superparamagnetic iron oxide nanoparticles for biomedical applications are usually coated with organic molecules to form a steric barrier against agglomeration. The stability of these coatings is well established in the synthesis medium but is more difficult to assess in physiological environment. To obtain a first theoretical estimate of their stability in such an environment, we perform density functional theory calculations of the adsorption of water, polyvinyl alcohol (PVA) and polyethylene glycol (PEG) coating molecules, as well as the monomer and dimer of glycine as a prototype short peptide, on the (110) surface of magnetite (Fe3O4) in vacuo. Our results show that PVA binds significantly stronger to the surface than both PEG and glycine, while the difference between the latter two is quite small. Depending on the coverage, the wateradsorption strength is intermediate between PVA and glycine. Due to its strongly interacting OH side groups, PVA is likely to remain bound to the surface in the presence of short peptides. This stability will have to be further assessed by molecular dynamics in the solvated state for which the present work forms the basis.

Molecular Machanics of Water Monomer-Maple Code

This is the Maple code to support the molecular dynamics of a water monomer molecule, allowing investigation of the classical vibrations of this molecule.

Molecular Mechanics of Water Monomer

Simulation of the classical molecular dynamics of a water molecule can be useful in explaining normal modes of motion, Fourier Transforms, and fundamental frequencies of vibration, as illustrated herein.

Targeting TRAF6 for Cancer Therapeutical Development

Tumor necrosis factor (TNF)-Receptor Associated Factors (TRAFs) are a family of signal transducer proteins. TRAF6 is a unique member of this family in that it is involved in not only the TNF superfamily, but the toll-like receptor (TLR)/IL-1R (TIR) superfamily. The formation of the complex consisting of Receptor Activator of Nuclear Factor κ B (RANK), with its ligand (RANKL) results in the recruitment of TRAF6, which activates NF-κB, JNK and MAP kinase pathways. TRAF6 is critical in signaling with leading to release of various growth factors in bone, and promotes osteoclastogenesis. TRAF6 has also been implicated as an oncogene in lung cancer and as a target in multiple myeloma. In the hopes of developing small molecule inhibitors of the TRAF6-RANK interaction, multiple steps were carried out. Computational prediction of hot spot residues on the protein-protein interaction of TRAF6 and RANK were examined. Three methods were used: Robetta, KFC2, and HotPoint, each of which uses a different methodology to determine if a residue is a hot spot. These hot spot predictions were considered the basis for resolving the binding site for in silico high-throughput screening using GOLD and the MyriaScreen database of drug/lead-like compounds. Computationally intensive molecular dynamics simulations highlighted the binding mechanism and TRAF6 structural changes upon hit binding. Compounds identified as hits were verified using a GST-pull down assay, comparing inhibition to a RANK decoy peptide. Since many drugs fail due to lack of efficacy and toxicity, predictive models for the evaluation of the LD50 and bioavailability of our TRAF6 hits, and these models can be used towards other drugs and small molecule therapeutics as well. Datasets of compounds and their corresponding bioavailability and LD50 values were curated based, and QSAR models were built using molecular descriptors of these compounds using the k-nearest neighbor (k-NN) method, and quality of these models were cross-validated.

CONFORMATIONAL DYNAMICS OF K-RAS AND H-RAS PROTEINS: IS THERE FUNCTIONAL SPECIFICITY AT THE CATALYTIC DOMAIN?

Ras proteins serve as crucial signaling modulators in cell proliferation through their ability to hydrolyze GTP and exist in a GTP “on” state and GTP “off” state. There are three different human Ras isoforms: H-ras, N-ras and K-ras (4A and 4B). Although their sequence identity is very high at the catalytic domain, these isoforms differ in their ability to activate different effectors and hence different signaling pathways. Much of the previous work on this topic has attributed this difference to the hyper variable region of Ras proteins, which contains most of the sequence variance among the isoforms and encodes specificity for differential distribution in the membrane. However, we hypothesize that sequence variation on lobe II of Ras catalytic domain alters dynamics and leads to differential preference for different effectors or modulators. In this work, we used all atom molecular dynamics to analyze the dynamics in the catalytic domain of H-ras and K-ras. We have also analyzed the dynamics of a transforming mutant of H-ras and K-ras and further studied the dynamics of an effectorselective mutant of H-ras. Collectively we have determined that wild type K-ras is more dynamic than H-ras and that the structure of the effector binding loop more closely resembles that of the T35S Raf-selective mutant, possibly giving us a new view and insight into the v mode of effector specificity. Furthermore we have determined that specific mutations at the same location perturb the conformational equilibrium differently in H-ras and K-ras and that an enhanced oncogenic potential may arise from different structural perturbations for each point mutation of a specific isoform.

Molecular studies of human high grade B cell non -Hodgkin's lymphomas

The non-Hodgkin's B cell lymphomas are a diverse group of neoplastic diseases. The incidence rate of the malignant tumors has been rising rapidly over the past twenty years in the United States and worldwide. The lack of insight to pathogenesis of the disease poses a significant problem in the early detection and effective treatment of the human malignancies. These studies attempted to investigate the molecular basis of pathogenesis of the human high grade B cell non-Hodgkin's lymphomas with a reverse genetic approach. The specific objective was to clone gene(s) which may play roles in development and progression of human high grade B cell non-Hodgkin's lymphomas.^ The messenger RNAs from two high grade B cell lymphoma lines, CJ and RR, were used for construction of cDNA libraries. Differential screening of the derived cDNA libraries yielded a 1.4 kb cDNA clone. The gene, designated as NHL-B1.4, was shown to be highly amplified and over-expressed in the high grade B cell lymphoma lines. It was not expressed in the peripheral blood lymphoid cells from normal donors. However, it was inducible in peripheral blood T lymphocytes by a T cell mitogen, PHA, but could not be activated in normal B cells by B cell mitogen PMA. Further molecular characterization revealed that the gene may have been rearranged in the RR and some other B cell lymphoma lines. The coding capacity of the cDNA has been confirmed by a rabbit reticulocyte lysate and wheat germ protein synthesis system. A recombinant protein with a molecular weight of approximate 30 kDa was visualized in autoradiogram. Polyclonal antisera have been generated by immunization of two rabbits with the NHL-B1.4 recombinant protein produced in the E. coli JM109. The derived antibody can recognize a natural protein with molecular weight of 49 kDa in cell lysate of activated peripheral T lymphocytes of normal donors and both the cell lysate and supernatant of RR B cell lymphoma lines. The possible biologic functions of the molecule has been tested preliminarily in a B lymphocyte proliferation assay. It was found that the Q-sepharose chromatograph purified supernatant of COS cell transfection could increase tritiated thymidine uptake by B lymphocytes but not by T lymphocytes. The B cell stimulatory activity of the supernatant of COS cell tranfection could be neutralized by the polyclonal antisera, indicating that the NHL-B1.4 gene product may be a molecule with BCGF-like activity.^ The expression profiles of NHL-B1.4 in normal and neoplastic lymphoid cells were consistent with the current B lymphocyte activation model and autocrine hypothesis of high grade B cell lymphomagenesis. These results suggested that the NHL-B1.4 cDNA may be a disease-related gene of human high grade B cell lymphomas, which may codes for a postulated B cell autocrine growth factor. ^

Influence of termal fluctuations in radiation damage cascade production and defect dynamics in tungsten

The detailed study of the deterioration suffered by the materials of the components of a nuclear facility, in particular those forming part of the reactor core, is a topic of great interest which importance derives in large technological and economic implications. Since changes in the atomic-structural properties of relevant components pose a risk to the smooth operation with clear consequences for security and life of the plant, controlling these factors is essential in any development of engineering design and implementation. In recent times, tungsten has been proposed as a structural material based on its good resistance to radiation, but still needs to be done an extensive study on the influence of temperature on the behavior of this material under radiation damage. This work aims to contribute in this regard. Molecular Dynamics (MD) simulations were carried out to determine the influence of temperature fluctuations on radiation damage production and evolution in Tungsten. We have particularly focused our study in the dynamics of defect creation, recombination, and diffusion properties. PKA energies were sampled in a range from 5 to 50 KeV. Three different temperature scenarios were analyzed, from very low temperatures (0-200K), up to high temperature conditions (300-500 K). We studied the creation of defects, vacancies and interstitials, recombination rates, diffusion properties, cluster formation, their size and evolution. Simulations were performed using Lammps and the Zhou EAM potential for W

Thermomechanical relaxation and different water states in cottonseed protein derived bioplastics

Thermomechanical relaxation events and different water states in cottonseed protein bioplastics are presented whilst investigating the effects of aldehyde cross-linking agents. Thermomechanical relaxation of cottonseed protein bioplastics associated with protein denaturation, moisture absorption and broad glass transitions (Tg) were observed from DSC and DMA measurements. It was shown that variation of the aldehyde influences the storage modulus at very low temperature (below Tg). From measurements of the water fusion point, enthalpy, vaporisation, and weight loss, three water states in the water-absorbed bioplastics are suggested; namely strongly-bound-to-polymer, weakly-bound-to-polymer and bulk-like water. The water content and unreacted cross-linking agents are influential factors in controlling formation of the different water states, whilst the selection of different aldehydes was found to be negligible. These results could be valuable for adjusting the thermomechanical relaxations of protein based bioplastics, and tailoring their properties in wet environments.

Measurements of attractive forces between proteins and end-grafted poly(ethylene glycol) chains

The surface force apparatus was used to measure directly the molecular forces between streptavidin and lipid bilayers displaying grafted Mr 2,000 poly(ethylene glycol) (PEG). These measurements provide direct evidence for the formation of relatively strong attractive forces between PEG and protein. At low compressive loads, the forces were repulsive, but they became attractive when the proteins were pressed into the polymer layer at higher loads. The adhesion was sufficiently robust that separation of the streptavidin and PEG uprooted anchored polymer from the supporting membrane. These interactions altered the properties of the grafted chains. After the onset of the attraction, the polymer continued to bind protein for several hours. The changes were not due to protein denaturation. These data demonstrate directly that the biological activity of PEG is not due solely to properties of simple polymers such as the excluded volume. It is also coupled to the competitive interactions between solvent and other materials such as proteins for the chain segments and to the ability of this material to adopt higher order intrachain structures.

A molecular level picture of the stabilization of A-DNA in mixed ethanol–water solutions

Advances in computer power, methodology, and empirical force fields now allow routine “stable” nanosecond-length molecular dynamics simulations of DNA in water. The accurate representation of environmental influences on structure remains a major, unresolved issue. In contrast to simulations of A-DNA in water (where an A-DNA to B-DNA transition is observed) and in pure ethanol (where disruption of the structure is observed), A-DNA in ≈85% ethanol solution remains in a canonical A-DNA geometry as expected. The stabilization of A-DNA by ethanol is likely due to disruption of the spine of hydration in the minor groove and the presence of ion-mediated interhelical bonds and extensive hydration across the major groove.

Crystal structure of Bacillus subtilis YabJ, a purine regulatory protein and member of the highly conserved YjgF family

The yabJ gene in Bacillus subtilis is required for adenine-mediated repression of purine biosynthetic genes in vivo and codes for an acid-soluble, 14-kDa protein. The molecular mechanism of YabJ is unknown. YabJ is a member of a large, widely distributed family of proteins of unknown biochemical function. The 1.7-Å crystal structure of YabJ reveals a trimeric organization with extensive buried hydrophobic surface and an internal water-filled cavity. The most important finding in the structure is a deep, narrow cleft between subunits lined with nine side chains that are invariant among the 25 most similar homologs. This conserved site is proposed to be a binding or catalytic site for a ligand or substrate that is common to YabJ and other members of the YER057c/YjgF/UK114 family of proteins.

Folding simulations of a three-stranded antiparallel β-sheet peptide

Protein folding is a grand challenge of the postgenomic era. In this paper, 58 folding events sampled during 47 molecular dynamics trajectories for a total simulation time of more than 4 μs provide an atomic detail picture of the folding of a 20-residue synthetic peptide with a stable three-stranded antiparallel β-sheet fold. The simulations successfully reproduce the NMR solution conformation, irrespective of the starting structure. The sampling of the conformational space is sufficient to determine the free energy surface and localize the minima and transition states. The statistically predominant folding pathway involves the formation of contacts between strands 2 and 3, starting with the side chains close to the turn, followed by association of the N-terminal strand onto the preformed 2–3 β-hairpin. The folding mechanism presented here, formation of a β-hairpin followed by consolidation, is in agreement with a computational study of the free energy surface of another synthetic three-stranded antiparallel β-sheet by Bursulaya and Brooks [(1999) J. Am. Chem. Soc. 121, 9947–9951]. Hence, it might hold in general for antiparallel β-sheets with short turns.

The importance of reactant positioning in enzyme catalysis: A hybrid quantum mechanics/molecular mechanics study of a haloalkane dehalogenase

Hybrid quantum mechanics/molecular mechanics calculations using Austin Model 1 system-specific parameters were performed to study the SN2 displacement reaction of chloride from 1,2-dichloroethane (DCE) by nucleophilic attack of the carboxylate of acetate in the gas phase and by Asp-124 in the active site of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10. The activation barrier for nucleophilic attack of acetate on DCE depends greatly on the reactants having a geometry resembling that in the enzyme or an optimized gas-phase structure. It was found in the gas-phase calculations that the activation barrier is 9 kcal/mol lower when dihedral constraints are used to restrict the carboxylate nucleophile geometry to that in the enzyme relative to the geometries for the reactants without dihedral constraints. The calculated quantum mechanics/molecular mechanics activation barriers for the enzymatic reaction are 16.2 and 19.4 kcal/mol when the geometry of the reactants is in a near attack conformer from molecular dynamics and in a conformer similar to the crystal structure (DCE is gauche), respectively. This haloalkane dehalogenase lowers the activation barrier for dehalogenation of DCE by 2–4 kcal/mol relative to the single point energies of the enzyme's quantum mechanics atoms in the gas phase. SN2 displacements of this sort in water are infinitely slower than in the gas phase. The modest lowering of the activation barrier by the enzyme relative to the reaction in the gas phase is consistent with mutation experiments.