HT-SuMD: making molecular dynamics simulations suitable for fragment-based screening. A comparative study with NMR.
Ontology highlight
ABSTRACT:
SUBMITTER: Ferrari F
PROVIDER: S-EPMC7598995 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Ribonucleic acid (RNA) plays a key regulatory role within the cell, cooperating with proteins to control the genome expression and several biological processes. Due to its characteristic structural features, this polymer can mold itself into different three-dimensional structures able to recognize target biomolecules with high affinity and specificity, thereby attracting the interest of drug developers and medicinal chemists. One successful example of the exploitation of RNA's structural and functional peculiarities is represented by aptamers, a class of therapeutic and diagnostic tools that can recognize and tightly bind several pharmaceutically relevant targets, ranging from small molecules to proteins, making use of the available structural and conformational freedom to maximize the complementarity with their interacting counterparts. In this scientific work, we present the first application of Supervised Molecular Dynamics (SuMD), an enhanced sampling Molecular Dynamics-based method for the study of receptor-ligand association processes in the nanoseconds timescale, to the study of recognition pathways between RNA aptamers and proteins, elucidating the main advantages and limitations of the technique while discussing its possible role in the rational design of RNA-based therapeutics.
| S-EPMC9706429 | biostudies-literature
Project description:Computational approaches supporting the early characterization of fragment molecular recognition mechanism represent a valuable complement to more expansive and low-throughput experimental techniques. In this retrospective study, we have investigated the geometric accuracy with which high-throughput supervised molecular dynamics simulations (HT-SuMD) can anticipate the experimental bound state for a set of 23 fragments targeting the SARS-CoV-2 main protease. Despite the encouraging results herein reported, in line with those previously described for other MD-based posing approaches, a high number of incorrect binding modes still complicate HT-SuMD routine application. To overcome this limitation, fragment pose stability has been investigated and integrated as part of our in-silico pipeline, allowing us to prioritize only the more reliable predictions.
| S-EPMC8250706 | biostudies-literature
Project description:Fragment-based drug discovery (FBDD) using NMR has become a central approach over the last twenty years for development of small molecule inhibitors against biological macromolecules, to control a variety of cellular processes. Yet, several considerations should be taken into account for obtaining a therapeutically relevant agent. In this review, we aim to list the considerations that make NMR fragment screening a successful process for yielding potent inhibitors. Factors that may govern the competence of NMR in fragment based drug discovery are discussed, as well as later steps that involve optimization of hits obtained by NMR-FBDD.
| S-EPMC6017141 | biostudies-literature
Project description:Fragment-based screening has evolved as a remarkable approach within the drug discovery process both in the industry and academia. Fragment screening has become a more structure-based approach to inhibitor development, but also towards development of pathway-specific clinical probes. However, it is often witnessed that the availability, immediate and long-term, of a high quality fragment-screening library is still beyond the reach of most academic laboratories. Within iNEXT (Infrastructure for NMR, EM and X-rays for Translational research), a EU-funded Horizon 2020 program, a collection of 782 fragments were assembled utilizing the concept of "poised fragments" with the aim to facilitate downstream synthesis of ligands with high affinity by fragment ligation. Herein, we describe the analytical procedure to assess the quality of this purchased and assembled fragment library by NMR spectroscopy. This quality assessment requires buffer solubility screening, comparison with LC/MS quality control and is supported by state-of-the-art software for high throughput data acquisition and on-the-fly data analysis. Results from the analysis of the library are presented as a prototype of fragment progression through the quality control process.
| S-EPMC7683495 | biostudies-literature
Project description:Methods of protein structure determination based on NMR chemical shifts are becoming increasingly common. The most widely used approaches adopt the molecular fragment replacement strategy, in which structural fragments are repeatedly reassembled into different complete conformations in molecular simulations. Although these approaches are effective in generating individual structures consistent with the chemical shift data, they do not enable the sampling of the conformational space of proteins with correct statistical weights. Here, we present a method of molecular fragment replacement that makes it possible to perform equilibrium simulations of proteins, and hence to determine their free energy landscapes. This strategy is based on the encoding of the chemical shift information in a probabilistic model in Markov chain Monte Carlo simulations. First, we demonstrate that with this approach it is possible to fold proteins to their native states starting from extended structures. Second, we show that the method satisfies the detailed balance condition and hence it can be used to carry out an equilibrium sampling from the Boltzmann distribution corresponding to the force field used in the simulations. Third, by comparing the results of simulations carried out with and without chemical shift restraints we describe quantitatively the effects that these restraints have on the free energy landscapes of proteins. Taken together, these results demonstrate that the molecular fragment replacement strategy can be used in combination with chemical shift information to characterize not only the native structures of proteins but also their conformational fluctuations.
| S-EPMC4183297 | biostudies-literature
Project description:Although proteins have represented the molecular target of choice in the development of new drug candidates, the pharmaceutical importance of ribonucleic acids has gradually been growing. The increasing availability of structural information has brought to light the existence of peculiar three-dimensional RNA arrangements, which can, contrary to initial expectations, be recognized and selectively modulated through small chemical entities or peptides. The application of classical computational methodologies, such as molecular docking, for the rational development of RNA-binding candidates is, however, complicated by the peculiarities characterizing these macromolecules, such as the marked conformational flexibility, the singular charges distribution, and the relevant role of solvent molecules. In this work, we have thus validated and extended the applicability domain of SuMD, an all-atoms molecular dynamics protocol that allows to accelerate the sampling of molecular recognition events on a nanosecond timescale, to ribonucleotide targets of pharmaceutical interest. In particular, we have proven the methodological ability by reproducing the binding mode of viral or prokaryotic ribonucleic complexes, as well as that of artificially engineered aptamers, with an impressive degree of accuracy.
| S-EPMC7057144 | biostudies-literature
Project description:Computational methods are powerful in capturing the results of experimental studies in terms of force fields that both explain and predict biological structures. Validation of molecular simulations requires comparison with experimental data to test and confirm computational predictions. Here we report a comprehensive database of NMR results for membrane phospholipids with interpretations intended to be accessible by non-NMR specialists. Experimental ¹³C-¹H and ²H NMR segmental order parameters (S(CH) or S(CD)) and spin-lattice (Zeeman) relaxation times (T(1Z)) are summarized in convenient tabular form for various saturated, unsaturated, and biological membrane phospholipids. Segmental order parameters give direct information about bilayer structural properties, including the area per lipid and volumetric hydrocarbon thickness. In addition, relaxation rates provide complementary information about molecular dynamics. Particular attention is paid to the magnetic field dependence (frequency dispersion) of the NMR relaxation rates in terms of various simplified power laws. Model-free reduction of the T(1Z) studies in terms of a power-law formalism shows that the relaxation rates for saturated phosphatidylcholines follow a single frequency-dispersive trend within the MHz regime. We show how analytical models can guide the continued development of atomistic and coarse-grained force fields. Our interpretation suggests that lipid diffusion and collective order fluctuations are implicitly governed by the viscoelastic nature of the liquid-crystalline ensemble. Collective bilayer excitations are emergent over mesoscopic length scales that fall between the molecular and bilayer dimensions, and are important for lipid organization and lipid-protein interactions. Future conceptual advances and theoretical reductions will foster understanding of biomembrane structural dynamics through a synergy of NMR measurements and molecular simulations.
| S-EPMC5176272 | biostudies-literature
Project description:The relative orientation of the two variable domains, VH and VL , influences the shape of the antigen binding site, that is, the paratope, and is essential to understand antigen specificity. ABangle characterizes the VH -VL orientation by using five angles and a distance and compares it to other known structures. Molecular dynamics simulations of antibody variable domains (Fvs) reveal fluctuations in the relative domain orientations. The observed dynamics between these domains are confirmed by NMR experiments on a single-chain variable fragment antibody (scFv) in complex with IL-1? and an antigen-binding fragment (Fab). The variability of these relative domain orientations can be interpreted as a structural feature of antibodies, which increases the antibody repertoire significantly and can enlarge the number of possible binding partners substantially. The movements of the VH and VL domains are well sampled with molecular dynamics simulations and are in agreement with the NMR ensemble. Fast Fourier transformation of the ABangle metrics allows to assign timescales of 0.1-10?GHz to the fastest collective interdomain movements. The results clearly show the necessity of dynamics to understand and characterize the favorable orientations of the VH and VL domains implying a considerable binding interface flexibility and reveal in all antibody fragments (Fab, scFv, and Fv) very similar VH -VL interdomain variations comparable to the distributions observed for known X-ray structures of antibodies. SIGNIFICANCE STATEMENT: Antibodies have become key players as therapeutic agents. The binding ability of antibodies is determined by the antigen-binding fragment (Fab), in particular the variable fragment region (Fv). Antigen-binding is mediated by the complementarity-determining regions consisting of six loops, each three of the heavy and light chain variable domain VH and VL . The relative orientation of the VH and VL domains influences the shape of the antigen-binding site and is a major objective in antibody design. In agreement with NMR experiments and molecular dynamics simulations, we show a considerable binding site flexibility in the low nanosecond timescale. Thus we suggest that this flexibility and its implications for binding and specificity should be considered when designing and optimizing therapeutic antibodies.
| S-EPMC7317758 | biostudies-literature
Project description:The TDP-43 is originally a nuclear protein but translocates to the cytoplasm in the pathological condition. TDP-43, as an RNA-binding protein, consists of two RNA Recognition Motifs (RRM1 and RRM2). RRMs are known to involve both protein-nucleotide and protein-protein interactions and mediate the formation of stress granules. Thus, they assist the entire TDP-43 protein with participating in neurodegenerative and cancer diseases. Consequently, they are potential therapeutic targets. Protein-observed and ligand-observed nuclear magnetic resonance (NMR) spectroscopy were used to uncover the small molecule inhibitors against the tandem RRM of TDP-43. We identified three hits weakly binding the tandem RRMs using the ligand-observed NMR fragment-based screening. The binding topology of these hits is then depicted by chemical shift perturbations (CSP) of the 15N-labeled tandem RRM and RRM2, respectively, and modeled by the CSP-guided High Ambiguity Driven biomolecular DOCKing (HADDOCK). These hits mainly bind to the RRM2 domain, which suggests the druggability of the RRM2 domain of TDP-43. These hits also facilitate further studies regarding the hit-to-lead evolution against the TDP-43 RRM domain.
| S-EPMC6651732 | biostudies-literature
Project description:PBRM1 is a subunit of the PBAF chromatin remodeling complex that uniquely contains six bromodomains. PBRM1 can operate as a tumor suppressor or tumor promoter. PBRM1 is a tumor promoter in prostate cancer, contributing to migratory and immunosuppressive phenotypes. Selective chemical probes targeting PBRM1 bromodomains are desired to elucidate the association between aberrant PBRM1 chromatin binding and cancer pathogenesis and the contributions of PBRM1 to immunotherapy. Previous PBRM1 inhibitors unselectively bind SMARCA2 and SMARCA4 bromodomains with nanomolar potency. We used our protein-detected NMR screening pipeline to screen 1968 fragments against the second PBRM1 bromodomain, identifying 17 hits with Kd values from 45 μM to >2 mM. Structure-activity relationship studies on the tightest-binding hit resulted in nanomolar inhibitors with selectivity for PBRM1 over SMARCA2 and SMARCA4. These chemical probes inhibit the association of full-length PBRM1 to acetylated histone peptides and selectively inhibit growth of a PBRM1-dependent prostate cancer cell line.
| S-EPMC9630929 | biostudies-literature