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: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:Determining the rate of forming the truly folded conformation of ultrafast folding proteins is an important issue for both experiments and simulations. The double-norleucine mutant of the 35-residue villin subdomain is the focus of recent computer simulations with atomistic molecular dynamics because it is currently the fastest folding protein. The folding kinetics of this protein have been measured in laser temperature-jump experiments using tryptophan fluorescence as a probe of overall folding. The conclusion from the simulations, however, is that the rate determined by fluorescence is significantly larger than the rate of overall folding. We have therefore employed an independent experimental method to determine the folding rate. The decay of the tryptophan triplet-state in photoselection experiments was used to monitor the change in the unfolded population for a sequence of the villin subdomain with one amino acid difference from that of the laser temperature-jump experiments, but with almost identical equilibrium properties. Folding times obtained in a two-state analysis of the results from the two methods at denaturant concentrations varying from 1.5-6.0 M guanidinium chloride are in excellent agreement, with an average difference of only 20%. Polynomial extrapolation of all the data to zero denaturant yields a folding time of 220 (+100,-70) ns at 283 K, suggesting that under these conditions the barrier between folded and unfolded states has effectively disappeared--the so-called "downhill scenario."
| S-EPMC3076883 | biostudies-literature
Project description:Cellular active small molecule inhibitors of Mycobacterium tuberculosis by NMR fragment screen.
2015-01-01 | GSE17424 | GEO