Site-specific incorporation of a fluorescent nucleobase analog enhances i-motif stability and allows monitoring of i-motif folding inside cells

The i-motif is an intriguing non-canonical DNA structure, whose role in the cell is still controversial. Development of methods to study i-motif formation under physiological conditions in living cells is necessary to study its potential biological functions. The cytosine analog 1,3-diaza-2-oxophenoxazine (tCO) is a fluorescent nucleobase able to form either hemiprotonated base pairs with cytosine residues, or neutral base pairs with guanines. We show here that when tCO is incorporated in the proximity of a G:C:G:C minor groove tetrad, it induces a strong thermal and pH stabilization, resulting in i-motifs with Tm of 39ºC at neutral pH. The structural determination by NMR methods reveals that the enhanced stability is due to a large stacking interaction between the guanines of the tetrad with the tCO nucleobase, which forms a tCO:C+ in the folded structure at unusually-high pHs, leading to an increased quenching in its fluorescence at neutral conditions. This quenching is much lower when tCO is base-paired to guanines and totally disappears when the oligonucleotide is unfolded. By taking profit of this property, we have been able to monitor i-motif folding in cells.

From crystallographic data to the solution structure of photoreceptors: the case of the AppA BLUF domain

Photoreceptor proteins bind a chromophore, which, upon light absorption, modifies its geometry or its interactions with the protein, finally inducing the structural change needed to switch the protein from an inactive to an active or signaling state. In the Blue Light-Using Flavin (BLUF) family of photoreceptors, the chromophore is a flavin and the changes have been connected with a rearrangement of the hydrogen bond network around it on the basis of spectroscopic changes measured for the dark-to-light conversion. However, the exact conformational change triggered by the photoexcitation is still elusive mainly because a clear consensus on the identity not only of the light activated state but also of the dark one has not been achieved. Here, we present an integrated investigation that combines microsecond MD simulations starting from the two conflicting crystal structures available for the AppA BLUF domain with calculations of NMR, IR and UV-Vis spectra using a polarizable QM/MM approach. Thanks to such a combined analysis of the three different spectroscopic responses, a robust characterization of the structure of the dark state in solution is given together with the uncovering of important flaws of the most popular molecular mechanisms present in the literature for the dark-to-light activation.

Ultrafast Transient Infrared Spectroscopy of Photoreceptors with Polarizable QM/MM Dynamics

Ultrafast transient infrared (TRIR) spectroscopy is widely used to measure the excitation-induced structural changes of protein-bound chromophores. Here, we design a novel and general strategy to compute TRIR spectra of photoreceptors by combining μs-long MM molecular dynamics with ps-long QM/AMOEBA Born-Oppenheimer molecular dynamics (BOMD) trajectories for both ground and excited electronic states. As a proof of concept, the strategy is here applied to AppA, a blue-light-utilizing flavin (BLUF) protein, found in bacteria. We first analyzed the short-time evolution of the embedded flavin upon excitation revealing that its dynamic Stokes shift is ultrafast and mainly driven by the internal reorganization of the chromophore. A different normal-mode representation was needed to describe ground- and excited-state IR spectra. In this way, we could assign all of the bands observed in the measured transient spectrum. In particular, we could characterize the flavin isoalloxazine-ring region of the spectrum, for which a full and clear description was missing.

The structural changes in the signaling mechanism of bacteriophytochromes in solution revealed by a multiscale computational investigation

Phytochromes are red-light sensing proteins, with important light-regulatory roles in different organisms, which are capturing an increasing interest in bioimaging and optogenetics. Upon absorption of light by the embedded bilin chromophore, they undergo structural changes that extend from the chromophore to the protein and finally drive the biological function. Up to now, the underlying mechanism still has to be characterized fully. Here we investigate the Pfr activated form of a bacterial phytochrome, by combining extensive molecular dynamics simulations with a polarizable QM/MM description of the spectroscopic properties, revealing a large structure relaxation in solution, compared to the crystal structure, both in the chromophore-binding pocket and in the overall structure of the phytochrome. Our results indicate that the final opening of the dimeric structure is preceded by an important internal reorganization of the phytochrome specific (PHY) domain involving a bend of the helical spine connecting the PHY domain with the chromophore-binding domain, opening the way to a new understanding of the activation pathway.

Role of Chemical Dynamics Simulations in Mass Spectrometry Studies of Collision-Induced Dissociation and Collisions of Biological Ions with Organic Surfaces

In this article, a perspective is given of chemical dynamics simulations of collisions of biological ions with surfaces and of collision-induced dissociation (CID) of ions. The simulations provide an atomic-level understanding of the collisions and, overall, are in quite good agreement with experiment. An integral component of ion/surface collisions is energy transfer to the internal degrees of freedom of both the ion and the surface. The simulations reveal how this energy transfer depends on the collision energy, incident angle, biological ion, and surface. With energy transfer to the ion's vibration fragmentation may occur, i.e. surface-induced dissociation (SID), and the simulations discovered a new fragmentation mechanism, called shattering, for which the ion fragments as it collides with the surface. The simulations also provide insight into the atomistic dynamics of soft-landing and reactive-landing of ions on surfaces. The CID simulations compared activation by multiple "soft" collisions, resulting in random excitation, versus high energy single collisions and nonrandom excitation. These two activation methods may result in different fragment ions. Simulations provide fragmentation products in agreement with experiments and, hence, can provide additional information regarding the reaction mechanisms taking place in experiment. Such studies paved the way on using simulations as an independent and predictive tool in increasing fundamental understanding of CID and related processes.

Elucidating the role of structural fluctuations, and intermolecular and vibronic interactions in the spectroscopic response of a bacteriophytochrome

Molecular dynamics and a multiscale polarizable QM/MM strategy allow reproducing absorption, circular dichroism, and resonance Raman spectra of a bacteriophytochrome.

Fragmentation Spectra Prediction and DNA Adducts Structural Determination

In this work, chemical dynamics simulations were optimized and used to predict fragmentation mass spectra for DNA adduct structural determination. O⁶-methylguanine (O⁶-Me-G) was used as a simple model adduct to calculate theoretical spectra for comparison with measured high-resolution fragmentation data. An automatic protocol was established to consider the different tautomers accessible at a given energy and obtain final theoretical spectra by insertion of an initial tautomer. In the work reported here, the most stable tautomer was chosen as the initial structure, but in general, any structure could be considered. Allowing for the formation of the various possible tautomers during simulation calculations was found to be important to getting a more complete fragmentation spectrum. The calculated theoretical results reproduce the experimental peaks such that it was possible to determine reaction pathways and product structures. The calculated tautomerization network was crucial to correctly identifying all the observed ion peaks, showing that a mobile proton model holds not only for peptide fragmentation but also for nucleobases. Finally, first principles results were compared to simple machine learning fragmentation models.

Chemical dynamics simulations of CID of peptide ions: comparisons between TIK(H+)2 and TLK(H+)2 fragmentation dynamics, and with thermal simulations

Gas phase unimolecular fragmentation of the two model doubly protonated tripeptides threonine-isoleucine-lysine (TIK) and threonine-leucine-lysine (TLK) is studied using chemical dynamics simulations. Attention is focused on different aspects of collision induced dissociation (CID): fragmentation pathways, energy transfer, theoretical mass spectra, fragmentation mechanisms, and the possibility of distinguishing isoleucine (I) and leucine (L). Furthermore, discussion is given regarding the differences between single collision CID activation, which results from a localized impact between the ions and a colliding molecule N₂, and previous thermal activation simulation results; Z. Homayoon, S. Pratihar, E. Dratz, R. Snider, R. Spezia, G. L. Barnes, V. Macaluso, A. Martin-Somer and W. L. Hase, J. Phys. Chem. A, 2016, 120, 8211-8227. Upon thermal activation unimolecular fragmentation is statistical and in accord with RRKM unimolecular rate theory. Simulations show that in collisional activation some non-statistical fragmentation occurs, including shattering, which is not present when the ions dissociate statistically. Products formed by non-statistical shattering mechanisms may be related to characteristic mass spectrometry peaks which distinguish the two isomers I and L.

Threshold for shattering fragmentation in collision-induced dissociation of the doubly protonated tripeptide TIK(H+)2

Determination of shattering threshold for unimolecular dissociation of a model tripeptide.

Unimolecular dissociation of peptides: statistical vs. non-statistical fragmentation mechanisms and time scales

In the present work we have investigated mechanisms of gas phase unimolecular dissociation of a relatively simple dipeptide, the di-proline anion, by means of chemical dynamics simulations, using the PM3 semi-empirical Hamiltonian. In particular, we have considered two activation processes that are representative limits of what occurs in collision induced dissociation experiments: (i) thermal activation, corresponding to several low energy collisions, in which the system is prepared with a microcanonical distribution of energy; (ii) collisional activation where a single shock of hundreds of kcal mol⁻¹ (300 kcal mol⁻¹ in the present case) can transfer sufficient energy to allow dissociation. From these two activation processes we obtained different product abundances, and for one particular fragmentation pathway a clear mechanistic difference for the two activation processes. This mechanism corresponds to the leaving of an OH⁻ group and subsequent formation of water by taking a proton from the remaining molecule. This last reaction is always observed in thermal activation while in collisional activation it is less favoured and the formation of OH⁻ as a final product is observed. More importantly, we show that while in thermal activation unimolecular dissociation follows exponential decay, in collision activation the initial population decays with non-exponential behaviour. Finally, from the thermal activation simulations it was possible to obtain rate constants as a function of temperature that show Arrhenius behaviour. Thus activation energies have also been extracted from these simulations.