A new paradigm for molecular dynamics databases: the COVID-19 database, the legacy of a titanic community effort.

Molecular dynamics (MD) simulations are keeping computers busy around the world, generating a huge amount of data that is typically not open to the scientific community. Pioneering efforts to ensure the safety and reusability of MD data have been based on the use of simple databases providing a limited set of standard analyses on single-short trajectories. Despite their value, these databases do not offer a true solution for the current community of MD users, who want a flexible analysis pipeline and the possibility to address huge non-Markovian ensembles of large systems. Here we present a new paradigm for MD databases, resilient to large systems and long trajectories, and designed to be compatible with modern MD simulations. The data are offered to the community through a web-based graphical user interface (GUI), implemented with state-of-the-art technology, which incorporates system-specific analysis designed by the trajectory providers. A REST API and associated Jupyter Notebooks are integrated into the platform, allowing fully customized meta-analysis by final users. The new technology is illustrated using a collection of trajectories obtained by the community in the context of the effort to fight the COVID-19 pandemic. The server is accessible at https://bioexcel-cv19.bsc.es/#/. It is free and open to all users and there are no login requirements. It is also integrated into the simulations section of the BioExcel-MolSSI COVID-19 Molecular Structure and Therapeutics Hub: https://covid.molssi.org/simulations/ and is part of the MDDB effort (https://mddbr.eu).

Using interactive Jupyter Notebooks and BioConda for FAIR and reproducible biomolecular simulation workflows.

Interactive Jupyter Notebooks in combination with Conda environments can be used to generate FAIR (Findable, Accessible, Interoperable and Reusable/Reproducible) biomolecular simulation workflows. The interactive programming code accompanied by documentation and the possibility to inspect intermediate results with versatile graphical charts and data visualization is very helpful, especially in iterative processes, where parameters might be adjusted to a particular system of interest. This work presents a collection of FAIR notebooks covering various areas of the biomolecular simulation field, such as molecular dynamics (MD), protein-ligand docking, molecular checking/modeling, molecular interactions, and free energy perturbations. Workflows can be launched with myBinder or easily installed in a local system. The collection of notebooks aims to provide a compilation of demonstration workflows, and it is continuously updated and expanded with examples using new methodologies and tools.

BioExcel Building Blocks Workflows (BioBB-Wfs), an integrated web-based platform for biomolecular simulations.

We present BioExcel Building Blocks Workflows, a web-based graphical user interface (GUI) offering access to a collection of transversal pre-configured biomolecular simulation workflows assembled with the BioExcel Building Blocks library. Available workflows include Molecular Dynamics setup, protein-ligand docking, trajectory analyses and small molecule parameterization. Workflows can be launched in the platform or downloaded to be run in the users' own premises. Remote launching of long executions to user's available High-Performance computers is possible, only requiring configuration of the appropriate access credentials. The web-based graphical user interface offers a high level of interactivity, with integration with the NGL viewer to visualize and check 3D structures, MDsrv to visualize trajectories, and Plotly to explore 2D plots. The server requires no login but is recommended to store the users' projects and manage sensitive information such as remote credentials. Private projects can be made public and shared with colleagues with a simple URL. The tool will help biomolecular simulation users with the most common and repetitive processes by means of a very intuitive and interactive graphical user interface. The server is accessible at https://mmb.irbbarcelona.org/biobb-wfs.

BioExcel Building Blocks REST API (BioBB REST API), programmatic access to interoperable biomolecular simulation tools.

MOTIVATION: The BioExcel Building Blocks (BioBB) library offers a broad collection of wrappers on top of common biomolecular simulation and bioinformatics tools. The possibility to access the library remotely and programmatically increases its usability, allowing individual and sporadic executions and enabling remote workflows. RESULTS: BioBB REST API extends and complements the BioBB library offering programmatic access to the collection of biomolecular simulation tools included in the BioExcel Building Blocks library. Molecular Dynamics setup, docking, structure modeling, free energy simulations and flexibility analyses are examples of functionalities included in the endpoints collection. All functionalities are accessible through standard REST API calls, voiding the need for tool installation. AVAILABILITY AND IMPLEMENTATION: All the information related to the BioBB REST API endpoints is accessible from https://mmb.irbbarcelona.org/biobb-api/. Links to extended documentation, including OpenAPI endpoints specification and examples, Read-The-Docs documentation and a complete workflow tutorial can be found in the Supplementary Table S1. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Sequence-Dependent Properties of the RNA Duplex.

Sequence-dependent properties of the DNA duplex have been accurately described using extensive molecular dynamics simulations. The RNA duplex meanwhile─which is typically represented as a sequence-averaged rigid rod─does not benefit from having equivalent molecular dynamics simulations. In this paper, we present a massive simulation effort using a set of ABC-optimized duplexes from which we derived tetramer-resolution properties of the RNA duplex and a simple mesoscopic model that can represent elastic properties of long RNA duplexes. Despite the extreme chemical similarity between DNA and RNA, the local and global elastic properties of the duplexes are very different. DNA duplexes show a complex and nonelastic pattern of flexibility, for instance, while RNA duplexes behave as an elastic system whose deformations can be represented by simple harmonic potentials. In RNA duplexes (RNA2), not only are intra- and interbase pair parameters (equilibrium and mechanical) different from those in the equivalent DNA duplex sequences (DNA2) but the correlations between movements also differ. Simple statements on the relative flexibility or stability of both polymers are meaningless and should be substituted by a more detailed description depending on the sequence and the type of deformation considered.

High-Throughput Prediction of the Impact of Genetic Variability on Drug Sensitivity and Resistance Patterns for Clinically Relevant Epidermal Growth Factor Receptor Mutations from Atomistic Simulations.

Mutations in the kinase domain of the epidermal growth factor receptor (EGFR) can be drivers of cancer and also trigger drug resistance in patients receiving chemotherapy treatment based on kinase inhibitors. A priori knowledge of the impact of EGFR variants on drug sensitivity would help to optimize chemotherapy and design new drugs that are effective against resistant variants before they emerge in clinical trials. To this end, we explored a variety of in silico methods, from sequence-based to "state-of-the-art" atomistic simulations. We did not find any sequence signal that can provide clues on when a drug-related mutation appears or the impact of such mutations on drug activity. Low-level simulation methods provide limited qualitative information on regions where mutations are likely to cause alterations in drug activity, and they can predict around 70% of the impact of mutations on drug efficiency. High-level simulations based on nonequilibrium alchemical free energy calculations show predictive power. The integration of these "state-of-the-art" methods into a workflow implementing an interface for parallel distribution of the calculations allows its automatic and high-throughput use, even for researchers with moderate experience in molecular simulations.

A multi-modal coarse grained model of DNA flexibility mappable to the atomistic level.

We present a new coarse grained method for the simulation of duplex DNA. The algorithm uses a generalized multi-harmonic model that can represent any multi-normal distribution of helical parameters, thus avoiding caveats of current mesoscopic models for DNA simulation and representing a breakthrough in the field. The method has been parameterized from accurate parmbsc1 atomistic molecular dynamics simulations of all unique tetranucleotide sequences of DNA embedded in long duplexes and takes advantage of the correlation between helical states and backbone configurations to derive atomistic representations of DNA. The algorithm, which is implemented in a simple web interface and in a standalone package reproduces with high computational efficiency the structural landscape of long segments of DNA untreatable by atomistic molecular dynamics simulations.

Oncogenic mutations at the EGFR ectodomain structurally converge to remove a steric hindrance on a kinase-coupled cryptic epitope.

Epidermal growth factor receptor (EGFR) signaling is initiated by a large ligand-favored conformational change of the extracellular domain (ECD) from a closed, self-inhibited tethered monomer, to an open untethered state, which exposes a loop required for strong dimerization and activation. In glioblastomas (GBMs), structurally heterogeneous missense and deletion mutations concentrate at the ECD for unclear reasons. We explore the conformational impact of GBM missense mutations, combining elastic network models (ENMs) with multiple molecular dynamics (MD) trajectories. Our simulations reveal that the main missense class, located at the I-II interface away from the self-inhibitory tether, can unexpectedly favor spontaneous untethering to a compact intermediate state, here validated by small-angle X-ray scattering (SAXS). Significantly, such intermediate is characterized by the rotation of a large ECD fragment (N-TR1), deleted in the most common GBM mutation, EGFRvIII, and that makes accessible a cryptic epitope characteristic of cancer cells. This observation suggested potential structural equivalence of missense and deletion ECD changes in GBMs. Corroborating this hypothesis, our FACS, in vitro, and in vivo data demonstrate that entirely different ECD variants all converge to remove N-TR1 steric hindrance from the 806-epitope, which we show is allosterically coupled to an intermediate kinase and hallmarks increased oncogenicity. Finally, the detected extraintracellular coupling allows for synergistic cotargeting of the intermediate with mAb806 and inhibitors, which is proved herein.

Molywood: streamlining the design and rendering of molecular movies.

MOTIVATION: High-quality dynamic visuals are needed at all levels of science communication, from the conference hall to the classroom. As scientific journals embrace new article formats, many key concepts-particularly, in structural biology-are also more easily conveyed as videos than still frames. Notwithstanding, the design and rendering of a complex molecular movie remain an arduous task. Here, we introduce Molywood, a robust and intuitive tool that builds on the capabilities of Visual Molecular Dynamics (VMD) to automate all stages of movie rendering. RESULTS: Molywood is a Python-based script that uses an integrated workflow to give maximal flexibility in movie design. It implements the basic concepts of actions, layers, grids and concurrency and requires no programming experience to run. AVAILABILITY AND IMPLEMENTATION: The script is freely available on GitLab (gitlab.com/KomBioMol/molywood) and PyPI (through pip), and features an extended documentation, tutorial and gallery hosted on mmb.irbbarcelona.org/molywood.

Modulation of the helical properties of DNA: next-to-nearest neighbour effects and beyond.

We used extensive molecular dynamics simulations to study the structural and dynamic properties of the central d(TpA) step in the highly polymorphic d(CpTpApG) tetranucleotide. Contrary to the assumption of the dinucleotide-model and its nearest neighbours (tetranucleotide-model), the properties of the central d(TpA) step change quite significantly dependent on the next-to-nearest (hexanucleotide) sequence context and in a few cases are modulated by even remote neighbours (beyond next-to-nearest from the central TpA). Our results highlight the existence of previously undescribed dynamical mechanisms for the transmission of structural information into the DNA and demonstrate the existence of certain sequences with special physical properties that can impact on the global DNA structure and dynamics.

Nucleosome Dynamics: a new tool for the dynamic analysis of nucleosome positioning.

We present Nucleosome Dynamics, a suite of programs integrated into a virtual research environment and created to define nucleosome architecture and dynamics from noisy experimental data. The package allows both the definition of nucleosome architectures and the detection of changes in nucleosomal organization due to changes in cellular conditions. Results are displayed in the context of genomic information thanks to different visualizers and browsers, allowing the user a holistic, multidimensional view of the genome/transcriptome. The package shows good performance for both locating equilibrium nucleosome architecture and nucleosome dynamics and provides abundant useful information in several test cases, where experimental data on nucleosome position (and for some cases expression level) have been collected for cells under different external conditions (cell cycle phase, yeast metabolic cycle progression, changes in nutrients or difference in MNase digestion level). Nucleosome Dynamics is a free software and is provided under several distribution models.

The static and dynamic structural heterogeneities of B-DNA: extending Calladine-Dickerson rules.

We present a multi-laboratory effort to describe the structural and dynamical properties of duplex B-DNA under physiological conditions. By processing a large amount of atomistic molecular dynamics simulations, we determine the sequence-dependent structural properties of DNA as expressed in the equilibrium distribution of its stochastic dynamics. Our analysis includes a study of first and second moments of the equilibrium distribution, which can be accurately captured by a harmonic model, but with nonlocal sequence-dependence. We characterize the sequence-dependent choreography of backbone and base movements modulating the non-Gaussian or anharmonic effects manifested in the higher moments of the dynamics of the duplex when sampling the equilibrium distribution. Contrary to prior assumptions, such anharmonic deformations are not rare in DNA and can play a significant role in determining DNA conformation within complexes. Polymorphisms in helical geometries are particularly prevalent for certain tetranucleotide sequence contexts and are always coupled to a complex network of coordinated changes in the backbone. The analysis of our simulations, which contain instances of all tetranucleotide sequences, allow us to extend Calladine-Dickerson rules used for decades to interpret the average geometry of DNA, leading to a set of rules with quantitative predictive power that encompass nonlocal sequence-dependence and anharmonic fluctuations.

Parmbsc1: a refined force field for DNA simulations.

We present parmbsc1, a force field for DNA atomistic simulation, which has been parameterized from high-level quantum mechanical data and tested for nearly 100 systems (representing a total simulation time of ∼ 140 µs) covering most of DNA structural space. Parmbsc1 provides high-quality results in diverse systems. Parameters and trajectories are available at http://mmb.irbbarcelona.org/ParmBSC1/.

How accurate are accurate force-fields for B-DNA?

Last generation of force-fields are raising expectations on the quality of molecular dynamics (MD) simulations of DNA, as well as to the belief that theoretical models can substitute experimental ones in several cases. However these claims are based on limited benchmarks, where MD simulations have shown the ability to reproduce already existing 'experimental models', which in turn, have an unclear accuracy to represent DNA conformation in solution. In this work we explore the ability of different force-fields to predict the structure of two new B-DNA dodecamers, determined herein by means of 1H nuclear magnetic resonance (NMR). The study allowed us to check directly for experimental NMR observables on duplexes previously not solved, and also to assess the reliability of 'experimental structures'. We observed that technical details in the annealing procedures can induce non-negligible local changes in the final structures. We also found that while not all theoretical simulations are equally reliable, those obtained using last generation of AMBER force-fields (BSC1 and BSC0OL15) show predictive power in the multi-microsecond timescale and can be safely used to reproduce global structure of DNA duplexes and fine sequence-dependent details.

DNA structure directs positioning of the mitochondrial genome packaging protein Abf2p.

The mitochondrial genome (mtDNA) is assembled into nucleo-protein structures termed nucleoids and maintained differently compared to nuclear DNA, the involved molecular basis remaining poorly understood. In yeast (Saccharomyces cerevisiae), mtDNA is a ∼80 kbp linear molecule and Abf2p, a double HMG-box protein, packages and maintains it. The protein binds DNA in a non-sequence-specific manner, but displays a distinct 'phased-binding' at specific DNA sequences containing poly-adenine tracts (A-tracts). We present here two crystal structures of Abf2p in complex with mtDNA-derived fragments bearing A-tracts. Each HMG-box of Abf2p induces a 90° bend in the contacted DNA, causing an overall U-turn. Together with previous data, this suggests that U-turn formation is the universal mechanism underlying mtDNA compaction induced by HMG-box proteins. Combining this structural information with mutational, biophysical and computational analyses, we reveal a unique DNA binding mechanism for Abf2p where a characteristic N-terminal flag and helix are crucial for mtDNA maintenance. Additionally, we provide the molecular basis for A-tract mediated exclusion of Abf2p binding. Due to high prevalence of A-tracts in yeast mtDNA, this has critical relevance for nucleoid architecture. Therefore, an unprecedented A-tract mediated protein positioning mechanism regulates DNA packaging proteins in the mitochondria, and in combination with DNA-bending and U-turn formation, governs mtDNA compaction.

Predicting the Limit of Intramolecular Hydrogen Bonding with Classical Molecular Dynamics.

The energetics of intramolecular recognition processes are governed by the balance of pre-organization and flexibility, which is often difficult to measure and hard to predict. Using classical MD simulations, we predict and quantify the effective strength of intramolecular hydrogen bonds between donor and acceptor sites separated by a variable alkyl linker in several solvents and crowded solutions. The balance of entropic and enthalpic contributions poses a solvent-dependent limit to the occurrence of intramolecular H-bonding. Still, free energies show a constant offset among different solvents with, for example, a 13 kJ mol-1 difference between water and chloroform. Molecular crowding shows little effect on the thermodynamic equilibrium, but induces variations on the H-bond kinetics. The results are in quantitative agreement with experiments in chloroform and showcase a general strategy to investigate molecular interactions in different environments, extending the limits of current experiments towards the prospective prediction of H-bond interactions in a variety of contexts.

Long-timescale dynamics of the Drew-Dickerson dodecamer.

We present a systematic study of the long-timescale dynamics of the Drew-Dickerson dodecamer (DDD: d(CGCGAATTGCGC)2) a prototypical B-DNA duplex. Using our newly parameterized PARMBSC1 force field, we describe the conformational landscape of DDD in a variety of ionic environments from minimal salt to 2 M Na(+)Cl(-) or K(+)Cl(-) The sensitivity of the simulations to the use of different solvent and ion models is analyzed in detail using multi-microsecond simulations. Finally, an extended (10 µs) simulation is used to characterize slow and infrequent conformational changes in DDD, leading to the identification of previously uncharacterized conformational states of this duplex which can explain biologically relevant conformational transitions. With a total of more than 43 µs of unrestrained molecular dynamics simulation, this study is the most extensive investigation of the dynamics of the most prototypical DNA duplex.

BIGNASim: a NoSQL database structure and analysis portal for nucleic acids simulation data.

Molecular dynamics simulation (MD) is, just behind genomics, the bioinformatics tool that generates the largest amounts of data, and that is using the largest amount of CPU time in supercomputing centres. MD trajectories are obtained after months of calculations, analysed in situ, and in practice forgotten. Several projects to generate stable trajectory databases have been developed for proteins, but no equivalence exists in the nucleic acids world. We present here a novel database system to store MD trajectories and analyses of nucleic acids. The initial data set available consists mainly of the benchmark of the new molecular dynamics force-field, parmBSC1. It contains 156 simulations, with over 120 µs of total simulation time. A deposition protocol is available to accept the submission of new trajectory data. The database is based on the combination of two NoSQL engines, Cassandra for storing trajectories and MongoDB to store analysis results and simulation metadata. The analyses available include backbone geometries, helical analysis, NMR observables and a variety of mechanical analyses. Individual trajectories and combined meta-trajectories can be downloaded from the portal. The system is accessible through http://mmb.irbbarcelona.org/BIGNASim/. Supplementary Material is also available on-line at http://mmb.irbbarcelona.org/BIGNASim/SuppMaterial/.

Small Details Matter: The 2'-Hydroxyl as a Conformational Switch in RNA.

While DNA is mostly a primary carrier of genetic information and displays a regular duplex structure, RNA can form very complicated and conserved 3D structures displaying a large variety of functions, such as being an intermediary carrier of the genetic information, translating such information into the protein machinery of the cell, or even acting as a chemical catalyst. At the base of such functional diversity is the subtle balance between different backbone, nucleobase, and ribose conformations, finely regulated by the combination of hydrogen bonds and stacking interactions. Although an apparently simple chemical modification, the presence of the 2'OH in RNA has a profound effect in the ribonucleotide conformational balance, adding an extra layer of complexity to the interactions network in RNA. In the present work, we have combined database analysis with extensive molecular dynamics, quantum mechanics, and hybrid QM/MM simulations to provide direct evidence on the dramatic impact of the 2'OH conformation on sugar puckering. Calculations provide evidence that proteins can modulate the 2'OH conformation to drive sugar repuckering, leading then to the formation of bioactive conformations. In summary, the 2'OH group seems to be a primary molecular switch contributing to specific protein-RNA recognition.

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.