Physicochemical stability study of protein–benzoic acid complexes using molecular dynamics simulations

Conformations of a highly expressed Z19 α-zein studied with AlphaFold2 and MD simulations

α-zeins are amphiphilic maize seed storage proteins with material properties suitable for a multitude of applications e.g., in renewable plastics, foods, therapeutics and additive manufacturing (3D-printing). To exploit their full potential, molecular-level insights are essential. The difficulties in experimental atomic-resolution characterization of α-zeins have resulted in a diversity of published molecular models. However, deep-learning α-zein models are largely unexplored. Therefore, this work studies an AlphaFold2 (AF2) model of a highly expressed α-zein using molecular dynamics (MD) simulations. The sequence of the α-zein cZ19C2 gave a loosely packed AF2 model with 7 α-helical segments connected by turns/loops. Compact tertiary structure was limited to a C-terminal bundle of three α-helices, each showing notable agreement with a published consensus sequence. Aiming to chart possible α-zein conformations in practically relevant solvents, rather than the native solid-state, the AF2 model was subjected to MD simulations in water/ethanol mixtures with varying ethanol concentrations. Despite giving structurally diverse endpoints, the simulations showed several patterns: In water and low ethanol concentrations, the model rapidly formed compact globular structures, largely preserving the C-terminal bundle. At ≥ 50 mol% ethanol, extended conformations prevailed, consistent with previous SAXS studies. Tertiary structure was partially stabilized in water and low ethanol concentrations, but was disrupted in ≥ 50 mol% ethanol. Aggregated results indicated minor increases in helicity with ethanol concentration. β-sheet content was consistently low (∼1%) across all conditions. Beyond structural dynamics, the rapid formation of branched α-zein aggregates in aqueous environments was highlighted. Furthermore, aqueous simulations revealed favorable interactions between the protein and the crosslinking agent glycidyl methacrylate (GMA). The proximity of GMA epoxide carbons and side chain hydroxyl oxygens simultaneously suggested accessible reactive sites in compact α-zein conformations and pre-reaction geometries for methacrylation. The findings may assist in expanding the applications of these technologically significant proteins, e.g., by guiding chemical modifications.

Campesterol and dithymoquinone as a potent inhibitors of SARS cov-2 main proteases-promising drug candidates for targeting its novel variants

The sudden outbreak of the COVID-19 pandemic has currently taken approximately 2.4 million lives, with no specific medication and fast-tracked tested vaccines for prevention. These vaccines have their own adverse effects, which have severely affected the global healthcare system. The discovery of the main protease structure of coronavirus (Mpro/Clpro) has resulted in the identification of compounds having antiviral potential, especially from the herbal system. In this study, the computer-associated drug design tools were utilised to analyze the reported phytoconstituents of Nigella sativa for their antiviral activity against the main protease. Fifty-eight compounds were subjected to pharmacological parameter analysis to determine their lead likeness in comparison to the standard drugs (chloroquine and nirmatrelvir) used in the treatment of SARS-CoV-2. Nearly 31 compounds were docked against five different SARS-CoV-2 main proteases, and all compounds showed better binding affinity and inhibition constant against the proteases. However, dithymoquinone and campesterol displayed the best binding scores and hence were further subjected to dynamics and MMPBSA study for 100 ns. The stability analysis shows that dithymoquinone and campesterol show less variation in fluctuation in residues compared to standard complexes. Moreover, dithymoquinone exhibited higher binding affinity and favorable interaction followed by campesterol as compared to the standard drug. The in silico computational analysis provides a promising hit for regulating the main proteases activity.Communicated by Ramaswamy H. Sarma.

Doping matters in carbon nanomaterial efficiency in environmental remediation

Carbon nanomaterials (CNMs) are rapidly emerging in materials science research due to their widespread environmental applications. They are useful for environmental pollutants' remediation through various methods. Heteroatom doping resulted in reliable approaches to overcome pristine CNMs challenges. The engineering of the dopants is believed to be a promising route to improve the efficiency of CNMs in environmental remediation. The idea of doping has been attractive since it allows the control of electronic properties due to the electron transfer between dopants and the host material and the dopants along with the bonding between analogous atoms and carbon atoms. This mini-review, through computational and experimental studies, puts special emphasis on the role of doping different CNMs as an efficient approach to enhance the environmental remediation.

Virtual Screening of Benzimidazole Derivatives as Potential Triose Phosphate Isomerase Inhibitors with Biological Activity against Leishmania mexicana

Leishmania mexicana (L. mexicana) is a causal agent of cutaneous leishmaniasis (CL), a “Neglected disease”, for which the search for new drugs is a priority. Benzimidazole is a scaffold used to develop antiparasitic drugs; therefore, it is interesting molecule against L. mexicana. In this work, a ligand-based virtual screening (LBVS) of the ZINC15 database was performed. Subsequently, molecular docking was used to predict the compounds with potential binding at the dimer interface of triosephosphate isomerase (TIM) of L. mexicana (LmTIM). Compounds were selected on binding patterns, cost, and commercial availability for in vitro assays against L. mexicana blood promastigotes. The compounds were analyzed by molecular dynamics simulation on LmTIM and its homologous human TIM. Finally, the physicochemical and pharmacokinetic properties were determined in silico. A total of 175 molecules with docking scores between −10.8 and −9.0 Kcal/mol were obtained. Compound E2 showed the best leishmanicidal activity (IC50 = 4.04 µM) with a value similar to the reference drug pentamidine (IC50 = 2.23 µM). Molecular dynamics analysis predicted low affinity for human TIM. Furthermore, the pharmacokinetic and toxicological properties of the compounds were suitable for developing new leishmanicidal agents.

Identification of lactoferrin-derived peptides as potential inhibitors against the main protease of SARS-CoV-2

COVID-19 is a global health emergency that causes serious concerns. A global effort is underway to identify drugs for the treatment of COVID-19. One possible solution to the present problem is to develop drugs that can inhibit SARS-CoV-2 main protease (M^pro), a coronavirus protein that been considered as one among many drug targets. In this work, lactoferrin from Bos taurus L. was in silico hydrolyzed. The bioactivity, water solubility, and ADMET properties of the generated peptides were predicted using various online tools. The molecular interactions between M^pro and the peptides were studied using molecular docking and molecular dynamic simulation. The results demonstrated that peptide GSRY was predicted to have better physicochemical properties, and the value of '-C DOCKER interaction energy' between peptide GSRY and M^pro was 80.8505 kcal/mol. The interaction between the peptide GSRY and the native ligand N3 co-crystallized with M^pro had overlapped amino acids, i.e., HIS163, GlY143, GLU166, GLN189 and MET165. Molecular dynamic simulation revealed that M^pro/GSRY complexes were stable. Collectively, the peptide GSRY may be a potential candidate drug against M^pro of SARS-CoV-2.

Solvent-Free Synthesis, In Vitro and In Silico Studies of Novel Potential 1,3,4-Thiadiazole-Based Molecules against Microbial Pathogens

A new series of 1,3,4-thiadiazoles was synthesized by the reaction of methyl 2-(4-hydroxy-3-methoxybenzylidene) hydrazine-1-carbodithioate (2) with selected derivatives of hydrazonoyl halide by grinding method at room temperature. The chemical structures of the newly synthesized derivatives were resolved from correct spectral and microanalytical data. Moreover, all synthesized compounds were screened for their antimicrobial activities using Escherichia coli, Pseudomonas aeruginosa, Proteus vulgaris, Bacillus subtilis, Staphylococcus aureus, and Candida albicans. However, compounds 3 and 5 showed significant antimicrobial activity against all tested microorganisms. The other prepared compounds exhibited either only antimicrobial activity against Gram-positive bacteria like compounds 4 and 6, or only antifungal activity like compound 7. A molecular docking study of the compounds was performed against two important microbial enzymes: tyrosyl-tRNA synthetase (TyrRS) and N-myristoyl transferase (Nmt). The tested compounds showed variety in binding poses and interactions. However, compound 3 showed the best interactions in terms of number of hydrogen bonds, and the lowest affinity binding energy (−8.4 and −9.1 kcal/mol, respectively). From the in vitro and in silico studies, compound 3 is a good candidate for the next steps of the drug development process as an antimicrobial drug.

Computational insights into binding mechanism of drugs as potential inhibitors against SARS-CoV-2 targets

Because of the scale of the novel coronavirus (COVID-19) pandemic and the swift transmission of this highly contagious respiratory virus, repurposing existing drugs has become an urgent treatment approach. The objective of our study is to unravel the binding mechanism of the Food and Drug Administration (FDA)-approved dexamethasone (Dex) and boceprevir (Boc) drugs with selected COVID-19 protein targets SARS-CoV-2 spike protein C-terminal domain (spike-CTD), main protease (M^pro), and interleukin-6 (IL-6). Another objective is to analyze the effects of binding Dex and Boc drugs on the interactions of viral spike protein to human angiotensin-converting enzyme 2 (hACE2). Molecular docking and one-microsecond-long molecular dynamics simulations of each of the six protein–drug complexes along with steered molecular dynamics (SMD) and umbrella sampling (US) methods have revealed the binding mode interactions and the physicochemical stability of the three targeted proteins with two drugs. Results have shown that both drugs bind strongly with the three protein targets through hydrogen bonding and hydrophobic interactions. A major finding from this study is how the binding of the drugs with viral spike protein affects its interactions at the binding interface with hACE2 protein. Simulations of drug-bound spike-CTD with hACE2 show that due to the presence of a drug at the binding interface of spike-CTD, hACE2 is being blocked from making putative interactions with viral protein at such interface. These important findings regarding the binding affinity and stability of the two FDA-approved drugs with the main targets of COVID-19 along with the effect of drugs on hACE2 interactions would contribute to COVID-19 drug discovery and development.

Lysozyme and Human Serum Albumin Proteins as Potential Nitric Oxide Cardiovascular Drug Carriers: Theoretical and Experimental Investigation

Nitric oxide-containing drugs present a critical remedy for cardiovascular diseases. Nitroglycerin (NG, O-NO) and S-nitrosoglutathione (SNG, S-NO) are the most common nitric oxide drugs for cardiovascular diseases. Insights regarding the binding affinity of NO drugs with lysozyme and human serum albumin (HSA) proteins and their dissociation mechanism will provide inquisitive information regarding the potential of the proteins as drug carriers. For the first time, the binding interactions and affinities are investigated using molecular docking, conventional molecular dynamics, steered molecular dynamics, and umbrella sampling to explore the ability of both proteins to act as nitric oxide drug carriers. The molecular dynamics simulation results showed higher stability of lysozyme-drug complexes compared to HSA. For lysozyme, cardiovascular drugs were bound in the protein cavity mainly by the electrostatic and hydrogen bond interactions with residues ASP53, GLN58, ILE59, ARG62, TRP64, ASP102, and TRP109. For HSA, key binding residues were ARG410, TYR411, LYS414, ARG485, GLU450, ARG486, and SER489. The free energy profiles produced from umbrella sampling also suggest that lysozyme-drug complexes had better binding affinity than HSA-drug. Binding characteristics of nitric oxide-containing drugs NG and SNG to lysozyme and HSA proteins were studied using fluorescence and UV-vis absorption spectroscopy. The relative change in the fluorescence intensity as a function of drug concentrations was analyzed using Stern-Volmer calculations. This was also confirmed by the change in the UV-vis spectra. Fluorescence quenching results of both proteins with the drugs, based on the binding constant values, demonstrated significantly weak binding affinity to NG and strong binding affinity to SNG. Both computational and experimental studies provided important data for understanding protein-drug interactions and will aid in developing potential drug carrier systems in cardiovascular diseases.