Search

Platelets are the critical target for preventing and treating pathological thrombus formation. However, despite current antiplatelet therapy, cardiovascular mortality remains high, and cardiovascular events continue in prescribed patients. In this study, first results were obtained with ortho-carbonyl hydroquinones as antiplatelet agents; we found that linking triphenylphosphonium cation to a bicyclic ortho-carbonyl hydroquinone moiety by a short alkyl chain significantly improved their antiplatelet effect by affecting the mitochondrial functioning. The mechanism of action involves uncoupling OXPHOS, which leads to an increase in mitochondrial ROS production and a decrease in the mitochondrial membrane potential and OCR. This alteration disrupts the energy production by mitochondrial function necessary for the platelet activation process. These effects are responsive to the complete structure of the compounds and not to isolated parts of the compounds tested. The results obtained in this research can be used as the basis for developing new antiplatelet agents that target mitochondria.

The papain-like protease (PLpro) found in coronaviruses that can be transmitted from animals to humans is a critical target in respiratory diseases linked to Severe Acute Respiratory Syndrome (SARS-CoV). Researchers have proposed designing PLpro inhibitors. In this study, a set of 89 compounds, including recently reported 2-phenylthiophenes with nanomolar inhibitory potency, were investigated as PLpro noncovalent inhibitors using advanced molecular modeling techniques. To develop the work with these inhibitors, multiple structures of the SARS-CoV-2 PLpro binding site were generated using a molecular sampling method. These structures were then clustered to select a group that represents the flexibility of the site. Subsequently, models of the protein-ligand complexes were created for the set of inhibitors within the chosen conformations. The quality of the complex models was assessed using LigRMSD software to verify similarities in the orientations of the congeneric series and interaction fingerprints to determine the recurrence of chemical interactions. With the multiple models constructed, a protocol was established to choose one per ligand, optimizing the correlation between the calculated docking energy values and the biological activities while incorporating the effect of the binding site’s flexibility. A strong correlation (R2 = 0.922) was found when employing this flexible docking protocol.

The design of TRPV1 antagonists and agonists has reached a new era since TRPV1 structures at near-atomic resolution are available. Today, the ligand-binding forms of several classical antagonists and agonists are known; therefore, the specific role of key TRPV1’s residues in binding of ligands can be elucidated. It is possible to place the well-defined pharmacophore of TRPV1 ligands, conformed by head, neck, and tail groups, in the right pocket regions of TRPV1. It will allow a more thorough use of molecular modelling methods to conduct more effective rational drug design protocols. In this work, important points about the interactions between TRPV1 and capsaicin-like compounds are spelled out, based on the known pharmacophore of the ligands and the already available TRPV1 structures. These points must be addressed to generate reliable poses of novel candidates and should be considered during the design of novel TRPV1 antagonists and agonists.

Multicomponent reactions were performed to develop novel α,β-unsaturated carbonyl depsipeptides and peptoids incorporating various chromophores such as cinnamic, coumarin, and quinolines. Thus, through the Passerini and Ugi multicomponent reactions (P-3CR and U-4CR), we obtained thirteen depsipeptides and peptoids in moderate to high yield following the established protocol and fundamentally varying the electron-rich carboxylic acid as reactants. UV/Vis spectroscopy was utilized to study the photophysical properties of the newly synthesized compounds. Differences between the carbonyl-substituted chromophores cause differences in electron delocalization that can be captured in the spectra. The near UV regions of all the compounds exhibited strong absorption bands. Compounds P2, P5, U2, U5, and U7 displayed absorption bands in the range of 250–350 nm, absorbing radiation in this broad region of the electromagnetic spectrum. A photostability study for U5 showed that its molecular structure does not change after exposure to UV radiation. Fluorescence analysis showed an incipient emission of U5, while U6 showed blue fluorescence under UV radiation. The photophysical properties and electronic structure were also determined by TD-DFT theoretical study.

The papain-like protease (PLpro) from zoonotic coronaviruses (CoVs) has been identified as a target with an essential role in viral respiratory diseases caused by Severe Acute Respiratory Syndrome-associated coronaviruses (SARS-CoVs). The design of PLpro inhibitors has been proposed as an alternative to developing potential drugs against this disease. In this work, 67 naphthalene-derived compounds as noncovalent PLpro inhibitors were studied using molecular modeling methods. Structural characteristics of the bioactive conformations of these inhibitors and their interactions at the SARS-CoV-1 PLpro binding site were reported here in detail, taking into account the flexibility of the protein residues. Firstly, a molecular docking protocol was used to obtain the orientations of the inhibitors. After this, the orientations were compared, and the recurrent interactions between the PLpro residues and ligand chemical groups were described (with LigRMSD and interaction fingerprints methods). In addition, efforts were made to find correlations between docking energy values and experimentally determined binding affinities. For this, the PLpro was sampled by using Gaussian Accelerated Molecular Dynamics (GaMD), generating multiple conformations of the binding site. Diverse protein conformations were selected and a cross-docking experiment was performed, yielding models of the 67 naphthalene-derived compounds adopting different binding modes. Representative complexes for each ligand were selected to obtain the highest correlation between docking energies and activities. A good correlation (R2 = 0.948) was found when this flexible docking protocol was performed.

Purpose: This study assessed the impact of different types of medial foot arch on postural stability and core center of gravity muscle activity among collegiate athletes. Methods: The study sample included 103 university-level athletes across various sports (soccer, rugby, basketball, volleyball, field tennis, table tennis, karate, and cheerleading) from the College of Magdalena (Colombia) who exhibited distinct types of medial foot arch: 32 high, 35 low, and 36 neutral arches. Surface electromyography (sEMG) was employed to assess conduction velocity, magnitude values, latency, and fatigue in focal muscles including the spinal erector (SE), internal oblique (IO), external oblique (EO), and rectus abdominis (AR), while measurements of static and dynamic postural control were also considered. Post hoc analysis was performed with Bonferroni correction for all electromyographically measured muscle groups, as well as for measurements of static and dynamic postural stability. Pearson’s or Spearman’s correlation tests were used to compare the different types of feet. Results: There were no substantial differences observed between the distinct types of feet in terms of focal muscle activity, static stability, or dynamics. Even though the mean values indicated higher muscle activity and stability among those with high foot arches and lower values among those with low arches compared to the neutral foot type, this observed difference was deemed statistically insignificant. We also observed a positive correlation between internal oblique muscle activity and the average power of dynamic postural stability, which remained consistent across all foot types. Our findings indicate that static instability is directly correlated with dynamic instability in the anteroposterior direction, while a clear inverse relationship was established in the lateral direction upon examining the variable correlations. Conclusions: The presence of high or low foot arches did not significantly impact the activity of the muscles responsible for maintaining the body’s center of gravity or postural stability among university-level athletes. This suggests the existence of neuromuscular compensation mechanisms that attempt to restore balance and compensate for any changes in postural stability caused by varying foot types. Through targeted training that emphasizes activation of the internal oblique muscle, athletes may see improved postural stability. Our findings indicate that static stabilization exercises can also prove beneficial in improving dynamic stability in the anteroposterior plane, while a more dynamic approach may be required to improve dynamic stability in the lateral plane.

Parkinson's disease (PD) is a neurodegenerative disorder that affects adult people whose treatment is palliative. Thus, we decided to test three dammarane triterpenes 1, 1a, 1b, and we determined that 1 and 1a inhibit β-aggregation through thioflavine T rather than 1b. Since compound 1 was most active, we determined the interaction between α-synuclein and 1 at 50 µM (Kd) through microscale thermophoresis. Also, we observed differences in height and diameter of aggregates, and α-synuclein remains unfolded in the presence of 1. Also, aggregates treated with 1 do not provoke neurites' retraction in N2a cells previously induced by retinoic acid. Finally, we studied the potential sites of interaction between 1 with α-synuclein fibrils using molecular modelling. Docking experiments suggest that 1 preferably interact with the site 2 of α-synuclein through hydrogen bonds with residues Y39 and T44.

Selected adenosine A2A receptor agonists (PSB-15826, PSB-12404, and PSB-16301) have been evaluated as new antiplatelet agents. In addition, radioligand-binding studies and receptor-docking experiments were performed in order to explain their differential biological effects on a molecular level. Among the tested adenosine derivatives, PSB-15826 was the most potent compound to inhibit platelet aggregation (EC50 0.32 ± 0.05 µmol/L) and platelet P-selectin cell-surface localization (EC50 0.062 ± 0.2 µmol/L), and to increase intraplatelets cAMP levels (EC50 0.24 ± 0.01 µmol/L). The compound was more active than CGS21680 (EC50 0.97±0.07 µmol/L) and equipotent to NECA (EC50 0.31 ± 0.05 µmol/L) in platelet aggregation induced by ADP. In contrast to the results from cAMP assays, Ki values determined in radioligand-binding studies were not predictive of the A2A agonists’ antiplatelet activity. Docking studies revealed the key molecular determinants of this new family of adenosine A2A receptor agonists: differences in activities are related to π-stacking interactions between the ligands and the residue His264 in the extracellular loop of the adenosine A2A receptor which may result in increased residence times. In conclusion, these results provide an improved understanding of the requirements of antiplatelet adenosine A2A receptor agonists.

Alzheimer’s disease is a common tauopathy where fibril formation and aggregates are the hallmark of the disease. Efforts targeting amyloid-β plaques have succeeded to remove plaques but failed in clinical trials to improve cognition; thus, the current therapeutic strategy is at preventing tau aggregation. Here, we demonstrated that four phenolic diterpenoids and rosmarinic acid inhibit fibrillization. Since, rosmarinic acid was the most active compound, we observe morphological changes in atomic force microscopy images after treatment. Hence, rosmarinic acid leads to a decrease in amide regions I and III, indicating that rosmarinic acid prevents β-sheet assembly. Molecular docking study inside the steric zipper model of the hexapeptide 306VQIVYK311 involved in fibrillization and β sheet formation, suggests that rosmarinic acid binds to the steric zipper with similar chemical interactions with respect to those observed for orange G, a known pharmacofore for amyloid.

The human Sigma1 receptor (S1R), which has been identified as a target with an important role in neuropsychological disorders, was first crystallized 3 years ago. Since S1R structure has no relation with another previous crystallized structures, the presence of the new crystal is an important hallmark for the design of agonists and antagonists against this important target. Some years ago, our group identified RC-33, a potent and selective S1R agonist, endowed with neuroprotective properties. In this work, drawing on new structural information, we studied the interactions of RC-33 and its analogs with the S1R binding site by using computational methods such as docking, interaction fingerprints, and receptor-guided alignment three dimensional quantitative structure–activity relationship (3D-QSAR). We found that RC-33 and its analogs adopted similar orientations within S1R binding site, with high similitude with orientations of the crystallized ligands; such information was used for identifying the residues involved in chemical interactions with ligands. Furthermore, the structure-activity relationship of the studied ligands was adequately described considering classical QSAR tests. All relevant aspects of the interactions between the studied compounds and S1R were covered here, through descriptions of orientations, binding interactions, and features that influence differential affinities. In this sense, the present results could be useful in the future design of novel S1R modulators.

Cyclin-dependent kinase 2 (CDK2) is an important member of the CDK family exerting its most important function in the regulation of the cell cycle. It catalyzes the transfer of the gamma phosphate group from an ATP (adenosine triphosphate) molecule to a Serine/Threonine residue of a peptide substrate. Due to the importance of this enzyme, and protein kinases in general, a detailed understanding of the reaction mechanism is desired. Thus, in this work the phosphoryl transfer reaction catalyzed by CDK2 was revisited and studied by means of hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. Our results suggest that the base-assisted mechanism is preferred over the substrate-assisted pathway when one Mg2+ is present in the active site, in agreement with a previous theoretical study. The base-assisted mechanism resulted to be dissociative, with a potential energy barrier of 14.3 kcal/mol, very close to the experimental derived value. An interesting feature of the mechanism is the proton transfer from Lys129 to the phosphoryl group at the second transition state, event that could be helping in neutralizing the charge on the phosphoryl group upon the absence of a second Mg2+ ion. Furthermore, important insights into the mechanisms in terms of bond order and charge analysis were provided. These descriptors helped to characterize the synchronicity of bond forming and breaking events, and to characterize charge transfer effects. Local interactions at the active site are key to modulate the charge distribution on the phosphoryl group and therefore alter its reactivity.

Several years ago, the crystallographic structures of the transient receptor potential vanilloid 1 (TRPV1) in the presence of agonists and antagonists were reported, providing structural information about its chemical activation and inactivation. TRPV1’s activation increases the transport of calcium and sodium ions, leading to the excitation of sensory neurons and the perception of pain. On the other hand, its antagonistic inactivation has been explored to design analgesic drugs. The interactions between the antagonists 5,5-diarylpentadienamides (DPDAs) and TRPV1 were studied here to explain why they inactivate TRPV1. The present work identified the structural features of TRPV1–DPDA complexes, starting with a consideration of the orientations of the ligands inside the TRPV1 binding site by using molecular docking. After this, a chemometrics analysis was performed (i) to compare the orientations of the antagonists (by using LigRMSD), (ii) to describe the recurrent interactions between the protein residues and ligand groups in the complexes (by using interaction fingerprints), and (iii) to describe the relationship between topological features of the ligands and their differential antagonistic activities (by using a quantitative structure–activity relationship (QSAR) with 2D autocorrelation descriptors). The interactions between the DPDA groups and the residues Y511, S512, T550, R557, and E570 (with a recognized role in the binding of classic ligands), and the occupancy of isoquinoline or 3-hydroxy-3,4-dihydroquinolin-2(1H)-one groups of the DPDAs in the vanilloid pocket of TRPV1 were clearly described. Based on the results, the structural features that explain why DPDAs inactivate TRPV1 were clearly exposed. These features can be considered for the design of novel TRPV1 antagonists.

The angiotensin-converting enzyme (ACE) is a two-domain dipeptidylcarboxypeptidase, which has a direct involvement in the control of blood pressure by performing the hydrolysis of angiotensin I to produce angiotensin II. At the same time, ACE hydrolyzes other substrates such as the vasodilator peptide bradykinin and the anti-inflammatory peptide N-acetyl-SDKP. In this sense, ACE inhibitors are bioactive substances with potential use as medicinal products for treatment or prevention of hypertension, heart failures, myocardial infarction, and other important diseases. This review examined the most recent literature reporting ACE inhibitors with the help of molecular modeling. The examples exposed here demonstrate that molecular modeling methods, including docking, molecular dynamics (MD) simulations, quantitative structure-activity relationship (QSAR), etc, are essential for a complete structural picture of the mode of action of ACE inhibitors, where molecular docking has a key role. Examples show that too many works identified ACE inhibitory activities of natural peptides and peptides obtained from hydrolysates. In addition, other works report non-peptide compounds extracted from natural sources and synthetic compounds. In all these cases, molecular docking was used to provide explanation of the chemical interactions between inhibitors and the ACE binding sites. For docking applications, most of the examples exposed here do not consider that: (i) ACE has two domains (nACE and cACE) with available X-ray structures, which are relevant for the design of selective inhibitors, and (ii) nACE and cACE binding sites have large dimensions, which leads to non-reliable solutions during docking calculations. In support of the solution of these problems, the structural information found in Protein Data Bank (PDB) was used to perform an interaction fingerprints (IFPs) analysis applied on both nACE and cACE domains. This analysis provides plots that identify the chemical interactions between ligands and both ACE binding sites, which can be used to guide docking experiments in the search of selective natural components or novel drugs. In addition, the use of hydrogen bond constraints in the S2 and S2′ subsites of nACE and cACE are suggested to guarantee that docking solutions are reliable.

PDE3s belong to the phosphodiesterases family, where the PDE3A isoform is the major subtype in platelets involved in the cAMP regulation pathway of platelet aggregation. PDE3A inhibitors have been designed as potential antiplatelet agents. In this work, a homology model of PDE3A was developed and used to obtain the binding modes of bicyclic heteroaromatic pyridazinones and pyrazolones. Most of the studied compounds adopted similar orientations within the PDE3A active site, establishing hydrogen bonds with catalytic amino acids. Besides, the structure-activity relationship of the studied inhibitors was described by using a field-based 3D-QSAR method. Different structure alignment strategies were employed, including template-based and docking-based alignments. Adequate correlation models were obtained according to internal and external validations. In general, QSAR models revealed that steric and hydrophobic fields describe the different inhibitory activities of the compounds, where the hydrogen bond donor and acceptor fields have minor contributions. It should be stressed that structural elements of PDE3A inhibitors are reported here, through descriptions of their binding interactions and their differential affinities. In this sense, the present results could be useful in the future design of more specific and potent PDE3A inhibitors that may be used for the treatment of cardiovascular diseases.

TASK-3 is a two-pore domain potassium (K2P) channel highly expressed in the hippocampus, cerebellum, and cortex. TASK-3 has been identified as an oncogenic potassium channel and it is overexpressed in different cancer types. For this reason, the development of new TASK-3 blockers could influence the pharmacological treatment of cancer and several neurological conditions. In the present work, we searched for novel TASK-3 blockers by using a virtual screening protocol that includes pharmacophore modeling, molecular docking, and free energy calculations. With this protocol, 19 potential TASK-3 blockers were identified. These molecules were tested in TASK-3 using patch clamp, and one blocker (DR16) was identified with an IC50 = 56.8 ± 3.9 μM. Using DR16 as a scaffold, we designed DR16.1, a novel TASK-3 inhibitor, with an IC50 = 14.2 ± 3.4 μM. Our finding takes on greater relevance considering that not many inhibitory TASK-3 modulators have been reported in the scientific literature until today. These two novel TASK-3 channel inhibitors (DR16 and DR16.1) are the first compounds found using a pharmacophore-based virtual screening and rational drug design protocol.

Ugi four component reaction (Ugi-4CR) isocyanide-based multicomponent reactions were used to synthesize diN-substituted glycyl-phenylalanine (diNsGF) derivatives. All of the synthesized compounds were characterized by spectroscopic and spectrometric techniques. In order to evaluate potential biological applications, the synthesized compounds were tested in computational models that predict the bioactivity of organic molecules by using only bi-dimensional molecular information. The diNsGF derivatives were predicted as cholinesterase inhibitors. Experimentally, all of the synthesized diNsGF derivatives showed moderate inhibitory activities against acetylcholinesterase (AChE) and poor activities against butyrylcholinesterase (BuChE). Compound 7a has significant activity and selectivity against AChE, which reveals that the diNsGF scaffold could be improved to reach novel candidates by combining other chemical components of the Ugi-4CR in a high-throughput combinatorial screening experiment. Molecular docking experiments of diNsGF derivatives inside AChE suggest that these compounds placed the phenylalanine group at the peripheral site of AChE. The orientations and chemical interactions of diNsGF derivatives were analyzed, and the changeable groups were identified for future exploration of novel candidates that could lead to the improvement of diNsGF derivative inhibitory activities.

In the last years, the interactions of flavonoids with protein kinases (PKs) have been described by using crystallographic experiments. Interestingly, different orientations have been found for one flavonoid inside different PKs and different chemical substitutions lead to different orientations of the flavonoid scaffold inside one PK. Accordingly, orientation predictions of novel analogues could help to the design of flavonoids with high PK inhibitory activities. With this in mind, we studied the binding modes of 37 flavonoids (flavones and chalcones) inside the cyclin-dependent PK CDK1 using docking experiments. We found that the compounds under study adopted two different orientations into the active site of CDK1 (orientations I and II in the manuscript). In addition, quantitative structure-activity relationship (QSAR) models using CoMFA and CoMSIA methodologies were constructed to explain the trend of the CDK1 inhibitory activities for the studied flavonoids. Template-based and docking-based alignments were used. Models developed starting from docking-based alignment were applied for describing the whole dataset and compounds with orientation I. Adequate R2 and Q2 values were obtained by each method; interestingly, only hydrophobic and hydrogen bond donor fields describe the differential potency of the flavonoids as CDK1 inhibitors for both defined alignments and subsets. Our current application of docking and QSAR together reveals important elements to be drawn for the design of novel flavonoids with increased PK inhibitory activities.

It is evident that 99mTc causes radical-mediated DNA damage due to Auger electrons, which were emitted simultaneously with the known γ-emission of 99mTc. We have synthesized a series of new 99mTc-labeled pyrene derivatives with varied distances between the pyrene moiety and the radionuclide. The pyrene motif is a common DNA intercalator and allowed us to test the influence of the radionuclide distance on damages of the DNA helix. In general, pUC 19 plasmid DNA enables the investigation of the unprotected interactions between the radiotracers and DNA that results in single-strand breaks (SSB) or double-strand breaks (DSB). The resulting DNA fragments were separated by gel electrophoresis and quantified by fluorescent staining. Direct DNA damage and radical-induced indirect DNA damage by radiolysis products of water were evaluated in the presence or absence of the radical scavenger DMSO. We demonstrated that Auger electrons directly induced both SSB and DSB in high efficiency when 99mTc was tightly bound to the plasmid DNA and this damage could not be completely prevented by DMSO, a free radical scavenger. For the first time, we were able to minimize this effect by increasing the carbon chain lengths between the pyrene moiety and the 99mTc nuclide. However, a critical distance between the 99mTc atom and the DNA helix could not be determined due to the significantly lowered DSB generation resulting from the interaction which is dependent on the type of the 99mTc binding motif. The effect of variable DNA damage caused by the different chain length between the pyrene residue and the Tc-core as well as the possible conformations of the applied Tc-complexes was supplemented with molecular dynamics (MD) calculations. The effectiveness of the DNA-binding 99mTc-labeled pyrene derivatives was demonstrated by comparison to non-DNA-binding 99mTcO4-, since nearly all DNA damage caused by 99mTcO4- was prevented by incubating with DMSO.

Aroma and flavor are important factors of fruit quality and consumer preference. The specific pattern of aroma is generated during ripening by the accumulation of volatiles compounds, which are mainly esters. Alcohol acyltransferase (AAT) (EC 2.3.1.84) catalyzes the esterification reaction of aliphatic and aromatic alcohols and acyl-CoA into esters in fruits and flowers. In Fragaria x ananassa, there are different volatiles compounds that are obtained from different alcohol precursors, where octanol and hexanol are the most abundant during fruit ripening. At present, there is not structural evidence about the mechanism used by the AAT to synthesize esters. Experimental data attribute the kinetic role of this enzyme to 2 amino acidic residues in a highly conserved motif (HXXXD) that is located in the middle of the protein. With the aim to understand the molecular and energetic aspects of volatiles compound production from F. x ananassa, we first studied the binding modes of a series of alcohols, and also different acyl-CoA substrates, in a molecular model of alcohol acyltransferase from Fragaria x ananassa (SAAT) using molecular docking. Afterwards, the dynamical behavior of both substrates, docked within the SAAT binding site, was studied using routine molecular dynamics (MD) simulations. In addition, in order to correlate the experimental and theoretical data obtained in our laboratories, binding free energy calculations were performed; which previous results suggested that octanol, followed by hexanol, presented the best affinity for SAAT. Finally, and concerning the SAAT molecular reaction mechanism, it is suggested from molecular dynamics simulations that the reaction mechanism may proceed through the formation of a ternary complex, in where the Histidine residue at the HXXXD motif deprotonates the alcohol substrates. Then, a nucleophilic attack occurs from alcohol charged oxygen atom to the carbon atom at carbonyl group of the acyl CoA. This mechanism is in agreement with previous results, obtained in our group, in alcohol acyltransferase from Vasconcellea pubescens (VpAAT1).

(E)-2-(benzo[d]thiazol-2-yl)-3-heteroarylacrylonitriles are described as a new class of selective inhibitors of acetylcholinesterase (AChE). The most potent compound in the series exhibited good AChE inhibitory activity (IC50 = 64 µM). Compound 7f was found to be more selective than galanthamine in inhibiting AChE and it showed a moderate selectivity index. Kinetic studies on AChE indicated that a competitive type of inhibition pattern exist for these acrylonitrile derivates. Molecular docking models of the ligand-AChE complexes suggest that compound 7g is located on the periphery of the AChE active site.

Human arginase I (hARGI) is an important enzyme involved in the urea cycle; its overexpression has been associated to cardiovascular and cerebrovascular diseases. In the last years, several congeneric sets of hARGI inhibitors have been reported with possible beneficial roles for the cardiovascular system. At the same time, crystallographic data have been reported including hARGI–inhibitor complexes, which can be considered for the design of novel inhibitors. In this work, the structure–activity relationship (SAR) of Cα substituted 2(S)-amino-6-boronohexanoic acid (ABH) derivatives as hARGI inhibitors was studied by using a three-dimensional quantitative structure–activity relationships (3D-QSAR) method. The predictivity of the obtained 3D-QSAR model was demonstrated by using internal and external validation experiments. The best model revealed that the differential hARGI inhibitory activities of the ABH derivatives can be described by using steric and electrostatic fields; the local effects of these fields in the activity are presented. In addition, binding modes of the above-mentioned compounds inside the hARGI binding site were obtained by using molecular docking. It was found that ABH derivatives adopted the same orientation reported for ABH within the hARGI active site, with the substituents at Cα exposed to the solvent with interactions with residues at the entrance of the binding site. The hARGI residues involved in chemical interactions with inhibitors were identified by using an interaction fingerprints (IFPs) analysis.

This study describes the in-silico design, synthesis, and evaluation of a cross-linked PVA hydrogel (CLPH) for the absorption of organophosphorus pesticide dimethoate from aqueous solutions. The crosslinking effectiveness of 14 dicarboxilic acids was evaluated through in-silico studies using semiempirical quantum mechanical calculations. According to the theoretical studies, the nanopore of PVA cross-linked with malic acid (CLPH-MA) showed the best interaction energy with dimethoate. Later, using all-atom molecular dynamics simulations, three hydrogels with different proportions of PVA:MA (10:2, 10:4, and 10:6) were used to evaluate their interactions with dimethoate. These results showed that the suitable crosslinking degree for improving the affinity for the pesticide was with 20% (W%) of the cross-linker. In the experimental absorption study, the synthesized CLPH-MA20 recovered 100% of dimethoate from aqueous solutions. Therefore, the theoretical data were correlated with the experimental studies. Surface morphology of CLPH-MA20 by Scanning Electron Microscopy (SEM) was analyzed. In conclusion, the ability of CLPH-MA20 to remove dimethoate could be used as a technological alternative for the treatment of contaminated water.

Most snake venom toxins are proteins, and participate to envenomation through a diverse array of bioactivities, such as bleeding, inflammation, and pain, cytotoxic, cardiotoxic or neurotoxic effects. The venom of a single snake species contains hundreds of toxins, and the venoms of the 725 species of venomous snakes represent a large pool of potentially bioactive proteins. Despite considerable discovery efforts, most of the snake venom toxins are still uncharacterized. Modern bioinformatics tools have been recently developed to mine snake venoms, helping focus experimental research on the most potentially interesting toxins. Some computational techniques predict toxin molecular targets, and the binding mode to these targets. This review gives an overview of current knowledge on the ~2200 sequences, and more than 400 three-dimensional structures of snake toxins deposited in public repositories, as well as of molecular modeling studies of the interaction between these toxins and their molecular targets. We also describe how modern bioinformatics have been used to study the snake venom protein phospholipase A2, the small basic myotoxin Crotamine, and the three-finger peptide Mambalgin.

Non-peptidic thrombin inhibitors (TIs; 177 compounds) with diverse groups at motifs P1 (such as oxyguanidine, amidinohydrazone, amidine, amidinopiperidine), P2 (such as cyanofluorophenylacetamide, 2-(2-chloro-6-fluorophenyl)acetamide), and P3 (such as phenylethyl, arylsulfonate groups) were studied using molecular modeling to analyze their interactions with S1, S2, and S3 subsites of the thrombin binding site. Firstly, a protocol combining docking and three dimensional quantitative structure-activity relationship was performed. We described the orientations and preferred active conformations of the studied inhibitors, and derived a predictive CoMSIA model including steric, donor hydrogen bond, and acceptor hydrogen bond fields. Secondly, the dynamic behaviors of some selected TIs (compounds 26, 133, 147, 149, 162, and 177 in this manuscript) that contain different molecular features and different activities were analyzed by creating the solvated models and using molecular dynamics (MD) simulations. We used the conformational structures derived from MD to accomplish binding free energetic calculations using MM-GBSA. With this analysis, we theorized about the effect of van der Waals contacts, electrostatic interactions and solvation in the potency of TIs. In general, the contents reported in this article help to understand the physical and chemical characteristics of thrombin-inhibitor complexes.

BackgroundChlorogenic acid is a potent phenolic antioxidant. However, its effect on platelet aggregation, a critical factor in arterial thrombosis, remains unclear. Consequently, chlorogenic acid-action mechanisms in preventing platelet activation and thrombus formation were examined.Methods and resultsChlorogenic acid in a dose-dependent manner (0.1 to 1 mmol/L) inhibited platelet secretion and aggregation induced by ADP, collagen, arachidonic acid and TRAP-6, and diminished platelet firm adhesion/aggregation and platelet-leukocyte interactions under flow conditions. At these concentrations chlorogenic acid significantly decreased platelet inflammatory mediators (sP-selectin, sCD40L, CCL5 and IL-1β) and increased intraplatelet cAMP levels/PKA activation. Interestingly, SQ22536 (an adenylate cyclase inhibitor) and ZM241385 (a potent A2A receptor antagonist) attenuated the antiplatelet effect of chlorogenic acid. Chlorogenic acid is compatible to the active site of the adenosine A2A receptor as revealed through molecular modeling. In addition, chlorogenic acid had a significantly lower effect on mouse bleeding time when compared to the same dose of aspirin.ConclusionsAntiplatelet and antithrombotic effects of chlorogenic acid are associated with the A2A receptor/adenylate cyclase/cAMP/PKA signaling pathway.

We have carried out a protocol in computational biochemistry including molecular dynamics (MD) simulations and MM/GBSA free energy calculations on the complex between the protein kinase A (PKA) and the specific peptide substrate Kemptide (LRRASLG). We made the same calculations on other PKA complexes that contain Kemptide derivatives (with mutations of the arginines, and with deletions of N and C-terminal amino acids). We predicted shifts in the free energy changes from the free PKA to PKA-substrate complex (ΔΔG(E→ES)) when Kemptide structure is modified (we consider that the calculated shifts correlate with the experimental shifts of the free energy changes from the free PKA to the transition states (ΔΔG(E→TS)) determined by the catalytic efficiency (k(cat)/K(M)) changes). Our results demonstrate that it is possible to predict the kinetic properties of protein kinases using simple computational biochemistry methods. As an additional benefit, these methods give detailed molecular information that permit the analysis of the atomic forces that contribute to the affinity between protein kinases and their substrates.

The inhibitory effect of adenosine on platelet aggregation is abrogated after the addition of adenosine-deaminase. Inosine is a naturally occurring nucleoside degraded from adenosine.The mechanisms of antiplatelet action of adenosine and inosine in vitro and in vivo, and their differential biological effects by molecular modeling were investigated.Adenosine (0.5, 1 and 2 mmol/L) inhibited phosphatidylserine exposure from 52±4% in the control group to 44±4 (p

Many protein kinase (PK) inhibitors have been reported in recent years, but only a few have been approved for clinical use. The understanding of the available molecular information using computational tools is an alternative to contribute to this process. With this in mind, we studied the binding modes of 77 maleimide derivates inside the PK glycogen synthase kinase 3 beta (GSK3β) using docking experiments. We found that the orientations that these compounds adopt inside GSK3β binding site prioritize the formation of hydrogen bond (HB) interactions between the maleimide group and the residues at the hinge region (residues Val135 and Asp133), and adopt propeller-like conformations (where the maleimide is the propeller axis and the heterocyclic substituents are two slanted blades). In addition, quantitative structure-activity relationship (QSAR) models using CoMSIA methodology were constructed to explain the trend of the GSK3β inhibitory activities for the studied compounds. We found a model to explain the structure-activity relationship of non-cyclic maleimide (NCM) derivatives (54 compounds). The best CoMSIA model (training set included 44 compounds) included steric, hydrophobic, and HB donor fields and had a good Q(2) value of 0.539. It also predicted adequately the most active compounds contained in the test set. Furthermore, the analysis of the plots of the steric CoMSIA field describes the elements involved in the differential potency of the inhibitors that can be considered for the selection of suitable inhibitors.

Liquid chromatography coupled with mass spectrometry is an outstanding methodology for fast analysis of phenolic compounds in biological samples. Twenty two compounds were quickly and accurately identified in the methanolic extract of the Antarctic lichen Ramalina terebrata for the first time using ultra high pressure liquid chromatography coupled with photodiode array detector and high resolution mass spectrometry (UHPLC-PDA-Q/Orbitrap/MS/MS). In addition, the extract and the four compounds isolated from this species were tested for the inhibitory activity of tau protein aggregation, which is a protein involved in Alzheimer’s disease (AD). All compounds showed null activity with the exception of parietin, which it was able to inhibit aggregation process of tau in a concentration range between 3 µg/mL (10 µM) to 28 µg/mL (100 µM). In addition, we show how parietin interact with tau 306VQIVYK311 hexapeptide inside of the microtubule binding domain (4R) with the help of molecular docking experiments. Finally, the constituents present in the methanolic extract could possibly contribute to the established anti-aggregation activity for this extract and this in-depth analysis of the chemical composition of R. terebrata could guide further research into its medicinal properties and potential uses.

Molecular docking is a computational chemistry method which has become essential for the rational drug design process. In this context, it has had great impact as a successful tool for the study of ligand–receptor interaction modes, and for the exploration of large chemical datasets through virtual screening experiments. Despite their unquestionable merits, docking methods are not reliable for predicting binding energies due to the simple scoring functions they use. However, comparisons between two or three complexes using the predicted binding energies as a criterion are commonly found in the literature. In the present work we tested how wise is it to trust the docking energies when two complexes between a target protein and enantiomer pairs are compared. For this purpose, a ligand library composed by 141 enantiomeric pairs was used, including compounds with biological activities reported against seven protein targets. Docking results using the software Glide (considering extra precision (XP), standard precision (SP), and high-throughput virtual screening (HTVS) modes) and AutoDock Vina were compared with the reported biological activities using a classification scheme. Our test failed for all modes and targets, demonstrating that an accurate prediction when binding energies of enantiomers are compared using docking may be due to chance. We also compared pairs of compounds with different molecular weights and found the same results.