Molecular modeling study for conformational changes of Sirtuin 2 due to substrate and inhibitor binding

Sirtuin is a member of NAD⁺-dependent deacetylase family. The structural details of Sirtuin 2 (SIRT2) complex will be very useful to discover the drug which might have beneficial effects on various diseases like cancer, diabetes, etc. Unfortunately, SIRT2 complex structure is not available yet, hence molecular docking was carried out to dock the substrate (NAD⁺ and acetylated lysine) and inhibitor (sirtinol) in the NAD⁺ binding site. The suitable binding orientation of substrate and inhibitor in the SIRT2 active site was selected and subjected to 5?ns molecular dynamics simulations to adjust the binding orientation of inhibitor and substrate as well as to identify the conformational changes in the active site. The result provides an insight about 3D SIRT2 structural details as well as the importance of F96 in deacetylation function. In addition, our simulations revealed the displacement of F96 upon substrate and inhibitor binding, inducing an extended conformation of loop3 and changing its interactions with the rest of SIRT2. We believe that our study could be helpful to gain a structural insight of SIRT2 and to design the receptor-based inhibitors.     

Evidence of salicylic acid pathway with EDS1 and PAD4 proteins by molecular dynamics simulation for grape improvement

Biotic stress is a major cause of heavy loss in grape productivity. In order to develop biotic stress-resistant grape varieties, the key defense genes along with its pathway have to be deciphered. In angiosperm plants, lipase-like protein phytoalexin deficient 4 (PAD4) is well known to be essential for systemic resistance against biotic stress. PAD4 functions together with its interacting partner protein enhanced disease susceptibility 1 (EDS1) to promote salicylic acid (SA)-dependent and SA-independent defense pathway. Existence and structure of key protein of systemic resistance EDS1 and PAD4 are not known in grapes. Before SA pathway studies are taken in grape, molecular evidence of EDS1: PAD4 complex is to be established. To establish this, EDS1 protein sequence was retrieved from NCBI and homologous PAD4 protein was generated using Arabidopsis thaliana as template and conserved domains were confirmed. In this study, computational methods were used to model EDS1 and PAD4 and simulated the interactions of EDS1 and PAD4. Since no structural details of the proteins were available, homology modeling was employed to construct three-dimensional structures. Further, molecular dynamic simulations were performed to study the dynamic behavior of the EDS1 and PAD4. The modeled proteins were validated and subjected to molecular docking analysis. Molecular evidence of stable complex of EDS1:PAD4 in grape supporting SA defense pathway in response to biotic stress is reported in this study. If SA defense pathway genes are explored, then markers of genes involved can play pivotal role in grape variety development especially against biotic stress leading to higher productivity.     

Asymmetric structure and domain binding interfaces of human tyrosyl‐tRNA synthetase studied by molecular dynamics simulations

Human tyrosyl‐tRNA synthetase (HsTyrRS) is composed of two structural modules: N‐terminal catalytic core and an EMAP II‐like C‐terminal domain. The structures of these modules are known, but no crystal structure of the full‐length HsTyrRS is currently available. An all‐atom model of the full‐length HsTyrRS was developed in this work. The structure, dynamics, and domain binding interfaces of HsTyrRS were investigated by extensive molecular dynamics (MD) simulations. Our data suggest that HsTyrRS in solution consists of a number of compact asymmetric conformations, which differ significantly by their rigidity, internal mobility, orientation of C‐terminal modules, and the strength of interdomain binding. Interfaces of domain binding obtained in MD simulations are in perfect agreement with our previous coarse‐grained hierarchical rotations technique simulations. Formation of the hydrogen bonds between R93 residue of the ELR cytokine motif and the residues A340 and E479 in the C‐module was observed. This observation supports the idea that the lack of cytokine activity in the full‐length HsTyrRS is explained by interactions between N‐modules and C‐modules, which block the ELR motif. Copyright © 2013 John Wiley & Sons, Ltd.     

Multiscale simulation studies of geometrical effects on solution transport through nanopores

Due to intrinsic properties, solid-state nanopores are widely used in nanopore technology. Different geometries (cylindrical (CY), hourglass (HG) and conical (CO)) of artificial nanopores have been fabricated and studied. Each was found to promote different transport abilities experimentally. To explore such pore effects, the combination of finite element (FE) and molecular dynamics (MD) simulations with applied electric filed (150 mV) were performed. The dimension of anion-selective protein pore was used as a nanopore template. Different pore geometries with a narrowest diameter ranging from 1.8 to 1.8 μm were studied here. Firstly, we found that the narrowest regions at a pore orifice in CO and constriction site in HG maximise water velocity and consequently control a water flow rate. Secondly, CY triggers the highest water flux, but low ion selectivity, whilst the funnel-like geometries (HG and CO) enhance the ion selectivity significantly. Both HG and CO show similar degrees of permeant flux and selectivity. The orifice and constriction site in CO and HG are the main player for selectivity and permeation control. Thirdly, the transport properties are tuneable by changing the flow direction in asymmetric CO pore. The tip-to-base flow in CO obviously promotes stronger anion selectivity than the base-to-tip one.     

Molecular dynamics simulation of chitinase I from Thermomyces lanuginosus SSBP to ensure optimal activity

Abstract

The fungal chitinase I obtained from Thermomyces lanuginosus SSBP, a thermophilic deuteromycete, has an optimum growth temperature and pH of 323.15 K and 6.5, respectively. This enzyme plays an important task in the defence mechanism of organisms against chitin-containing parasites by hydrolysing β-1, 4-linkages in chitin. It acts as both anti-fungal and biofouling agents, with some being thermostable and suitable for the industrial applications. Three-dimensional model of chitinase I enzyme was predicted and analysed using various bioinformatics tools. The structure of chitinase I exhibited a well-defined TIM barrel topology with an eight-stranded α/β domain. Structural analysis and folding studies at temperatures ranging from 300 to 375 K using 10 ns molecular dynamics simulations clearly showed the stability of the protein was evenly distributed even at higher temperatures, in accordance with the experimental results. We also carried out a number of 20 ns constant pH molecular dynamics simulations of chitinase I at a pH range 2–6 in a solvent. This work was aimed at establishing the optimum activity and stability profiles of chitinase I. We observed a strong conformational pH dependence of chitinase I and the enzyme retained their characteristic TIM barrel topology at low pH.     

Wing 1 of protein HOP2 is as important as helix 3 in DNA binding by MD simulation

The repair of programmed DNA double-strand breaks through recombination is required for proper association and disjunction of the meiotic homologous chromosomes. Meiosis-specific protein HOP2 plays essential roles in recombination by promoting recombinase activities. The N-terminal domain of HOP2 interacts with DNA through helix 3 (H3) and wing 1 (W1). Mutations in wing 1 (Y65A/K67A/Q68A) slightly weakened the binding but mutations in helices 2 and 3 (Q30A/K44A/K49A) nearly abolished the binding. To better understand such differential effects at atomic level, molecular dynamics simulations were employed. Despite losing some hydrogen bonds, the W1-mutant DNA complex was rescued by stronger hydrophobic interactions. For the wild type and W1-mutant, the protein was found to slide along the DNA grooves as the DNA rolls along its double-helix axis. This motion could be functionally important to facilitate the precise positioning of the single-stranded DNA with the homologous double-stranded DNA. The sliding motion was reduced in the W1-mutant. The H-mutant nearly lost all intermolecular interactions. Moreover, an additional mutation in wing 1 (Y65A/K67A/Q68A/K69A) also caused complete complex dissociation. Therefore, both wing 1 and helix 3 make important contribution to the DNA binding, which could be important to the strand invasion function of HOP2 homodimer and HOP2-MND1 heterodimer. Similar to cocking a medieval crossbow with the archer’s foot placed in the stirrup, wing 1 may push the minor groove to cause distortion while helix 3 grabs the major groove.     

Molecular dynamics simulations of human and dog gastric lipases: insights into domain movements

Mammalian gastric lipases are stable and active under acidic conditions and also in the duodenal lumen. There has been considerable interest in acid stable lipases owing to their potential application in the treatment of pancreatic exocrine insufficiency. In order to gain insights into the domain movements of these enzymes, molecular dynamics simulations of human gastric lipase was performed at an acidic pH and under neutral conditions. For comparative studies, simulation of dog gastric lipase was also performed at an acidic pH. Analyses show, that in addition to the lid region, there is another region of high mobility in these lipases. The potential role of this novel region is discussed.     

Structural and Dynamical Studies of Humanin in Water and TFE/Water Mixture: A Molecular Dynamics Simulation

Summary

Proline-rich peptides are known to adopt preferentially the extended polyproline II (PPII) helical conformation, which is involved in several protein-protein recognition events. By resorting to molecular modelling techniques, we wished to investigate the extent to which PPII helices could be used for the formation of isohelical peptide-DNA complexes leading to the selective recognition of the major groove of B-DNA. For that purpose, we have grafted to a cationic intercalator, 9-amino-acridine, an oligopeptide having the sequence: Pro-Arg-Pro-Pro-Arg-Pro-Pro-Arg-Pro-Pro-Asp-Pro-Pro. Each residue in the sequence was set in the D configuration, to prevent enzymatic hydrolysis, and each Arg residue was designed to target O6/N7 of a guanine base following the intercalation site. The Asp residue was designed to target a cytosine base, whilst simultaneously forming a bidentate complex with the Arg three residues upstream. Energy-minimization, using the JUMNA procedure, led to the following conclusions: 1) major groove binding is favoured over minor groove or exclusive binding to the phosphates by large energy differences, of over 50 and 90 kcal/mole, respectively; 2) the two best bound sequences are those having three successive guanine bases on the same DNA strand, immediately adjacent to the intercalation site. Sequence d(CGGGC G), encountered in the Primer Binding Site of the HIV retrovirus, thus ranks amongst the best-bound sequences; 3) replacement of an individual guanine amongst the three ones upstream of the intercalation site, by an adenine base, weakens by > 6 kcal/mole the binding energetics; 4) the conformational rigidity of the DNA-bound PPII helix should enable for a modulation of the base sequence selectivity, by appropriate replacements of the Arg and Asp residues. Thus sequence CGGCAAG, also encountered in the HIV genome, could be targeted by an oligopeptide having the sequence Pro-Arg-Pro-Pro-Asp-Pro-Pro- Asn-Pro-Pro-Asn-Pro-Pro-Arg-Ala.