Computing Ligands Bound to Proteins Using MELD-Accelerated MD.
Ontology highlight
ABSTRACT:
SUBMITTER: Liu C
PROVIDER: S-EPMC7572789 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:It is challenging to predict the docked conformations of two proteins. Current methods are susceptible to errors from treating proteins as rigid bodies and from an inability to compute relative Boltzmann populations of different docked conformations. Here, we show that by using the ClusPro server as a front end to generate possible protein-protein contacts, and using Modeling Employing Limited Data (MELD) accelerated molecular dynamics (MELD × MD) as a back end for atomistic simulations, we can find 16/20 native dimer structures of small proteins as those having the lowest free energy, starting from good-bound-backbone structures. We show that atomistic MD free energies can be used to identify native protein dimer structures.
| S-EPMC6690486 | biostudies-literature
Project description:A current challenge is to compute the native structures of proteins from their amino acid sequences. A main approach of bioinformatics is threading, in which a protein to be predicted is computationally threaded onto protein fragments of similar sequence having an already known structure. However, ∼15% of proteins cannot be folded in this way; this has been called the glass ceiling, and the proteins are called nonthreadables. For these, physical molecular dynamics (MD) modeling is promising because it does not require templates. We find that MD, when used with an accelerator called MELD, can fold many nonthreadables. For 41 nonthreadable proteins with fewer than 125 residues, MELD-accelerated MD (MELD × MD) folds 20 of them to better than 4 Å error. In 10 cases, MELD × MD succeeds even when the force field does not properly encode the native state. In 11 cases, MELD × MD foretells its own success; seeing large Boltzmann populations in the simulations predicts it has converged to the correct native state. MELD × MD acceleration can be applied to a broad physical protein modeling range.
| S-EPMC6705390 | biostudies-literature
Project description:MotivationFast and accurate prediction of protein-ligand binding structures is indispensable for structure-based drug design and accurate estimation of binding free energy of drug candidate molecules in drug discovery. Recently, accurate pose prediction methods based on short Molecular Dynamics (MD) simulations, such as MM-PBSA and MM-GBSA, among generated docking poses have been used. Since molecular structures obtained from MD simulation depend on the initial condition, taking the average over different initial conditions leads to better accuracy. Prediction accuracy of protein-ligand binding poses can be improved with multiple runs at different initial velocity.ResultsThis paper shows that a machine learning method, called Best Arm Identification, can optimally control the number of MD runs for each binding pose. It allows us to identify a correct binding pose with a minimum number of total runs. Our experiment using three proteins and eight inhibitors showed that the computational cost can be reduced substantially without sacrificing accuracy. This method can be applied for controlling all kinds of molecular simulations to obtain best results under restricted computational resources.Availability and implementationCode and data are available on GitHub at https://github.com/tsudalab/bpbi.Contactterayama@cbms.k.u-tokyo.ac.jp or tsuda@k.u-tokyo.ac.jp.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC6030886 | biostudies-literature
Project description:It is challenging to determine the structures of protein fibrils such as amyloids. In principle, Molecular Dynamics (MD) modeling can aid experiments, but normal MD has been impractical for these large multi-molecules. Here, we show that MELD accelerated MD (MELD x MD) can give amyloid structures from limited data. Five long-chain fibril structures are accurately predicted from NMR and Solid State NMR (SSNMR) data. Ten short-chain fibril structures are accurately predicted from more limited restraints information derived from the knowledge of strand directions. Although the present study only tests against structure predictions - which are the most detailed form of validation currently available - the main promise of this physical approach is ultimately in going beyond structures to also give mechanical properties, conformational ensembles, and relative stabilities.
| S-EPMC8342454 | biostudies-literature
Project description:PurposeTo develop generic optimization strategies for image reconstruction using graphical processing units (GPUs) in magnetic resonance imaging (MRI) and to exemplarily report on our experience with a highly accelerated implementation of the nonlinear inversion (NLINV) algorithm for dynamic MRI with high frame rates.MethodsThe NLINV algorithm is optimized and ported to run on a multi-GPU single-node server. The algorithm is mapped to multiple GPUs by decomposing the data domain along the channel dimension. Furthermore, the algorithm is decomposed along the temporal domain by relaxing a temporal regularization constraint, allowing the algorithm to work on multiple frames in parallel. Finally, an autotuning method is presented that is capable of combining different decomposition variants to achieve optimal algorithm performance in different imaging scenarios.ResultsThe algorithm is successfully ported to a multi-GPU system and allows online image reconstruction with high frame rates. Real-time reconstruction with low latency and frame rates up to 30 frames per second is demonstrated.ConclusionNovel parallel decomposition methods are presented which are applicable to many iterative algorithms for dynamic MRI. Using these methods to parallelize the NLINV algorithm on multiple GPUs, it is possible to achieve online image reconstruction with high frame rates.
| S-EPMC5804376 | biostudies-literature
Project description:The ability to compute from first principles the infrared spectrum of a protein in solution phase representing a biological system would provide a useful connection to atomistic models of protein structure and dynamics. Indeed, such calculations are a vital complement to 2DIR experimental measurements, allowing the observed signals to be interpreted in terms of detailed structural and dynamical information. In this article, we have studied nine structurally and spectroscopically well-characterized proteins, representing a range of structural types. We have simulated the equilibrium conformational dynamics in an explicit point charge water model. Using the resulting trajectories based on MD simulations, we have computed the one and two dimensional infrared spectra in the Amide I region, using an exciton approach, in which a local mode basis of carbonyl stretches is considered. The role of solvent in shifting the Amide I band (by 30 to 50 cm-1 ) is clearly evident. Similarly, the conformational dynamics contribute to the broadening of peaks in the spectrum. The inhomogeneous broadening in both the 1D and 2D spectra reflects the significant conformational diversity observed in the simulations. Through the computed 2D cross-peak spectra, we show how different pulse schemes can provide additional information on the coupled vibrations. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
| S-EPMC5434914 | biostudies-literature
Project description:The study of bound-state conformations of ligands interacting with proteins is important to the understanding of protein function and the design of drugs that alter function. Traditionally, transferred nuclear Overhauser effects (trNOEs), measured from NMR spectra of ligands in rapid exchange between bound and free states, have been used in these studies, owing to the inherent heavy weighting of bound-state data in the averaged ligand signals. In principle, residual dipolar couplings (RDCs) provide a useful complement to NOE data in that they provide orientational constraints as opposed to distance constraints, but use in ligand-binding applications has been limited due to the absence of heavy weighting of bound-state data. A widely applicable approach to increasing the weighting of bound-state data in averaged RDCs measured on ligands is presented. The approach rests on association of a His-tagged protein with a nickel-chelate-carrying lipid inserted into the lipid bilayer-like alignment media used in the acquisition of RDCs. The approach is validated through the observation of bound-state RDCs for the disaccharide, lactose, bound to the carbohydrate recognition domain of the mammalian lectin, galectin-3.
| S-EPMC2542485 | biostudies-literature
Project description:This paper describes the use of magnetic levitation (MagLev) to measure the association of proteins and ligands. The method starts with diamagnetic gel beads that are functionalized covalently with small molecules (putative ligands). Binding of protein to the ligands within the bead causes a change in the density of the bead. When these beads are suspended in a paramagnetic aqueous buffer and placed between the poles of two NbFeB magnets with like poles facing, the changes in the density of the bead on binding of protein result in changes in the levitation height of the bead that can be used to quantify the amount of protein bound. This paper uses a reaction-diffusion model to examine the physical principles that determine the values of rate and equilibrium constants measured by this system, using the well-defined model system of carbonic anhydrase and aryl sulfonamides. By tuning the experimental protocol, the method is capable of quantifying either the concentration of protein in a solution, or the binding affinities of a protein to several resin-bound small molecules simultaneously. Since this method requires no electricity and only a single piece of inexpensive equipment, it may find use in situations where portability and low cost are important, such as in bioanalysis in resource-limited settings, point-of-care diagnosis, veterinary medicine, and plant pathology. It still has several practical disadvantages. Most notably, the method requires relatively long assay times and cannot be applied to large proteins (>70 kDa), including antibodies. The design and synthesis of beads with improved characteristics (e.g., larger pore size) has the potential to resolve these problems.
| S-EPMC3319098 | biostudies-literature
Project description:Here, we demonstrate a method to capture local dynamics on a time scale 3 orders of magnitude beyond state-of-the-art simulation approaches. We apply accelerated molecular dynamics simulations for conformational sampling and extract reweighted backbone dihedral distributions. Local dynamics are characterized by torsional probabilities, resulting in residue-wise dihedral entropies. Our approach is successfully validated for three different protein systems of increasing size: alanine dipeptide, bovine pancreatic trypsin inhibitor (BPTI), and the major birch pollen allergen Bet v 1a. We demonstrate excellent agreement of flexibility profiles with both large-scale computer simulations and NMR experiments. Thus, our method provides efficient access to local protein dynamics on extended time scales of high biological relevance.
| S-EPMC4980685 | biostudies-literature
Project description:Ligands bound to protein assemblies provide critical information for function, yet are often difficult to capture and define. Here we develop a top-down method, 'nativeomics', unifying 'omics' (lipidomics, proteomics, metabolomics) analysis with native mass spectrometry to identify ligands bound to membrane protein assemblies. By maintaining the link between proteins and ligands, we define the lipidome/metabolome in contact with membrane porins and a mitochondrial translocator to discover potential regulators of protein function.
| S-EPMC7332344 | biostudies-literature