New Insights into the LANCL2- Binding Mode towards the Evaluation of New LANCL Agonists

The lanthionine synthetase C-like (LANCL) proteins include LANCL2, which is expressed in the central nervous system (CNS) and in peripheral tissues. LANCL2 exhibits glutathionylation activity and is involved in the neutralization of reactive electrophiles. Several studies explored LANCL2 activation as a validated pharmacological target for diabetes and inflammatory bowel disease. In this context, LANCL2 was found to bind the natural product abscisic acid (ABA), whose pre-clinical effectiveness in different inflammatory diseases was reported in the literature. More recently, LANCL2 attracted more attention as a valuable resource in the field of neurodegenerative disorders. ABA was found to regulate neuro-inflammation and synaptic plasticity to enhance learning and memory, exhibiting promising neuroprotective effects. Up until now, a limited number of LANCL2 ligands are known; among them, BT-11 is the only compound patented and investigated for its anti-inflammatory properties. To guide the design of novel putative LANCL2 agonists, a computational study including molecular docking and long molecular dynamic (MD) simulations of both ABA and BT-11 was carried out. The results pointed out the main LANCL2 ligand chemical features towards the following virtual screening of a novel putative LANCL2 agonist (AR-42). Biochemical assays on rat H9c2 cardiomyocytes showed a similar, LANCL2-mediated stimulation by BT-11 and by AR-42 of the mitochondrial proton gradient and of the transcriptional activation of the AMPK/PGC-1α/Sirt1 axis, the master regulator of mitochondrial function, effects that are previously observed with ABA. These results may allow the development of LANCL2 agonists for the treatment of mitochondrial dysfunction, a common feature of chronic and degenerative diseases.

Molecular Dynamics and Machine Learning Give Insights on the Flexibility-Activity Relationships in Tyrosine Kinome

Tyrosine kinases are a subfamily of kinases with critical roles in cellular machinery. Dysregulation of their active or inactive forms is associated with diseases like cancer. This study aimed to holistically understand their flexibility-activity relationships, focusing on pockets and fluctuations. We studied 43 different tyrosine kinases by collecting 120 μs of molecular dynamics simulations, pocket and residue fluctuation analysis, and a complementary machine learning approach. We found that the inactive forms often have increased flexibility, particularly at the DFG motif level. Noteworthy, thanks to these long simulations combined with a decision tree, we identified a semiquantitative fluctuation threshold of the DGF+3 residue over which the kinase has a higher probability to be in the inactive form.

Probing Interplays between Human XBP1u Translational Arrest Peptide and 80S Ribosome

The ribosome stalling mechanism is a crucial biological process, yet its atomistic underpinning is still elusive. In this framework, the human XBP1u translational arrest peptide (AP) plays a central role in regulating the unfolded protein response (UPR) in eukaryotic cells. Here, we report multimicrosecond all-atom molecular dynamics simulations designed to probe the interactions between the XBP1u AP and the mammalian ribosome exit tunnel, both for the wild type AP and for four mutant variants of different arrest potencies. Enhanced sampling simulations allow investigating the AP release process of the different variants, shedding light on this complex mechanism. The present outcomes are in qualitative/quantitative agreement with available experimental data. In conclusion, we provide an unprecedented atomistic picture of this biological process and clear-cut insights into the key AP–ribosome interactions.

Relevance of Ebola virus VP35 homo-dimerization on the type I interferon cascade inhibition

Ebola virus high lethality relies on its ability to efficiently bypass the host innate antiviral response, which senses the viral dsRNA through the RIG-I receptor and induces type I interferon α/β production. In the bypassing action, the Ebola virus protein VP35 plays a pivotal role at multiple levels of the RIG-I cascade, masking the viral 5′-triphosphorylated dsRNA from RIG-I, and interacting with other cascade components. The VP35 type I interferon inhibition is exerted by the C-terminal domain, while the N-terminal domain, containing a coiled-coil region, is primarily required for oligomerization. However, mutations at key VP35 residues L90/93/107A (VP35-3m) in the coiled-coil region were reported to affect oligomerization and reduce type I interferon antagonism, indicating a possible but unclear role of homo-oligomerization on VP35 interaction with the RIG-I pathway components. In this work, we investigated the VP35 dimerization thermodynamics and its contribution to type I interferon antagonism by computational and biological methods. Focusing on the coiled-coil region, we combined coarse-grained and all-atom simulations on wild type VP35 and VP35-3m homo-dimerization. According to our results, wild type VP35 coiled-coil is able to self-assemble into dimers, while VP35-3m coiled-coil shows poor propensity to even dimerize. Free-energy calculations confirmed the key role of L90, L93 and L107 in stabilizing the coiled-coil homo-dimeric structure. In vitro type I interferon antagonism studies, using full-length wild type VP35 and VP35-3m, revealed that VP35 homo-dimerization is an essential preliminary step for dsRNA binding, which appears to be the main factor of the VP35 RIG-I cascade inhibition, while it is not essential to block the other steps.

Dynamics behind affinity maturation of an anti-HCMV antibody family influencing antigen binding

The investigation of antibody affinity maturation and its effects on antigen binding is important with respect to understanding the regulation of the immune response. To shed light on this crucial process, we analyzed two Igs neutralizing the human cytomegalovirus: the primary germline antibody M2J1 and its related mature antibody 8F9. Both antibodies target the AD-2S1 epitope of the gB envelope protein and are considered to establish similar interactions with the cognate antigen. We used molecular dynamics simulations to understand the effect of mutations on the antibody-antigen interactions. The results provide a qualitative explanation for the increased 8F9 peptide affinity compared with that of M2J1. The emerging atomistic-detailed description of these complexes reveals the molecular effects of the somatic hypermutations occurring during affinity maturation.

Kissing loop interaction in adenine riboswitch: insights from umbrella sampling simulations

INTRODUCTION

Riboswitches are cis-acting regulatory RNA elements prevalently located in the leader sequences of bacterial mRNA. An adenine sensing riboswitch cis-regulates adeninosine deaminase gene (add) in Vibrio vulnificus. The structural mechanism regulating its conformational changes upon ligand binding mostly remains to be elucidated. In this open framework it has been suggested that the ligand stabilizes the interaction of the distal "kissing loop" complex. Using accurate full-atom molecular dynamics with explicit solvent in combination with enhanced sampling techniques and advanced analysis methods it could be possible to provide a more detailed perspective on the formation of these tertiary contacts.

METHODS

In this work, we used umbrella sampling simulations to study the thermodynamics of the kissing loop complex in the presence and in the absence of the cognate ligand. We enforced the breaking/formation of the loop-loop interaction restraining the distance between the two loops. We also assessed the convergence of the results by using two alternative initialization protocols. A structural analysis was performed using a novel approach to analyze base contacts.

RESULTS

Contacts between the two loops were progressively lost when larger inter-loop distances were enforced. Inter-loop Watson-Crick contacts survived at larger separation when compared with non-canonical pairing and stacking interactions. Intra-loop stacking contacts remained formed upon loop undocking. Our simulations qualitatively indicated that the ligand could stabilize the kissing loop complex. We also compared with previously published simulation studies.

DISCUSSION AND CONCLUSIONS

Kissing complex stabilization given by the ligand was compatible with available experimental data. However, the dependence of its value on the initialization protocol of the umbrella sampling simulations posed some questions on the quantitative interpretation of the results and called for better converged enhanced sampling simulations.

Using reweighted pulling simulations to characterize conformational changes in riboswitches

Riboswitches are RNA sequences located in noncoding portions of mRNA that can sense specific ligands and subsequently control gene expression. The ligand-binding event induces conformational changes in the riboswitch that are then transmitted to the gene expression apparatus. Probing the mechanisms of such a fine regulation at atomic resolution is very difficult experimentally and molecular dynamics (MD) could be used to quantify the ligand-dependent behavior of a riboswitch. However, since the accessible time scale of fully atomistic simulations is limited, this can only be done using enhanced sampling techniques. Here, we discuss the application of steered MD to the characterization of the ligand-dependent stability of the aptamer terminal helix in the add adenine-sensing riboswitch. The employed techniques are discussed in detail and sample input files are provided. We show that with a limited computational effort it is possible to quantify, in terms of free energy, the stacking interaction between the ligand and the terminal helix, obtaining results in agreement with thermodynamic experiments.

The role of nucleobase interactions in RNA structure and dynamics

The intricate network of interactions observed in RNA three-dimensional structures is often described in terms of a multitude of geometrical properties, including helical parameters, base pairing/stacking, hydrogen bonding and backbone conformation. We show that a simple molecular representation consisting in one oriented bead per nucleotide can account for the fundamental structural properties of RNA. In this framework, canonical Watson-Crick, non-Watson-Crick base-pairing and base-stacking interactions can be unambiguously identified within a well-defined interaction shell. We validate this representation by performing two independent, complementary tests. First, we use it to construct a sequence-independent, knowledge-based scoring function for RNA structural prediction, which compares favorably to fully atomistic, state-of-the-art techniques. Second, we define a metric to measure deviation between RNA structures that directly reports on the differences in the base–base interaction network. The effectiveness of this metric is tested with respect to the ability to discriminate between structurally and kinetically distant RNA conformations, performing better compared to standard techniques. Taken together, our results suggest that this minimalist, nucleobase-centric representation captures the main interactions that are relevant for describing RNA structure and dynamics.

Ligand-induced stabilization of the aptamer terminal helix in the add adenine riboswitch

Riboswitches are structured mRNA elements that modulate gene expression. They undergo conformational changes triggered by highly specific interactions with sensed metabolites. Among the structural rearrangements engaged by riboswitches, the forming and melting of the aptamer terminal helix, the so-called P1 stem, is essential for genetic control. The structural mechanisms by which this conformational change is modulated upon ligand binding mostly remain to be elucidated. Here, we used pulling molecular dynamics simulations to study the thermodynamics of the P1 stem in the add adenine riboswitch. The P1 ligand-dependent stabilization was quantified in terms of free energy and compared with thermodynamic data. This comparison suggests a model for the aptamer folding in which direct P1-ligand interactions play a minor role on the conformational switch when compared with those related to the ligand-induced aptamer preorganization.