Metadynamics simulations of ligands binding to protein surfaces: a novel tool for rational drug design

Toward the Prediction of Binding Events in Very Flexible, Allosteric, Multidomain Proteins

Knowledge of the structures formed by proteins and small molecules is key to understand the molecular principles of chemotherapy and for designing new and more effective drugs. During the early stage of a drug discovery program, it is customary to predict ligand-protein complexes in silico, particularly when screening large compound databases. While virtual screening based on molecular docking is widely used for this purpose, it generally fails in mimicking binding events associated with large conformational changes in the protein, particularly when the latter involve multiple domains. In this work, we describe a new methodology to generate bound-like conformations of very flexible and allosteric proteins bearing multiple binding sites by exploiting only information on the unbound structure and the putative binding sites. The protocol is validated on the paradigm enzyme adenylate kinase, for which we generated a significant fraction of bound-like structures. A fraction of these conformations, employed in ensemble-docking calculations, allowed to find native-like poses of substrates and inhibitors (binding to the active form of the enzyme), as well as catalytically incompetent analogs (binding the inactive form). Our protocol provides a general framework for the generation of bound-like conformations of challenging drug targets that are suitable to host different ligands, demonstrating high sensitivity to the fine chemical details that regulate protein's activity. We foresee applications in virtual screening, in the prediction of the impact of amino acid mutations on structure and dynamics, and in protein engineering.

Advances in methods and concepts provide new insight into antibiotic fluxes across the bacterial membrane

The sophisticated envelope of Gram-negative bacteria modulates the uptake of small molecules in a side-chain-sensitive manner. Despite intensive theoretical and experimental investigations, a general set of pathways underpinning antibiotic uptake has not been identified. This manuscript discusses the passive influx versus active efflux of antibiotics, considering the responsible membrane proteins and the transported molecules. Recent methods have analyzed drug transport across the bacterial membrane in order to understand their activity. The combination of in vitro, in cellulo and in silico methods shed light on the key, mainly electrostatic, interactions between the molecule surface, porins and transporters during permeation. A key factor is the relationship between the dose of an active compound near its target and its antibacterial activity during the critical early window. Today, methodology breakthroughs provide fruitful tools to precisely dissect drug transport, identify key steps in drug resistance associated with membrane impermeability and efflux, and highlight key parameters to generate more effective drugs.

Predicting binding events in very flexible, allosteric, multi-domain proteins

Knowledge of the structures formed by proteins and small molecules is key to understand the molecular principles of chemotherapy and for designing new and more effective drugs. During the early stage of a drug discovery program, it is customary to predict ligand-protein complexes in silico, particularly when screening large compound databases. While virtual screening based on molecular docking is widely used for this purpose, it generally fails in mimicking binding events associated with large conformational changes in the protein, particularly when the latter involve multiple domains. In this work, we describe a new methodology to generate bound-like conformations of very flexible and allosteric proteins bearing multiple binding sites by exploiting only information on the unbound structure and the putative binding sites. The protocol is validated on the paradigm enzyme adenylate kinase, for which we generated a significant fraction of bound-like structures. A fraction of these conformations, employed in ensemble-docking calculations, allowed to find native-like poses of substrates and inhibitors (binding to the active form of the enzyme), as well as catalytically incompetent analogs (binding the inactive form). Our protocol provides a general framework for the generation of bound-like conformations of challenging drug targets that are suitable to host different ligands, demonstrating high sensitivity to the fine chemical details that regulate protein's activity. We foresee applications in virtual screening, in the prediction of the impact of amino acid mutations on structure and dynamics, and in protein engineering.

Collective Variable-Based Enhanced Sampling: From Human Learning to Machine Learning

Enhanced-sampling algorithms relying on collective variables (CVs) are extensively employed to study complex (bio)chemical processes that are not amenable to brute-force molecular simulations. The selection of appropriate CVs characterizing the slow movement modes is of paramount importance for reliable and efficient enhanced-sampling simulations. In this Perspective, we first review the application and limitations of CVs obtained from chemical and geometrical intuition. We also introduce path-sampling algorithms, which can identify path-like CVs in a high-dimensional free-energy space. Machine-learning algorithms offer a viable approach to finding suitable CVs by analyzing trajectories from preliminary simulations. We discuss both the performance of machine-learning-derived CVs in enhanced-sampling simulations of experimental models and the challenges involved in applying these CVs to realistic, complex molecular assemblies. Moreover, we provide a prospective view of the potential advancements of machine-learning algorithms for the development of CVs in the field of enhanced-sampling simulations.

Paramagnetic NMR to study iron sulfur proteins: 13C detected experiments illuminate the vicinity of the metal center

The robustness of NMR coherence transfer in proximity of a paramagnetic center depends on the relaxation properties of the nuclei involved. In the case of Iron-Sulfur Proteins, different pulse schemes or different parameter sets often provide complementary results. Tailored versions of HCACO and CACO experiments significantly increase the number of observed Cα/C' connectivities in highly paramagnetic systems, by recovering many resonances that were lost due to paramagnetic relaxation. Optimized 13C direct detected experiments can significantly extend the available assignments, improving the overall knowledge of these systems. The different relaxation properties of Cα and C' nuclei are exploited in CACO vs COCA experiments and the complementarity of the two experiments is used to obtain structural information. The two [Fe2S2]+ clusters containing NEET protein CISD3 and the one [Fe4S4]2+ cluster containing HiPIP protein PioC have been taken as model systems. We show that tailored experiments contribute to decrease the blind sphere around the cluster, to extend resonance assignment of cluster bound cysteine residues and to retrieve details on the topology of the iron-bound ligand residues.