Machine learning (ML) methods offer a promising route to the construction of universal molecular potentials with high accuracy and low computational cost. It is becoming evident that integrating physical principles into these models, or utilizing them in a Δ-ML scheme, significantly enhances their robustness and transferability. This paper introduces PM6-ML, a Δ-ML method that synergizes the semiempirical quantum-mechanical (SQM) method PM6 with a state-of-the-art ML potential applied as a universal correction. The method demonstrates superior performance over standalone SQM and ML approaches and covers a broader chemical space than its predecessors. It is scalable to systems with thousands of atoms, which makes it applicable to large biomolecular systems. Extensive benchmarking confirms PM6-ML's accuracy and robustness. Its practical application is facilitated by a direct interface to MOPAC. The code and parameters are available at https://github.com/Honza-R/mopac-ml.

• Publikační typ
• časopisecké články MeSH

We report the generation of gas-phase riboguanosine radicals that were tagged at ribose with a fixed-charge 6-(trimethylammonium)hexane-1-aminocarbonyl group. The radical generation relied on electron transfer from fluoranthene anion to noncovalent dibenzocrown-ether dication complexes which formed nucleoside cation radicals upon one-electron reduction and crown-ether ligand loss. The cation radicals were characterized by collision-induced dissociation (CID), photodissociation (UVPD), and UV-vis action spectroscopy. Identification of charge-tagged guanosine radicals was challenging because of spontaneous dissociations by loss of a hydrogen atom and guanine that occurred upon storing the ions in the ion trap without further excitation. The loss of H proceeded from an exchangeable position on N-7 in guanine that was established by deuterium labeling and was the lowest energy dissociation of the guanosine radicals according to transition-state energy calculations. Rate constant measurements revealed an inverse isotope effect on the loss of either hydrogen or deuterium with rate constants kH = 0.25-0.26 s-1 and kD = 0.39-0.54 s-1. We used time-dependent density functional theory calculations, including thermal vibronic effects, to predict the absorption spectra of several protomeric radical isomers. The calculated spectra of low-energy N-7-H guanine-radical tautomers closely matched the action spectra. Transition-state-theory calculations of the rate constants for the loss of H-7 and guanine agreed with the experimental rate constants for a narrow range of ion effective temperatures. Our calculations suggest that the observed inverse isotope effect does not arise from the isotope-dependent differences in the transition-state energies. Instead, it may be caused by the dynamics of post-transition-state complexes preceding the product separation.

• Publikační typ
• časopisecké články MeSH

We report a gas-phase UV photodissociation study investigating non-covalent interactions between neutral hydrophobic pentapeptides and peptide ions incorporating a diazirine-tagged photoleucine residue. Phenylalanine (Phe) and proline (Pro) were chosen as the conformation-affecting residues that were incorporated into a small library of neutral pentapeptides. Gas-phase ion-molecule complexes of these peptides with photo-labeled pentapeptides were subjected to photodissociation. Selective photocleavage of the diazirine ring at 355 nm formed short-lived carbene intermediates that underwent cross-linking by insertion into H-X bonds of the target peptide. The cross-link positions were established from collision-induced dissociation tandem mass spectra (CID-MS3) providing sequence information on the covalent adducts. Effects of the amino acid residue (Pro or Phe) and its position in the target peptide sequence were evaluated. For proline-containing peptides, interactions resulting in covalent cross-links in these complexes became more prominent as proline was moved towards the C-terminus of the target peptide sequence. The photocross-linking yields of phenylalanine-containing peptides depended on the position of both phenylalanine and photoleucine. Density functional theory calculations were used to assign structures of low-energy conformers of the (GLPMG + GLL*LK + H)+ complex. Born-Oppenheimer molecular dynamics trajectory calculations were used to capture the thermal motion in the complexes within 100 ps and determine close contacts between the incipient carbene and the H-X bonds in the target peptide. This provided atomic-level resolution of potential cross-links that aided spectra interpretation and was in agreement with experimental data. Graphical Abstract ᅟ.

• Klíčová slova
• Born-Oppenheimer molecular dynamics, Cross-linking, Diazirine tags, Peptide-peptide complexes, Photodissociation,
• Publikační typ
• časopisecké články MeSH

Noncovalent complexes of hydrophobic peptides GLLLG and GLLLK with photoleucine (L*) tagged peptides G(L* n L m )K (n = 1,3, m = 2,0) were generated as singly charged ions in the gas phase and probed by photodissociation at 355 nm. Carbene intermediates produced by photodissociative loss of N2 from the L* diazirine rings underwent insertion into X-H bonds of the target peptide moiety, forming covalent adducts with yields reaching 30%. Gas-phase sequencing of the covalent adducts revealed preferred bond formation at the C-terminal residue of the target peptide. Site-selective carbene insertion was achieved by placing the L* residue in different positions along the photopeptide chain, and the residues in the target peptide undergoing carbene insertion were identified by gas-phase ion sequencing that was aided by specific (13)C labeling. Density functional theory calculations indicated that noncovalent binding to GL*L*L*K resulted in substantial changes of the (GLLLK + H)(+) ground state conformation. The peptide moieties in [GL*L*LK + GLLLK + H](+) ion complexes were held together by hydrogen bonds, whereas dispersion interactions of the nonpolar groups were only secondary in ground-state 0 K structures. Born-Oppenheimer molecular dynamics for 100 ps trajectories of several different conformers at the 310 K laboratory temperature showed that noncovalent complexes developed multiple, residue-specific contacts between the diazirine carbons and GLLLK residues. The calculations pointed to the substantial fluidity of the nonpolar side chains in the complexes. Diazirine photochemistry in combination with Born-Oppenheimer molecular dynamics is a promising tool for investigations of peptide-peptide ion interactions in the gas phase. Graphical Abstract ᅟ.

• Klíčová slova
• Born-Oppenheimer molecular dynamics, Diazirine chromophores, Laser photodissociation, Peptide–peptide ion complexes, Photoleucine,
• MeSH
• diazomethan chemie   MeSH
• fotolýza MeSH
• ionty chemie   MeSH
• leucin chemie   MeSH
• molekulární modely MeSH
• peptidy chemie   MeSH
• plyny chemie   MeSH
• simulace molekulární dynamiky MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Názvy látek
• diazomethan MeSH
• ionty MeSH
• leucin MeSH
• peptidy MeSH
• plyny MeSH

Human aldo-keto reductase 1B15 (AKR1B15) is a newly discovered enzyme which shares 92% amino acid sequence identity with AKR1B10. While AKR1B10 is a well characterized enzyme with high retinaldehyde reductase activity, involved in the development of several cancer types, the enzymatic activity and physiological role of AKR1B15 are still poorly known. Here, the purified recombinant enzyme has been subjected to substrate specificity characterization, kinetic analysis and inhibitor screening, combined with structural modeling. AKR1B15 is active towards a variety of carbonyl substrates, including retinoids, with lower kcat and Km values than AKR1B10. In contrast to AKR1B10, which strongly prefers all-trans-retinaldehyde, AKR1B15 exhibits superior catalytic efficiency with 9-cis-retinaldehyde, the best substrate found for this enzyme. With ketone and dicarbonyl substrates, AKR1B15 also shows higher catalytic activity than AKR1B10. Several typical AKR inhibitors do not significantly affect AKR1B15 activity. Amino acid substitutions clustered in loops A and C result in a smaller, more hydrophobic and more rigid active site in AKR1B15 compared with the AKR1B10 pocket, consistent with distinct substrate specificity and narrower inhibitor selectivity for AKR1B15.

• MeSH
• aldehydreduktasa antagonisté a inhibitory   genetika   metabolismus   MeSH
• aldo-keto reduktasy MeSH
• diterpeny MeSH
• inhibitory enzymů farmakologie   MeSH
• katalytická doména genetika   MeSH
• kinetika MeSH
• konformace proteinů MeSH
• lidé MeSH
• molekulární modely MeSH
• oxidoreduktasy působící na aldehydy nebo donory oxo-skupin antagonisté a inhibitory   genetika   metabolismus   MeSH
• rekombinantní proteiny genetika   metabolismus   MeSH
• retinaldehyd metabolismus   MeSH
• sekvence aminokyselin MeSH
• strukturní homologie proteinů MeSH
• substituce aminokyselin MeSH
• substrátová specifita MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• 9-cis-retinal MeSH Prohlížeč
• AKR1B10 protein, human MeSH Prohlížeč
• AKR1B15 protein, human MeSH Prohlížeč
• aldehydreduktasa MeSH
• aldo-keto reduktasy MeSH
• diterpeny MeSH
• inhibitory enzymů MeSH
• oxidoreduktasy působící na aldehydy nebo donory oxo-skupin MeSH
• rekombinantní proteiny MeSH
• retinaldehyd MeSH

Tuberculosis, the second leading infectious disease killer after HIV, remains a top public health priority. The causative agent of tuberculosis, Mycobacterium tuberculosis (Mtb), which can cause both acute and clinically latent infections, reprograms metabolism in response to the host niche. Phosphoenolpyruvate carboxykinase (Pck) is the enzyme at the center of the phosphoenolpyruvate-pyruvate-oxaloacetate node, which is involved in regulating the carbon flow distribution to catabolism, anabolism, or respiration in different states of Mtb infection. Under standard growth conditions, Mtb Pck is associated with gluconeogenesis and catalyzes the metal-dependent formation of phosphoenolpyruvate. In non-replicating Mtb, Pck can catalyze anaplerotic biosynthesis of oxaloacetate. Here, we present insights into the regulation of Mtb Pck activity by divalent cations. Through analysis of the X-ray structure of Pck-GDP and Pck-GDP-Mn2+ complexes, mutational analysis of the GDP binding site, and quantum mechanical (QM)-based analysis, we explored the structural determinants of efficient Mtb Pck catalysis. We demonstrate that Mtb Pck requires presence of Mn2+ and Mg2+ cations for efficient catalysis of gluconeogenic and anaplerotic reactions. The anaplerotic reaction, which preferably functions in reducing conditions that are characteristic for slowed or stopped Mtb replication, is also effectively activated by Fe2+ in the presence of Mn2+ or Mg2+ cations. In contrast, simultaneous presence of Fe2+ and Mn2+ or Mg2+ inhibits the gluconeogenic reaction. These results suggest that inorganic ions can contribute to regulation of central carbon metabolism by influencing the activity of Pck. Furthermore, the X-ray structure determination, biochemical characterization, and QM analysis of Pck mutants confirmed the important role of the Phe triad for proper binding of the GDP-Mn2+ complex in the nucleotide binding site and efficient catalysis of the anaplerotic reaction.

• MeSH
• aktivace enzymů MeSH
• fosfoenolpyruvátkarboxykinasa (závislá na ATP) chemie   genetika   metabolismus   MeSH
• glukoneogeneze MeSH
• katalýza MeSH
• kationty dvojmocné MeSH
• kineze MeSH
• konformace proteinů MeSH
• lidé MeSH
• molekulární modely MeSH
• molekulární sekvence - údaje MeSH
• multimerizace proteinu MeSH
• mutace MeSH
• Mycobacterium tuberculosis enzymologie   genetika   MeSH
• nukleotidy metabolismus   MeSH
• sekvence aminokyselin MeSH
• sekvenční seřazení MeSH
• vazebná místa MeSH
• vodíková vazba MeSH
• vztahy mezi strukturou a aktivitou MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• fosfoenolpyruvátkarboxykinasa (závislá na ATP) MeSH
• kationty dvojmocné MeSH
• nukleotidy MeSH

The structural characterization of the insulin-insulin receptor (IR) interaction still lacks the conformation of the crucial B21-B30 insulin region, which must be different from that in its storage forms to ensure effective receptor binding. Here, it is shown that insulin analogues modified by natural amino acids at the TyrB26 site can represent an active form of this hormone. In particular, [AsnB26]-insulin and [GlyB26]-insulin attain a B26-turn-like conformation that differs from that in all known structures of the native hormone. It also matches the receptor interface, avoiding substantial steric clashes. This indicates that insulin may attain a B26-turn-like conformation upon IR binding. Moreover, there is an unexpected, but significant, binding specificity of the AsnB26 mutant for predominantly the metabolic B isoform of the receptor. As it is correlated with the B26 bend of the B-chain of the hormone, the structures of AsnB26 analogues may provide the first structural insight into the structural origins of differential insulin signalling through insulin receptor A and B isoforms.

• Klíčová slova
• active conformation, complex, insulin, insulin receptor, isothermal titration microcalorimetry, molecular dynamics,
• MeSH
• fenylalanin MeSH
• fibroblasty metabolismus   MeSH
• inzulin analogy a deriváty   chemie   genetika   metabolismus   MeSH
• konformace proteinů MeSH
• krystalografie rentgenová MeSH
• kultivované buňky MeSH
• lidé MeSH
• lymfocyty metabolismus   MeSH
• molekulární modely MeSH
• mutace MeSH
• myši knockoutované MeSH
• myši MeSH
• potkani Wistar MeSH
• receptor inzulinu chemie   metabolismus   MeSH
• substituce aminokyselin MeSH
• zvířata MeSH
• Check Tag
• lidé MeSH
• mužské pohlaví MeSH
• myši MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• fenylalanin MeSH
• inzulin MeSH
• receptor inzulinu MeSH

We evaluate the performance of the most widely used wavefunction, density functional theory, and semiempirical methods for the description of noncovalent interactions in a set of larger, mostly dispersion-stabilized noncovalent complexes (the L7 data set). The methods tested include MP2, MP3, SCS-MP2, SCS(MI)-MP2, MP2.5, MP2.X, MP2C, DFT-D, DFT-D3 (B3-LYP-D3, B-LYP-D3, TPSS-D3, PW6B95-D3, M06-2X-D3) and M06-2X, and semiempirical methods augmented with dispersion and hydrogen bonding corrections: SCC-DFTB-D, PM6-D, PM6-DH2 and PM6-D3H4. The test complexes are the octadecane dimer, the guanine trimer, the circumcoronene…adenine dimer, the coronene dimer, the guanine-cytosine dimer, the circumcoronene…guanine-cytosine dimer, and an amyloid fragment trimer containing phenylalanine residues. The best performing method is MP2.5 with relative root mean square deviation (rRMSD) of 4 %. It can thus be recommended as an alternative to the CCSD(T)/CBS (alternatively QCISD(T)/CBS) benchmark for molecular systems which exceed current computational capacity. The second best non-DFT method is MP2C with rRMSD of 8 %. A method with the most favorable "accuracy/cost" ratio belongs to the DFT family: BLYP-D3, with an rRMSD of 8 %. Semiempirical methods deliver less accurate results (the rRMSD exceeds 25 %). Nevertheless, their absolute errors are close to some much more expensive methods such as M06-2X, MP2 or SCS(MI)-MP2, and thus their price/performance ratio is excellent.

• Publikační typ
• časopisecké články MeSH

Schistosomiasis caused by a parasitic blood fluke of the genus Schistosoma afflicts over 200 million people worldwide. Schistosoma mansoni cathepsin B1 (SmCB1) is a gut-associated peptidase that digests host blood proteins as a source of nutrients. It is under investigation as a drug target. To further this goal, we report three crystal structures of SmCB1 complexed with peptidomimetic inhibitors as follows: the epoxide CA074 at 1.3 Å resolution and the vinyl sulfones K11017 and K11777 at 1.8 and 2.5 Å resolutions, respectively. Interactions of the inhibitors with the subsites of the active-site cleft were evaluated by quantum chemical calculations. These data and inhibition profiling with a panel of vinyl sulfone derivatives identify key binding interactions and provide insight into the specificity of SmCB1 inhibition. Furthermore, hydrolysis profiling of SmCB1 using synthetic peptides and the natural substrate hemoglobin revealed that carboxydipeptidase activity predominates over endopeptidolysis, thereby demonstrating the contribution of the occluding loop that restricts access to the active-site cleft. Critically, the severity of phenotypes induced in the parasite by vinyl sulfone inhibitors correlated with enzyme inhibition, providing support that SmCB1 is a valuable drug target. The present structure and inhibitor interaction data provide a footing for the rational design of anti-schistosomal inhibitors.

• MeSH
• hemoglobiny chemie   MeSH
• hydrolýza účinky léků   MeSH
• inhibitory proteas chemie   MeSH
• kathepsin D antagonisté a inhibitory   chemie   genetika   MeSH
• krystalografie rentgenová MeSH
• lékové transportní systémy * MeSH
• lidé MeSH
• peptidomimetika chemie   MeSH
• peptidy chemie   MeSH
• proteiny červů antagonisté a inhibitory   chemie   genetika   MeSH
• Schistosoma mansoni enzymologie   genetika   MeSH
• schistosomiasis mansoni farmakoterapie   enzymologie   genetika   MeSH
• vztahy mezi strukturou a aktivitou MeSH
• zvířata MeSH
• Check Tag
• lidé MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Názvy látek
• hemoglobiny MeSH
• inhibitory proteas MeSH
• kathepsin D MeSH
• peptidomimetika MeSH
• peptidy MeSH
• proteiny červů 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
• Názvy látek
• CDK2 protein, human MeSH Prohlížeč
• cyklin-dependentní kinasa 2 MeSH
• inhibitory proteinkinas MeSH
• ligandy MeSH