We describe a novel approach of reaction representation as a combination of two mixtures: a mixture of reactants and a mixture of products. In turn, each mixture can be encoded using an earlier reported approach involving simplex descriptors (SiRMS). The feature vector representing these two mixtures results from either concatenated product and reactant descriptors or the difference between descriptors of products and reactants. This reaction representation doesn't need an explicit labeling of a reaction center. The rigorous "product-out" cross-validation (CV) strategy has been suggested. Unlike the naïve "reaction-out" CV approach based on a random selection of items, the proposed one provides with more realistic estimation of prediction accuracy for reactions resulting in novel products. The new methodology has been applied to model rate constants of E2 reactions. It has been demonstrated that the use of the fragment control domain applicability approach significantly increases prediction accuracy of the models. The models obtained with new "mixture" approach performed better than those required either explicit (Condensed Graph of Reaction) or implicit (reaction fingerprints) reaction center labeling.
This article focuses on designing mutations of the PA-IIL lectin from Pseudomonas aeruginosa that lead to change in specificity. Following the previous results revealing the importance of the amino acid triad 22-23-24 (so-called specificity-binding loop), saturation in silico mutagenesis was performed, with the intent of finding mutations that increase the lectin's affinity and modify its specificity. For that purpose, a combination of docking, molecular dynamics and binding free energy calculation was used. The combination of methods revealed mutations that changed the performance of the wild-type lectin and its mutants to their preferred partners. The mutation at position 22 resulted in 85% in inactivation of the binding site, and the mutation at 23 did not have strong effects thanks to the side chain being pointed away from the binding site. Molecular dynamics simulations followed by binding free energy calculation were performed on mutants with promising results from docking, and also at those where the amino acid at position 24 was replaced for bulkier or longer polar chain. The key mutants were also prepared in vitro and their binding properties determined by isothermal titration calorimetry. Combination of the used methods proved to be able to predict changes in the lectin performance and helped in explaining the data observed experimentally.
- MeSH
- bakteriální adheziny chemie genetika metabolismus MeSH
- design s pomocí počítače MeSH
- lektiny chemie genetika metabolismus MeSH
- metabolismus sacharidů MeSH
- mutace MeSH
- mutageneze * MeSH
- počítačová simulace MeSH
- Pseudomonas aeruginosa chemie genetika metabolismus MeSH
- simulace molekulární dynamiky MeSH
- simulace molekulového dockingu MeSH
- termodynamika MeSH
- vazebná místa MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Twelve homology models of the human M2 muscarinic receptor using different sets of templates have been designed using the Prime program or the modeller program and compared to crystallographic structure (PDB:3UON). The best models were obtained using single template of the closest published structure, the M3 muscarinic receptor (PDB:4DAJ). Adding more (structurally distant) templates led to worse models. Data document a key role of the template in homology modeling. The models differ substantially. The quality checks built into the programs do not correlate with the RMSDs to the crystallographic structure and cannot be used to select the best model. Re-docking of the antagonists present in crystallographic structure and relative binding energy estimation by calculating MM/GBSA in Prime and the binding energy function in YASARA suggested it could be possible to evaluate the quality of the orthosteric binding site based on the prediction of relative binding energies. Although estimation of relative binding energies distinguishes between relatively good and bad models it does not indicate the best one. On the other hand, visual inspection of the models for known features and knowledge-based analysis of the intramolecular interactions allows an experimenter to select overall best models manually.
- MeSH
- konformace proteinů * MeSH
- krystalografie rentgenová * MeSH
- lidé MeSH
- molekulární modely MeSH
- receptor muskarinový M2 chemie MeSH
- sekvence aminokyselin MeSH
- sekvenční homologie aminokyselin MeSH
- simulace molekulového dockingu MeSH
- vazebná místa MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Molecular dynamics simulations of complexes between Norwalk virus RNA dependent RNA polymerase and its natural CTP and 2dCTP (both containing the O5'-C5'-C4'-O4' sequence of atoms bridging the triphosphate and sugar moiety) or modified coCTP (C5'-O5'-C4'-O4'), cocCTP (C5'-O5'-C4'-C4'') substrates were produced by means of CUDA programmable graphical processing units and the ACEMD software package. It enabled us to gain microsecond MD trajectories clearly showing that similar nucleoside triphosphates can bind surprisingly differently into the active site of the Norwalk virus RNA dependent RNA polymerase. It corresponds to their different modes of action (CTP-substrate, 2dCTP-poor substrate, coCTP-chain terminator, cocCTP-inhibitor). Moreover, extremely rare events-as repetitive pervasion of Arg182 into a potentially reaction promoting arrangement-were captured.
- MeSH
- cytidintrifosfát analogy a deriváty metabolismus MeSH
- infekce viry z čeledi Caliciviridae virologie MeSH
- lidé MeSH
- Norovirus enzymologie metabolismus MeSH
- RNA-dependentní RNA-polymerasa metabolismus MeSH
- simulace molekulární dynamiky MeSH
- simulace molekulového dockingu MeSH
- substrátová specifita MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
A semiempirical quantum mechanical PM6-DH2 method accurately covering the dispersion interaction and H-bonding was used to score fifteen structurally diverse CDK2 inhibitors. The geometries of all the complexes were taken from the X-ray structures and were reoptimised by the PM6-DH2 method in continuum water. The total scoring function was constructed as an estimate of the binding free energy, i.e., as a sum of the interaction enthalpy, interaction entropy and the corrections for the inhibitor desolvation and deformation energies. The applied scoring function contains a clear thermodynamical terms and does not involve any adjustable empirical parameter. The best correlations with the experimental inhibition constants (ln K (i)) were found for bare interaction enthalpy (r (2) = 0.87) and interaction enthalpy corrected for ligand desolvation and deformation energies (r (2) = 0.77); when the entropic term was considered, however, the correlation becomes worse but still acceptable (r (2) = 0.52). The resulting correlation based on the PM6-DH2 scoring function is better than previously published function based on various docking/scoring, SAR studies or advanced QM/MM approach, however, the robustness is limited by number of available experimental data used in the correlation. Since a very similar correlation between the experimental and theoretical results was found also for a different system of the HIV-1 protease, the suggested scoring function based on the PM6-DH2 method seems to be applicable in drug design, even if diverse protein-ligand complexes have to be ranked.
- MeSH
- cyklin-dependentní kinasa 2 antagonisté a inhibitory metabolismus MeSH
- inhibitory proteinkinas chemie farmakologie MeSH
- kvantová teorie MeSH
- lidé MeSH
- ligandy MeSH
- molekulární modely MeSH
- racionální návrh léčiv MeSH
- termodynamika MeSH
- vazba proteinů MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
