GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit

Predicting the Thermodynamics and Kinetics of Helix Formation in a Cyclic Peptide Model

The peptide Ac-(cyclo-2,6)-R[KAAAD]-NH2 (cyc-RKAAAD) is a short cyclic peptide known to adopt a remarkably stable single turn α-helix in water. Due to its simplicity and the availability of thermodynamic and kinetic experimental data, cyc-RKAAAD poses as an ideal model for evaluating the aptness of current molecular dynamics (MD) simulation methodologies to accurately sample conformations that reproduce experimentally observed properties. In this work, we extensively sample the conformational space of cyc-RKAAAD using microsecond-timescale MD simulations. We characterize the peptide conformational preferences in terms of secondary structure propensities and, using Cartesian-coordinate principal component analysis (cPCA), construct its free energy landscape, thus obtaining a detailed weighted discrimination between the helical and nonhelical subensembles. The cPCA state discrimination, together with a Markov model built from it, allowed us to estimate the free energy of unfolding (-0.57 kJ/mol) and the relaxation time (∼0.435 μs) at 298.15 K, which are in excellent agreement with the experimentally reported values (-0.22 kJ/mol and 0.42 μs, Serrano, A. L.; Tucker, M. J.; Gai, F. J. Phys. Chem. B, 2011, 115, 7472-7478.). Additionally, we present simulations conducted using two enhanced sampling methods: replica-exchange molecular dynamics (REMD) and bias-exchange metadynamics (BE-MetaD). We compare the free energy landscape obtained by these two methods with the results from MD simulations and discuss the sampling and computational gains achieved. Overall, the results obtained attest to the suitability of modern simulation methods to explore the conformational behavior of peptide systems with a high level of realism.

CO2 migration pathways in oxalate decarboxylase and clues about its active site

Oxalate decarboxylase catalyzes the decarboxylation of oxalate to formate and CO2 in the presence of molecular oxygen. This enzyme has two domains, each containing a Mn(II) ion coordinated with three histidine residues. The specific domain in which the decarboxylation process takes place is still a matter of investigation. Herein, the transport of the product, i.e., CO2, from the reaction center to the surface of the enzyme is studied using atomistic molecular dynamics simulations. The specific pathway for the migration of the molecule as well as its microscopic interactions with the amino acid residues lining the path is delineated. Further, the transport of CO2 is shown to occur in a facile manner from only domain I and not from domain II, indicating that the former is likely to be the active site of the enzyme.

Biochemical properties of nematode O-acetylserine(thiol)lyase paralogs imply their distinct roles in hydrogen sulfide homeostasis

O-Acetylserine(thiol)lyases (OAS-TLs) play a pivotal role in a sulfur assimilation pathway incorporating sulfide into amino acids in microorganisms and plants, however, these enzymes have not been found in the animal kingdom. Interestingly, the genome of the roundworm Caenorhabditis elegans contains three expressed genes predicted to encode OAS-TL orthologs (cysl-1-cysl-3), and a related pseudogene (cysl-4); these genes play different roles in resistance to hypoxia, hydrogen sulfide and cyanide. To get an insight into the underlying molecular mechanisms we purified the three recombinant worm OAS-TL proteins, and we determined their enzymatic activities, substrate binding affinities, quaternary structures and the conformations of their active site shapes. We show that the nematode OAS-TL orthologs can bind O-acetylserine and catalyze the canonical reaction although this ligand may more likely serve as a competitive inhibitor to natural substrates instead of being a substrate for sulfur assimilation. In addition, we propose that S-sulfocysteine may be a novel endogenous substrate for these proteins. However, we observed that the three OAS-TL proteins are conformationally different and exhibit distinct substrate specificity. Based on the available evidences we propose the following model: CYSL-1 interacts with EGL-9 and activates HIF-1 that upregulates expression of genes detoxifying sulfide and cyanide, the CYSL-2 acts as a cyanoalanine synthase in the cyanide detoxification pathway and simultaneously produces hydrogen sulfide, while the role of CYSL-3 remains unclear although it exhibits sulfhydrylase activity in vitro. All these data indicate that C. elegans OAS-TL paralogs have distinct cellular functions and may play different roles in maintaining hydrogen sulfide homeostasis.

Structural insights into RIP3-mediated necroptotic signaling

RIP3 is an essential upstream kinase in necroptosis. The pseudokinase MLKL functions as a substrate of RIP3 to mediate downstream signaling. The molecular mechanism by which RIP3 recognizes and phosphorylates MLKL remains unknown. Here, we report the crystal structures of the mouse RIP3 kinase domain, the MLKL kinase-like domain, and a binary complex between the two. Both RIP3 and MLKL adopt the canonical kinase fold. Free RIP3 exists in an active conformation, whereas MLKL-bound RIP3 is stabilized by AMP-PNP to adopt an inactive conformation. The formation of the RIP3-MLKL complex, involving their respective N- and C-lobes, is accompanied by pronounced conformational changes of the αC helix and activation loop in RIP3 and the corresponding structural elements in MLKL. RIP3-mediated MLKL phosphorylation, though important for downstream signaling, is dispensable for stable complex formation between RIP3 and MLKL. Our study serves as a framework for mechanistic understanding of RIP3-mediated necroptotic signaling.

Deciphering the synergism of endogenous glycoside hydrolase families 1 and 9 from Coptotermes gestroi

Termites can degrade up to 90% of the lignocellulose they ingest using a repertoire of endogenous and symbiotic degrading enzymes. Termites have been shown to secrete two main glycoside hydrolases, which are GH1 (EC 3.2.1.21) and GH9 (EC 3.2.1.4) members. However, the molecular mechanism for lignocellulose degradation by these enzymes remains poorly understood. The present study was conducted to understand the synergistic relationship between GH9 (CgEG1) and GH1 (CgBG1) from Coptotermes gestroi, which is considered the major urban pest of São Paulo State in Brazil. The goal of this work was to decipher the mode of operation of CgEG1 and CgBG1 through a comprehensive biochemical analysis and molecular docking studies. There was outstanding degree of synergy in degrading glucose polymers for the production of glucose as a result of the endo-β-1,4-glucosidase and exo-β-1,4-glucosidase degradation capability of CgEG1 in concert with the high catalytic performance of CgBG1, which rapidly converts the oligomers into glucose. Our data not only provide an increased comprehension regarding the synergistic mechanism of these two enzymes for cellulose saccharification but also give insight about the role of these two enzymes in termite biology, which can provide the foundation for the development of a number of important applied research topics, such as the control of termites as pests as well as the development of technologies for lignocellulose-to-bioproduct applications.

Structural intermediates and folding events in the early assembly of the ribosomal small subunit

Using all-atom explicit solvent molecular dynamics (MD) simulations, we investigated the early structural intermediates of the 5' domain of the 16S rRNA in Escherichia coli upon the removal of the primary binding r-proteins S4, S17, and S20 and the secondary binding r-protein S16. Removal of each r-protein corresponded to the disappearance of subdomains with correlated dynamics. Correlation-based network analysis of the MD trajectories of the naked rRNA showed that the different subdomains are connected via multiple pathways with high betweenness. These pathways cross at the internal loop of helix 17 (h17) in the five-way junction (5WJ). The structure of the internal loop is disrupted by the binding of S17 and rescued by the addition of S16, suggesting an important function of the secondary binding protein in biasing the rRNA folding landscape toward the native basin. Using structure-based Gō simulations, we investigated the folding barriers of the lower four-way junction (4WJ) with h6, which is the primary binding site of S20 and the first to be transcribed. The folding of the 4WJ is consistent with the protection patterns observed in hydroxyl radical footprinting. Results from the all-atom simulations show that the fluctuations in the 5WJ are independent of the fluctuations in the 4WJ, suggesting that the subdomains fold independently and are stabilized by primary r-proteins.

A simple method to measure protein side-chain mobility using NMR chemical shifts

Protein side-chain motions are involved in many important biological processes including enzymatic catalysis, allosteric regulation, and the mediation of protein-protein, protein-DNA, protein-RNA, and protein-cofactor interactions. NMR spectroscopy has long been used to provide insights into the motions of side-chain groups. Currently, the method of choice for studying side-chain dynamics by NMR is the measurement of methyl (2)H autorelaxation. Methyl (2)H autorelaxation exhibits simple relaxation mechanisms and can be straightforwardly converted to meaningful dynamic parameters. However, methyl groups can only be found in 6 of 19 side-chain bearing amino acids. Consequently, only a sparse assessment of protein side-chain dynamics is possible. Therefore, there is a significant interest in developing novel methods of studying side-chain motions that can be applied to all types of side-chains. Here, we show how side-chain chemical shifts can be used to determine the magnitude of fast side-chain motions in proteins. The chemical shift method is applicable to all side-chain bearing residues and does not require any additional measurements beyond standard NMR experiments for backbone and side-chain assignments.

Molecular dynamics simulations on the Tre1 G protein-coupled receptor: exploring the role of the arginine of the NRY motif in Tre1 structure

BACKGROUND

The arginine of the D/E/NRY motif in Rhodopsin family G protein-coupled receptors (GPCRs) is conserved in 96% of these proteins. In some GPCRs, this arginine in transmembrane 3 can form a salt bridge with an aspartic acid or glutamic acid in transmembrane 6. The Drosophila melanogaster GPCR Trapped in endoderm-1 (Tre1) is required for normal primordial germ cell migration. In a mutant form of the protein, Tre1sctt, eight amino acids RYILIACH are missing, resulting in a severe disruption of primordial germ cell development. The impact of the loss of these amino acids on Tre1 structure is unknown. Since the missing amino acids in Tre1sctt include the arginine that is part of the D/E/NRY motif in Tre1, molecular dynamics simulations were performed to explore the hypothesis that these amino acids are involved in salt bridge formation and help maintain Tre1 structure.

RESULTS

Structural predictions of wild type Tre1 (Tre1+) and Tre1sctt were subjected to over 250 ns of molecular dynamics simulations. The ability of the model systems to form a salt bridge between the arginine of the D/E/NRY motif and an aspartic acid residue in transmembrane 6 was analyzed. The results indicate that a stable salt bridge can form in the Tre1+ systems and a weak salt bridge or no salt bridge, using an alternative arginine, is likely in the Tre1sctt systems.

CONCLUSIONS

The weak salt bridge or lack of a salt bridge in the Tre1sctt systems could be one possible explanation for the disrupted function of Tre1sctt in primordial germ cell migration. These results provide a framework for studying the importance of the arginine of the D/E/NRY motif in the structure and function of other GPCRs that are involved in cell migration, such as CXCR4 in the mouse, zebrafish, and chicken.

Molecular dynamics simulation of SDS and CTAB micellization and prediction of partition equilibria with COSMOmic

Molecular dynamics (MD) simulations of the self-assembly of different ionic surfactants have been performed in order to obtain representative micellar structures. Subsequently, these structures were used to predict the partition behavior of various solutes in these micelles with COSMOmic, an extension of COSMO-RS. This paper includes multiple self-assembled micelles of SDS (sodium dodecyl sulfate, anionic surfactant) and CTAB (cetyltrimethylammoniumbromide, cationic surfactant) at different concentrations. Micellar size, density profiles, and shape (eccentricity) have been investigated. However, the size strongly depends on the functional definition of a micelle. For this reason, we present a method based on the free monomer concentration in aqueous solution as an optimization criterion for the micelle definition. The combination of MD with COSMOmic has the benefit of combining detailed atomistic information from MD with fast calculations of COSMOmic. For the first time the influence of micelle structure on pratition equilibria, predicted with COSMOmic, were investigated. In case of SDS more than 4600 and for CTAB more than 800 single micelles have been studied. The predictions of the partition coefficients with COSMOmic are in good agreement with experimental data. Additionally, the most favorable locations of selected molecules in the micelles as well as probable energy barriers are determined even for complex solutes such as toluene, propanolol, ephedrine, acetone, phenol, lidocaine, syringic acid, coumarin, isovanillin, ferulic acid, and vanillic acid. This method can therefore be applied as a potential screening tool for solutes (e.g., drugs) to find the optimal solute-surfactant combination.

Structural and dynamic properties of thymopoietin mimetics

We propose a hypothesis that the T-cell receptor is a possible target of thymic hormones. We modelled the conformational dynamics of thymopentin and its structural variants in solution, as well as the interactions of these short peptides with the proposed molecular target. Thymopentin is a five-amino-acid fragment of the thymic hormone thymopoietin (residues 32 to 36) that reproduces the immunomodulatory activity of the complete hormone. Using molecular dynamics and flexible docking methods, we demonstrated high-affinity binding of thymopentin and its prospective mimetics with the T-cell receptor. The calculated biological activity spectra of thymopentin and its two promising modifications can be used in immunomodulatory activity screenings with live systems.

Role of N-terminal extension of Bacillus stearothermophilus RNase H2 and C-terminal extension of Thermotoga maritima RNase H2

Bacillus stearothermophilus RNase H2 (BstRNH2) and Thermotoga maritima RNase H2 (TmaRNH2) have N-terminal and C-terminal extensions, respectively, as compared with Aquifex aeolicus RNase H2 (AaeRNH2). To analyze the role of these extensions, BstRNH2 and TmaRNH2 without these extensions were constructed, and their biochemical properties were compared with those of their intact partners and AaeRNH2. The far-UV CD spectra of all proteins were similar, suggesting that the protein structure is not significantly altered by removal of these extensions. However, both the junction ribonuclease and RNase H activities of BstRNH2 and TmaRNH2, as well as their substrate-binding affinities, were considerably decreased by removal of these extensions. The stability of BstRNH2 and TmaRNH2 was also decreased by removal of these extensions. The activity, substrate binding affinity and stability of TmaRNH2 without the C-terminal 46 residues were partly restored by the attachment of the N-terminal extension of BstRNH2. These results suggest that the N-terminal extension of BstRNH2 functions as a substrate-binding domain and stabilizes the RNase H domain. Because the C-terminal extension of TmaRNH2 assumes a helix hairpin structure and does not make direct contact with the substrate, this extension is probably required to make the conformation of the substrate-binding site functional. AaeRNH2 showed comparable junction ribonuclease activity to those of BstRNH2 and TmaRNH2, and was more stable than these proteins, indicating that bacterial RNases H2 do not always require an N-terminal or C-terminal extension to increase activity, substrate-binding affinity, and/or stability.

Molecular dynamics of the full-length p53 monomer

The p53 protein is frequently mutated in a very large proportion of human tumors, where it seems to acquire gain-of-function activity that facilitates tumor onset and progression. A possible mechanism is the ability of mutant p53 proteins to physically interact with other proteins, including members of the same family, namely p63 and p73, inactivating their function. Assuming that this interaction might occurs at the level of the monomer, to investigate the molecular basis for this interaction, here, we sample the structural flexibility of the wild-type p53 monomeric protein. The results show a strong stability up to 850 ns in the DNA binding domain, with major flexibility in the N-terminal transactivations domains (TAD1 and TAD2) as well as in the C-terminal region (tetramerization domain). Several stable hydrogen bonds have been detected between N-terminal or C-terminal and DNA binding domain, and also between N-terminal and C-terminal. Essential dynamics analysis highlights strongly correlated movements involving TAD1 and the proline-rich region in the N-terminal domain, the tetramerization region in the C-terminal domain; Lys120 in the DNA binding region. The herein presented model is a starting point for further investigation of the whole protein tetramer as well as of its mutants.

Assessing the effect of loop mutations in the folding space of β2-microglobulin with molecular dynamics simulations

We use molecular dynamics simulations of a full atomistic Gō model to explore the impact of selected DE-loop mutations (D59P and W60C) on the folding space of protein human β₂-microglobulin (Hβ₂m), the causing agent of dialysis-related amyloidosis, a conformational disorder characterized by the deposition of insoluble amyloid fibrils in the osteoarticular system. Our simulations replicate the effect of mutations on the thermal stability that is observed in experiments in vitro. Furthermore, they predict the population of a partially folded state, with 60% of native internal free energy, which is akin to a molten globule. In the intermediate state, the solvent accessible surface area increases up to 40 times relative to the native state in 38% of the hydrophobic core residues, indicating that the identified species has aggregation potential. The intermediate state preserves the disulfide bond established between residue Cys25 and residue Cys80, which helps maintain the integrity of the core region, and is characterized by having two unstructured termini. The movements of the termini dominate the essential modes of the intermediate state, and exhibit the largest displacements in the D59P mutant, which is the most aggregation prone variant. PROPKA predictions of pK_a suggest that the population of the intermediate state may be enhanced at acidic pH explaining the larger amyloidogenic potential observed in vitro at low pH for the WT protein and mutant forms.

Convergence of QM/MM and Cluster Models for the Spectroscopic Properties of the Oxygen-Evolving Complex in Photosystem II

The latest crystal structure of photosystem II at 1.9 Å resolution, which resolves the topology of the Mn4CaO5 oxygen evolving complex (OEC) at atomistic detail, enables a better correlation between structural features and spectroscopic properties than ever before. Building on the refined crystallographic model of the OEC and the protein, we present combined quantum mechanical/molecular mechanical (QM/MM) studies of the spectroscopic properties of the natural catalyst embedded in the protein matrix. Focusing on the S2 state of the catalytic cycle, we examine the convergence of not only structural parameters but also of the intracluster magnetic interactions in terms of exchange coupling constants and of experimentally relevant (55)Mn, (17)O, and (14)N hyperfine coupling constants with respect to QM/MM partitioning using five QM regions of increasing size. This enables us to assess the performance of the method and to probe second sphere effects by identifying amino acid residues that principally affect the spectroscopic properties of the OEC. Comparison between QM-only and QM/MM treatments reveals that whereas QM/MM models converge quickly to stable values, the QM cluster models need to incorporate significantly larger parts of the second coordination sphere and surrounding water molecules to achieve convergence for certain properties. This is mainly due to the sensitivity of the QM-only models to fluctuations in the hydrogen bonding network and ligand acidity. Additionally, a hydrogen bond that is typically omitted in QM-only treatments is shown to determine the hyperfine coupling tensor of the unique Mn(III) ion by regulating the rotation plane of the ligated D1-His332 imidazole ring, the only N-donor ligand of the OEC.

Blind prediction of charged ligand binding affinities in a model binding site

Predicting absolute protein-ligand binding affinities remains a frontier challenge in ligand discovery and design. This becomes more difficult when ionic interactions are involved because of the large opposing solvation and electrostatic attraction energies. In a blind test, we examined whether alchemical free-energy calculations could predict binding affinities of 14 charged and 5 neutral compounds previously untested as ligands for a cavity binding site in cytochrome c peroxidase. In this simplified site, polar and cationic ligands compete with solvent to interact with a buried aspartate. Predictions were tested by calorimetry, spectroscopy, and crystallography. Of the 15 compounds predicted to bind, 13 were experimentally confirmed, while 4 compounds were false negative predictions. Predictions had a root-mean-square error of 1.95 kcal/mol to the experimental affinities, and predicted poses had an average RMSD of 1.7Å to the crystallographic poses. This test serves as a benchmark for these thermodynamically rigorous calculations at predicting binding affinities for charged compounds and gives insights into the existing sources of error, which are primarily electrostatic interactions inside proteins. Our experiments also provide a useful set of ionic binding affinities in a simplified system for testing new affinity prediction methods.

Mapping of Functional Groups in Metal-Organic Frameworks

We determined the heterogeneous mesoscale spatial apportionment of functional groups in a series of multivariate metal-organic frameworks (MTV-MOF-5) containing BDC (1,4-benzenedicarboxylate) linkers with different functional groups—B (BDC-NH2), E (BDC-NO2), F [(BDC-(CH3)2], H [BDC-(OC3H5)2], and I [BDC-(OC7H7)2]—using solid-state nuclear magnetic resonance measurements combined with molecular simulations. Our analysis reveals that these methods discern between random (EF), alternating (EI and EHI), and various cluster (BF) forms of functional group apportionments. This combined synthetic, characterization, and computational approach predicts the adsorptive properties of crystalline MTV-MOF systems. This methodology, developed in the context of ordered frameworks, is a first step in resolving the more general problem of spatial disorder in other ordered materials, including mesoporous materials, functionalized polymers, and defect distributions within crystalline solids.

Calculating the sensitivity and robustness of binding free energy calculations to force field parameters

Binding free energy calculations offer a thermodynamically rigorous method to compute protein-ligand binding, and they depend on empirical force fields with hundreds of parameters. We examined the sensitivity of computed binding free energies to the ligand's electrostatic and van der Waals parameters. Dielectric screening and cancellation of effects between ligand-protein and ligand-solvent interactions reduce the parameter sensitivity of binding affinity by 65%, compared with interaction strengths computed in the gas-phase. However, multiple changes to parameters combine additively on average, which can lead to large changes in overall affinity from many small changes to parameters. Using these results, we estimate that random, uncorrelated errors in force field nonbonded parameters must be smaller than 0.02 e per charge, 0.06 Å per radius, and 0.01 kcal/mol per well depth in order to obtain 68% (one standard deviation) confidence that a computed affinity for a moderately-sized lead compound will fall within 1 kcal/mol of the true affinity, if these are the only sources of error considered.

Insights on the solubility of CO2 in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide from the microscopic point of view

Emissions of greenhouse gases due to human activities have been well documented as well as the effects on global warming resulting from it. Efforts to reduce greenhouse gases at the source are crucial to curb climate change, but due to insignificant economic incentives to reduce usage of fossil fuels, not a lot of progress has been made by this route. This necessitates additional measures to reduce the occurrence of greenhouse gases in the atmosphere. Here we used theoretical methods to study the solubility of carbon dioxide in ionic liquids (ILs) since sequestration of CO2 in ILs has been proposed as a possible technology for reducing the emissions of CO2 to the atmosphere. Ionic liquids form a class of solvents with melting temperatures below 100 °C and, due to very low vapor pressures, which are not volatile. We have performed molecular dynamics (MD) simulations of 1-ethyl-3-methylimidazolium (C2mim) bis(trifluoromethylsulfonyl)imide (Tf2N) and its mixtures with carbon dioxide in order to investigate the CO2 concentration effect on the CO2-cation and CO2-anion interactions. A systematic investigation of CO2 concentration effects on resulting equilibrium liquid structure, and the local environment of the ions is provided. The Quantum Theory of Atoms in Molecules (QTAIM) was used to determine the interaction energy for CO2-cation and CO2-anion complexes from uncorrelated structures derived from MD simulations. A spatial distribution function analysis demonstrates the specific interactions between CO2 and the ionic liquid. Our findings indicate that the total volume of the system increases with the CO2 concentration, with a molar volume of CO2 of about 0.038 L/mol, corresponding to liquid CO2 under a pressure of 100 bar. In other words, the IL effectively pressurizes the CO2 inside its matrix. The thermodynamics of CO2 solvation in C2 min-Tf2N were computed using free energy techniques, and the solubility of CO2 is found to be higher in this IL (-3.7 ± 1 kcal/mol) than in water (+0.2 kJ/mol), predominantly due to anion-CO2 interactions.

The Role of Aromatic-Aromatic Interactions in Strand-Strand Stabilization of β-Sheets

Aromatic-aromatic interactions have long been believed to play key roles in protein structure, folding, and binding functions. However, we still lack full understanding of the contributions of aromatic-aromatic interactions to protein stability and the timing of their formation during folding. Here, using an aromatic ladder in the β-barrel protein, cellular retinoic acid-binding protein 1 (CRABP1), as a case study, we find that aromatic π stacking plays a greater role in the Phe65-Phe71 cross-strand pair, while in another pair, Phe50-Phe65, hydrophobic interactions are dominant. The Phe65-Phe71 pair spans β-strands 4 and 5 in the β-barrel, which lack interstrand hydrogen bonding, and we speculate that it compensates energetically for the absence of strand-strand backbone interactions. Using perturbation analysis, we find that both aromatic-aromatic pairs form after the transition state for folding of CRABP1, thus playing a role in the final stabilization of the β-sheet rather than in its nucleation as had been earlier proposed. The aromatic interaction between strands 4 and 5 in CRABP1 is highly conserved in the intracellular lipid-binding protein (iLBP) family, and several lines of evidence combine to support a model wherein it acts to maintain barrel structure while allowing the dynamic opening that is necessary for ligand entry. Lastly, we carried out a bioinformatics analysis and found 51 examples of aromatic-aromatic interactions across non-hydrogen-bonded β-strands outside the iLBPs, arguing for the generality of the role played by this structural motif.

Molecular dynamics study of carbon nanotube as a potential dual-functional inhibitor of HIV-1 integrase

HIV-1 integrase (IN) plays an important role in integrating viral DNA into human genome, which has been considered as the drug target for anti-AIDS therapy. The appearance of drug-resistance mutants urgently requires novel inhibitors that act on non-active site of HIV-1 IN. Nanoparticles have such unique geometrical and chemical properties, which inspires us that nanoparticles like nanotubes may serve as better HIV-1 IN inhibitors than the conventional inhibitors. To test this hypothesis, we performed molecular dynamics (MD) simulation to study the binding of a carbon nanotube (CNT) to a full-length HIV-1 IN. The results showed that the CNT could stably bind to the C-terminal domain (CTD) of HIV-1 IN. The CNT also induced a domain-shift which disrupted the binding channel for viral DNA. Further MD simulation showed that a HIV-1 IN inhibitor, 5ClTEP was successfully sealed inside the uncapped CNT. These results indicate that the CNT may serve as a potential dual-functional HIV-1 IN inhibitor, not only inducing conformation change as an allosteric inhibitor but also carrying small-molecular inhibitors as a drug delivery system.

GSATools: analysis of allosteric communication and functional local motions using a structural alphabet

Motivation: GSATools is a free software package to analyze conformational ensembles and to detect functional motions in proteins by means of a structural alphabet. The software integrates with the widely used GROMACS simulation package and can generate a range of graphical outputs. Three applications can be supported: (i) investigation of the conformational variability of local structures; (ii) detection of allosteric communication; and (iii) identification of local regions that are critical for global functional motions. These analyses provide insights into the dynamics of proteins and allow for targeted design of functional mutants in theoretical and experimental studies.

Availability: The C source code of the GSATools, along with a set of pre-compiled binaries, is freely available under GNU General Public License from http://mathbio.nimr.mrc.ac.uk/wiki/GSATools.

Contact:alessandro.pandini@kcl.ac.uk or jkleinj@nimr.mrc.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

The influence of positional isomerism on G-quadruplex binding and anti-proliferative activity of tetra-substituted naphthalene diimide compounds

The synthesis together with biophysical and biological evaluation of a series of tetra-substituted naphthalene diimide (ND) compounds, are presented. These compounds are positional isomers of a recently-described series of quadruplex-binding ND derivatives, in which the two N-methyl-piperidine-alkyl side-chains have now been interchanged with the positions of side-chains bearing a range of end-groups. Molecular dynamics simulations of a pair of positional isomers are in accord with the quadruplex stabilization and biological data for these compounds. Analysis of structure-activity data indicates that for compounds where the side-chains are not of equivalent length then the positional isomers described here tend to have improved cell proliferation potency and in some instances, superior quadruplex stabilization ability.

Progress in developing Poisson-Boltzmann equation solvers

This review outlines the recent progress made in developing more accurate and efficient solutions to model electrostatics in systems comprised of bio-macromolecules and nano-objects, the last one referring to objects that do not have biological function themselves but nowadays are frequently used in biophysical and medical approaches in conjunction with bio-macromolecules. The problem of modeling macromolecular electrostatics is reviewed from two different angles: as a mathematical task provided the specific definition of the system to be modeled and as a physical problem aiming to better capture the phenomena occurring in the real experiments. In addition, specific attention is paid to methods to extend the capabilities of the existing solvers to model large systems toward applications of calculations of the electrostatic potential and energies in molecular motors, mitochondria complex, photosynthetic machinery and systems involving large nano-objects.

Tubulin tyrosine nitration regulates microtubule organization in plant cells

During last years, selective tyrosine nitration of plant proteins gains importance as well-recognized pathway of direct nitric oxide (NO) signal transduction. Plant microtubules are one of the intracellular signaling targets for NO, however, the molecular mechanisms of NO signal transduction with the involvement of cytoskeletal proteins remain to be elucidated. Since biochemical evidence of plant α-tubulin tyrosine nitration has been obtained recently, potential role of this posttranslational modification in regulation of microtubules organization in plant cell is estimated in current paper. It was shown that 3-nitrotyrosine (3-NO2-Tyr) induced partially reversible Arabidopsis primary root growth inhibition, alterations of root hairs morphology and organization of microtubules in root cells. It was also revealed that 3-NO2-Tyr intensively decorates such highly dynamic microtubular arrays as preprophase bands, mitotic spindles and phragmoplasts of Nicotiana tabacum Bright Yellow-2 (BY-2) cells under physiological conditions. Moreover, 3D models of the mitotic kinesin-8 complexes with the tail of detyrosinated, tyrosinated and tyrosine nitrated α-tubulin (on C-terminal Tyr 450 residue) from Arabidopsis were reconstructed in silico to investigate the potential influence of tubulin nitrotyrosination on the molecular dynamics of α-tubulin and kinesin-8 interaction. Generally, presented data suggest that plant α-tubulin tyrosine nitration can be considered as its common posttranslational modification, the direct mechanism of NO signal transduction with the participation of microtubules under physiological conditions and one of the hallmarks of the increased microtubule dynamics.