Rajeev Prabhakar

Parameter calibration of transition-metal elements for the spin-polarized self-consistent-charge density-functional tight-binding (DFTB) method : Sc, Ti, Fe, Co, and Ni

Recently developed parameters for five first-row transition-metal elements (M = Sc, Ti, Fe, Co, and Ni) in combination with H, C, N, and O as well as the same metal (M-M) for the spin-polarized self-consistent-charge density-functional tight-binding (DFTB) method have been calibrated. To test their performance a couple sets of compounds have been selected to represent a variety of interactions and bonding schemes that occur frequently in transition-metal containing systems. The results show that the DFTB method with the present parameters in most cases reproduces structural properties very well, but the bond energies and the relative energies of different spin states only qualitatively compared to the B3LYP/SDD+6-31G(d) density functional (DFT) results. An application to the ONIOM(DFT:DFTB) indicates that DFTB works well as the low level method for the ONIOM calculation.

Comparative Molecular Dynamics Studies of Wild-Type and Oxidized Forms of Full-Length Alzheimer Amyloid β-Peptides Aβ(1−40) and Aβ(1−42)

In this study, all-atom 50 ns molecular dynamics simulations are performed on the full-length amyloid beta (Abeta) monomers (WT-Abeta(1-40) and WT-Abeta(1-42)) and their oxidized forms (Met35(O)-Abeta(1-40) and Met35(O)-Abeta(1-42)) in aqueous solution. The effects of the oxidation state of Met35 and the presence of dipeptide (Ile41-Ala42) on the secondary structures of the three distinct regions (the central hydrophobic core region 17-21 (LVFFA), the loop 23-28 (DVGSNK), and the second hydrophobic domain 29-35 (GAIIGLM)) of all monomers have been analyzed in detail, and results are compared with the available experimental information. Our simulations indicate that the WT-Abeta(1-40) monomer adopts an overall beta-hairpin-like structure, which is promoted by the turn region (24-27). This turn region is stabilized through salt-bridge formation between the Asp23 and Lys28 residues. In contrast, the overall structure of the oxidized (Met35(O)-Abeta(1-40)) monomer can be divided into three well-defined bend regions separated by coil segments. These structural differences may be critical for the measured decrease in the rate of aggregation of Met35(O)-Abeta(1-40) peptide. In the WT-Abeta(1-42) monomer, in comparison to the WT-Abeta(1-40), the Asp23-Lys28 salt bridge is absent, and consequently, the turn in the middle (24-27) region has a smaller curvature. The observed difference in the aggregation rates of these two peptides may be related to the opening of the turn (24-27) stabilized by the Asp23-Lys28 salt bridge. For WT-Abeta(1-42), in the absence of this salt bridge, the unfolding and aggregation events may be more favorable than for WT-Abeta(1-40).

Cotinine reduces amyloid-β aggregation and improves memory in Alzheimer's disease mice

Alzheimers disease (AD) affects millions of people world-wide and new effective and safe therapies are needed. Cotinine, the main metabolite of nicotine, has a long half-life and does not have cardiovascular or addictive side effects in humans. We studied the effect of cotinine on amyloid-β (Aβ) aggregation as well as addressed its impact on working and reference memories. Cotinine reduced Aβ deposition, improved working and reference memories, and inhibited Aβ oligomerization in the brains of transgenic (Tg) 6799 AD mice. In vitro studies confirmed the inhibitory effect of cotinine on Aβ1-42 aggregation. Cotinine stimulated Akt signaling, including the inhibition of glycogen synthase kinase 3β (GSK3β), which promotes neuronal survival and the synaptic plasticity processes underlying learning and memory in the hippocampus and cortex of wild type and Tg6799 AD mice. Simulation of the cotinine-Aβ1-42 complex using molecular dynamics showed that cotinine may interact with key histidine residues of Aβ1-42, altering its structure and inhibiting its aggregation. The good safety profile in humans and its beneficial effects suggest that cotinine may be an excellent therapeutic candidate for the treatment of AD.

Crystallographic, spectroscopic, and computational analysis of a flavin C4a-oxygen adduct in choline oxidase.

Flavin C4a-OO(H) and C4a-OH adducts are critical intermediates proposed in many flavoenzyme reaction mechanisms, but they are rarely detected even by rapid transient kinetics methods. We observe a trapped flavin C4a-OH or C4a-OO(H) adduct by single-crystal spectroscopic methods and in the 1.86 A resolution X-ray crystal structure of choline oxidase. The microspectrophotometry results show that the adduct forms rapidly in situ at 100 K upon exposure to X-rays. Density functional theory calculations establish the electronic structures for the flavin C4a-OH and C4a-OO(H) adducts and estimate the stabilization energy of several active site hydrogen bonds deduced from the crystal structure. We propose that the enzyme-bound FAD is reduced in the X-ray beam. The aerobic crystals then form either a C4a-OH or C4a-OO(H) adduct, but an insufficient proton inventory prevents their decay at cryogenic temperatures.

Unraveling the Photoluminescence Response of Light-Switching Ruthenium(II) Complexes Bound to Amyloid-β

Photoluminescent molecules are widely used for real-time monitoring of peptide aggregation. In this Article, we detail both experimental and computational modeling to elucidate the interaction between [Ru(bpy)2dppz](2+) and amyloid-β (Aβ(1-40)) aggregates. The transition from monomeric to fibrillar Aβ is of interest in the study of Alzheimers disease. Concentration-dependent experiments allowed the determination of a dissociation constant of 2.1 μM, while Job plots provided a binding stoichiometry of 2.6 Aβ monomers per [Ru(bpy)2dppz](2+). Our computational approach that combines molecular docking (both rigid and flexible) and all-atom molecular dynamics (MD) simulations predicts that the hydrophobic cleft between Val18 and Phe20 is a plausible binding site, which could also explain the increase in photoluminescence of [Ru(bpy)2dppz](2+) upon binding. This binding site is parallel to the fibril axis, in marked contrast to the binding site of these complexes in DNA (perpendicular to the DNA axis). Other binding sites may exist at the edges of the Aβ fibril, but they are actually of low abundance in an Aβ fibril several micrometers long. The assignment of the binding site was confirmed by binding studies in an Aβ fragment (Aβ(25-35)) that lacked the amino acids necessary to form the binding site. The agreement between the experimental and computational work is remarkable and provides a general model that can be used for studying the interaction of amyloid-binding molecules to Aβ.

Dimerization of the full-length Alzheimer amyloid β-peptide (Aβ42) in explicit aqueous solution: a molecular dynamics study.

In this study, the mechanism of dimerization of the full-length Alzheimer amyloid beta (Aβ42) peptide and structural properties of the three most stable dimers have been elucidated through 0.8 μs classical molecular dynamics (MD) simulations. The Aβ42 dimer has been reported to be the smallest neurotoxic species that adversely affects both memory and synaptic plasticity. On the basis of interactions between the distinct regions of the Aβ42 monomer, 10 different starting configurations were developed from their native folded structures. However, only six of them were found to form dimers and among them the three most stable (X(P), C-C(AP), and N-N(P)) were chosen for the detailed analysis. The structural properties of these dimers were compared with the available experimental and theoretical data. The MD simulations show that hydrophobic regions of both monomers play critical roles in the dimerization process. The high content of the α-helical structure in all the dimers is in line with its experimentally proposed role in the oligomerization. The formation of a zipper-like structure in X(P) is also in accordance with its existence in the aggregates of several short amyloidogenic peptides. The computed values of translational (D(T)) and rotational (D(R)) diffusion constants of 0.63 × 10(-6) cm(2)/s and 0.035 ns(-1), respectively, for this dimer are supported by the corresponding values of the Aβ42 monomer. These simulations have also elucidated several other key structural properties of these peptides. This information will be very useful to design small molecules for the inhibition and disruption of the critical Aβ42 dimers.

Computational modeling of substrate specificity and catalysis of the β-secretase (BACE1) enzyme.

In this combined MD simulation and DFT study, interactions of the wild-type (WT) amyloid precursor protein (APP) and its Swedish variant (SW), Lys670 → Asn and Met671 → Leu, with the beta-secretase (BACE1) enzyme and their cleavage mechanisms have been investigated. BACE1 catalyzes the rate-limiting step in the generation of 40-42 amino acid long Alzheimer amyloid beta (Aβ) peptides. All key structural parameters such as position of the flap, volume of the active site, electrostatic binding energy, structures, and positions of the inserts A, D, and F and 10s loop obtained from the MD simulations show that, in comparison to the WT-substrate, BACE1 exhibits greater affinity for the SW-substrate and orients it in a more reactive conformation. The enzyme-substrate models derived from the MD simulations were further utilized to investigate the general acid/base mechanism used by BACE1 to hydrolytically cleave these substrates. This mechanism proceeds through the following two steps: (1) formation of the gem-diol intermediate and (2) cleavage of the peptide bond. For the WT-substrate, the overall barrier of 22.4 kcal/mol for formation of the gem-diol intermediate is 3.3 kcal/mol higher than for the SW-substrate (19.1 kcal/mol). This process is found to be the rate-limiting in the entire mechanism. The computed barrier is in agreement with the measured barrier of ca. 18.00 kcal/mol for the WT-substrate and supports the experimental observation that the cleavage of the SW-substrate is 60 times more efficient than the WT-substrate.

Disulfide bridges remain intact while native insulin converts into amyloid fibrils

Amyloid fibrils are β-sheet-rich protein aggregates commonly found in the organs and tissues of patients with various amyloid-associated diseases. Understanding the structural organization of amyloid fibrils can be beneficial for the search of drugs to successfully treat diseases associated with protein misfolding. The structure of insulin fibrils was characterized by deep ultraviolet resonance Raman (DUVRR) and Nuclear Magnetic Resonance (NMR) spectroscopy combined with hydrogen-deuterium exchange. The compositions of the fibril core and unordered parts were determined at single amino acid residue resolution. All three disulfide bonds of native insulin remained intact during the aggregation process, withstanding scrambling. Three out of four tyrosine residues were packed into the fibril core, and another aromatic amino acid, phenylalanine, was located in the unordered parts of insulin fibrils. In addition, using all-atom MD simulations, the disulfide bonds were confirmed to remain intact in the insulin dimer, which mimics the fibrillar form of insulin.

Molecular dynamics study to investigate the effect of chemical substitutions of methionine 35 on the secondary structure of the amyloid β(Aβ(1-42)) monomer in aqueous solution

In this study, all-atom molecular dynamics simulations in the explicit water solvent are performed to investigate conformational changes in the secondary structure of the Abeta(1-42) monomer associated with the substitution of the Cgamma-methylene position of the Met35 amino acid residue by sulfoxide (Met35(O)), sulfone (Met35(O2)), and norleucine (Met35(CH2)). The effects of these substitutions on the structural changes that occur in three distinct regions (the central hydrophobic core (CHC) region 17-21 (LVFFA), stable turn segment 24-27 (VGSN), and second hydrophobic region 29-35 (GAIIGLM)) of all monomers have been analyzed in detail, and results are compared with experiments. Our 20 ns simulations indicate that the most significant changes take place in the second hydrophobic region of the Met35(O) and Met35(O2) monomers. However, for the Met35(CH2) monomer, this region does not exhibit significant deviations. In comparison to the wild-type (WT)-Abeta(1-42) monomer, for Met35(O) the second hydrophobic region is characterized by the formation of internal beta-sheets separated by stable turns, whereas for Met35(O2) it exhibits a more helical conformation. These substantial changes in the secondary structure can be explained in terms of an increase in the computed dipole moment and solvent accessible surface area (SASA) per residue of these substituents. These structural modifications can affect interaction between monomers, which in turn may influence the oligomerization process involved in Alzheimers disease (AD).

Protonation states of the catalytic dyad of β-secretase (BACE1) in the presence of chemically diverse inhibitors: a molecular docking study.

In this molecular docking study, the protonation states of the catalytic Asp dyad of the beta-secretase (BACE1) enzyme in the presence of eight chemically diverse inhibitors have been predicted. BACE1 catalyzes the rate-determining step in the generation of Alzheimer amyloid beta peptides and is widely considered as a promising therapeutic target. All the inhibitors were redocked into their corresponding X-ray structures using a combination of eight different protonation states of the Asp dyad for each inhibitor. Five inhibitors were primarily found to favor two different monoprotonated states, and the remaining three favor a dideprotonated state. In addition, five of them exhibited secondary preference for a diprotonated state. These results show that the knowledge of a single protonation state of the Asp dyad is not sufficient to search for the novel inhibitors of BACE1 and the most plausible state for each inhibitor must be determined prior to conducting in-silico screening.