Bringing LTL Model Checking to Biologists

The BioModelAnalyzer (BMA) is a web based tool for the development of discrete models of biological systems. Through a graphical user interface, it allows rapid development of complex models of gene and protein interaction networks and stability analysis without requiring users to be proficient computer programmers. Whilst stability is a useful specification for testing many systems, testing temporal specifications in BMA presently requires the user to perform simulations. Here we describe the LTL module, which includes a graphical and natural language interfaces to testing LTL queries. The graphical interface allows for graphical construction of the queries and presents results visually in keeping with the current style of BMA. The Natural language interface complements the graphical interface by allowing a gentler introduction to formal logic and exposing educational resources

A two-domain elevator mechanism for sodium/proton antiport

Sodium/proton (Na+/H+) antiporters, located at the plasma membrane in every cell, are vital for cell homeostasis1. In humans, their dysfunction has been linked to diseases, such as hypertension, heart failure and epilepsy, and they are well-established drug targets2. The best understood model system for Na+/H+ antiport is NhaA from Escherichia coli1, 3, for which both electron microscopy and crystal structures are available4, 5, 6. NhaA is made up of two distinct domains: a core domain and a dimerization domain. In the NhaA crystal structure a cavity is located between the two domains, providing access to the ion-binding site from the inward-facing surface of the protein1, 4. Like many Na+/H+ antiporters, the activity of NhaA is regulated by pH, only becoming active above pH 6.5, at which point a conformational change is thought to occur7. The only reported NhaA crystal structure so far is of the low pH inactivated form4. Here we describe the active-state structure of a Na+/H+ antiporter, NapA from Thermus thermophilus, at 3 Å resolution, solved from crystals grown at pH 7.8. In the NapA structure, the core and dimerization domains are in different positions to those seen in NhaA, and a negatively charged cavity has now opened to the outside. The extracellular cavity allows access to a strictly conserved aspartate residue thought to coordinate ion binding1, 8, 9 directly, a role supported here by molecular dynamics simulations. To alternate access to this ion-binding site, however, requires a surprisingly large rotation of the core domain, some 20° against the dimerization interface. We conclude that despite their fast transport rates of up to 1,500 ions per second3, Na+/H+ antiporters operate by a two-domain rocking bundle model, revealing themes relevant to secondary-active transporters in general

Helical unwinding and side-chain unlocking unravel the outward open conformation of the melibiose transporter

Molecular dynamics simulations have been used to study the alternate access mechanism of the melibiose transporter from Escherichia coli. Starting from the outward-facing partially occluded form, 2 out of 12 simulations produced an outward full open form and one partially open, whereas the rest yielded fully or partially occluded forms. The shape of the outward-open form resembles other outward-open conformations of secondary transporters. During the transporter opening, conformational changes in some loops are followed by changes in the periplasm region of transmembrane helix 7. Helical curvature relaxation and unlocking of hydrophobic and ionic locks promote the outward opening of the transporter making accessible the substrate binding site. In particular, FRET studies on mutants of conserved aromatic residues of extracellular loop 4 showed lack of substrate binding, emphasizing the importance of this loop for making crucial interactions that control the opening of the periplasmic side. This study indicates that the alternate access mechanism for the melibiose transporter fits better into a flexible gating mechanism rather than the archetypical helical rigid-body rocker-switch mechanism

Membrane cholesterol access into a G-protein-coupled receptor

Cholesterol is a key component of cell membranes with a proven modulatory role on the function and ligand-binding properties of G-protein-coupled receptors (GPCRs). Crystal structures of prototypical GPCRs such as the adenosine A2A receptor (A2AR) have confirmed that cholesterol finds stable binding sites at the receptor surface suggesting an allosteric role of this lipid. Here we combine experimental and computational approaches to show that cholesterol can spontaneously enter the A2AR-binding pocket from the membrane milieu using the same portal gate previously suggested for opsin ligands. We confirm the presence of cholesterol inside the receptor by chemical modification of the A2AR interior in a biotinylation assay. Overall, we show that cholesterol’s impact on A2AR-binding affinity goes beyond pure allosteric modulation and unveils a new interaction mode between cholesterol and the A2AR that could potentially apply to other GPCRs