Visualizing protein-ligand binding with chemical energy-wise decomposition (CHEWD): application to ligand binding in the kallikrein-8 S1 Site

Anti-leishmanial activity of Eleutherine plicata Herb. and predictions of isoeleutherin and its analogues

Introduction: Leishmaniasis is caused by protozoa of the genus Leishmania, classified as tegumentary and visceral. The disease treatment is still a serious problem, due to the toxic effects of available drugs, the costly treatment and reports of parasitic resistance, making the search for therapeutic alternatives urgent. This study assessed the in vitro anti-leishmanial potential of the extract, fractions, and isoeleutherin from Eleutherine plicata, as well as the in silico interactions of isoeleutherin and its analogs with Trypanothione Reductase (TR), in addition to predicting pharmacokinetic parameters. Methods: From the ethanolic extract of E. plicata (EEEp) the dichloromethane fraction (FDEp) was obtained, and isoeleutherin isolated. All samples were tested against promastigotes, and parasite viability was evaluated. Isoeleutherin analogues were selected based on similarity in databases (ZINK and eMolecules) to verify the impact on structural change. Results and Discussion: The extract and its fractions were not active against the promastigote form (IC50 > 200 μg/mL), while isoeleutherin was active (IC50 = 25 μg/mL). All analogues have high intestinal absorption (HIA), cell permeability was moderate in Caco2 and low to moderate in MDCK. Structural changes interfered with plasma protein binding and blood-brain barrier permeability. Regarding metabolism, all molecules appear to be CYP3A4 metabolized and inhibited 2-3 CYPs. Molecular docking and molecular dynamics assessed the interactions between the most stable configurations of isoeleutherin, analogue compound 17, and quinacrine (control drug). Molecular dynamics simulations demonstrated stability and favorable interactions with TR. In summary, fractionation contributed to antileishmanial activity and isoleutherin seems to be promising. Structural alterations did not contribute to improve pharmacokinetic aspects and analogue 17 proved to be more promising than isoeleutherin, presenting better stabilization in TR.

(-)-5-Demethoxygrandisin B a New Lignan from Virola surinamensis (Rol.) Warb. Leaves: Evaluation of the Leishmanicidal Activity by In Vitro and In Silico Approaches

Leishmaniasis is a complex disease caused by infection with different Leishmania parasites. The number of medications used for its treatment is still limited and the discovery of new drugs is a valuable approach. In this context, here we describe the in vitro leishmanicidal activity and the in silico interaction between trypanothione reductase (TryR) and (-)-5-demethoxygrandisin B from the leaves of Virola surinamensis (Rol.) Warb. The compound (-)-5-demethoxygrandisin B was isolated from V. surinamensis leaves, a plant found in the Brazilian Amazon, and it was characterized as (7R,8S,7'R,8'S)-3,4,5,3',4'-pentamethoxy-7,7'-epoxylignan. In vitro antileishmanial activity was examined against Leishmania amazonensis, covering both promastigote and intracellular amastigote phases. Cytotoxicity and nitrite production were gauged using BALB/c peritoneal macrophages. Moreover, transmission electron microscopy was applied to probe ultrastructural alterations, and flow cytometry assessed the shifts in the mitochondrial membrane potential. In silico methods such as molecular docking and molecular dynamics assessed the interaction between the most stable configuration of (-)-5-demethoxygrandisin B and TryR from L. infantum (PDB ID 2JK6). As a result, the (-)-5-demethoxygrandisin B was active against promastigote (IC50 7.0 µM) and intracellular amastigote (IC50 26.04 µM) forms of L. amazonensis, with acceptable selectivity indexes. (-)-5-demethoxygrandisin B caused ultrastructural changes in promastigotes, including mitochondrial swelling, altered kDNA patterns, vacuoles, vesicular structures, autophagosomes, and enlarged flagellar pockets. It reduced the mitochondria membrane potential and formed bonds with important residues in the TryR enzyme. The molecular dynamics simulations showed stability and favorable interaction with TryR. The compound targets L. amazonensis mitochondria via TryR enzyme inhibition.

An in silico molecular modeling approach of halolactone derivatives as potential inhibitors for human immunodeficiency virus type-1 reverse transcriptase enzyme

Acquired Immune Deficiency Syndrome (AIDS) is an infectious disease caused by Human Immunodeficiency Virus (HIV) infection and its replication requires the Reverse Transcriptase (RT) enzyme. RT plays a key role in the HIV life cycle, making it one of the most important targets for designing new drugs. Thus, in order to increase therapeutic options against AIDS, halolactone derivatives (D-halolactone) that have been showed as potential non-nucleoside inhibitors of the RT enzyme were studied. In the present work, a series of D-halolactone were investigated by molecular modeling studies, combining Three-dimensional Quantitative Structure-Activity Relationship (3 D-QSAR), molecular docking and Molecular Dynamics (MD) techniques, to understand the molecular characteristics that promote biological activity. The internal and external validation parameters indicated that the 3 D-QSAR model has good predictive capacity and statistical significance. Contour maps provided useful information on the structural characteristics of compounds for anti-HIV-1 activity. The docking results showed that D-halolactone present good complementarity by the RT allosteric site. In MD simulations it was observed that the formation of enzyme-ligand complexes were favorable, and from the free energy decomposition it was found that Leu100, Val106, Tyr181, Try188, and Trp229 are key residues for stabilization in the enzymatic site. Thus, the results showed that the proposed models can be used to design promising HIV-1 RT inhibitors. Communicated by Ramaswamy H. Sarma.

Exploring the druggability of the binding site of aurovertin, an exogenous allosteric inhibitor of FOF1-ATP synthase

In addition to playing a central role in the mitochondria as the main producer of ATP, F_OF₁-ATP synthase performs diverse key regulatory functions in the cell membrane. Its malfunction has been linked to a growing number of human diseases, including hypertension, atherosclerosis, cancer, and some neurodegenerative, autoimmune, and aging diseases. Furthermore, inhibition of this enzyme jeopardizes the survival of several bacterial pathogens of public health concern. Therefore, F_OF₁-ATP synthase has emerged as a novel drug target both to treat human diseases and to combat antibiotic resistance. In this work, we carried out a computational characterization of the binding sites of the fungal antibiotic aurovertin in the bovine F₁ subcomplex, which shares a large identity with the human enzyme. Molecular dynamics simulations showed that although the binding sites can be described as preformed, the inhibitor hinders inter-subunit communications and exerts long-range effects on the dynamics of the catalytic site residues. End-point binding free energy calculations revealed hot spot residues for aurovertin recognition. These residues were also relevant to stabilize solvent sites determined from mixed-solvent molecular dynamics, which mimic the interaction between aurovertin and the enzyme, and could be used as pharmacophore constraints in virtual screening campaigns. To explore the possibility of finding species-specific inhibitors targeting the aurovertin binding site, we performed free energy calculations for two bacterial enzymes with experimentally solved 3D structures. Finally, an analysis of bacterial sequences was carried out to determine conservation of the aurovertin binding site. Taken together, our results constitute a first step in paving the way for structure-based development of new allosteric drugs targeting F_OF₁-ATP synthase sites of exogenous inhibitors.

Pharmacophore-based virtual screening from phytocannabinoids as antagonist r-CB1

Search for new pharmacological alternatives for obesity is based on the design and development of compounds that can aid in weight loss so that they can be used safely and effectively over a long period while maintaining their function. The endocannabinoid system is related to obesity by increasing orexigenic signals and reducing satiety signals. Cannabis sativa is a medicinal plant of polypharmaceutical potential that has been widely studied for various medicinal purposes. The in silico evaluation of their natural cannabinoids (also called phytocannabinoids) for anti-obesity purpose stems from the existence of synthetic cannabinoid compounds that have already presented this result, but which did not guarantee patient safety. In order to find new molecules from C. sativa phytocannabinoids, with the potential to interact peripherally with the pharmacological target cannabinoid receptor 1, a pharmacophore-based virtual screening was performed, including the evaluation of physicochemical, pharmacokinetic, toxicological predictions and molecular docking. The results obtained from the ZINC¹² database pointed to Zinc 69 (ZINC33053402) and Zinc 70 (ZINC19084698) molecules as promising anti-obesity agents. Molecular dynamics (MD) studies disclose that both complexes were stable by analyzing the RMSD (root mean square deviation) values, and the binding free energy values demonstrate that the selected structures can interact and inhibit their catalytic activity.

Docking and Molecular Dynamic of Microalgae Compounds as Potential Inhibitors of Beta-Lactamase

Bacterial resistance is responsible for a wide variety of health problems, both in children and adults. The persistence of symptoms and infections are mainly treated with β-lactam antibiotics. The increasing resistance to those antibiotics by bacterial pathogens generated the emergence of extended-spectrum β-lactamases (ESBLs), an actual public health problem. This is due to rapid mutations of bacteria when exposed to antibiotics. In this case, β-lactamases are enzymes used by bacteria to hydrolyze the beta-lactam rings present in the antibiotics. Therefore, it was necessary to explore novel molecules as potential β-lactamases inhibitors to find antibacterial compounds against infection caused by ESBLs. A computational methodology based on molecular docking and molecular dynamic simulations was used to find new microalgae metabolites inhibitors of β-lactamase. Six 3D β-lactamase proteins were selected, and the molecular docking revealed that the metabolites belonging to the same structural families, such as phenylacridine (4-Ph), quercetin (Qn), and cryptophycin (Cryp), exhibit a better binding score and binding energy than commercial clinical medicine β-lactamase inhibitors, such as clavulanic acid, sulbactam, and tazobactam. These results indicate that 4-Ph, Qn, and Cryp molecules, homologous from microalgae metabolites, could be used, likely as novel β-lactamase inhibitors or as structural templates for new in-silico pharmaceutical designs, with the possibility of combatting β-lactam resistance.

Potential implications of vouacapan compounds for insecticidal activity: an in silico study

BACKGROUND

From the fruits and seeds of the species of Pterodon, it is possible to obtain two main products: the essential oil and oleoresin. In oleoresin, numerous vouacapan compounds have been demonstrated to have biological potential, including insecticidal activity.

OBJECTIVE

In silico studies were performed to identify potential candidates for natural insecticides among the vouacapans present in the genus Pterodon.

MATERIALS AND METHODS

Molecular docking and molecular dynamics studies were performed to analyze the interaction of vouacapan compounds with acetylcholinesterase of Drosophila melanogaster. Pharmacokinetic parameters regarding physicochemical properties, plasma protein binding, and activity in the central nervous system were evaluated. The toxicological properties of the selected molecules were predicted using Malathion as the reference compound.

RESULTS

6α,7β-dimethoxivouacapan-17-ene (15) showed a high number of interactions and scores in molecular docking studies. This result suggests that this compound exhibits an inhibitory activity of the enzyme acetylcholinesterase. Regarding physicochemical properties, this compound showed the best results, besides presenting low cutaneous permeability values, suggesting null absorption. Molecular dynamics studies demonstrated few conformational changes in the structure of the complex formed by compound 4 and acetylcholinesterase enzyme throughout the simulation time.

CONCLUSION

It was determined that compound 4 (vouacapan 6α,7β,17β,19-tetraol) is an excellent candidate for usage as a natural insecticide.

Exploration of the structural requirements of Aurora Kinase B inhibitors by a combined QSAR, modelling and molecular simulation approach

Aurora kinase B plays an important role in the cell cycle to orchestrate the mitotic process. The amplification and overexpression of this kinase have been implicated in several human malignancies. Therefore, Aurora kinase B is a potential drug target for anticancer therapies. Here, we combine atom-based 3D-QSAR analysis and pharmacophore model generation to identify the principal structural features of acylureidoindolin derivatives that could potentially be responsible for the inhibition of Aurora kinase B. The selected CoMFA and CoMSIA model showed significant results with cross-validation values (q²) of 0.68, 0.641 and linear regression values (r²) of 0.971, 0.933 respectively. These values support the statistical reliability of our model. A pharmacophore model was also generated, incorporating features of reported crystal complex structures of Aurora kinase B. The pharmacophore model was used to screen commercial databases to retrieve potential lead candidates. The resulting hits were analyzed at each stage for diversity based on the pharmacophore model, followed by molecular docking and filtering based on their interaction with active site residues and 3D-QSAR predictions. Subsequently, MD simulations and binding free energy calculations were performed to test the predictions and to characterize interactions at the molecular level. The results suggested that the identified compounds retained the interactions with binding residues. Binding energy decomposition identified residues Glu155, Trp156 and Ala157 of site B and Leu83 and Leu207 of site C as major contributors to binding affinity, complementary to 3D-QSAR results. To best of our knowledge, this is the first comparison of WaterSwap field and 3D-QSAR maps. Overall, this integrated strategy provides a basis for the development of new and potential AK-B inhibitors and is applicable to other protein targets.

Drug repurposing and computational modeling for discovery of inhibitors of the main protease (Mpro) of SARS-CoV-2

The main protease (M^pro or 3CL^pro) is a conserved cysteine protease from the coronaviruses and started to be considered an important drug target for developing antivirals, as it produced a deadly outbreak of COVID-19. Herein, we used a combination of drug reposition and computational modeling approaches including molecular docking, molecular dynamics (MD) simulations, and the calculated binding free energy to evaluate a set of drugs in complex with the M^pro enzyme. Particularly, our results show that darunavir and triptorelin drugs have favorable binding free energy (−63.70 and −77.28 kcal mol⁻¹, respectively) in complex with the M^pro enzyme. Based on the results, the structural and energetic features that explain why some drugs can be repositioned to inhibit M^pro from SARS-CoV-2 were exposed. These features should be considered for the design of novel M^pro inhibitors.

Structural and energetic features explain why some drugs can be repositioned to inhibit M^pro from SARS-CoV-2.

A potential interaction between the SARS-CoV-2 spike protein and nicotinic acetylcholine receptors

Changeux et al. (Changeux et al. C. R. Biol. 343:33-39.) recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs) and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs and of the spike to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4β2 and α7 subtypes and the muscle-like αβγδ receptor from Tetronarce californica. Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behavior of the bound Y674-R685 is highly dependent on the receptor subtype; it adopts extended conformations in the α4β2 and α7 complexes but is more compact when bound to the muscle-like receptor. In the α4β2 and αβγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket in which it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1, and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a nonparallel arrangement to one another.

Simulations support the interaction of the SARS-CoV-2 spike protein with nicotinic acetylcholine receptors

Changeux et al. recently suggested that the SARS-CoV-2 spike (S) protein may interact with nicotinic acetylcholine receptors (nAChRs). Such interactions may be involved in pathology and infectivity. Here, we use molecular simulations of validated atomically detailed structures of nAChRs, and of the S protein, to investigate this 'nicotinic hypothesis'. We examine the binding of the Y674-R685 loop of the S protein to three nAChRs, namely the human α4β2 and α7 subtypes and the muscle-like αβγδ receptor from Tetronarce californica. Our results indicate that Y674-R685 has affinity for nAChRs and the region responsible for binding contains the PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. In particular, R682 has a key role in the stabilisation of the complexes as it forms interactions with loops A, B and C in the receptor's binding pocket. The conformational behaviour of the bound Y674-R685 region is highly dependent on the receptor subtype, adopting extended conformations in the α4β2 and α7 complexes and more compact ones when bound to the muscle-like receptor. In the α4β2 and αβγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket where it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1 and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of the simulations of the full-length S protein show that the Y674-R685 region is accessible for binding, and suggest a potential binding orientation of the S protein with nAChRs.

Computational insight into the selective allosteric inhibition for PTP1B versus TCPTP: a molecular modelling study

All over the world, diabetes mellitus type 2 has spread as a problematic pandemic. Despite currently available treatments, approved drugs still show undesirable side effects and loss of efficacy or target symptoms instead of causes. Protein tyrosine phosphatase 1B (PTP1B), since its discovery, has emerged as a very promising target against this disease. Although the information regarding the enzyme is immense, little is known about the selectivity between this enzyme and its closest homologue, lymphocyte T tyrosine phosphatase (TCPTP), which is responsible for complicated side effects. In this study, on the basis of different computational approaches, we are able to highlight the importance of a phenylalanine residue located in PTP1B, but not in TCPTP, as a crucial hotspot that causes selectivity and stability for the whole ligand bound system. These results not only allow to explain the selectivity determinants of PTP1B but also provide a useful guide for the design of new allosteric inhibitors. Communicated by Ramaswamy H. Sarma.

Molecular Dynamics Simulation Framework to Probe the Binding Hypothesis of CYP3A4 Inhibitors

The Cytochrome P450 family of heme-containing proteins plays a major role in catalyzing phase I metabolic reactions, and the CYP3A4 subtype is responsible for the metabolism of many currently marketed drugs. Additionally, CYP3A4 has an inherent affinity for a broad spectrum of structurally diverse chemical entities, often leading to drug-drug interactions mediated by the inhibition or induction of the metabolic enzyme. The current study explores the binding of selected highly efficient CYP3A4 inhibitors by docking and molecular dynamics (MD) simulation protocols and their binding free energy calculated using the WaterSwap method. The results indicate the importance of binding pocket residues including Phe57, Arg105, Arg106, Ser119, Arg212, Phe213, Thr309, Ser312, Ala370, Arg372, Glu374, Gly481 and Leu483 for interaction with CYP3A4 inhibitors. The residue-wise decomposition of the binding free energy from the WaterSwap method revealed the importance of binding site residues Arg106 and Arg372 in the stabilization of all the selected CYP3A4-inhibitor complexes. The WaterSwap binding energies were further complemented with the MM(GB/PB)SA results and it was observed that the binding energies calculated by both methods do not differ significantly. Overall, our results could guide towards the use of multiple computational approaches to achieve a better understanding of CYP3A4 inhibition, subsequently leading to the design of highly specific and efficient new chemical entities with suitable ADMETox properties and reduced side effects.

A Visualization Tool for Cryo-EM Protein Validation with an Unsupervised Machine Learning Model in Chimera Platform

Background: Cryo-electron microscopy (cryo-EM) has become a major technique for protein structure determination. However, due to the low quality of cryo-EM density maps, many protein structures derived from cryo-EM contain outliers introduced during the modeling process. The current protein model validation system lacks identification features for cryo-EM proteins making it not enough to identify outliers in cryo-EM proteins. Methods: This study introduces an efficient unsupervised outlier detection model for validating protein models built from cryo-EM technique. The current model uses a high-resolution X-ray dataset (<1.5 Å) as the reference dataset. The distal block distance, side-chain length, phi, psi, and first chi angle of the residues in the reference dataset are collected and saved as a database of the histogram-based outlier score (HBOS). The HBOS value of the residues in target cryo-EM proteins can be read from this HBOS database. Results: Protein residues with a HBOS value greater than ten are labeled as outliers by default. Four datasets containing proteins derived from cryo-EM density maps were tested with this probabilistic anomaly detection model. Conclusions: According to the proposed model, a visualization assistant tool was designed for Chimera, a protein visualization platform.