Molecular dynamics simulations of acylpeptide hydrolase bound to chlorpyrifosmethyl oxon and dichlorvos

Synthesis, In Vivo Anticonvulsant Activity Evaluation and In Silico Studies of Some Quinazolin-4(3H)-One Derivatives

Two series, "a" and "b", each consisting of nine chemical compounds, with 2,3-disubstituted quinazolin-4(3H)-one scaffold, were synthesized and evaluated for their anticonvulsant activity. They were investigated as dual potential positive allosteric modulators of the GABAA receptor at the benzodiazepine binding site and inhibitors of carbonic anhydrase II. Quinazolin-4(3H)-one derivatives were evaluated in vivo (D1-3 = 50, 100, 150 mg/kg, administered intraperitoneally) using the pentylenetetrazole (PTZ)-induced seizure model in mice, with phenobarbital and diazepam, as reference anticonvulsant agents. The in silico studies suggested the compounds act as anticonvulsants by binding on the allosteric site of GABAA receptor and not by inhibiting the carbonic anhydrase II, because the ligands-carbonic anhydrase II predicted complexes were unstable in the molecular dynamics simulations. The mechanism targeting GABAA receptor was confirmed through the in vivo flumazenil antagonism assay. The pentylenetetrazole experimental anticonvulsant model indicated that the tested compounds, 1a-9a and 1b-9b, present a potential anticonvulsant activity. The evaluation, considering the percentage of protection against PTZ, latency until the onset of the first seizure, and reduction in the number of seizures, revealed more favorable results for the "b" series, particularly for compound 8b.

Discovery of A Novel Series of Quinazoline-Thiazole Hybrids as Potential Antiproliferative and Anti-Angiogenic Agents

Considering the pivotal role of angiogenesis in solid tumor progression, we developed a novel series of quinazoline-thiazole hybrids (SA01-SA07) as antiproliferative and anti-angiogenic agents. Four out of the seven compounds displayed superior antiproliferative activity (IC50 =1.83-4.24 µM) on HepG2 cells compared to sorafenib (IC50 = 6.28 µM). The affinity towards the VEGFR2 kinase domain was assessed through in silico prediction by molecular docking, molecular dynamics studies, and MM-PBSA. The series displayed a high degree of similarity to sorafenib regarding the binding pose within the active site of VEGFR2, with a different orientation of the 4-substituted-thiazole moieties in the allosteric pocket. Molecular dynamics and MM-PBSA evaluations identified SA05 as the hybrid forming the most stable complex with VEGFR2 compared to sorafenib. The impact of the compounds on vascular cell proliferation was assessed on EA.hy926 cells. Six compounds (SA01-SA05, SA07) displayed superior anti-proliferative activity (IC50 = 0.79-5.85 µM) compared to sorafenib (IC50 = 6.62 µM). The toxicity was evaluated on BJ cells. Further studies of the anti-angiogenic effect of the most promising compounds, SA04 and SA05, through the assessment of impact on EA.hy296 motility using a wound healing assay and in ovo potential in a CAM assay compared to sorafenib, led to the confirmation of the anti-angiogenic potential.

HPLC Analysis and In Vitro and In Silico Evaluation of the Biological Activity of Polyphenolic Components Separated with Solvents of Various Polarities from Helichrysum italicum

Helichrysum italicum has piqued the interest of many researchers in recent years, mostly for its essential oil, but increasingly for its polyphenolic content as well. In the current study, we examine the polyphenolic composition of H. italicum grown in Bulgaria. The polyphenolic complex was fractionated with solvents of various polarities, including hexane, chloroform, ethyl acetate, and butanol, in order to assess the biological impact of the components. HPLC-PDA and UHPLC-MS/MS were used to examine all fractions. The green coffee fingerprint profile was employed as a "surrogate standard" in the polyphenolic components detection approach. From the UHPLC-MS/MS analysis, we identified 60 components of the polyphenolic complex such as quercetin 3-O-glucuronide, quercetin acetyl-glycoside, isorhamnetin acetyl-glycoside, isorhamnetin caffeoyl-glycoside, quercetin caffeoyl-malonyl-glycoside, isorhamnetin coumaroyl-glycoside, coumaroyl-caffeoylquinic acid, and diCQA-acetyl-derivative were first reported in the composition of H. italicum. The biological activity of the fractions was evaluated in vitro and in silico, which included the fight against oxidative stress (hydrogen peroxide scavenging activity (HPSA), hydroxyl radical scavenging activity (HRSA), metal-chelating activity (MChA)) and nitrosative (nitric oxide scavenging activity) (NOSA)), in vitro anti-inflammatory, and anti-arthritic activity. Results are presented as IC50 ± SD μg/mL. The analysis showed that the EtOAc fraction was characterized by highest HPSA (57.12 ± 1.14 μg/mL), HRSA (92.23 ± 1.10 μg/mL), MChA (5.60 ± 0.17 μg/mL), and NOSA (89.81 ± 2.09 μg/mL), while the hexane and chloroform fractions showed significantly higher in vitro anti-inflammatory activity (30.48 ± 2.33 μg/mL, 62.50 ± 1.69 μg/mL) compared to the standard ibuprofen. All three fractions showed potential anti-arthritic activity (102.93 ± 8.62 μg/mL, 108.92 ± 4.42 μg/mL, 84.19 ± 3.89 μg/mL).

Studies on the Selectivity Mechanism of Wild-Type E. coli Thioesterase ‘TesA and Its Mutants for Medium- and Long-Chain Acyl Substrates

E. coli thioesterase ‘TesA is an important enzyme in fatty acid production. Medium-chain fatty acids (MCFAs, C6-C10) are of great interest due to their similar physicochemical properties to petroleum-based oleo-chemicals. It has been shown that wild-type ‘TesA had better selectivity for long-chain acyl substrates (≥C16), while the two mutants ‘TesAE142D/Y145G and ‘TesAM141L/E142D/Y145G had better selectivity for medium-chain acyl substrates. However, it is difficult to obtain the selectivity mechanism of substrates for proteins by traditional experimental methods. In this study, in order to obtain more MCFAs, we analyzed the binding mode of proteins (‘TesA, ‘TesAE142D/Y145G and ‘TesAM141L/E142D/Y145G) and substrates (C16/C8-N-acetylcysteamine analogs, C16/C8-SNAC), the key residues and catalytic mechanisms through molecular docking, molecular dynamics simulations and the molecular mechanics Poisson–Boltzmann surface area (MM/PBSA). The results showed that several main residues related to catalysis, including Ser10, Asn73 and His157, had a strong hydrogen bond interaction with the substrates. The mutant region (Met141-Tyr146) and loop107–113 were mainly dominated by Van der Waals contributions to the substrates. For C16-SNAC, except for ‘TesAM141L/E142D/Y145G with large conformational changes, there were strong interactions at both head and tail ends that distorted the substrate into a more favorable high-energy conformation for the catalytic reaction. For C8-SNAC, the head and tail found it difficult to bind to the enzyme at the same time due to insufficient chain length, which made the substrate binding sites more variable, so ‘TesAM141L/E142D/Y145G with better binding sites had the strongest activity, and ‘TesA had the weakest activity, conversely. In short, the matching substrate chain and binding pocket length are the key factors affecting selectivity. This will be helpful for the further improvement of thioesterases.

Computational Study on the Effect of Inactivating/Activating Mutations on the Inhibition of MEK1 by Trametinib

Activation of the mitogen-activated protein kinase (MAPK) signaling pathway regulated by human MAP kinase 1 (MEK1) is associated with the carcinogenesis and progression of numerous cancers. In addition, two active mutations (P124S and E203K) have been reported to enhance the activity of MEK1, thereby eventually leading to the tumorigenesis of cancer. Trametinib is an MEK1 inhibitor for treating EML4-ALK-positive, EGFR-activated, and KRAS-mutant lung cancers. Therefore, in this study, molecular docking and molecular dynamic (MD) simulations were performed to explore the effects of inactive/active mutations (A52V/P124S and E203K) on the conformational changes of MEK1 and the changes in the interaction of MEK1 with trametinib. Moreover, steered molecular dynamic (SMD) simulations were further utilized to compare the dissociation processes of trametinib from the wild-type (WT) MEK1 and two active mutants (P124S and E203K). As a result, trametinib had stronger interactions with the non-active MEK1 (WT and A52V mutant) than the two active mutants (P124S and E203K). Moreover, two active mutants may make the allosteric channel of MEK1 wider and shorter than that of the non-active types (WT and A52V mutant). Hence, trametinib could dissociate from the active mutants (P124S and E203K) more easily compared with the WT MEK1. In summary, our theoretical results demonstrated that the active mutations may attenuate the inhibitory effects of MEK inhibitor (trametinib) on MEK1, which could be crucial clues for future anti-cancer treatment.

Molecular dynamics investigations of oligosaccharides recognized by family 16 and 22 carbohydrate binding modules

As a non-catalytic domain, carbohydrate binding modules (CBMs) are often considered to play some key roles in the degradation and recognition of polysaccharides catalyzed by cellulases. In this work, we investigated the recognition dynamics of cello- or xylo-saccharides by two typical CBMs (CBM16-1 and CBM22-2), which are grouped into Type B CBMs. By combining extensive molecular dynamics, principle component analysis, and binding free energy calculations, we constructed several complex models of the two CBMs in both complex cello- and xylo-oligosaccharides. The corresponding substrate recognition affinity and critical residues having significant contributions were systematically investigated. The residues containing aromatic side chain groups were shown to contribute significantly to substrate binding. The calculated binding free energies were in fairly good agreement with the experimental measurements with the absolute mean error of 0.69 kcal mol^-1. The overall electrostatic interactions were shown to have negative effects on substrate recognition. Further metadynamics simulations revealed the substrate dissociation process.

Effects of Tyr555 and Trp678 on the processivity of cellobiohydrolase A from Ruminiclostridium thermocellum: A simulation study

Cellobiohydrolase A from Ruminiclostridium thermocellum (Cbh9A) is a processive exoglucanase from family 9 and is an important cellobiohydrolase that hydrolyzes cello-oligosaccharide into cellobiose. Residues Tyr555 and Trp678 considerably affect catalytic activity, but their mechanisms are still unknown. To investigate how the Tyr555 and Trp678 affect the processivity of Cbh9A, conventional molecular dynamics, steered molecular dynamics, and free energy calculation were performed to simulate the processive process of wild type (WT)-Cbh9A, Y555S mutant, and W678G mutant. Analysis of simulation results suggests that the binding free energies between the substrate and WT-Cbh9A are lower than those of Y555S and W678G mutants. The pull forces and energy barrier in Y555S and W678G mutants also reduced significantly during the steered molecular dynamics (SMD) simulation compared with that of the WT-Cbh9A. And the potential mean force calculations showed that the pulling energy barrier of Y555S and W678G mutants is much lower than that of WT-Cbh9A.

Probing inhibition mechanisms of adenosine deaminase by using molecular dynamics simulations

Adenosine deaminase (ADA) catalyzes the deamination of adenosine, which is important in purine metabolism. ADA is ubiquitous to almost all human tissues, and ADA abnormalities have been reported in various diseases, including rheumatoid arthritis. ADA can be divided into two conformations based on the inhibitor that it binds to: open and closed forms. Here, we chose three ligands, namely, FR117016 (FR0), FR221647 (FR2) (open form), and HDPR (PRH, closed form), to investigate the inhibition mechanism of ADA and its effect on ADA through molecular dynamics simulations. In open forms, Egap and electrostatic potential (ESP) indicated that electron transfer might occur more easily in FR0 than in FR2. Binding free energy and hydrogen bond occupation revealed that the ADA-FR0 complex had a more stable structure than ADA-FR2. The probability of residues Pro159 to Lys171 of ADA-FR0 and ADA-FR2 to form a helix moderately increased compared with that in nonligated ADA. In comparison with FR0 and FR2 PRH could maintain ADA in a closed form to inhibit the function of ADA. The α7 helix (residues Thr57 to Ala73) of ADA in the closed form was mostly unfastened because of the effect of PRH. The number of H bonds and the relative superiority of the binding free energy indicated that the binding strength of PRH to ADA was significantly lower than that of an open inhibitor, thereby supporting the comparison of the inhibitory activities of the three ligands. Alanine scanning results showed that His17, Gly184, Asp295, and Asp296 exerted the greatest effects on protein energy, suggesting that they played crucial roles in binding to inhibitors. This study served as a theoretical basis for the development of new ADA inhibitors.

Understanding the interactions of different substrates with wild-type and mutant acylaminoacyl peptidase using molecular dynamics simulations

Acylaminoacylpeptidase (AAP) belongs to peptidase protein family, which can degrade amyloid β-peptide forms in the brains of patients, and hence leads to Alzheimer's disease. And so, AAP is considered to be a novel target in the design of drugs against Alzheimer's disease. In this investigation, six molecular dynamics simulations were used to find that the interaction between the wild-type and R526V AAP with two different substrates (p-nitrophenylcaprylate and Ac-Leu-p-nitroanilide). Our results were as follows: firstly, Ac-Leu-p-nitroanilide bound to R526V AAP to form a more disordered loop (residues 552-562) in the α/β-hydrolase fold like of AAP, which caused an open and inactive AAP domain form, secondly, binding p-nitrophenylcaprylate and Ac-Leu-p-nitroanilide to AAP can decrease the flexibility of residues 225-250, 260-270, and 425-450, in which the ordered secondary structures may contain the suitable geometrical structure and so it is useful to serine attack. Our theoretical results showed that the binding of the two substrates can induce specific conformational changes responsible for the diverse AAP catalytic specificity. These theoretical substrate-induced structural diversities can help explain the abilities of AAPs to recognize and hydrolyze extremely different substrates.

Exploration of binding and inhibition mechanism of a small molecule inhibitor of influenza virus H1N1 hemagglutinin by molecular dynamics simulation

Influenza viruses are a major public health threat worldwide. The influenza hemagglutinin (HA) plays an essential role in the virus life cycle. Due to the high conservation of the HA stem region, it has become an especially attractive target for inhibitors for therapeutics. In this study, molecular simulation was applied to study the mechanism of a small molecule inhibitor (MBX2329) of influenza HA. Behaviors of the small molecule under neutral and acidic conditions were investigated, and an interesting dynamic binding mechanism was found. The results suggested that the binding of the inhibitor with HA under neutral conditions facilitates only its intake, while it interacts with HA under acidic conditions using a different mechanism at a new binding site. After a series of experiments, we believe that binding of the inhibitor can prevent the release of HA1 from HA2, further maintaining the rigidity of the HA2 loop and stabilizing the distance between the long helix and short helices. The investigated residues in the new binding site show high conservation, implying that the new binding pocket has the potential to be an effective drug target. The results of this study will provide a theoretical basis for the mechanism of new influenza virus inhibitors.

Understanding the Phosphorylation Mechanism by Using Quantum Chemical Calculations and Molecular Dynamics Simulations

Phosphorylation is one of the most frequent post-translational modifications on proteins. It regulates many cellular processes by modulation of phosphorylation on protein structure and dynamics. However, the mechanism of phosphorylation-induced conformational changes of proteins is still poorly understood. Here, we report a computational study of three representative groups of tyrosine in ADP-ribosylhydrolase 1, serine in BTG2, and serine in Sp100C by using six molecular dynamics (MD) simulations and quantum chemical calculations. Added phosphorylation was found to disrupt hydrogen bond, and increase new weak interactions (hydrogen bond and hydrophobic interaction) during MD simulations, leading to conformational changes. Quantum chemical calculations further indicate that the phosphorylation on tyrosine, threonine, and serine could decrease the optical band gap energy (Egap), which can trigger electronic transitions to form or disrupt interactions easily. Our results provide an atomic and electronic description of how phosphorylation facilitates conformational and dynamic changes in proteins, which may be useful for studying protein function and protein design.

Models for the binding channel of wild type and mutant transthyretin with glabridin

Our results indicate that additional high-occupancy hydrogen bonds were observed at the binding interface between the two dimers in V30A TTR, while stabilisation hydrophobic interactions between residues in the mutant AB loop decreased.

Exploring the different ligand escape pathways in acylaminoacyl peptidase by random acceleration and steered molecular dynamics simulations

Acylaminoacyl peptidase (APH, EC 3.4.19.1) is a novel class of serine-type protease belonging to the prolyl oligopeptidase (POP) family.