Road-blocker HSP disease mutation disrupts pre-organization for ATP hydrolysis in kinesin through a second sphere control.

Kinesin motor proteins perform several essential cellular functions powered by the adenosine triphosphate (ATP) hydrolysis reaction. Several single-point mutations in the kinesin motor protein KIF5A have been implicated to hereditary spastic paraplegia disease (HSP), a lethal neurodegenerative disease in humans. In earlier studies, we have shown that a series of HSP-related mutations can impair the kinesin's long-distance displacement or processivity by modulating the order-disorder transition of the linker connecting the heads to the coiled coil. On the other hand, the reduction of kinesin's ATP hydrolysis reaction rate by a distal asparagine-to-serine mutation is also known to cause HSP disease. However, the molecular mechanism of the ATP hydrolysis reaction in kinesin by this distal mutation is still not fully understood. Using classical molecular dynamics simulations combined with quantum mechanics/molecular mechanics calculations, the pre-organization geometry required for optimal hydrolysis in kinesin motor bound to α/ß-tubulin is determined. This optimal geometry has only a single salt-bridge (of the possible two) between Arg203-Glu236, putting a reactive water molecule at a perfect position for hydrolysis. Such geometry is also needed to create the appropriate configuration for proton translocation during ATP hydrolysis. The distal asparagine-to-serine mutation is found to disrupt this optimal geometry. Therefore, the current study along with our previous one demonstrates how two different effects on kinesin dynamics (processivity and ATP hydrolysis), caused by a different set of genotypes, can give rise to the same phenotype leading to HSP disease.

The Role of Copper Salts and O2 in the Mechanism of C≡N Bond Activation for Facilitating Nitrogen Transfer Reactions.

C≡N bond scission can be a potential avenue for the functionalization of chemical bonds. We have conducted a computational study, using density functional theory (DFT) and ab initio multireference CASSCF methods, to unravel the intricate mechanistic pathways traversed in the copper-promoted, dioxygen-assisted reaction for the formation of aryl isocyanate species from aryl aldehyde. This aryl isocyanate species acts as an active species for C≡N bond cleavage of coordinated cyanide anion enabling nitrogen transfer to various aldehydes. Electronic structure analysis revealed that under all the reaction conditions radical-based pathways are operative, which is in agreement with the experimental findings. The major driving force is a CuII/I redox cycle initiated by single-electron transfer from the carbon center of the nitrile moiety. Our study reveals that the copper salts act as the "electron pool" in this unique nitrogen transfer reaction forming an aryl isocyanate species from aryl aldehydes.

Externally Regulated Specific Molecular Recognition Driven Pathway Selectivity in Supramolecular Polymerization.

This article reveals 4-dimethylaminopyridine (DMAP) regulated pathway selectivity in the supramolecular polymerization of a naphthalene-diimide derivative (NDI-1), appended with a carboxylic acid group. In decane, NDI-1 produces ill-defined aggregate (Agg-1) due to different H-bonding motifs of the -COOH group. With one mole equivalent DMAP, the NDI-1/DMAP complex introduces new nucleation condition and exhibits a cooperative supramolecular polymerization producing J-aggregated fibrillar nanostructure (Agg-2). With 10 % DMAP and fast cooling (10 K/min), similar nucleation and open chain H-bonding with the free monomer in an anti-parallel arrangement produces identical J-aggregate (Agg-2a). With 2.5 % DMAP and slow cooling (1 K/min), a distinct nucleation and supramolecular polymerization pathway emerge leading to the thermodynamically controlled Agg-3 with face-to-face stacking and 2D-morphology. Slow cooling with 5-10 % DMAP produces a mixture of Agg-2a and Agg-3. Computational modelling studies provide valuable insights into the internal order and the pathway complexity.

Mechanistic study of the ATP hydrolysis reaction in dynein motor protein.

Dynein, a large and complex motor protein, harnesses energy from adenosine triphosphate (ATP) hydrolysis to regulate essential cellular activities. The ATP hydrolysis mechanism for the dynein motor is still shrouded in mystery. Herein, molecular dynamics simulations of a dynein motor disclosed that two water molecules are present close to the γ-phosphate of ATP and Glu1742 at the AAA1 site of dynein. We have proposed three possible mechanisms for the ATP hydrolysis. We divulge by using a quantum mechanics/molecular mechanics (QM/MM) study that two water molecules and Glu1742 are crucial for facilitating the ATP hydrolysis reaction in dynein. Moreover, the ATP hydrolysis step is initiated by the activation of lytic water (W1) by Glu1742 through relay proton transfers with the help of auxiliary water (W2) yielding HPO42- and ADP, as a product. In the next step, a proton is shifted back from Glu1742 to generate inorganic phosphate (H2PO4-) via another relay proton transfer event. The overall activation barrier for the Glu1742 assisted ATP hydrolysis is found to be the most favourable pathway compared to other plausible pathways. We also unearthed that ATP hydrolysis in dynein follows a so-called associative-like pathway in its rate-limiting step. Our study ascertained the important indirect roles of the two amino acids (such as Arg2109, Asn1792) and Mg2+ ion in the ATP hydrolysis of dynein. Additionally, multiple sequence alignment of the different organisms of dynein motors has conveyed the evolutionary importance of the Glu1742, Asn1742, and Arg2109 residues, respectively. As similar mechanisms are also prevalent in other motors, and GTPase and ATPase enzymes, the present finding spells out the definitive requirement of a proton relay process through an extended water-chain as one of the key components in an enzymatic ATP hydrolysis reaction.

A computational study of the dechlorination of ß-hexachlorocyclohexane (ß-HCH) catalyzed by the haloalkane dehalogenase LinB.

LinB, a haloalkane dehalogenase from Sphingomonas paucimobilis UT26, is known to metabolize halohydrocarbons to halide ions and the respective alcohols. Its broad substrate specificity allowed its consideration for bioremediation. Herein, we have shown its catalytic action toward ß-hexachlorocyclohexane (ß-HCH) - an example of large-size substrates that can be accommodated in its active site. We have analyzed the capability of combined QM/MM schemes to describe in detail the SN2 dechlorination reaction between ß-HCH and Asp108 in the active site of LinB. Free energy surfaces have been calculated using one and two dimensional potentials of mean force (PMF) obtained at the PM3/MM (MM=amberff99SB, TIP3P) level of theory. The overestimated energetic barriers by the PM3 Hamiltonian were corrected using a DFT functional (M06-2X). The resulted activation energies (16 and 19 kcal mol(-1) from 1D and 2D-PMF profiles, respectively) for the dechlorination reaction of ß-HCH in the active site of LinB enzyme are in qualitative agreement with the experimentally determined value of 17 kcal mol(-1). The binding of ß-HCH to the active site of LinB has been compared to the binding of smaller 1-chlorobutane (1-CB) and larger δ-hexabromocyclododecane (δ-HBCD).

Conformational Effects on the Absorption Spectra of Phytochromes.

Bacteriophytochrome is a photoreceptor protein that contains the biliverdin (BV) chromophore as its active component. The spectra of BV upon mutation remain remarkably unchanged, as far as spectral positions are concerned. This points toward the minimal effect of electrostatic effects on the electronic structure of the chromophore. However, the relative intensities of the Q and Soret bands of the chromophore change dramatically upon mutation. In this work, we delve into the molecular origin of this unusual intensity modulation. Using extensive classical MD and QM/MM calculations, we show that due to mutation, the conformational population of the chromophore changes significantly. The noncovalent interactions, especially the stacking interactions, lead to extra stabilization of the cyclic form in the D207H mutated species as opposed to the open form in the wild-type BV. Thus, unlike the commonly observed direct electrostatic effect on the spectral shift, in the case of BV the difference observed is in varying intensities, and this in turn is driven by a conformational shift due to enhanced stacking interaction.

Impact of the Hydrogen-Bonding Functional Group on Hydrogelation of Amphiphilic Naphthalene-diimide Derivatives and Nonspecific Protein Adsorption.

This manuscript reports the effect of hydrogen-bonding functionality on the supramolecular assembly of naphthalene-diimide (NDI)-derived amphiphilic building blocks in water. All the molecules contain a central NDI chromophore, functionalized with a hydrophilic oligo-oxyethylene (OE) wedge in one arm and a phenyl group on the opposite arm. They differ by a single H-bonding functionality, which links the NDI chromophore and the phenyl moiety. The H-bonding functionalities are amide, thioamide, urea, and urethane in NDI-A, NDI-TA, NDI-U, and NDI-UT, respectively. All of these molecules exhibit π-stacking in water, as evident from their distinct UV/vis absorption spectra when compared to that of the monomeric dye in THF. However, among these four, only NDI-A and NDI-TA show hydrogelation, while the other two precipitate out of the medium. The NDI-A hydrogel also exhibits transient stability and leads to a crystalline precipitate within ∼5 h. Only NDI-TA produces stable transparent hydrogel with the entangled fibrillar morphology that is typical for gelators. Both NDI-A and NDI-TA showed a thermoresponsive property with a lower critical solution temperature of about 41-42 °C. Powder XRD studies show a parallel orientation for NDI-A and an antiparallel orientation for NDI-TA. Computational studies support this experimental observation and indicate that the NDI-A assembly is highly stabilized by strong H-bonding among the amide groups and π-stacking interaction in the parallel orientation. On the other hand, due to weak H-bonding among the thioamide groups, the binding energy of the parallelly oriented NDI-TA was significantly lower and the optimized structure was disordered. Instead, its antiparallel orientation was more stable, with criss-cross aligned H-bonding interactions and π-π interactions between adjacent aromatic rings. The NDI-TA hydrogel with less ordered OE chains on the surface showed prominent adsorption of serum protein BSA. In sharp contrast, NDI-A did not exhibit any notable interaction with BSA, as evident from the ITC studies.

Photoelectrochemical Water Oxidation over Novel Semiconducting Zinc-Based Metal-Thiolate Framework.

Designing an efficient catalyst for a sustainable photoelectrochemical water oxidation reaction is very challenging in the context of renewable energy research. Here, we have introduced a new semiconducting porous zinc-thiolate framework via successful stitching of an "N" donor linker with a triazine-based tristhiolate secondary building unit in the overall architecture. The introduction of both linker and tristhiolate ligand synergistically modifies the architecture by making it a rigid, crystalline, three-dimensional, thermally stable, and porous framework. Our novel zinc-thiolate framework is used as an n-type semiconductor as revealed from the solid-state UV-vis DRS spectroscopic analysis, ac and dc conductivity analysis, and Mott-Schottky plot. This n-type semiconductor-based zinc-thiolate framework is utilized in the photoelectrochemical water oxidation reaction. It displayed a very high efficiency for a visible-light-driven oxygen evolution reaction (OER) in a KOH medium using standard Ag/AgCl as the reference electrode. The superiority of this material was further revealed from the low onset potential (0.822 mV vs RHE), high photocurrent density (0.204 mA cm-2), good stability, and high O2 evolution rate (77 µmol g-1 of oxygen evolution within 2 h), and a good efficiency (ABPE 0.42%, IPCE 29.6% and APCE 34.5%). Furthermore, the porosity in the overall framework seems to be a blessing to the photoelectrochemical performance due to better mass diffusion of the electrolyte. A detailed mechanism for the OER reaction was analyzed through density functional theory analysis suggesting the potential future of this Zn-thiolate framework for achieving a high efficiency in the sustainable water oxidation reaction.

Seeking the Source of Catalytic Efficiency of Lindane Dehydrochlorinase, LinA.

Herein we present the results of an in-depth simulation study of LinA and its two variants. In our analysis, we combined the exploration of protein conformational dynamics with and without bound substrates (hexachlorocyclohexane (HCH) isomers) performed using molecular dynamics simulation followed by the extraction of the most frequently visited conformations and their characteristics with a detailed description of the interactions taking place in the active site between the respective HCH molecule and the first shell residues by using symmetry-adapted perturbation theory (SAPT) calculations. A detailed investigation of the conformational space of LinA substates has been accompanied by description of enzymatic catalytic steps carried out using a hybrid quantum mechanics/molecular mechanics (QM/MM) potential along with the computation of the potential of mean force (PMF) to estimate the free energy barriers for the studied transformations: dehydrochlorination of γ-, (-)-α-, and (+)-α-HCH by LinA-type I and -type II variants. The applied combination of computational techniques allowed us not only to characterize two LinA types but also to point to the most important differences between them and link their features to catalytic efficiency each of them possesses toward the respective ligand. More importantly it has been demonstrated that type I protein is more mobile, its active site has a larger volume, and the dehydrochlorination products are stabilized more strongly than in the case of type II enzyme, due to differences in the residues present in the active sites. Additionally, interaction energy calculations revealed very interesting patterns not predicted before but having the potential to be utilized in any attempts of improving LinA catalytic efficiency. On the basis of all these observations, LinA-type I protein seems to be more preorganized for the dehydrochlorination reaction it catalyzes than the type II variant.

Carbon-bromine bond cleavage - A perspective from bromine and carbon kinetic isotope effects on model debromination reactions.

In this work, we explore the effect of solvent on 13C and 81Br kinetic isotope effects (KIEs) during elimination of bromine substituent from brominated organic compounds promoted by hydroxyl anion. In the present study, we investigate HBr elimination from 2-bromoethylbenzene in three different polar media (water, ethanol, and acetonitrile) as well as bromide ion elimination from 1,2-dibromoethane upon nucleophilic substitution by the hydroxyl ion in aqueous solution using carbon and bromine isotope analysis as mechanistic tools. We reconsider the hypothesis that the magnitude of leaving group halide KIE should visibly depend on the solvent and bond-breaking in a protic solvent should be accompanied by hydrogen bonding which would result in less zero-point energy loss than in an aprotic solvent. Modeling the elimination reaction using the available popular theoretical methods along with different approaches for including environment effects we demonstrate in the presented study no interpretable effect of the solvent on the transition state structure and hence on the theoretically predicted KIEs. The comparison of the magnitudes of carbon and bromine kinetic isotope effects for two different mechanistic pathways (elimination vs substitution) is also discussed.

Modulation of the Conformational Dynamics of Apo-Adenylate Kinase through a π-Cation Interaction.

Large-scale conformational transition from open to closed state of adenylate kinase (ADK) is essential for its catalytic cycle. Apo-ADK undergoes conformational transition in a way that closely resembles an open-to-closed conformational transition. Here, equilibrium simulations, free-energy simulations, and quantum mechanics/molecular mechanics (QM/MM) calculations in combination with several bioinformatics approaches have been used to explore the molecular origin of this conformational transition in apo-ADK. In addition to its conventional open state, Escherichia coli apo-ADK adopts conformations that resemble a closed-like intermediate, the "half-open-half-closed" (HOHC) state, and a π-cation interaction can account for the stability of this HOHC state. Energetics and the electronic properties of this π-cation interaction have been explored using QM/MM calculations. Upon rescinding the π-cation interaction, the conformational landscape of the apo-ADK changes completely. The apo-ADK population is shifted completely toward the open state. This π-cation interaction is highly conserved in bacterial ADK; the cationic guanidinium moiety of a conserved ARG interacts with the delocalized π-electron cloud of either PHE or TYR. Interestingly, this study demonstrates the modulation of a principal protein dynamics by a conserved specific π-cation interaction across different organisms.

Measurement and Prediction of Chlorine Kinetic Isotope Effects in Enzymatic Systems.

Approaches to determine chlorine kinetic isotope effects (Cl-KIEs) on enzymatic dehalogenations are discussed and illustrated by representative examples. Three aspects are considered. First methodology for experimental measurement of Cl-KIEs, with stress being on FAB-IRMS technique developed in our laboratory, is described. Subsequently, we concentrate our discussion on the consequences of reaction complexity in the interpretation of experimental values, a problem especially important in cases of polychlorinated reactants. The most fruitful studies of enzymatic dehalogenations by Cl-KIEs require their theoretical evaluation, hence the computational focus of the second part of this chapter.

Destabilization of Hydrophobic Core of Chicken Villin Headpiece in Guanidinium Chloride Induced Denaturation: Hint of π-Cation Interaction.

Despite their routine use as protein denaturants, the comprehensive understanding of the molecular mechanisms by which urea and guanidinium chloride (GdmCl) disrupts proteins' structure is still lacking. Here, we use steered molecular dynamics simulations along with the umbrella sampling technique to elucidate the mechanism of unfolding of chicken villin headpiece (HP-36) in these two denaturants. We find that while urea denatures protein predominantly by forming hydrogen bonds with the protein backbone, GdmCl commences unfolding by weakening of the hydrophobic interactions present in the core. The potential of mean force calculation indicates the reduction of hydrophobic interactions between two benzene moieties in 6 M GdmCl as compared to 6 M urea. We observe a near parallel orientation between the guanidinium cation and aromatic side chains of the HP-36 suggesting π-cation type stacking interactions which play a crucial role in weakening of the hydrophobic interaction. We use QM/MM optimization calculations to estimate the energetics of this π-cation interaction. Additionally, the consistency of the unfolding paths between high temperature (400 K) unfolding simulations and steered molecular dynamics simulations strengthens the proposed molecular mechanism of unfolding further.

Dehydrochlorination of Hexachlorocyclohexanes Catalyzed by the LinA Dehydrohalogenase. A QM/MM Study.

The elucidation of the catalytic role of LinA dehydrohalogenase in the degradation processes of hexachlorocyclohexane (HCH) isomers is extremely important to further studies on the bioremediation of HCH polluted areas. Herein, QM/MM free energy simulations are employed to provide the details of the dehydrochlorination reaction of two HCH isomers (γ and ß). In particular, the role of the protonation state of one of the catalytic residues-His73-is explored. Based on our calculations, two distinct minimum free energy pathways (concerted and stepwise) were found for γ-HCH and ß-HCH. The choice of the reaction channel for the dehydrochlorination reactions of γ- and ß-HCH was shown to depend on the initial mutual orientations of the reacting species in the active site and the protonation form of His73. The sequential pathway comprises the transfer of the proton (Hδ1) between His73 and Asp25 and subsequently the H1/Cl2 pair elimination from the substrate molecule. Within a concerted mechanism, the dehydrochlorination reaction of γ-/ß-HCH is initiated with neutral His73 and the Hδ1 proton is transferred upon final product formation. We found that the concerted pathway for ß-HCH results in significantly higher free energy of activation than the stepwise route and therefore can be disregarded as not a feasible mechanism. On the other hand, the reaction that occurs with much lower energetic barrier requires a stronger base (i.e., anionic His73) to abstract the proton (H1) from the substrate molecule. The presence of such transient form of His results in higher energy than the respective Michaelis complex and was observed only in the stepwise pathway for both isomers. Furthermore, we have concluded that both pathways (concerted and stepwise) are feasible for the dehydrochlorination reaction of γ-HCH. The activation free energies obtained from the M05-2X/6-31+G(d,p) corrected path coordinate PMF profiles for the dehydrochlorination reactions of the γ-/ß-HCH are in good agreement with the experimental values.