In SilicoStructural and Functional Analysis of Fragments of the Ankyrin Repeat Protein p18INK4c

Molecular cloning and characterization of NPR1 gene from Arachis hypogaea

The NPR1 gene was an important regulator for a plant disease resistance. The cDNA of NPR1 gene was cloned from peanut cultivar Ri Hua 1 by rapid amplification of cDNA ends-polymerase chain reaction (RACE-PCR). The full length cDNA of Arachis hypogaea NPR1 consisted of 2,078 base pairs with a 1,446 bp open-reading frame encoding 481 amino acids. The predicted NPR1 contained the highly conserved functional domains (BTB/POZ domain from M1 to D116), protein-protein interaction domains (three ankyrin repeats from K158 to L186; N187 to L217 and R221 to D250) and one NPR1-like domain (C262 to S469). The DNA sequence of the NPR1 gene was 2,332 or 2,223 bp. Both two sequences contained three introns and four exons. The NPR1 transcripts were expressed mainly in roots and leaves, while fewer signals were detected in the stems. Amount of the NPR1 transcript was significantly increased 1 h after salicylic acid challenge and was eventually 5.3 times greater than that in the control group. Both the DNA sequence and the coding sequence were obtained from eight cultivars and nine wild species of Arachis. Maximum likelihood analyses of d N/d S ratios for 25 sequences from different species showed that different selection pressures may have acted on different branches.

Potential human cholesterol esterase inhibitor design: benefits from the molecular dynamics simulations and pharmacophore modeling studies

Human pancreatic cholesterol esterase (hCEase) is one of the lipases found to involve in the digestion of large and broad spectrum of substrates including triglycerides, phospholipids, cholesteryl esters, etc. The presence of bile salts is found to be very important for the activation of hCEase. Molecular dynamic simulations were performed for the apoform and bile salt complexed form of hCEase using the co-ordinates of two bile salts from bovine CEase. The stability of the systems throughout the simulation time was checked and two representative structures from the highly populated regions were selected using cluster analysis. These two representative structures were used in pharmacophore model generation. The generated pharmacophore models were validated and used in database screening. The screened hits were refined for their drug-like properties based on Lipinski's rule of five and ADMET properties. The drug-like compounds were further refined by molecular docking simulation using GOLD program based on the GOLD fitness score, mode of binding, and molecular interactions with the active site amino acids. Finally, three hits of novel scaffolds were selected as potential leads to be used in novel and potent hCEase inhibitor design. The stability of binding modes and molecular interactions of these final hits were re-assured by molecular dynamics simulations.

The pectin lyases in Arabidopsis thaliana: evolution, selection and expression profiles

Pectin lyases are a group of enzymes that are thought to contribute to many biological processes, such as the degradation of pectin. However, until this study, no comprehensive study incorporating phylogeny, chromosomal location, gene duplication, gene organization, functional divergence, adaptive evolution, expression profiling and functional networks has been reported for Arabidopsis. Sixty-seven pectin lyase genes have been identified, and most of them possess signal sequences targeting the secretory pathway. Phylogenetic analyses identified five gene groups with considerable conservation among groups. Pectin lyase genes were non-randomly distributed across chromosomes and clustering was evident. Functional divergence and adaptive evolution analyses suggested that purifying selection was the main force driving pectin lyase evolution, although some critical sites responsible for functional divergence might be the consequence of positive selection. A stigma- and receptacle-specific expression promoter was identified, and it had increased expression in response to wounding. Two hundred and eighty-eight interactions were identified by functional network analyses, and most of these were involved in cellular metabolism, cellular transport and localization, and stimulus responses. This investigation contributes to an improved understanding of the complexity of the Arabidopsis pectin lyase gene family.

Insights into the thermal stabilization and conformational transitions of DNA by hyperthermophile protein Sso7d: molecular dynamics simulations and MM-PBSA analysis

In the assembly of DNA-protein complex, the DNA kinking plays an important role in nucleoprotein structures and gene regulation. Molecular dynamics (MD) simulations were performed on specific protein-DNA complexes in this study to investigate the stability and structural transitions of DNA depending on temperature. Furthermore, we introduced the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) approach to analyze the interactions between DNA and protein in hyperthermophile. Focused on two specific Sso7d-DNA complexes (PDB codes: 1BNZ and 1BF4), we performed MD simulations at four temperatures (300, 360, 420, and 480 K) and MM-PBSA at 300 and 360 K to illustrate detailed information on the changes of DNA. Our results show that Sso7d stabilizes DNA duplex over a certain temperature range and DNA molecules undergo B-like to A-like form transitions in the binary complex with the temperature increasing, which are consistent with the experimental data. Our work will contribute to a better understanding of protein-DNA interaction.

The effect of fucoidan on tyrosinase: computational molecular dynamics integrating inhibition kinetics

Fucoidan is a complex sulfated polysaccharide extracted from brown seaweed and has a wide variety of biological activities. In this study, we investigated the inhibitory effect of fucoidan on tyrosinase via a combination of inhibition kinetics and computational simulations. Fucoidan reversibly inhibited tyrosinase in a mixed-type manner. Time-interval kinetics showed that the inhibition was processed as first order with biphasic processes. For further insight, we simulated dockings with various sizes of molecular models (monomer to decamer) of fucoidan and showed that the best binding energy change results were obtained from the pentamer (-1.89 kcal/mol) and the hexamer (-1.97 kcal/mol) models of AutoDock Vina. The molecular dynamics simulation confirmed the binding mechanisms between tyrosinase and fucoidan and suggested that fucoidan mostly interacts with several residues including copper ions located in the active site. Our study suggests that fucoidan might be a potential natural antipigment agent.

Mechanism of antisense oligonucleotide interaction with natural RNAs

Oligonucleotides find several numbers of applications: as diagnostic probes, RT and PCR primers and antisense agents due to their ability of forming specific interactions with complementary nucleotide sequences within nucleic acids. These interactions are strongly affected by accessibility of the target sequence in the RNA structure. In the present work the mechanism of invasion of RNA structure by oligonucleotide was investigated using a model system: yeast tRNA(Phe) and oligonucleotides complementary to the 3'-part of this molecule. Kinetics of interaction of oligonucleotides with in vitro transcript of yeast tRNAPhe was studied using stopped-flow technique with fluorescence quenching detection, 5'-DABCYL labeled oligonucleotide was hybridized with 3'-fluorescein labeled tRNA(Phe). The results evidence for a four-step invasion process of the oligonucleotide-RNA complex formation. The process is initiated by formation of transition complexes with nucleotides in the T-loop and ACCA sequence. This complex formation is followed by RNA unfolding and formation of an extended heteroduplex with the oligonucleotide via strand displacement process. Computer modeling of oligonucleotide-tRNA(Phe) interaction revealed potential factors that could favor transition complexes formation and confirmed the proposed mechanism, showing the oligonucleotide to be a molecular "wedge". Our data evidence that oligonucleotide invasion into structured RNA is initiated by loop-single strand interactions, similar to the initial step of the antisense RNA-RNA interactions. The obtained results can be used for choosing efficient oligonucleotide probes.

Geometric parameters defining the structure of proteins--relation to early-stage folding step

Two geometrical parameters describing the structure of a polypeptide: V-dihedral angle between two sequential peptide bond planes and R-radius of curvature are used for structural classification of polypeptide structure in proteins. The relation between these two parameters was the basis for the definition of the conformational sub-space for early-stage structural forms. The cluster analysis of V and lnR, applied to the selected proteins of well-defined secondary structure (according to DSSP classification) and to proteins without any introductory classified analysis, revealed that several of the discriminated groups of proteins agree with the assumed model of early-stage conformational sub-space. This analysis shows that protein structures may be represented in VR space instead of Phi, Psi angles space, thus lowering the conformational space dimensionality. The VR model allows classification of traditional secondary structure elements as well as different Random Coil motifs, which broadens the range of recognized structural categories (compared to standard secondary structure elements).

In search of allosteric modulators of a7-nAChR by solvent density guided virtual screening

Nicotinic acetylcholine receptors (nAChR) are pentameric ligand gated ion channels whose activity can be modulated by endogenous neurotransmitters as well as by synthetic ligands that bind the same or distinct sites from the natural ligand. The subtype of α7 nAChR has been considered as a potenial therapeutic target for Alzheimer's disease, schizophrenia and other neurological and psychiatric disorders. Here we have developed a homology model of α7 nAChR based on two high resolution crystal structures with Brookhaven Protein Data Bank (PDB) codes 2QC1 and 2WN9 for threading on one monomer and then for building a pentamer, respectively. A number of small molecule binding sites are identified using Pocket Finder (J. An, M. Tortov, and R. Abagyan, Molecular & Cellular Proteomics, 4.6, 752-761 (2005)) of Internal Coordinate Mechanics (ICM). Remarkably, these computer-identified sites match perfectly with ordered solvent densities found in the high-resolution crystal structure of α1 nAChR, suggesting that the surface cavities in the α7 nAChR model are likely binding sites of small molecules. A high throughput virtual screening by flexible ligand docking of 5008 small molecule compounds was performed at three potential allosteric modulator (AM) binding sites of α7 nAChR using Molsoft ICM software (R. Abagyan, M. Tortov and D. Kuznetsov, J Comput Chem 15, 488-506, (1994)). Some experimentally verified allosteric modulators of α7 like CCMI comp-6, LY 7082101, 5-HI, TQS, PNU-120596, genistein, and NS-1738 ranked among top 100 compounds, while the rest of the compounds in the list could guide further search for new allosteric modulators.

Using computer-aided drug design and medicinal chemistry strategies in the fight against diabetes

The aim of this work is to present a simple, practical and efficient protocol for drug design, in particular Diabetes, which includes selection of the illness, good choice of a target as well as a bioactive ligand and then usage of various computer aided drug design and medicinal chemistry tools to design novel potential drug candidates in different diseases. We have selected the validated target dipeptidyl peptidase IV (DPP-IV), whose inhibition contributes to reduce glucose levels in type 2 diabetes patients. The most active inhibitor with complex X-ray structure reported was initially extracted from the BindingDB database. By using molecular modification strategies widely used in medicinal chemistry, besides current state-of-the-art tools in drug design (including flexible docking, virtual screening, molecular interaction fields, molecular dynamics, ADME and toxicity predictions), we have proposed 4 novel potential DPP-IV inhibitors with drug properties for Diabetes control, which have been supported and validated by all the computational tools used herewith.

Screening from the world's largest TCM database against H1N1 virus

The swine influenza virus (H1N1) 2009 pandemic highlights the importance of having effective anti-viral strategies. Recently, oseltamivir (Tamiflu) resistant influenza viruses are identified; which further emphasizes the urgency in developing new antiviral agents. In influenza virus replication cycle, viral surface glycoprotein, hemagglutinin, is responsible for viral entry into host cells. Hence, a potentially effective antiviral strategy is to inhibit viral entry mechanism. To develop novel antiviral agent that inhibits viral entry, we analyzed 20,000 traditional Chinese medicine (TCM) ingredients in hemagglutinin subtype H1 sialic acid binding site found on H1N1 virus. We then performed molecular dynamics simulations to investigate receptor-ligand interaction of the candidates obtained from docking. Here, we report three TCM derivatives that have high binding affinities to H1 sialic acid binding site residues based on structure-based calculations. The top three derivatives, xylopine_2, rosmaricine_14 and rosmaricine_15, all have an amine group that interact with Glu83 and a pyridinium group that interact with Asp103. Molecular dynamics simulations show that these derivatives form strong hydrogen bonding with Glu83 but interact transiently with Asp103. We therefore suggest that an enhanced hemagglutinin inhibitor, based on our scaffold, should be designed to bind both Glu83 and Asp103 with high affinity.

Thermal stability and unfolding pathways of Sso7d and its mutant F31A: insight from molecular dynamics simulation

The thermo-stability and unfolding behaviors of a small hyperthermophilic protein Sso7d as well as its single-point mutation F31A are studied by molecular dynamics simulation at temperatures of 300 K, 371 K and 500 K. Simulations at 300 K show that the F31A mutant displays a much larger flexibility than the wild type, which implies that the mutation obviously decreases the protein's stability. In the simulations at 371 K, although larger fluctuations were observed, both of these two maintain their stable conformations. High temperature simulations at 500 K suggest that the unfolding of these two proteins evolves along different pathways. For the wild-type protein, the C-terminal alpha-helix is melted at the early unfolding stage, whereas it is destroyed much later in the unfolding process of the F31A mutant. The results also show that the mutant unfolds much faster than its parent protein. The deeply buried aromatic cluster in the F31A mutant dissociates quickly relative to the wild-type protein at high temperature. Besides, it is found that the triple-stranded antiparallel β-sheet in the wild-type protein plays an important role in maintaining the stability of the entire structure.

Microsecond scale replica exchange molecular dynamic simulation of villin headpiece: an insight into the folding landscape

Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 µs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 Å was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins.

A novel OmpY porin from Yersinia pseudotuberculosis: structure, channel-forming activity and trimer thermal stability

A novel OmpY porin was predicted based on the Yersinia pseudotuberculosis genome analysis. Whereas it has the different genomic annotation such as "outer membrane protein N" (ABS46310.1) in str. IP 31758 or "outer membrane protein C2, porin" (YP_070481.1) in str. IP32953, it might be warranted to rename the OmpN/OmpC2 to OmpY, "outer membrane protein Y", where letter "Y" pertained to Yersinia. Both phylogenetic analysis and genomic localization clearly support that the OmpY porin belongs to a new group of general bacterial porins. The recombinant OmpY protein with its signal sequence was overexpressed in porin-deficient Escherichia coli strain. The mature rOmpY was shown to insert into outer membrane as a trimer. The OmpY porin, isolated from the outer membrane, was studied employing spectroscopic, electrophoretic and bilayer lipid membranes techniques. The far UV CD spectrum of rOmpY was essentially identical to that of Y. pseudotuberculosis OmpF. The near UV CD spectrum of rOmpY was weaker and smoother than that of OmpF. The rOmpY single-channel conductance was 180 ± 20 pS in 0.1 M NaCl and was lower than that of the OmpF porin. As was shown by electrophoretic and bilayer lipid membrane experiments, the rOmpY trimers were less thermostable than the OmpF trimers. The porins differed in the trimer-monomer transition temperature by about 20°C. The three-dimensional structural models of the Y. pseudotuberculosis OmpY and OmpF trimers were generated and the intra- and intermonomeric interactions stabilizing the porins were investigated. The difference in the thermal stability of OmpY and OmpF trimers was established to correlate with the difference in intermonomeric polar contacts.

Explicit Water Models Affect the Specific Solvation and Dynamics of Unfolded Peptides While the Conformational Behavior and Flexibility of Folded Peptides Remain Intact

Conventional molecular dynamics simulations on 50 ns to 1 μs time scales were used to study the effects of explicit solvent models on the conformational behavior and solvation of two oligopeptide solutes: α-helical EK-peptide (14 amino acids) and a β-hairpin chignolin (10 amino acids). The widely used AMBER force fields (ff99, ff99SB, and ff03) were combined with four of the most commonly used explicit solvent models (TIP3P, TIP4P, TIP5P, and SPC/E). Significant differences in the specific solvation of chignolin among the studied water models were identified. Chignolin was highly solvated in TIP5P, whereas reduced specific solvation was found in the TIP4P, SPC/E, and TIP3P models for kinetic, thermodynamic, and both kinetic and thermodynamic reasons, respectively. The differences in specific solvation did not influence the dynamics of structured parts of the folded peptide. However, substantial differences between TIP5P and the other models were observed in the dynamics of unfolded chignolin, stability of salt bridges, and specific solvation of the backbone carbonyls of EK-peptide. Thus, we conclude that the choice of water model may affect the dynamics of flexible parts of proteins that are solvent-exposed. On the other hand, all water models should perform similarly for well-structured folded protein regions. The merits of the TIP3P model include its high and overestimated mobility, which accelerates simulation processes and thus effectively increases sampling.