pH-Dependent Capping Interactions Induce Large-Scale Structural Transitions in i-Motifs.

We study here a DNA oligonucleotide having the ability to form two different i-motif structures whose relative stability depends on pH and temperature. The major species at neutral pH is stabilized by two C:C+ base pairs capped by two minor groove G:C:G:C tetrads. The high pH and thermal stability of this structure are mainly due to the favorable effect of the minor groove tetrads on their adjacent positively charged C:C+ base pairs. At pH 5, we observe a more elongated i-motif structure consisting of four C:C+ base pairs capped by two G:T:G:T tetrads. Molecular dynamics calculations show that the conformational transition between the two structures is driven by the protonation state of key cytosines. In spite of large conformational differences, the transition between the acidic and neutral structures can occur without unfolding of the i-motif. These results represent the first case of a conformational switch between two different i-motif structures and illustrate the dramatic pH-dependent plasticity of this fascinating DNA motif.

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.

Asymmetric base-pair opening drives helicase unwinding dynamics.

The opening of a Watson-Crick double helix is required for crucial cellular processes, including replication, repair, and transcription. It has long been assumed that RNA or DNA base pairs are broken by the concerted symmetric movement of complementary nucleobases. By analyzing thousands of base-pair opening and closing events from molecular simulations, here, we uncover a systematic stepwise process driven by the asymmetric flipping-out probability of paired nucleobases. We demonstrate experimentally that such asymmetry strongly biases the unwinding efficiency of DNA helicases toward substrates that bear highly dynamic nucleobases, such as pyrimidines, on the displaced strand. Duplex substrates with identical thermodynamic stability are thus shown to be more easily unwound from one side than the other, in a quantifiable and predictable manner. Our results indicate a possible layer of gene regulation coded in the direction-dependent unwindability of the double helix.

Plasticity in oligomerization, operator architecture, and DNA binding in the mode of action of a bacterial B12-based photoreceptor.

Newly discovered bacterial photoreceptors called CarH sense light by using 5'-deoxyadenosylcobalamin (AdoCbl). They repress their own expression and that of genes for carotenoid synthesis by binding in the dark to operator DNA as AdoCbl-bound tetramers, whose light-induced disassembly relieves repression. High-resolution structures of Thermus thermophilus CarHTt have provided snapshots of the dark and light states and have revealed a unique DNA-binding mode whereby only three of four DNA-binding domains contact an operator comprising three tandem direct repeats. To gain further insights into CarH photoreceptors and employing biochemical, spectroscopic, mutational, and computational analyses, here we investigated CarHBm from Bacillus megaterium We found that apoCarHBm, unlike monomeric apoCarHTt, is an oligomeric molten globule that forms DNA-binding tetramers in the dark only upon AdoCbl binding, which requires a conserved W-X9-EH motif. Light relieved DNA binding by disrupting CarHBm tetramers to dimers, rather than to monomers as with CarHTt CarHBm operators resembled that of CarHTt, but were larger by one repeat and overlapped with the -35 or -10 promoter elements. This design persisted in a six-repeat, multipartite operator we discovered upstream of a gene encoding an Spx global redox-response regulator whose photoregulated expression links photooxidative and general redox responses in B. megaterium Interestingly, CarHBm recognized the smaller CarHTt operator, revealing an adaptability possibly related to the linker bridging the DNA- and AdoCbl-binding domains. Our findings highlight a remarkable plasticity in the mode of action of B12-based CarH photoreceptors, important for their biological functions and development as optogenetic tools.

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.

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.

Characterization of p38α autophosphorylation inhibitors that target the non-canonical activation pathway.

p38α is a versatile protein kinase that can control numerous processes and plays important roles in the cellular responses to stress. Dysregulation of p38α signaling has been linked to several diseases including inflammation, immune disorders and cancer, suggesting that targeting p38α could be therapeutically beneficial. Over the last two decades, numerous p38α inhibitors have been developed, which showed promising effects in pre-clinical studies but results from clinical trials have been disappointing, fueling the interest in the generation of alternative mechanisms of p38α modulation. Here, we report the in silico identification of compounds that we refer to as non-canonical p38α inhibitors (NC-p38i). By combining biochemical and structural analyses, we show that NC-p38i efficiently inhibit p38α autophosphorylation but weakly affect the activity of the canonical pathway. Our results demonstrate how the structural plasticity of p38α can be leveraged to develop therapeutic opportunities targeting a subset of the functions regulated by this pathway.

RNA unwinding from reweighted pulling simulations.

The forming and melting of complementary base pairs in RNA duplexes are conformational transitions required to accomplish a plethora of biological functions. Yet the dynamic steps of these transitions have not been quantitatively characterized at the molecular level. In this work, the base opening process was first enforced by atomistic pulling simulations and then analyzed with a novel reweighting scheme, which allowed the free-energy profile along any suitable reaction coordinate, for example, solvation, to be reconstructed. The systematic application of such approach to different base-pair combinations provides a molecular motion picture of helix opening, which is validated by comparison with an extensive set of experimental observations and links them to the enzyme-dependent unwinding mechanism. The RNA intrinsic dynamics disclosed in this work could rationalize the directionality observed in RNA-processing molecular machineries.

Probing allosteric regulations with coevolution-driven molecular simulations.

Protein-mediated allosteric regulations are essential in biology, but their quantitative characterization continues to posit formidable challenges for both experiments and computations. Here, we combine coevolutionary information, multiscale molecular simulations, and free-energy methods to interrogate and quantify the allosteric regulation of functional changes in protein complexes. We apply this approach to investigate the regulation of adenylyl cyclase (AC) by stimulatory and inhibitory G proteins­a prototypical allosteric system that has long escaped from in-depth molecular characterization. We reveal a surprisingly simple ON/OFF regulation of AC functional dynamics through multiple pathways of information transfer. The binding of G proteins reshapes the free-energy landscape of AC following the classical population-shift paradigm. The model agrees with structural and biochemical data and reveals previously unknown experimentally consistent intermediates. Our approach showcases a general strategy to explore uncharted functional space in complex biomolecular regulations.

Single-molecule pulling simulations can discern active from inactive enzyme inhibitors.

Understanding ligand-protein recognition and interaction processes is of primary importance for structure-based drug design. Traditionally, several approaches combining docking and molecular dynamics (MD) simulations have been exploited to investigate the physicochemical properties of complexes of pharmaceutical interest. Even if the geometric properties of a modeled protein-ligand complex can be well predicted by computational methods, it is challenging to rank a series of analogues in a consistent fashion with biological data. In the unique beta-hydroxyacyl-ACP dehydratase of Plasmodium falciparum (PfFabZ), the application of standard molecular docking and MD simulations was partially sufficient to shed light on the activity of previously discovered inhibitors. Complementing docking results with atomistic simulations in the steered molecular dynamics (SMD) framework, we devised an in silico approach to study molecular interactions and to compare the binding characteristics of ligand analogues. We hypothesized an interaction model that both explained the biological activity of known ligands, and provided insight into designing novel enzyme inhibitors. Mimicking single-molecule pulling experiments, we used SMD-derived force profiles to discern active from inactive compounds for the first time. A new compound was designed and its biological activity toward the PfFabZ enzyme predicted. Finally, the computational predictions were experimentally confirmed, highlighting the robustness of the drug design approach presented herein.

Exploring the Conformational Landscape of Bioactive Small Molecules.

By using a combination of classical Hamiltonian replica exchange with high-level quantum mechanical calculations on more than one hundred drug-like molecules, we explored here the energy cost associated with binding of drug-like molecules to target macromolecules. We found that, in general, the drug-like molecules present bound to proteins in the Protein Data Bank (PDB) can access easily the bioactive conformation and in fact for 73% of the studied molecules the "bioactive" conformation is within 3kBT from the most-stable conformation in solution as determined by DFT/SCRF calculations. Cases with large differences between the most-stable and the bioactive conformations appear in ligands recognized by ionic contacts, or very large structures establishing many favorable interactions with the protein. There are also a few cases where we observed a non-negligible uncertainty related to the experimental structure deposited in PDB. Remarkably, the rough automatic force field used here provides reasonable estimates of the conformational ensemble of drugs in solution. The outlined protocol can be used to better estimate the cost of adopting the bioactive conformation.

DFFR: A New Method for High-Throughput Recalibration of Automatic Force-Fields for Drugs.

We present drug force-field recalibration (DFFR), a new method for refining of automatic force-fields used to represent small drugs in docking and molecular dynamics simulations. The method is based on fine-tuning of torsional terms to obtain ensembles that reproduce observables derived from reference data. DFFR is fast and flexible and can be easily automatized for a high-throughput regime, making it useful in drug-design projects. We tested the performance of the method in a few model systems and also in a variety of druglike molecules using reference data derived from: (i) density functional theory coupled to a self-consistent reaction field (DFT/SCRF) calculations on highly populated conformers and (ii) enhanced sampling quantum mechanical/molecular mechanics (QM/MM) where the drug is reproduced at the QM level, while the solvent is represented by classical force-fields. Extension of the method to include other sources of reference data is discussed.

Bioactive Conformational Ensemble Server and Database. A Public Framework to Speed Up In Silico Drug Discovery.

Modern high-throughput structure-based drug discovery algorithms consider ligand flexibility, but typically with low accuracy, which results in a loss of performance in the derived models. Here we present the bioactive conformational ensemble (BCE) server and its associated database. The server creates conformational ensembles of drug-like ligands and stores them in the BCE database, where a variety of analyses are offered to the user. The workflow implemented in the BCE server combines enhanced sampling molecular dynamics with self-consistent reaction field quantum mechanics (SCRF/QM) calculations. The server automatizes all of the steps to transform one-dimensional (1D) or 2D representation of drugs into 3D molecules, which are then titrated, parametrized, hydrated, and optimized before being subjected to Hamiltonian replica-exchange (HREX) molecular dynamics simulations. Ensembles are collected and subjected to a clustering procedure to derive representative conformers, which are then analyzed at the SCRF/QM level of theory. All structural data are organized in a noSQL database accessible through a graphical interface and in a programmatic manner through a REST API. The server allows the user to define a private workspace and offers a deposition protocol as well as input files for "in house" calculations in those cases where confidentiality is a must. The database and the associated server are available at https://mmb.irbbarcelona.org/BCE.

Multi-target-directed coumarin derivatives: hAChE and BACE1 inhibitors as potential anti-Alzheimer compounds.

The complex etiology of Alzheimer's disease (AD) prompts scientists to develop multifunctional compounds to combat causes and symptoms of such neurodegeneration. To this aim we designed, synthesized, and tested a series of compounds by introducing halophenylalkylamidic functions on the scaffold of AP2238, which is a dual binding site acetylcholinesterase inhibitor. The inhibitory activity was successfully extended to the beta-site amyloid precursor protein cleavage enzyme, leading to the discovery of a potent inhibitor of this enzyme (3) and affording multifunctional compounds (2, 6, 8) for the treatment of AD.

Atomic-Level Characterization of the Chain-Flipping Mechanism in Fatty-Acids Biosynthesis.

During fatty acids biosynthesis the elongating acyl chain is sequestered within the core of the highly conserved acyl carrier protein (ACP). At each catalytic step, the acyl intermediates are transiently delivered from ACP to the active site of the enzymatic counterparts and, at the same time, are protected from the solvent to prevent nonselective reactivity. Yet, the molecular determinants of such a universal transition-termed chain flipping-remain poorly understood. Here we capture the atomic-level details of the chain-flipping mechanism by using metadynamics simulations. We observe the fatty-acid chain gliding through the protein-protein interface with barely 30% of its surface exposed to water molecules. The small ACP's helix III acts as gatekeeper of the process, and we find its conformational plasticity critical for a successful substrate transfer. The results are in agreement with a wide range of experimental observations and provide unprecedented insight on the molecular determinants and driving forces of the chain-flipping process.

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.

Probing riboswitch binding sites with molecular docking, focused libraries, and in-line probing assays.

Molecular docking calculations combined with chemically focused libraries can bring insight in the exploration of the structure-activity relationships for a series of related compounds against an RNA target. Yet, the in silico engine must be fueled by experimental observations to drive the research into a more effective ligand-discovery path. Here we show how molecular docking predictions can be coupled with in-line probing assays to explore the available chemical and configurational space in a riboswitch binding pocket.

Design, synthesis, and biological and crystallographic evaluation of novel inhibitors of Plasmodium falciparum enoyl-ACP-reductase (PfFabI).

Malaria, a disease of worldwide significance, is responsible for over one million deaths annually. The liver-stage of Plasmodium's life cycle is the first, obligatory, but clinically silent step in malaria infection. The P. falciparum type II fatty acid biosynthesis pathway (PfFAS-II) has been found to be essential for complete liver-stage development and has been regarded as a potential antimalarial target for the development of drugs for malaria prophylaxis and liver-stage eradication. In this paper, new coumarin-based triclosan analogues are reported and their biological profile is explored in terms of inhibitory potency against enzymes of the PfFAS-II pathway. Among the tested compounds, 7 and 8 showed the highest inhibitory potency against Pf enoyl-ACP-reductase (PfFabI), followed by 15 and 3. Finally, we determined the crystal structures of compounds 7 and 11 in complex with PfFabI to identify their mode of binding and to confirm outcomes of docking simulations.

Automated docking screens: a feasibility study.

Molecular docking is the most practical approach to leverage protein structure for ligand discovery, but the technique retains important liabilities that make it challenging to deploy on a large scale. We have therefore created an expert system, DOCK Blaster, to investigate the feasibility of full automation. The method requires a PDB code, sometimes with a ligand structure, and from that alone can launch a full screen of large libraries. A critical feature is self-assessment, which estimates the anticipated reliability of the automated screening results using pose fidelity and enrichment. Against common benchmarks, DOCK Blaster recapitulates the crystal ligand pose within 2 A rmsd 50-60% of the time; inferior to an expert, but respectrable. Half the time the ligand also ranked among the top 5% of 100 physically matched decoys chosen on the fly. Further tests were undertaken culminating in a study of 7755 eligible PDB structures. In 1398 cases, the redocked ligand ranked in the top 5% of 100 property-matched decoys while also posing within 2 A rmsd, suggesting that unsupervised prospective docking is viable. DOCK Blaster is available at http://blaster.docking.org .

Mechanical features of Plasmodium falciparum acyl carrier protein in the delivery of substrates.

The Acyl Carrier Protein (ACP) is a key element in the biosynthesis of fatty acids being responsible for the acyl group shuttling and delivery within a series of related enzymes. The molecular mechanism of the delivery process is poorly known, and its characterization is essential for in-depth understanding the biosynthetic machinery. A steered molecular dynamics approach has been applied to shed light on the putative delivery pathway, suggesting the small alpha3-helix act as gatekeeper for the transfer process. Preventing the delivery mechanism would be an innovative strategy for the development of pathway-based antimalarial compounds.