Exploring the concept of interaction computing through the discrete algebraic analysis of the Belousov–Zhabotinsky reaction

Interaction computing (IC) aims to map the properties of integrable low-dimensional non-linear dynamical systems to the discrete domain of finite-state automata in an attempt to reproduce in software the self-organizing and dynamically stable properties of sub-cellular biochemical systems. As the work reported in this paper is still at the early stages of theory development it focuses on the analysis of a particularly simple chemical oscillator, the Belousov–Zhabotinsky (BZ) reaction. After retracing the rationale for IC developed over the past several years from the physical, biological, mathematical, and computer science points of view, the paper presents an elementary discussion of the Krohn–Rhodes decomposition of finite-state automata, including the holonomy decomposition of a simple automaton, and of its interpretation as an abstract positional number system. The method is then applied to the analysis of the algebraic properties of discrete finite-state automata derived from a simplified Petri net model of the BZ reaction. In the simplest possible and symmetrical case the corresponding automaton is, not surprisingly, found to contain exclusively cyclic groups. In a second, asymmetrical case, the decomposition is much more complex and includes five different simple non-abelian groups whose potential relevance arises from their ability to encode functionally complete algebras. The possible computational relevance of these findings is discussed and possible conclusions are drawn.     

3D NMR structure of a complex between the amyloid beta peptide (1–40) and the polyphenol ε-viniferin glucoside: Implications in Alzheimer's disease

Background

Alzheimer's disease (AD) is a progressive neurodegenerative disorder. There is a consensus that Aβ is a pathologic agent and that its toxic effects, which are at present incompletely understood, may occur through several potential mechanisms. Polyphenols are known to have wide-ranging properties with regard to health and for helping to prevent various diseases like neurodegenerative disorders. Thus inhibiting the formation of toxic Aβ assemblies is a reasonable hypothesis to prevent and perhaps treat AD

Methods

Solution NMR and molecular modeling were used to obtain more information about the interaction between the Aβ_1–40 and the polyphenol ε-viniferin glucoside (EVG) and particularly the Aβ residues involved in the complex.

Results

The study demonstrates the formation of a complex between two EVG molecules and Aβ_1–40 in peptide characteristic regions that could be in agreement with the inhibition of aggregation. Indeed, in previous studies, we reported that EVG strongly inhibited in vitro the fibril formation of the full length peptides Aβ_1–40 and Aβ_1–42, and had a strong protective effect against PC12 cell death induced by these peptides.

Conclusion

For the full length peptide Aβ_1–40, the binding sites observed could explain the EVG inhibitory effect on fibrillization and thus prevent amyloidogenic neurotoxicity.

General significance

Even though this interaction might be important at the biological level to explain the protective effect of polyphenols in neurodegenerative diseases, caution is required when extrapolating this in vitro model to human physiology.     

Binding of the plant alkaloid aristololactam-β-d-glucoside and antitumor antibiotic daunomycin to single stranded polyribonucleotides

Background

Interaction of the plant alkaloid aristololactam-β-d-glucoside and the antitumor drug daunomycin with single stranded RNAs poly(G), poly(I), poly(C) and poly(U) has been investigated.

Methods

Biophysical techniques of absorption, fluorescence, competition dialysis, circular dichroism, and microcalorimetry have been used.

Results

Absorption and fluorescence studies have revealed noncooperative binding of ADG and DAN to the single stranded RNAs. The binding affinity of ADG varied as poly(G) > poly(I) > > poly(C) > poly(U). The affinity of DAN was one order higher than that of ADG and varied as poly(G) > poly(I) > poly(U) > poly(C). This binding preference was further confirmed by competition dialysis assay. The thermodynamics of the binding was characterised to be favourable entropy and enthalpic terms but their contributions were different for different systems. The major non-polyelectrolytic contribution to the binding revealed from salt dependent data appears to be arising mostly from stacking of DAN and ADG molecules with the bases leading to partial intercalation to single stranded RNA structures. Small negative heat capacity values have been observed in all the four cases.

Conclusions

This study presents the comparative structural and thermodynamic profiles of the binding of aristololactam-β-d-glucoside and daunomycin to single stranded polyribonucleotides.

General significance

These results suggest strong, specific but differential binding of these drug molecules to the single stranded RNAs and highlight the role of their structural differences in the interaction profile.     

Incorporation of impurity anions into DSP: insights into structure and stability from computer modelling

DSP is an important by-product of alumina production via the Bayer process. Under Western Australian processing conditions, the DSP has a sodalite-type structure that can incorporate anions within its framework. This is particularly useful for removal of impurity anions from liquor recycled in the circuit. As a first step to gaining a fundamental understanding of the incorporation process, we have undertaken molecular mechanics calculations to examine the interaction energy between a series of anions and the sodalite framework, as a measure of the affinity of the anions for the sodalite cage. Our calculations predict that the ions have an increased affinity for the cage along the series aluminate, chloride, carbonate, sulfate and hydroxide. These calculations have successfully predicted the trends that we observe from competitive-uptake experiments in our laboratory.     

Molecular dynamics simulation of fluorination effect for solvation of trifluoromethylbenzoic acid isomers in supercritical carbon dioxide

A molecular dynamics (MD) simulation was applied to carbon dioxide+trifluoromethylbenzoic acid isomer and carbon dioxide+methylbenzoic acid isomer systems to investigate the interactions between carbon dioxide and the solutes. The pair correlation functions between the carbon dioxide and trifluoromethyl group or methyl group in the solutes were calculated to study the fluorination effect of solvation. As a result, it was found that the interactions between carbon dioxide and trifluoromethyl group in trifluoromethylbenzoic acid isomers were stronger than those between carbon dioxide and the methyl group in methylbenzoic acid isomers. The simulation results had the same tendency as the experimental solubility enhancements and coincided with the trend of the interaction parameters of the Peng-Robinson equation of state that were determined from the solubility data.     

Alpha7 Helix Plays an Important Role in the Conformational Stability of PTP1B

Abstract

The C-terminus of Protein Tyrosine Phosphatase 1B (PTP1B) includes an α-helix (α7), which forms an allosteric binding site 20 Å away from the active site. This helix is specific to PTP1B and its truncation decreases the catalytic activity significantly. Here, molecular dynamics (MD) simulations in the presence and absence of α7 were performed to investigate the role played by α7. The highly mobile α7 was found to maintain its contacts with loop 11 (L11)- α3 helix throughout the simulations. The interactions of Tyr152 on L11, Tyr176, Thr177 on the catalytically important WPD loop and Ser190 on α3 are important for the conformational stability and the concerted motions of the regions surrounding the WPD loop. In the absence of α7, L11 and WPD loop move away from their crystal structure conformations, resulting in the loss of the interactions in this region, and a decrease in the residue displacement correlations in the vicinity of WPD loop. Therefore, we suggest that one of the functionally important roles of α7 may be to limit the L11 and α3 motions, and, facilitate the WPD loop motions. Truncation of α7 in PTP1B is found to affect distant regions as well, such as the substrate recognition site and the phosphate binding-loop (P-loop), changing the conformations of these regions significantly. Our results show that the PTP1B specific α7 is important for the conformation and dynamics of the WPD loop, and also may play a role in ligand binding.     

Structural Basis for the Different Stability and Activity between the Cdk5 Complexes with p35 and p39 Activators

Cyclin-dependent kinase 5 (Cdk5) is a brain-specific membrane-bound protein kinase that is activated by binding to the p35 or p39 activator. Previous studies have focused on p35-Cdk5, and little is known regarding p39-Cdk5. The lack of functional understanding of p39-Cdk5 is due, in part, to the labile property of p39-Cdk5, which dissociates and loses kinase activity in nonionic detergent conditions. Here we investigated the structural basis for the instability of p39-Cdk5. p39 and p35 contain N-terminal p10 regions and C-terminal Cdk5 activation domains (AD). Although p35 and p39 show higher homology in the C-terminal AD than the N-terminal region, the difference in stability is derived from the C-terminal AD. Based on the crystal structures of the p25 (p35 C-terminal region including AD)-Cdk5 complex, we simulated the three-dimensional structure of the p39 AD-Cdk5 complex and found differences in the hydrogen bond network between Cdk5 and its activators. Three amino acids of p35, Asp-259, Asn-266, and Ser-270, which are involved in hydrogen bond formation with Cdk5, are changed to Gln, Gln, and Pro in p39. Because these three amino acids in p39 do not participate in hydrogen bond formation, we predicted that the number of hydrogen bonds between p39 and Cdk5 was reduced compared with p35 and Cdk5. Using substitution mutants, we experimentally validated that the difference in the hydrogen bond network contributes to the different properties between Cdk5 and its activators.     

A Single Divergent Exon Inhibits Ankyrin-B Association with the Plasma Membrane

Vertebrate ankyrin-B and ankyrin-G exhibit divergent subcellular localization and function despite their high sequence and structural similarity and common origin from a single ancestral gene at the onset of chordate evolution. Previous studies of ankyrin family diversity have focused on the C-terminal regulatory domain. Here, we identify an ankyrin-B-specific linker peptide connecting the ankyrin repeat domain to the ZU5₂-UPA module that inhibits binding of ankyrin-B to membrane protein partners E-cadherin and neurofascin 186 and prevents association of ankyrin-B with epithelial lateral membranes as well as neuronal plasma membranes. The residues of the ankyrin-B linker required for autoinhibition are encoded by a small exon that is highly divergent between ankyrin family members but conserved in the ankyrin-B lineage. We show that the ankyrin-B linker suppresses activity of the ANK repeat domain through an intramolecular interaction, likely with a groove on the surface of the ANK repeat solenoid, thereby regulating the affinities between ankyrin-B and its binding partners. These results provide a simple evolutionary explanation for how ankyrin-B and ankyrin-G have acquired striking differences in their plasma membrane association while maintaining overall high levels of sequence similarity.     

Roles of Long-range Electrostatic Domain Interactions and K+ in Phosphoenzyme Transition of Ca2+-ATPase

Sarcoplasmic reticulum Ca²⁺-ATPase couples the motions and rearrangements of three cytoplasmic domains (A, P, and N) with Ca²⁺ transport. We explored the role of electrostatic force in the domain dynamics in a rate-limiting phosphoenzyme (EP) transition by a systematic approach combining electrostatic screening with salts, computer analysis of electric fields in crystal structures, and mutations. Low KCl concentration activated and increasing salt above 0.1 m inhibited the EP transition. A plot of the logarithm of the transition rate versus the square of the mean activity coefficient of the protein gave a linear relationship allowing division of the activation energy into an electrostatic component and a non-electrostatic component in which the screenable electrostatic forces are shielded by salt. Results show that the structural change in the transition is sterically restricted, but that strong electrostatic forces, when K⁺ is specifically bound at the P domain, come into play to accelerate the reaction. Electric field analysis revealed long-range electrostatic interactions between the N and P domains around their hinge. Mutations of the residues directly involved and other charged residues at the hinge disrupted in parallel the electric field and the structural transition. Favorable electrostatics evidently provides a low energy path for the critical N domain motion toward the P domain, overcoming steric restriction. The systematic approach employed here is, in general, a powerful tool for understanding the structural mechanisms of enzymes.