Ignacio Faustino

Martini Coarse-Grained Force Field: Extension to DNA

We systematically parameterized a coarse-grained (CG) model for DNA that is compatible with the Martini force field. The model maps each nucleotide into six to seven CG beads and is parameterized following the Martini philosophy. The CG nonbonded interactions are based on partitioning of the nucleobases between polar and nonpolar solvents as well as base-base potential of mean force calculations. The bonded interactions are fit to single-stranded DNA (ssDNA) atomistic simulations and an elastic network is used to retain double-stranded DNA (dsDNA) and other specific DNA conformations. We present the implementation of the Martini DNA model and demonstrate the properties of individual bases, ssDNA as well as dsDNA, and DNA-protein complexes. The model opens up large-scale simulations of DNA interacting with a wide range of other (bio)molecules that are available within the Martini framework.

Exploring polymorphisms in B-DNA helical conformations

The traditional mesoscopic paradigm represents DNA as a series of base-pair steps whose energy response to equilibrium perturbations is elastic, with harmonic oscillations (defining local stiffness) around a single equilibrium conformation. In addition, base sequence effects are often analysed as a succession of independent XpY base-pair steps (i.e. a nearest-neighbour (NN) model with only 10 unique cases). Unfortunately, recent massive simulations carried out by the ABC consortium suggest that the real picture of DNA flexibility may be much more complex. The paradigm of DNA flexibility therefore needs to be revisited. In this article, we explore in detail one of the most obvious violations of the elastic NN model of flexibility: the bimodal distributions of some helical parameters. We perform here an in-depth statistical analysis of a very large set of MD trajectories and also of experimental structures, which lead to very solid evidence of bimodality. We then suggest ways to improve mesoscopic models to account for this deviation from the elastic regime.

Coupled binding mechanism of three sodium ions and aspartate in the glutamate transporter homologue GltTk

Glutamate transporters catalyse the thermodynamically unfavourable transport of anionic amino acids across the cell membrane by coupling it to the downhill transport of cations. This coupling mechanism is still poorly understood, in part because the available crystal structures of these transporters are of relatively low resolution. Here we solve crystal structures of the archaeal transporter GltTk in the presence and absence of aspartate and use molecular dynamics simulations and binding assays to show how strict coupling between the binding of three sodium ions and aspartate takes place.

NAFlex: a web server for the study of nucleic acid flexibility

We present NAFlex, a new web tool to study the flexibility of nucleic acids, either isolated or bound to other molecules. The server allows the user to incorporate structures from protein data banks, completing gaps and removing structural inconsistencies. It is also possible to define canonical (average or sequence-adapted) nucleic acid structures using a variety of predefined internal libraries, as well to create specific nucleic acid conformations from the sequence. The server offers a variety of methods to explore nucleic acid flexibility, such as a colorless wormlike-chain model, a base-pair resolution mesoscopic model and atomistic molecular dynamics simulations with a wide variety of protocols and force fields. The trajectories obtained by simulations, or imported externally, can be visualized and analyzed using a large number of tools, including standard Cartesian analysis, essential dynamics, helical analysis, local and global stiffness, energy decomposition, principal components and in silico NMR spectra. The server is accessible free of charge from the mmb.irbbarcelona.org/NAFlex webpage.

Improved nucleic acid descriptors for siRNA efficacy prediction

Although considerable progress has been made recently in understanding how gene silencing is mediated by the RNAi pathway, the rational design of effective sequences is still a challenging task. In this article, we demonstrate that including three-dimensional descriptors improved the discrimination between active and inactive small interfering RNAs (siRNAs) in a statistical model. Five descriptor types were used: (i) nucleotide position along the siRNA sequence, (ii) nucleotide composition in terms of presence/absence of specific combinations of di- and trinucleotides, (iii) nucleotide interactions by means of a modified auto- and cross-covariance function, (iv) nucleotide thermodynamic stability derived by the nearest neighbor model representation and (v) nucleic acid structure flexibility. The duplex flexibility descriptors are derived from extended molecular dynamics simulations, which are able to describe the sequence-dependent elastic properties of RNA duplexes, even for non-standard oligonucleotides. The matrix of descriptors was analysed using three statistical packages in R (partial least squares, random forest, and support vector machine), and the most predictive model was implemented in a modeling tool we have made publicly available through SourceForge. Our implementation of new RNA descriptors coupled with appropriate statistical algorithms resulted in improved model performance for the selection of siRNA candidates when compared with publicly available siRNA prediction tools and previously published test sets. Additional validation studies based on in-house RNA interference projects confirmed the robustness of the scoring procedure in prospective studies.

Functionalization of the 3′-Ends of DNA and RNA Strands with N-ethyl-N-coupled Nucleosides: A Promising Approach To Avoid 3′-Exonuclease-Catalyzed Hydrolysis of Therapeutic Oligonucleotides

The development of nucleic acid derivatives to generate novel medical treatments has become increasingly popular, but the high vulnerability of oligonucleotides to nucleases limits their practical use. We explored the possibility of increasing the stability against 3′‐exonucleases by replacing the two 3′‐terminal nucleotides by N‐ethyl‐N‐coupled nucleosides. Molecular dynamics simulations of 3′‐N‐ethyl‐N‐modified DNA:Klenow fragment complexes suggested that this kind of alteration has negative effects on the correct positioning of the adjacent scissile phosphodiester bond at the active site of the enzyme, and accordingly was expected to protect the oligonucleotide from degradation. We verified that these modifications conferred complete resistance to 3′‐exonucleases. Furthermore, cellular RNAi experiments with 3′‐N‐ethyl‐N‐modified siRNAs showed that these modifications were compatible with the RNAi machinery. Overall, our experimental and theoretical studies strongly suggest that these modified oligonucleotides could be valuable for therapeutic applications.

Photoswitching of DNA Hybridization Using a Molecular Motor

Reversible control over the functionality of biological systems via external triggers may be used in future medicine to reduce the need for invasive procedures. Additionally, externally regulated biomacromolecules are now considered as particularly attractive tools in nanoscience and the design of smart materials, due to their highly programmable nature and complex functionality. Incorporation of photoswitches into biomolecules, such as peptides, antibiotics, and nucleic acids, has generated exciting results in the past few years. Molecular motors offer the potential for new and more precise methods of photoregulation, due to their multistate switching cycle, unidirectionality of rotation, and helicity inversion during the rotational steps. Aided by computational studies, we designed and synthesized a photoswitchable DNA hairpin, in which a molecular motor serves as the bridgehead unit. After it was determined that motor function was not affected by the rigid arms of the linker, solid-phase synthesis was employed to incorporate the motor into an 8-base-pair self-complementary DNA strand. With the photoswitchable bridgehead in place, hairpin formation was unimpaired, while the motor part of this advanced biohybrid system retains excellent photochemical properties. Rotation of the motor generates large changes in structure, and as a consequence the duplex stability of the oligonucleotide could be regulated by UV light irradiation. Additionally, Molecular Dynamics computations were employed to rationalize the observed behavior of the motor–DNA hybrid. The results presented herein establish molecular motors as powerful multistate switches for application in biological environments.

cgHeliParm: analysis of dsDNA helical parameters for coarse-grained MARTINI molecular dynamics simulations

Summary We introduce cgHeliParm, a python program that provides the conformational analysis of Martini-based coarse-grained double strand DNA molecules. The software calculates the helical parameters such as base, base pair and base pair step parameters. cgHeliParm can be used for the analysis of coarse grain Martini molecular dynamics trajectories without transformation into atomistic models. Availability and implementation This package works with Python 2.7 on MacOS and Linux. The program is freely available for download from https://github.com/ifaust83/cgheliparm. Together with the main script, the base reference files CG_X_std.lib, a number of examples and R scripts are also available from the same website. A tutorial on the use and application is also available at http://cgmartini.nl/index.php/tutorials-general-introduction/tutorial-martini-dna. Contact i.faustino@rug.nl. Supplementary information Supplementary data are available at Bioinformatics online.

Targeting RNA structure in SMN2 reverses spinal muscular atrophy molecular phenotypes

Modification of SMN2 exon 7 (E7) splicing is a validated therapeutic strategy against spinal muscular atrophy (SMA). However, a target-based approach to identify small-molecule E7 splicing modifiers has not been attempted, which could reveal novel therapies with improved mechanistic insight. Here, we chose as a target the stem-loop RNA structure TSL2, which overlaps with the 5′ splicing site of E7. A small-molecule TSL2-binding compound, homocarbonyltopsentin (PK4C9), was identified that increases E7 splicing to therapeutic levels and rescues downstream molecular alterations in SMA cells. High-resolution NMR combined with molecular modelling revealed that PK4C9 binds to pentaloop conformations of TSL2 and promotes a shift to triloop conformations that display enhanced E7 splicing. Collectively, our study validates TSL2 as a target for small-molecule drug discovery in SMA, identifies a novel mechanism of action for an E7 splicing modifier, and sets a precedent for other splicing-mediated diseases where RNA structure could be similarly targeted.Spinal muscular atrophy (SMA) is an autosomal recessive disorder with no present cure. Here the authors perform an in vitro screening leading to the identification of a small molecule that alters the conformational dynamics of the TSL2 RNA structure and acts as a modulator of SMN exon 7 splicing.

Molecular Mechanism of Lipid Nanodisk Formation by Styrene-Maleic Acid Copolymers

Experimental characterization of membrane proteins often requires solubilization. A recent approach is to use styrene-maleic acid (SMA) copolymers to isolate membrane proteins in nanometer-sized membrane disks, or so-called SMA lipid particles (SMALPs). The approach has the advantage of allowing direct extraction of proteins, keeping their native lipid environment. Despite the growing popularity of using SMALPs, the molecular mechanism behind the process remains poorly understood. Here, we unravel the molecular details of the nanodisk formation by using coarse-grained molecular dynamics simulations. We show how SMA copolymers bind to the lipid bilayer interface, driven by the hydrophobic effect. Due to the concerted action of multiple adsorbed copolymers, large membrane defects appear, including small, water-filled pores. The copolymers can stabilize the rim of these pores, leading to pore growth and membrane disruption. Although complete solubilization is not seen on the timescale of our simulations, self-assembly experiments show that small nanodisks are the thermodynamically preferred end state. Our findings shed light on the mechanism of SMALP formation and on their molecular structure. This can be an important step toward the design of optimized extraction tools for membrane protein research.