In silico analysis of SARS-CoV-2 papain-like protease potential inhibitors

The emergent outbreak caused by severe acute respiratory syndrome coronavirus 2 continues spreading and causing huge social and economic disruption. Papain-like protease (PLpro) has a crucial role in the cleavage of viral polyproteins, and disruption of host responses. PLpro is considered an important goal for the development of SARS-CoV-2 inhibitors. ZINC101291108 (lead 1) and ZINC16449029 (lead 2) were identified as potent SARS-CoV-2 PLpro inhibitors with IC₅₀ values of 0.085 μM and 0.063 μM, respectively. Molecular dynamics simulations (MD) were carried out for lead 1, 2 and several reported SARS-CoV-2 inhibitors. Analysis results of the simulations confirmed the stability of both compounds and showed that they adopted two confirmations along the simulation period. The per-residue decomposition results revealed that the key residues involved in inhibitor binding were E167, P247, P248, Y264, Y268 and Q269. H-bond analyses showed H-bonds with G266 and N267 and salt bridges with G209 and Y273, which are essential for strengthening the substrate-binding pocket. Both inhibitors showed hydrophobic interactions with the S4 site and BL2 loop residues. The RMSD of the BL2 loop with the two inhibitors was investigated, and the results showed that the Y268 and Q269 BL2 loop residues moved outward to accommodate the large size of lead 2. The van der Waals interaction was the main energy contribution that stabilized lead 2, while van der Waals and electrostatic interactions were the main energy contributions stabilizing lead 1. Rational design strategies were suggested to replace the 2-(2-hydroxybenzylidene) hydrazine moiety with naphthalene or nitrobenzene at the P4 position of lead 2 and introduce polar substituents as aniline and benzoate groups at position P1 to enhance hydrophobic interactions and H-bonds, respectively.

The emergent outbreak caused by severe acute respiratory syndrome coronavirus 2 continues spreading and causing huge social and economic disruption.     

In silico study on identification of novel MALT1 allosteric inhibitors

Mucosa-associated lymphoid tissue lymphoma translocation protein 1 (MALT1), which plays a crucial role in the nuclear factor-kappa B (NF-κB) activation signaling pathway as a paracaspase, is a new target for immunomodulatory and antitumor drugs. Here, novel inhibitors that target MALT1 allosteric sites were identified by virtual screening FDA-approved drug databases. Paliperidone, a compound that binds to the allosteric site of MALT1, is investigated. An in vitro study found that the proteolytic activity of MALT1 substrate cleavage was blocked by paliperidone. Meanwhile, the MALT1 proteolytic activity was reversible, as demonstrated by the partial recovery of the MALT1 substrate cleavage following compound wash out. The docking analysis of the interaction of MALT1 and paliperidone suggested that two hydrogen bonds formed in the allosteric pocket of MALT1. MALT1 and paliperidone achieved a good equilibrium, as demonstrated by 100 ns molecular dynamic (MD) simulations conducted with the program Gromacs. However, the catalytically active site of the MALT1 complex with paliperidone remained in an inactive conformation. Thus, paliperidone, a noncompetitive and allosteric inhibitor, was screened through in silico and in vitro methods. This study will be of significance for the development of effective and selective drugs that can treat MALT1-driven cancer or autoimmune diseases.

Paliperidone was screened as an effective and selective drug that can treat MALT1-driven cancer or autoimmune diseases.     

In silico investigation of the role of vitamins in cancer therapy through inhibition of MCM7 oncoprotein

An overabundance of MCM7 protein, a component of the minichromosome maintenance complex that normally initiates DNA replication, has been reported to cause different types of cancers with aggressive malignancy. Inhibition of MCM7 may lead to a significant reduction in cancer-associated cell proliferation. Despite such significance of MCM7 in cancer, the protein structure is yet to be resolved experimentally. This significantly halts the structure-guided ligand designing for cancer therapy targeting the MCM7. The present study aims to resolve the tertiary structure of MCM7 and repurpose the FDA-approved clinically used drugs for cancer therapy by targeting MCM7 protein. The secondary and 3D structures of MCM7 were generated using multiple bioinformatics tools, including the Self-Optimized Prediction Method with Alignment (SOPMA), SWISS-MODEL, and I-TASSER. The reliability of the modeled structure was assessed using PROCHECK. Initially, a structure-guided virtual screening was performed on the approved drug library to identify potential hits against MCM7. The detailed molecular mechanism of receptor interactions of the identified hits was evaluated using extensive molecular dynamics simulation. The results from this study reveal an intriguing discovery of the potential of ergocalciferol (vitamin D2), cholecalciferol (vitamin D3), ergosterol (precursor of vitamin D2) and menaquinone (vitamin K2) as oncoprotein inhibitors for cancer therapy via inhibition of MCM7.

An in silico drug repurposing strategy combining docking and molecular dynamics simulation identifies the anticancer potential of vitamins targeting the MCM7 protein.     

In vitro and in silico determination of glutaminyl cyclase inhibitors

Alzheimer''s disease (AD) is the most common form of neurodegenerative disease currently. It is widely accepted that AD is characterized by the self-assembly of amyloid beta (Aβ) peptides. The human glutaminyl cyclase (hQC) enzyme is characterized by association with Aβ peptide generation. The development of hQC inhibitors could prevent the self-aggregation of Aβ peptides, resulting in impeding AD. Utilizing structural knowledge of the hQC substrates and known hQC inhibitors, new heterocyclic and peptidomimetic derivatives were synthesized and were able to inhibit the hQC enzyme. The inhibiting abilities of these compounds were evaluated using a fluorometric assay. The binding mechanism at the atomic level was estimated using molecular docking, free energy perturbation, and quantum chemical calculation methods. The predicted log(BBB) and human intestinal absorption values indicated that these compounds are able to permeate the blood–brain barrier and be well-absorbed through the gastrointestinal tract. Overall, 5,6-dimethoxy-N-(3-(5-methyl-1H-imidazol-1-yl)propyl)-1H-benzo[d]imidazol-2-amine (1_2) was indicated as a potential drug for AD treatment.

Rational design of new hQC inhibitors.     

In silico study of inhibitory capacity of sacubitril/valsartan toward neprilysin and angiotensin receptor

Heart failure (HF) is a life-threatening condition that occurs when the heart cannot pump enough blood and oxygen to meet the body''s needs. It affects mostly the elderly, commonly from the male population, especially those with obesity, diabetes, or some other chronic condition. It can be treated with different medications, and promising results were shown by a relatively new medicament called Entresto. Results obtained from molecular docking and molecular dynamics simulations to examine the inhibitory capacity of Entresto are presented in this study. Parameters obtained by the molecular docking simulations show that both parts of Entresto (sacubitril (SAC) and valsartan (VAL)) interact with targeted proteins, and inhibit their physiological function. Simulations of molecular dynamics revealed some interesting inhibitory patterns. SAC was discovered to produce structural alterations in neprilysin by binding to it, reducing neprilysin''s physiological activity. In addition to blocking the active site, SAC binding causes the enzyme''s structure to become less compact over time, causing changes in its biochemical characteristics and preventing the enzyme from performing its biological function. Similar to SAC, VAL also causes deviations in the structure of angiotensin receptors. The angiotensin receptor GPCR (G-protein-coupled receptors) is immersed in the lipid bilayer, and changes in the tertiary structure are only visible through RMSD and RMSF, not by examining R_g. In this regard, MD simulations validated the results of molecular docking simulations, demonstrating that both SAC and VAL had inhibitory potential towards the neprilysin and angiotensin receptors, respectively.

The article presents results obtained from molecular docking and molecular dynamic simulations to examine the inhibitory capacity of Entresto.     

In silico study of natural compounds from sesame against COVID-19 by targeting Mpro,PLpro and RdRp

Natural products and traditional medicine products with known safety profiles are a promising source for the discovery of new drug leads. Natural products as sesame were reported to exhibit potential to protect from COVID-19 disease. In our study, the total methanolic extract of Sesamum indicum L. seeds (sesame) were led to isolation of seven known compounds, five lignan; sesamin 1, sesamolin 2, pinoresinol 3, hydroxymatairesinol 6, spicatolignan 7, together with two simple phenolic compounds; ferulic acid 4 and vanillic acid 5. All isolated compounds were evaluated in silico against three important SARS-CoV-2 protein targets; main protease (M^pro), papain-like protease (PL^pro) and RNA-dependent RNA polymerase (RdRp) which possessed crucial role in replication and proliferation of the virus inside the human cell. The results revealed that compound 6 has the high affinity against the three main proteins, specially towards the SARS-CoV-2 M^pro that exceeded the currently used SARS-CoV-2 M^pro inhibitor darunavir as well as, exhibiting a similar binding energy at SARS CoV-2 PLpro when compared with the co-crystallized ligand. This activity continued to include the RdRp as it displayed a comparable docking score with remdesivir. Inferiorly, compounds 1 and 2 showed also similar triple inhibitory effect against the three main proteins while compound 7 exhibited a dual inhibitory effect against SARS CoV-2 PL^Pro and RdRp. Further molecular dynamic simulation experiments were performed to validate these docking experiments and to calculate their binding free energies (ΔGs). Compounds 1, 2, 3, 6, and 7 showed comparable binding stability inside the active site of each enzyme with ΔG values ranged from −4.9 to −8.8 kcal mol⁻¹. All the compounds were investigated for their ADME and drug likeness properties, which showed acceptable ADME properties and obeying Lipinski''s rule of five parameters. It can be concluded that the isolated compounds from sesame lignans could be an alternative source for the development of new natural leads against COVID-19.

Natural products and traditional medicine products with known safety profiles are a promising source for the discovery of new drug leads.     

In silico prediction and interaction of resveratrol on methyl-CpG binding proteins by molecular docking and MD simulations study

Resveratrol enhances the BRCA1 gene expression and the MBD family of proteins bind to the promoter region of the BRCA1 gene. However, the molecular interaction is not yet reported. Here we have analyzed the binding affinity of resveratrol with MBD proteins. Our results suggest that resveratrol binds to the MBD proteins with higher binding affinity toward MeCP2 protein (ΔG = −6.5) by sharing four hydrogen bonds as predicted by molecular docking studies. Further, the molecular dynamics simulations outcomes showed that the backbones of all three protein–ligand complexes are stabilized after the period of 75 ns, constantly fluctuating around the deviations of 0.4 Å, 0.5 Å and 0.7 Å for MBD1, MBD2 and MeCP2, respectively. The inter-molecular hydrogen bonding trajectory analysis for protein–ligand complexes also support the strong binding between MeCP2–resveratrol complex. Further, binding free energy calculations showed binding energy of −94.764 kJ mol⁻¹, −53.826 kJ mol⁻¹ and −36.735 kJ mol⁻¹ for MeCP2–resveratrol, MBD2–resveratrol and MBD1–resveratrol complexes, respectively, which also supported our docking results. Our study also highlighted that the MBD family of proteins forms a binding interaction with other signaling proteins that are involved in various cancer initiation pathways.

Resveratrol enhances the BRCA1 gene expression and the MBD family of proteins bind to the promoter region of the BRCA1 gene.     

In silico design of mimosine containing peptides as new efficient chelators of aluminum

The design of new and efficient chelators that can remove aluminium(iii), a metal with increasing recognition as a potential toxic agent, from biological systems is an area of high therapeutic relevance. In the present paper, we present an extensive computational study of a new promising type of these chelators based on mimosine containing peptides. The reason to choose mimosine is that the sidechain of this residue is similar to deferiprone, a ligand known to tightly interact with highly-valent metals, and in particular with Al(iii). In this article we analyze systematically, using a combination of methods that include QM/MM MD simulations, how the size and sequence of the polypeptides can alter the fundamental binding patterns to aluminum, in comparison with the binding to deferiprone. Particular attention is given towards the identification of the smallest peptide that interacts efficiently with aluminum, since polypeptide size is a fundamental factor to allow a given polypeptide to efficiently cross the cell membrane. The results indicate that the longest peptides, with 8 or 9 amino acids, show no difficulties interacting with Al(iii) in an optimum arrangement. By contrast, when the peptide contains five or six amino acids Al(iii) is pentacoordinated, reducing the stability of the resultant complex. In summary, our study demonstrates that the mimosine containing peptides can efficiently coordinate highly valent metals such as Al(iii), with a subtle dependence of the binding on the specific chain-lengths of the polypeptide. We believe that the present study sheds light on the adequacy of this new type of chelator towards aluminum binding.

A novel chelator of aluminum is presented, a peptide containing three mimosine residues.     

In vitro and in silico studies of SARS-CoV-2 main protease Mpro inhibitors isolated from Helichrysum bracteatum

Discovering SARS-CoV-2 inhibitors from natural sources is still a target that has captured the interest of many researchers. In this study, the compounds (1–18) present in the methanolic extract of Helichrysum bracteatum were isolated, identified, and their in vitro inhibitory activities against SARS-CoV-2 main protease (M^pro) was evaluated using fluorescence resonance energy transfer assay (FRET-based assay). Based on 1D and 2D spectroscopic techniques, compounds (1–18) were identified as 24-β-ethyl-cholesta-5(6),22(23),25(26)-triene-3-ol (1), α-amyrin (2), linoleic acid (3), 24-β-ethyl-cholesta-5(6),22(23),25(26)-triene-3-O-β-d-glucoside (4), 1,3-propanediol-2-amino-1-(3′,4′-methylenedioxyphenyl) (5), (−)-(7R,8R,8′R)-acuminatolide (6), (+)-piperitol (7), 5,7,4′-trihydroxy-8,3′-dimethoxy flavanone (8), 5,7,4′-trihydroxy-6-methoxy flavanone (9), 4′,5-dihydroxy-3′,7,8-trimethoxyflavone (10), 5,7-dihydroxy-3′,4′,5′,8-tetramethoxy flavone (11), 1,3-propanediol-2-amino-1-(4′-hydroxy-3′-methoxyphenyl) (12), 3′,5′,5,7-tetrahydroxy-6-methoxyflavanone (13), simplexoside (piperitol-O-β-d-glucoside) (14), pinoresinol monomethyl ether-β-d-glucoside (15), orientin (16), luteolin-3′-O-β-d-glucoside (17), and 3,5-dicaffeoylquinic acid (18). Compounds 6, 12, and 14 showed comparable inhibitory activities against SARS-CoV-2 M^pro with IC₅₀ values of 0.917 ± 0.05, 0.476 ± 0.02, and 0.610 ± 0.03 μM, respectively, compared with the control lopinavir with an IC₅₀ value of 0.225 ± 0.01 μM. The other tested compounds showed considerable inhibitory activities. The molecular docking study for the tested compounds was carried out to correlate their binding modes and affinities for the SARS-CoV-2 M^pro enzyme with the in vitro results. Analyzing the results of the in vitro assay together with the obtained in silico results led to the conclusion that phenylpropanoids, lignans, and flavonoids could be considered suitable drug leads for developing anti-COVID-19 therapeutics. Moreover, the phenylpropanoid skeleton oxygenated at C3, C4 of the phenyl moiety and at C1, C3 of the propane parts constitute an essential core of the SARS-CoV-2 M^pro inhibitors, and thus could be proposed as a scaffold for the design of new anti-COVID-19 drugs.

Compounds isolated and identified from Helichrysum bracteatum leaves showed promising in vitro inhibitory activities against SARS-CoV-2 main protease (M^pro). Thus, could be considered suitable drug leads for developing anti-COVID-19 therapeutics.     

In silico functional and tumor suppressor role of hypothetical protein PCNXL2 with regulation of the Notch signaling pathway

Since the last decade, various genome sequencing projects have led to the accumulation of an enormous set of genomic data; however, numerous protein-coding genes still need to be functionally characterized. These gene products are called “hypothetical proteins”. The hypothetical protein pecanex-like protein 2 Homo sapiens (PCNXL2) is found to be mutated in colorectal carcinoma with microsatellite instability; therefore, annotation of the function of PCNXL2 in tumorigenesis is very important. In the present study, bioinformatics analysis of PCNXL2 was performed at the molecular level to assess its role in the progression of cancer for designing new anti-cancer drugs. The retrieved sequence of PCNXL2 was functionally and structurally characterized through the web tools Pfam, Batch CD (conserved domain) search, ExPASy, COACH and I-TASSER directed for pathway analysis and design to explore the intercellular interactions of PCNXL2 involved in cancer development. The present study has shown that PCNXL2 encodes multi-pass transmembrane proteins whose tumor suppressor function may involve regulating Notch signaling by transporting protons across the membrane to provide suitable membrane potential for γ secretase function, which may liberate the Notch intracellular domain NICD from the receptor to inside the cell. Furthermore, domain A of PCNXL2 may exhibit nuclear transport activity of NICD from the cytoplasm to the nucleus through interaction with a nuclear localization signal that may act as an activator for Notch signaling in the nucleus. Conclusively, the tumor suppressor role of PCNXL2 by regulation of the Notch signaling pathway and its functional and structural characteristics are important findings. However, further studies are required to validate the putative role of PCNXL2 as a cancer biomarker in cancer development.

Since the last decade, various genome sequencing projects have led to the accumulation of an enormous set of genomic data; however, numerous protein-coding genes still need to be functionally characterized.     

In silico exploration of lignin peroxidase for unraveling the degradation mechanism employing lignin model compounds

Lignin peroxidase is a heme-containing biocatalyst, well-known for its diverse applications in the fields from environmental chemistry to biotechnology. LiP-mediated oxidative catalysis is H₂O₂-dependent, and can oxidize phenolic, and non-phenolic substrates by oxidative cleavage of the C–C and C–O bonds of lignin. In contrast to fungi-derived LiP, the binding affinity of bacterial-derived LiP to lignin at the molecular level is poorly known to date. Tremendous wet-lab studies have been unveiled that provide degradation and biotransformation information on kraft lignin, whilst studies on the completely transformed compounds and the degradation of each transformed compounds simultaneously during degradation are scarce. To gain an understanding of the degradation process using docking, and MDS based studies, we assessed the binding affinity of selected lignin model compounds with bacterial origin LiP and validated such docked complexes exploiting 30 ns molecular dynamics simulations. We selected and picked a total of 12 lignin model compounds for molecular modeling analysis, namely two chlorinated lignin model compounds (monomer) (2-chlorosyringaldehyde and 5-chlorovanillin), eight standard lignin model compounds (veratryl alcohol, syringyl alcohol, sinapyl alcohol, methyl hydroquinone, guaiacol, coniferyl alcohol, catechol, and 4-methoxy phenol), while, two 4-O-5, and β-O-4 linkage-based multimeric model compounds (dimer: 2-methoxy-6-(2-methoxy-4-methylphenoxy)-4-methylphenol; trimer: syringyl β-O-4 syringyl β-O-4 sinapyl alcohol). Far more specific binding residues were observed from XP-Glide docking, as TYR, HIP (protonated histidine), PHE, VAL, ASP, THR, LYS and GLN. The binding affinity was confirmed by the Gibbs free energy or binding energy (ΔG) score; furthermore, it is found that the maximum binding energy seems to be observed for 4-methoxyphenol with a Glide score of −3.438 with Pi–Pi stacking and H-bond type bonding interactions, whilst the lowest XP Gscore as −8.136 with Pi–Pi stacking and H-bond (side chain) type bonding interactions were found for the trimer model compound. The docked complexes were further evaluated for deep rigorous structural and functional fluctuation analyses through high-performance molecular dynamics simulations-DESMOND, after a post simulation run of 30 ns. The RMSD trajectory analyses of the protein-ligands were found to be in the equilibrium state at the end of simulation run for multimeric lignin model compounds. In addition, ionic ligand–protein interaction occurs among chlorinated compounds, while hydrophobic and H-bond contacts have frequently been observed in all lignin-model compounds. The findings herein demonstrate that bacterial LiP can effectively catalyze multiple lignin model compounds, and it might further be used as an effective tool for sustainable mitigation of diverse environmental contaminants.

The findings herein demonstrate that bacterial LiP can effectively catalyze multiple lignin model compounds, and it might further be used as an effective tool for sustainable mitigation of diverse environmental contaminants.     

In silico identification of potential SARS COV-2 2′-O-methyltransferase inhibitor: fragment-based screening approach and MM-PBSA calculations

In the present era, there are many efforts trying to face the emerging and successive waves of the COVID-19 pandemic. This has led to considering new and unusual targets for SARS CoV-2. 2′-O-Methyltransferase (nsp16) is a key and attractive target in the SARS CoV-2 life cycle since it is responsible for the viral RNA protection via a cap formation process. In this study, we propose a new potential inhibitor for SARS COV-2 2′-O-methyltransferase (nsp16). A fragment library was screened against the co-crystal structure of the SARS COV-2 2′-O-methyltransferase complexed with Sinefungin (nsp16 – PDB ID: 6WKQ), and consequently the best proposed fragments were linked via a de novo approach to build molecule AP-20. Molecule AP-20 displayed a superior docking score to Sinefungin and reproduced the key interactions in the binding site of 2′-O-methyltransferase. Three molecular dynamic simulations of the 2′-O-methyltransferase apo structure and its complexed forms with AP-20 and Sinefungin were performed for 150 nano-seconds to provide insights on the dynamic nature of such setups and to assess the stability of the proposed AP-20/enzyme complex. AP-20/enzyme complex demonstrated better stability for the ligand–enzyme complex compared to Sinefungin in a respective setup. Furthermore, MM-PBSA binding free energy calculations showed a better profile for AP-20/enzyme complex compared to Sinefungin/enzyme complex emphasizing the potential inhibitory effect of AP-20 on SARS COV-2 2′-O-methyltransferase. We endorse our designed molecule AP-20 to be further explored via experimental evaluations to confront the spread of the emerging COVID-19. Also, in silico ADME profiling has ascribed to AP-20 an excellent safety and metabolic stability profile.

The identification of AP-20 as a potential SARS COV-2 2′-O-methyltransferase inhibitor: fragment-based screening approach and MM-PBSA calculations.     

In silico identification of SARS-CoV-2 spike (S) protein–ACE2 complex inhibitors from eight Tecoma species and cultivars analyzed by LC-MS

Coronavirus (CoV) is a positive RNA genome virus causing a global panic nowadays. Tecoma is a medicinally-valuable genus in the Bignoniaceae family, with some of its species exhibiting anti-HIV activity. This encouraged us to conduct an in silico exploration of some phytocompounds in Tecoma species cultivated in Egypt, namely Tecoma capensis and its four varieties i.e. yellow, harmony, pink and red, T. grandiflora Loisel., T. radicans L., and one hybrid i.e. Tecoma × smithii W. Watson. LC/MS-based metabolite profiling of the studied Tecoma plants resulted in the dereplication of 12 compounds (1–12) belonging to different phytochemical classes viz. alkaloids, iridoids, flavonoids and fatty acid esters. The in silico inhibitory action of these compounds against SARS-CoV-2 spike protein C-terminal domain in complex with human ACE2 was assessed via molecular docking. Succinic acid decyl-3-oxobut-2-yl ester (10), a fatty acid ester, possessed the best binding affinity (−6.77 kcal mol⁻¹), as compared to hesperidin (13) (−7.10 kcal mol⁻¹).

In silico exploration of 12 Tecoma phytocompounds that could serve as potential inhibitors of SARS-CoV entry to host cells.     

In silico and in vitro design of cordycepin encapsulation in liposomes for colon cancer treatment

Cordycepin or 3′-deoxyadenosine is an interesting anti-cancer drug candidate that is found in abundance in the fungus Cordyceps militaris. It inhibits cellular growth of many cancers including lung carcinoma, melanoma, bladder cancer, and colon cancer by inducing apoptosis, anti-proliferation, anti-metastasis and by arresting the cell cycle. Cordycepin has, however, poor stability and low solubility in water, resulting in loss of its bioactivity. Liposomes can be used to overcome these obstacles. Our aim is to improve cordycepin''s anti-colon cancer activity by liposome encapsulation. Cordycepin-encapsulated liposomes were designed and fabricated based on a combination of theoretical and experimental studies. Molecular dynamics (MD) simulations and free energy calculations suggest that phosphatidylcholine (PC) lipid environment is favorable for cordycepin adsorption. Cordycepin passively permeates into PC lipid bilayers without membrane damage and strongly binds to the lipids'' polar groups by flipping its deoxyribose sugar toward the bilayer center. Our fabricated liposomes containing 10 : 1 molar ratio of egg yolk PC : cholesterol showed encapsulation efficiency (%EE) of 99% using microfluidic hydrodynamic focusing (MHF) methods. In our in vitro study using the HT-29 colon cancer cell line, cordycepin was able to inhibit growth by induction of apoptosis. Cell viability was significantly decreased below 50% at 125 μg mL⁻¹ dosage after 48 h treatment with non-encapsulated and encapsulated cordycepin. Importantly, encapsulation provided (1) a 2-fold improvement in the inhibition of cancer cell growth at 125 μg mL⁻¹ dosage and (2) 4-fold increase in release time. These in silico and in vitro studies indicate that cordycepin-encapsulated liposomes could be a potent drug candidate for colon cancer therapy.

Cordycepin-encapsulated liposomes could be a potent drug candidate for cancer therapy.     

In silico decryption of serotonin–receptor binding: local non-covalent interactions and long-range conformational changes

Serotonin–receptor binding is the key step in the process behind serotonin functionality, including several psychological and physiological behaviours. This study is focused on identifying the main non-covalent interactions controlling the stability of serotonin–receptor complexes as well as the main conformational changes in the receptor due to serotonin–receptor binding using classical molecular dynamics simulations and quantum chemical calculations. A qualitative analysis based on two order parameters ((i) the centre of mass distance and (ii) the angle between the surface normals of each aromatic residue and serotonin in the binding site) on the serotonin–receptor complex trajectory suggests that the T-type stacking interaction is predominant in the binding site. Quantum chemical calculations of the stacking interaction energy provide the quantitative contributions of important aromatic residues to the stabilization of the complex. Furthermore, a three body stacking interaction (named ‘L’-type) was observed and likely contributes to the stability of the complex. Direct and water-mediated hydrogen bonding between the residues in the binding site and serotonin contributes to the complex stability. Principal component analysis of the molecular dynamics simulation trajectory of the serotonin–receptor complex and the apo-receptor in water indicates that the whole receptor is significantly stabilized due to serotonin binding. An analysis based on the dynamic cross correlation function reflects the strong correlation between trans-membrane (TM)3, TM5, TM6 (containing residues responsible for the stacking interaction and hydrogen bonding) and mini-G₀ which may participate in signal transduction leading to the functionality of serotonin.

This study is focused on identifying the main non-covalent interactions controlling the stability of serotonin–receptor complexes as well as the main conformational changes in the receptor due to serotonin–receptor binding.     

In silico approach: biological prediction of nordentatin derivatives as anticancer agent inhibitors in the cAMP pathway

A combination of computational techniques has been carried out to predict the binding of nordentatin derivatives based on pyranocoumarin semi-synthesis with the target protein from the expression of the PDE4B gene. The inhibition of the cAMP pathway is the main target of anti-cancer drugs, which is responsible for uncontrolled cell division in cancer. Modeling was done using a combination of semi-empirical methods and the density functional theory (PM3-DFT/6-31G*/B3LYP) to obtain the optimal structure of a small ligand that could be modeled. Studies on the interaction of the ligands and amino acid residues on protein targets were carried out using a combination of molecular docking and molecular dynamic simulation. Molecular docking based on functional grid scores showed a very good native ligand pose with an RMSD of 0.93 Å in determining the initial coordinates of the ligand–receptor interactions. Furthermore, the amino acid residues responsible for interaction through H-bonds were Tyr103, His104, His177, Met217, and Gln313. The binding free energy (kcal mol⁻¹) results of the candidates were PS-1 (−36.84 ± 0.31), PS-2 (−35.34 ± 0.28), PS-3 (−26.65 ± 0.30), PS-5 (−42.66 ± 0.26), PS-7 (−35.33 ± 0.23), and PS-9 (−32.57 ± 0.20), which are smaller than that of the native ligand Z72 (−24.20 ± 0.19), and thus these have good potential as drugs that can inhibit the cAMP pathway. These results provide theoretical information for the efficient inhibition of the cAMP pathway in the future.

A combination of computational techniques has been carried out to predict the binding of nordentatin derivatives based on pyranocoumarin semi-synthesis with the target protein from the expression of the PDE4B gene.     

In silico screening for oligopeptides useful as capture and reporting probes for interleukin-6 biosensing

IL-6 is an important interleukin associated with inflammation and several diseases such as cancer. Evaluation of its levels in human blood sera is a critical step for an accurate diagnosis of the diseases. Our goal is to design peptides that can selectively bind in different poses with good affinities to IL-6. For this purpose, we started from the crystal structures of different IL-6/protein complexes available in the Protein Data Bank (PDB) to select short peptides in the interaction zones, in which we intentionally introduced point mutations to increase their stability and affinity. To examine their usefulness as capture and reporting probes for the IL-6 biosensing, the five peptides and their interaction with IL-6 were studied in saline aqueous solution. Molecular docking, MD, and MM-PBSA were used to investigate the affinity and stability of these complexes. The conformational changes, the distance between the mass centers, the gyration radii, and the numbers of hydrogen bonds were analyzed to select the most suitable candidates. Three peptides, namely CTE17, CAY15 and CSE25, have the highest affinities presenting significant numbers of residues that have contact frequencies greater than 50% of simulation run time and are the most promising candidates. CTE17 and CSE25 showed they can form a stable sandwich with the target protein. For sake of comparison, we examined the previously known peptides (FND20, INL19 and CEK17) having affinity to IL-6 and the affinity of the lead i.e. CSE25 to two other interleukin family members (IL-4 and to IL-10).

In silico design by docking and molecular dynamics of short peptides that can selectively recognize IL-6 for biosensing purposes.     

An association study of ABCG2 rs2231142 on the concentrations of allopurinol and its metabolites

ABCG2 is a gene that codes for the human breast cancer resistance protein (BCRP). It is established that rs2231142 G>T, a single nucleotide polymorphism of the ABCG2 gene, is associated with gout and poor response to allopurinol, a uric acid‐lowering agent used to treat this condition. It has also been suggested that oxypurinol, the primary active metabolite of allopurinol, is a substrate of the BCRP. We thus hypothesized that carrying the rs2231142 variant would be associated with decreased oxypurinol concentrations, which would explain the lower reduction in uric acid. We performed a cross‐sectional study to investigate the association between the ABCG2 rs2231142 variant and oxypurinol, allopurinol, and allopurinol riboside concentrations in 459 participants from the Montreal Heart Institute Hospital Cohort. Age, sex, weight, use of diuretics, and estimated glomerular filtration rate were all significantly associated with oxypurinol plasma concentration. No association was found between rs2231142 and oxypurinol, allopurinol and allopurinol riboside plasma concentrations. Rs2231142 was not significantly associated with daily allopurinol dose in the overall population, but an association was observed in men, with T carriers receiving higher doses. Our results do not support a major role of ABCG2 in the pharmacokinetics of allopurinol or its metabolites. The underlying mechanism of the association between rs2231142 and allopurinol efficacy requires further investigation.

Study Highlights

• WHAT IS THE CURRENT KNOWLEDGE ON THE TOPIC?

ABCG2 is a gene that codes for the human BCRP. It is established that rs2231142 G>T, a single nucleotide polymorphism of the ABCG2 gene, is associated with gout and poor response to allopurinol, a uric acid‐lowering agent used to treat this condition. It is not clear if pharmacokinetic changes are involved in the underlying mechanism of those observations.

• WHAT QUESTION DID THIS STUDY ADDRESS?

This study addressed whether the rs2231142 loss‐of‐function variant is associated with decreased oxypurinol concentrations.

• WHAT DOES THIS STUDY ADD TO OUR KNOWLEDGE?

Our results do not support a major role of ABCG2 in the pharmacokinetics of allopurinol or its metabolites.

• HOW MIGHT THIS CHANGE CLINICAL PHARMACOLOGY OR TRANSLATIONAL SCIENCE?

Further investigations of the links between ABCG2 rs2231142 and the pharmacodynamics of allopurinol are needed.     

In vitro studies of maleidride-forming enzymes

In vitro assays of enzymes involved in the biosynthesis of maleidrides from polyketides in fungi were performed. The results show that the enzymes are closely related to primary metabolism enzymes of the citric acid cycle in terms of stereochemical preferences, but with an expanded substrate selectivity. A key citrate synthase can react both saturated and unsaturated acyl CoA substrates to give solely anti substituted citrates. This undergoes anti-dehydration to afford an unsaturated precursor which is cyclised in vitro by ketosteroid-isomerase-like enzymes to give byssochlamic acid.

In vitro synthesis of byssochkamic acid 12 was achieved from hexenoyl CoA 14dvia anhydride 1.