Free Energy Landscape Analysis Research Articles

High throughput virtual screening and validation of Plant-Based EGFR L858R kinase inhibitors against Non-Small cell lung Cancer: An integrated approach Utilizing GC–MS, network Pharmacology, Docking, and molecular dynamics

Lung cancer ranks as the 2nd most common cancer globally. It’s the most prevalent cancer in men and the 2nd most common in women. The prominent events in EGFR-mutated non-small-cell lung cancer (NSCLC) include the emergence of the L858R mutation within EGFR exon 21. Despite the promising efficacy of EGFR inhibitors in managing lung cancer, the development of acquired resistance poses a significant hurdle. In the current investigation, we focused on the screening of two phytochemicals, namely Dehydrocostus lactone and Mokkolactone, derived from the Saussurea lappa plant, as potential inhibitors targeting EGFR L858R mutant lung cancer. The chloroform and ethanol extract of the plant demonstrated anti-proliferative activity through the Resazurin chemosensitivity assay, exhibiting an IC50 value of 37.90 ± 0.29 µg/ml with selectivity index 2.4. Through a GC–MS study, we identified 11 phytochemicals for further insilico analysis. These compounds underwent ADMET assessment followed by drug likeliness analysis before being subjected to molecular docking against EGFR L858R, identified through protein–protein interaction network analysis. All phytochemicals exhibited binding energy scores ranging from −6.9 to −8.1 kcal/mol. Dehydrocostus lactone and Mokkolactone were specifically identified for their binding profile. Findings from 100 ns molecular dynamics simulations demonstrated their enhanced stability compared to the reference ligand DJK. This was evident in the root mean square deviation (RMSD) values, ranging from 0.23 ± 0.01 nm to 0.30 ± 0.05 nm, the radius of gyration values, from 1.71 ± 0.01 nm to 1.72 ± 0.01 nm, and the solvent accessible surface area values, from 155.39 ± 2.40 nm2 to 159.32 ± 2.14 nm2. Additionally, favourable characteristics were observed in terms of hydrogen bonding, principal component analysis, and free energy landscape analysis. Examination of their electronic structure via density functional theory revealed efficient properties, with the highest occupied molecular orbital-least unoccupied molecular orbital energy gap values ranging from −3.984 eV to −6.547 eV. Further, in vivo analysis is required to gain a more comprehensive understanding and efficacy of these identified phytochemicals against lung cancer.

The C-terminal activating domain promotes Panx1 channel opening.

Pannexin 1 (Panx1) constitutes a large pore channel responsible for the release of ATP from apoptotic cells. Strong evidence indicates that caspase-mediated cleavage of the C-terminus promotes the opening of the Panx1 channel by unplugging the pore. However, this simple pore- plugging mechanism alone cannot account for the observation that a Panx1 construct ending before the caspase cleavage site remains closed. Here, we show that a helical region located immediately before the caspase cleavage site, referred to as the "C-terminal activating domain (CAD)," plays a pivotal role in facilitating Panx1 activation. Electrophysiology and mutagenesis studies uncovered that two conserved leucine residues within the CAD plays a pivotal role. Cryo- EM analysis of the construct ending before reaching the CAD demonstrated that the N-terminus extends into an intracellular pocket. In contrast, the construct including the CAD revealed that this domain occupies the intracellular pocket, causing the N-terminus to flip upward within the pore. Analysis of electrostatic free energy landscape in the closed conformation indicated that the intracellular side of the ion permeation pore may be occupied by anions like ATP, creating an electrostatic barrier for anions attempting to permeate the pore. When the N-terminus flips up, it diminishes the positively charged surface, thereby reducing the drive to accumulate anions inside the pore. This dynamic change in the electrostatic landscape likely contributes to the selection of permeant ions. Collectively, these experiments put forth a novel mechanism in which C-terminal cleavage liberates the CAD, causing the repositioning of the N-terminus to promote Panx1 channel opening.

Mechanism of folding and stability of Met80Gly mutant of cytochrome-c

The present study deals with identifying conformational ensembles underlying the folding pathway of a model protein, the wild type (WT) horse cytochrome-c (cyt-c) and its mutant (M80G) in water and 8 M urea at different temperatures (300–450 K). Outcomes of this study revealed that the structural dynamics of both proteins remain restricted to the marginal changes in water and urea at 300 and 350 K. The dynamic equilibrium shifted towards unfolding at 400 K, and structures of both proteins are lost completely at 450 K, however, it occurred distinctly in water and urea. Free energy landscape analyses suggested that the unfolding of the WT and M80G proteins happens through different conformational states, the native (N), meta-stable intermediate (IN), intermediate (I), partially unfolded (IU) and completely unfolded (U) states. In urea at 400 K, IN state of M80G mutant is relatively long-lived as compared to the WT protein. However, the structure of M80G protein is lost more frequently in urea at 450 K. Our in vitro studies of the urea-induced denaturation of the WT cyt-c and its M80G mutant at pH 6.0 and 25 °C showed that the mutant is more stable than the WT protein. This conclusion is reached from the analysis of the denaturation curve (mean residue ellipticity at 222 nm versus urea concentration plot).

An integrated investigation combining network pharmacology and computer simulation dissects the potential mechanism of anti-obesity of Monascus pigments (MPs)

Monascus pigments (MPs), the secondary metabolites of Monascus-fermented foods, have been widely studied for their anti-obesity effects. However, the specific targets and potential mechanism of MPs for anti-obesity remain unclear. In this study, we applied network pharmacology and computer simulation to investigate the anti-obesity potential and action mechanism of MPs. Based on public databases, 117 potential therapeutic obesity targets of MPs were predicted. Five key targets (STAT3, ESR1, BCL2, PPARG, and MMP9) were obtained by protein-protein interaction (PPI) analysis and microarray data validation. The enrichment results showed that MPs exerted anti-obesity effects through multiple targets and multiple pathways, especially insulin resistance. Molecular docking demonstrated the binding of five key targets to MPs ligands through hydrophobic interactions and hydrogen bonding. Molecular dynamics simulation and binding free energy analysis revealed a higher stability in the binding of ESR1, BCL2, and PPARG targets with MPs ligands. In addition, free energy landscape (FEL) and independent gradient model (IGM) analysis showed that van der Waals forces were also the main driving force in stabilizing the binding of key targets and MPs ligands. This study systematically illustrated the potential anti-obesity mechanism of MPs and established the foundation for the development and application of MPs as anti-obesity functional additives.

The two-step strategy for enhancing the specific activity and thermostability of alginate lyase AlyG2 with mechanism for improved thermostability

To overcome the trade-off challenge encountered in the engineering of alginate lyase AlyG2 from Seonamhaeicola algicola Gy8T and to expand its potential industrial applications, we devised a two-step strategy encompassing activity enhancement followed by thermal stability engineering. To enhance the specific activity of efficient AlyG2, we strategically substituted residues with bulky steric hindrance proximal to the active pocket with glycine or alanine. This led to the generation of three promising positive mutants, with particular emphasis on the T91S mutant, exhibiting a 1.91-fold specific activity compared to the wild type. To mitigate the poor thermal stability of T91S, mutants with negative ΔΔG values in the thermal flexibility region were screened out. Notably, the S72Ya mutant not only displayed 17.96 % further increase in specific activity but also exhibited improved stability compared to T91S, manifesting as a remarkable 30.97 % increase in relative activity following a 1-hour incubation at 42 °C. Furthermore, enhanced kinetic stability was observed. To gain deeper insights into the mechanism underlying the enhanced thermostability of the S72Ya mutant, we conducted molecular dynamics simulations, principal component analysis (PCA), dynamic cross-correlation map (DCCM), and free energy landscape (FEL) analysis. The results unveiled a reduction in the flexibility of the surface loop, a stronger correlation dynamic and a narrower motion subspace in S72Ya system, along with the formation of more stable hydrogen bonds. Collectively, our findings suggest amino acids substitutions resulting in smaller side chains proximate to the active site can positively impact enzyme activity, while reducing the flexibility of surface loops emerges as a pivotal factor in conferring thermal stability. These insights offer valuable guidance and a framework for the engineering of other enzyme types.

Cheminformatics approach to identify andrographolide derivatives as dual inhibitors of methyltransferases (nsp14 and nsp16) of SARS-CoV-2

The Covid-19 pandemic outbreak has accelerated tremendous efforts to discover a therapeutic strategy that targets severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to control viral infection. Various viral proteins have been identified as potential drug targets, however, to date, no specific therapeutic cure is available against the SARS-CoV-2. To address this issue, the present work reports a systematic cheminformatic approach to identify the potent andrographolide derivatives that can target methyltransferases of SARS-CoV-2, i.e. nsp14 and nsp16 which are crucial for the replication of the virus and host immune evasion. A consensus of cheminformatics methodologies including virtual screening, molecular docking, ADMET profiling, molecular dynamics simulations, free-energy landscape analysis, molecular mechanics generalized born surface area (MM-GBSA), and density functional theory (DFT) was utilized. Our study reveals two new andrographolide derivatives (PubChem CID: 2734589 and 138968421) as natural bioactive molecules that can form stable complexes with both proteins via hydrophobic interactions, hydrogen bonds and electrostatic interactions. The toxicity analysis predicts class four toxicity for both compounds with LD50 value in the range of 500–700 mg/kg. MD simulation reveals the stable formation of the complex for both the compounds and their average trajectory values were found to be lower than the control inhibitor and protein alone. MMGBSA analysis corroborates the MD simulation result and showed the lowest energy for the compounds 2734589 and 138968421. The DFT and MEP analysis also predicts the better reactivity and stability of both the hit compounds. Overall, both andrographolide derivatives exhibit good potential as potent inhibitors for both nsp14 and nsp16 proteins, however, in-vitro and in vivo assessment would be required to prove their efficacy and safety in clinical settings. Moreover, the drug discovery strategy aiming at the dual target approach might serve as a useful model for inventing novel drug molecules for various other diseases.

Investigating the role of functional mutations in leucine binding to Sestrin2 in aging and age-associated degenerative pathologies using structural and molecular simulation approaches

Leucine is the native known ligand of Sestrin2 (Sesn2) and its interaction with Sesn2 is particularly important, as it influences the activity of mTOR in aging and its associated pathologies. It is important to find out how leucine interacts with Sesn2 and how mutations in the binding pocket of leucine affect the binding of leucine. Therefore, this study was committed to investigating the impact of non-synonymous mutations by incorporating a broad spectrum of simulation techniques, from molecular dynamics to free energy calculations. Our study was designed to model the atomic-scale interactions between leucine and mutant forms of Sesn2. Our results demonstrated that the interaction paradigm for the mutants has been altered thus showing a significant decline in the hydrogen bonding network. Moreover, these mutations compromised the dynamic stability by altering the conformational flexibility, sampling time, and leucine-induced structural constraints that consequently caused variation in the binding and structural stability. Molecular dynamics-based flexibility analysis revealed that the regions 217–339 and 371–380 demonstrated a higher fluctuation. Noteworthy, these regions correspond to a linker (217–339) and a loop (371–380) that cover the leucine binding cavity that is critical for the ‘latch’ mechanism in the N-terminal, which is essential for leucine binding. Further validation of reduced binding and modified internal motions caused by the mutants was obtained through binding free energy calculations, principal components analysis (PCA), and free energy landscape (FEL) analysis. By unraveling the molecular intricacies of Sesn2-leucine interactions and their mutations, we hope to pave the way for innovative strategies to combat the inevitable tide of aging and its associated diseases. Communicated by Ramaswamy H. Sarma

Long-timescale atomistic simulations uncover loss-of-function mechanisms of uncharacterized Angiogenin mutants associated with ALS

Amyotrophic Lateral Sclerosis (ALS) is a devastating neurodegenerative disease characterized by progressive degeneration of motor neurons, resulting in respiratory failure and mortality within 3–5 years. Mutations in the Angiogenin (ANG) cause loss of ribonucleolytic and nuclear translocation activities, contributing to ALS pathogenesis. This study focused on investigating two uncharacterized ANG mutations, T11S and R122H, newly identified in the Project Mine consortium. Using extensive computational analysis, including structural modeling and microsecond-timescale molecular dynamics (MD) simulations, we observed conformational changes in the catalytic residue His114 of ANG induced by T11S and R122H mutations. These alterations impaired ribonucleolytic activity, as inferred through molecular docking and binding free energy calculations. Gibbs free energy landscape and residue-residue interaction network analysis further supported our findings, revealing the energetic states and allosteric pathway from the mutated site to His114. Additionally, we assessed the binding of NCI-65828, an inhibitor of ribonucleolytic activity of ANG, and found reduced effectiveness in binding to T11S and R122H mutants when His114 assumed a non-native conformation. This highlights the crucial role of His114 and its association with ALS. Elucidating the relationship between physical structure and functional dynamics of frequently mutated ANG mutants is essential for understanding ALS pathogenesis and developing more effective therapeutic interventions.

Testing the feasibility of targeting a conserved region on the S2 domain of the SARS-CoV-2 spike protein

The efficacy of vaccines against the SARS-CoV-2 virus significantly declines with the emergence of mutant strains, prompting investigation into the feasibility of targeting highly conserved but often cryptic regions on the S2 domain of spike protein. Using tools from molecular dynamics, we find that exposure of a conserved S2 epitope located in the central helices below the receptor binding domains would require large-scale motion beyond receptor binding domain up-down motion, but, along the reaction coordinates we explored, it is unlikely to be exposed by such large-scale dynamic fluctuations of the S1 domain without any external facilitating factors, despite some previous computational evidence suggesting transient exposure of this region. Furthermore, glycans, particularly those on N165 and N234, hinder S2-exposing opening dynamics, and thus stabilize spike in addition to immunologically shielding the protein surface. Although the S2 epitope region examined here is central to large-scale conformational changes during viral entry, free energy landscape analysis obtained using the path coordinate formalism reveals no inherent “loaded spring” effect, suggesting that a vaccine immunogen would tend to present the epitope in a prefusion-like conformation and may be effective in neutralization. These findings contribute to a deeper understanding of the dynamic origins of the function of the spike protein, as well as further characterizing the feasibility of the S2 epitope as a therapeutic target.

In silico validation of hyaluronic acid – drug conjugates based targeted drug delivery for the treatment of COVID-19

The impact of COVID-19 urges scientists to develop targeted drug delivery to manage Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) viral infections with a fast recovery rate. The aim of the study is to develop Hyaluronic Acid (HA) drug conjugates of viral drugs to target two important enzymes (Mpro and PLpro) of SARS-CoV-2. Three antiviral drugs, namely Dexamethasone (DEX), Favipiravir (FAV), and Remdesivir (REM), were chosen for HA conjugation due to their reactive functional groups. Free forms of drugs (DEX, FAV, REM) and HA drug conjugates (HA-DEX, HA-FAV, HA-REM, HA-RHA, HA-RHE) were validated against Mpro (PDB ID 6LU7) and PLpro (PDB 7LLZ), which play an essential role in the replication and reproduction of the SARS-CoV-2 virus. The results of the present study revealed that HA-drug conjugates possess higher binding affinity and the best docking score towards the Mpro and PLpro target proteins of SARS-CoV-2 than their free forms of drugs. ADMET screening resulted that HA-drug conjugates exhibited better pharmacokinetic profiles than their pure forms of drugs. Further, molecular dynamic simulation studies, essential dynamics and free energy landscape analyses show that HA antiviral drug conjugates possess good trajectories and energy status, with the PLpro target protein (PDB 7LLZ) of SARS-CoV-2 through long-distance (500 ns) simulation screening. The research work recorded the best drug candidate for Cell-Targeted Drug Delivery (CTDD) for SARS-CoV-2-infected cells through hyaluronic acid conjugates of antiviral drugs. Communicated by Ramaswamy H. Sarma

Structure-guided engineering and molecular simulations to design a potent monoclonal antibody to target aP2 antigen for adaptive immune response instigation against type 2 diabetes.

Diabetes mellitus (DM) is recognized as one of the oldest chronic diseases and has become a significant public health issue, necessitating innovative therapeutic strategies to enhance patient outcomes. Traditional treatments have provided limited success, highlighting the need for novel approaches in managing this complex disease. In our study, we employed graph signature-based methodologies in conjunction with molecular simulation and free energy calculations. The objective was to engineer the CA33 monoclonal antibody for effective targeting of the aP2 antigen, aiming to elicit a potent immune response. This approach involved screening a mutational landscape comprising 57 mutants to identify modifications that yield significant enhancements in binding efficacy and stability. Analysis of the mutational landscape revealed that only five substitutions resulted in noteworthy improvements. Among these, mutations T94M, A96E, A96Q, and T94W were identified through molecular docking experiments to exhibit higher docking scores compared to the wild-type. Further validation was provided by calculating the dissociation constant (KD), which showed a similar trend in favor of these mutations. Molecular simulation analyses highlighted T94M as the most stable complex, with reduced internal fluctuations upon binding. Principal components analysis (PCA) indicated that both the wild-type and T94M mutant displayed similar patterns of constrained and restricted motion across principal components. The free energy landscape analysis underscored a single metastable state for all complexes, indicating limited structural variability and potential for high therapeutic efficacy against aP2. Total binding free energy (TBE) calculations further supported the superior performance of the T94M mutation, with TBE values demonstrating the enhanced binding affinity of selected mutants over the wild-type. Our findings suggest that the T94M substitution, along with other identified mutations, significantly enhances the therapeutic potential of the CA33 antibody against DM by improving its binding affinity and stability. These results not only contribute to a deeper understanding of antibody-antigen interactions in the context of DM but also provide a valuable framework for the rational design of antibodies aimed at targeting this disease more effectively.

Zingiberaceae-derived phytomolecules inhibit Japanese encephalitis virus RNA dependent RNA polymerase: a molecular dynamics study

The Japanese encephalitis virus, (JEV), is a flavivirus mostly transmitted by Culex mosquitoes mostly present in Southeast Asia and the Western Pacific region. Ardeid-wading birds are the natural reservoir of JEV; nonetheless, pigs are frequently a key amplifying host during epidemics in human populations. Although more domestic animals and wildlife are JEV hosts, it is unclear how these animals fit into the ecology and epidemiology of the virus. Even though there is no specific therapy, vaccines are available to prevent this infection. However, current vaccinations do not work against every clinical isolate and can cause neurological problems in certain people. In this study, we have screened 501 phytochemical compounds from various plants from the Zingeberaceae family against the RdRp protein of JEV. Based on this, the top five compounds (IMPHY014466, IMPHY004928, IMPHY007097, IMPHY014179 and IMPHY005010) were selected based on the obtained docking scores, which was above −8.0 Kcal/mol. Further, the binding affinity of these selected ligands was also analysed using molecular interaction, and the presence of interactions like hydrogen bonds, hydrophobic bonds and polar bonds with respective active residues were identified and studied elaborately. Furthermore, the dynamic stability of the docked RdRp protein with these selected phytochemicals was studied using Molecular dynamic simulation and essential dynamics. The free energy landscape analysis also provided information about the energy transition responsible stability of the complex. The results obtained advocated phytochemical compounds from the zingeberaceae family for future experimental validation, as these compounds exhibited significant potential as JEV antagonists. Communicated by Ramaswamy H. Sarma

Multi-targeted therapeutic potential of stigmasterol from the Euphorbia ammak plant in treating lung and breast cancer

Cancer is the most prevalent disease globally, which presents a significant challenge to the healthcare industry, with breast and lung cancer being predominant malignancies. This study used RNA-seq data from the TCGA database to identify potential biomarkers for lung and breast cancer. Tumor Necrosis Factor (TNFAIP8) and Sulfite Oxidase (SUOX) showed significant expression variation and were selected for further study using structure-based drug discovery (SBDD). Compounds derived from the Euphorbia ammak plant were selected for in-silico study with both TNFAIP8 and SUOX. Stigmasterol had the greatest binding scores (normalized scores of −8.53 kcal/mol and −9.69 kcal/mol) with both proteins, indicating strong stability in their binding pockets throughout the molecular dynamics’ simulation. Although Stigmasterol first changed its initial conformation (RMSD = 0.5 nm with the starting conformation) in SUOX, it eventually reached a stable conformation (RMSD of 1.5 nm). The compound on TNFAIP8 showed a persistent shape (RMSD of 0.35 nm), indicating strong protein stability. The binding free energy of the complex was calculated using the MM/GBSA technique; TNFAIP8 had a ΔGTOTAL of −24.98 kcal/mol, with TYR160 being the most significant residue, contributing −2.52 kcal/mol. On the other hand, the SUOX complex had a binding free energy of −16.87 kcal/mol, with LEU151 being the primary contributor (-1.17 kcal/mol). Analysis of the complexes' free energy landscape unveiled several states with minimum free energy, indicating robust interactions between the protein and ligand. In its conclusion, this work emphasises the favourable ability of Stigmasterol to bind with prospective targets for lung and breast cancer, indicating the need for more experimental study.

Understanding the basis of thermostability for enzyme “Nanoluc” towards designing industry-competent engineered variants

As a leading contender in the study of luminescence, nanoluciferase has recently attracted attention and proven effective in a wide variety of research areas. Although numerous attempts have been made to improve activity, there has yet to be a thorough exploration of further possibilities to improve thermostability. In this study, protein engineering in tandem with molecular dynamics simulation at various temperatures (300 K, 400 K, 450 K and 500 K) was used to improve our understanding of nanoluciferase dynamics and identification of factors that could significantly enhance the thermostability. Based on these, three novel mutations have been narrowed down, which were hypothesised to improve thermostability. Root mean square deviation and root mean square fluctuation studies confirmed higher stability of mutant at high temperature. Solvent-accessible surface area and protein unfolding studies revealed a decreased tendency of mutant to unfold at higher temperatures. Further free energy landscape and principal component analysis was adapted to get deeper insights into the thermodynamic and structural behavior of these proteins at elevated temperature. Thus, this study provides a deeper insight into the dynamic factors for thermostability and introduces a novel, enhanced nanoluciferase candidate with potential use in industry. Communicated by Ramaswamy H. Sarma

Targeting PDE4A for therapeutic potential: exploiting drug repurposing approach through virtual screening and molecular dynamics

cAMP-specific 3′,5′-cyclic phosphodiesterase 4 A (PDE4A) holds a pivotal role in modulating intracellular levels of cyclic adenosine monophosphate (cAMP). Targeting PDE4A with novel therapeutic agents shows promise in addressing neurological disorders (e.g. Alzheimer’s and Parkinson’s diseases), mood disorders (depression, anxiety), inflammatory conditions (asthma, chronic obstructive pulmonary disease), and even cancer. In this study, we present a comprehensive approach that integrates virtual screening and molecular dynamics (MD) simulations to identify potential inhibitors of PDE4A from the existing pool of FDA-approved drugs. The initial compound selection was conducted focusing on binding affinity scores, which led to the identification of several high-affinity compounds with potential PDE4A binding properties. From the refined selection process, two promising compounds, Fluspirilene and Dihydroergocristine, emerged as strong candidates, displaying substantial affinity and specificity for the PDE4A binding site. Interaction analysis provided robust evidence of their binding capabilities. To gain deeper insights into the dynamic behavior of Fluspirilene and Dihydroergocristine in complex with PDE4A, we conducted 300 ns MD simulations, principal components analysis (PCA), and free energy landscape (FEL) analysis. These analyses revealed that Fluspirilene and Dihydroergocristine binding stabilized the PDE4A structure and induced minimal conformational changes, highlighting their potential as potent binders. In conclusion, our study systematically explores repurposing existing FDA-approved drugs as PDE4A inhibitors through a comprehensive virtual screening pipeline. The identified compounds, Fluspirilene and Dihydroergocristine, exhibit a strong affinity for PDE4A, displaying characteristics that support their suitability for further development as potential therapeutic agents for conditions associated with PDE4A dysfunction. Communicated by Ramaswamy H. Sarma

Variational analysis and AI algorithm implementation of free energy landscapes of molecular system

Accurate description of the free energy landscape (FES) is the basis for understanding complex molecular systems, and for further realizing molecular design, manufacture and industrialization. Major challenges include multiple metastable states, which usually are separated by high potential barriers and are not linearly separable, and may exist at multiple levels of time and spatial scales. Consequently FES is not suitable for analytical analysis and brute force simulation. To address these challenges, many enhanced sampling methods have been developed. However, utility of them usually involves many empirical choices, which hinders research advancement, and also makes error control very unimportant. Although variational calculus has been widely applied and achieved great success in physics, engineering and statistics, its application in complex molecular systems has just begun with the development of neural networks. This brief review is to summarize the background, major developments, current limitations, and prospects of applying variation in this field. It is hoped to facilitate the AI algorithm development for complex molecular systems in general, and to promote the further methodological development in this line of research in particular.

Elucidating the pH-dependent conformational transitions of eryngin: Insights from molecular dynamics simulations

Eryngin is a novel antifungal peptide derived from the Pleurotus eryngii mushroom with potential drugproperties. This study aims to explore the stability and structural integrity of eryngin at different pH levels found within the human body. The molecular dynamic simulation method showed that eryngin maintains its overall structure within pH ranging from 1.5 to 8. The free energy landscape analysis showed that the free energy was distributed in the range of 12.80–19.30 kJ/mol. This investigation suggests that eryngin could potentially function effectively as a treatment for various diseases at different body sites with varying pH environments, owing to its ability to retain its secondary structure. Additional research is necessary to fully comprehend the biological activity of eryngin and its potential as a therapeutic agent.

Exploring antibiofilm potential of some new imidazole analogs against C. albicans: synthesis, antifungal activity, molecular docking and molecular dynamics studies

A series of 1, 2, 4, 5-tetrasubstituted imidazole derivatives were synthesized and their antibiofilm potential against Candida albicans was evaluated in vitro. Two of the synthesized derivatives 5e (IC50 = 25 µg/mL) and 5m (IC50 = 6 µg/mL),displayed better antifungal and antibiofilm potential than the standard drug Fluconazole (IC50 = 40 µg/mL) against C. albicans. Based on the in vitro results, we escalated the real time polymerase chain reaction (RT-PCR) analysis to gain knowledge of the enzymes expressed in the generation and maintenance of biofilms and the mechanism of biofilm inhibition by the synthesized analogues. We then investigated the possible interactions of the synthesized compounds in inhibiting agglutinin-like proteins, namely Als3, Als4 and Als6 were prominently down-regulated using in-silico molecular docking analysis against the previously available crystal structure of Als3 and constructed structure of Als4 and Als6 using the SWISS-MODEL server. The stability and energy of the agglutinin-like proteins–ligand complexes were evaluated using molecular dynamics simulations (MDS). According to the 100 ns MDS, all the compounds remained stable, formed a maximum of 3, and on average 2 hydrogen bonds, and Gibb’s free energy landscape analysis suggested greater affinity of the compounds 5e and 5m toward Als4 protein. Communicated by Ramaswamy H. Sarma

Identifying Phosphodiesterase-5 Inhibitors with Drug Repurposing Approach: Implications in Vasodysfunctional Disorders.

Phosphodiesterase type 5 (PDE5) is a multidomain protein that plays a crucial role in regulating cellular cyclic guanosine monophosphate (cGMP), a key signaling molecule involved in various physiological processes. Dysregulation of PDE5 and cGMP signaling is associated with a range of vasodysfunctional disorders, necessitating the development of effective therapeutic interventions. This study adopts comprehensive approach, combining virtual screening and molecular dynamics (MD) simulations, to repurpose FDA-approved drugs as potential PDE5 inhibitors. The initial focus involves selecting compounds based on their binding affinity. Shortlisted compounds undergo a meticulous analysis for their drug profiling and biological significance, followed by the activity evaluation and interaction analysis. Notably, based on binding potential and drug profiling, two molecules, Dutasteride and Spironolactone, demonstrate strong potential as PDE5 inhibitors. Furthermore, all atom MD simulations were employed (500 ns) to explore dynamic behavior of Dutasteride and Spironolactone in complexes with PDE5. Principal components analysis (PCA) and free energy landscape (FEL) analyses are further leveraged to decipher that the binding of Dutasteride and Spironolactone stabilizes the structure of PDE5 with minimal conformational changes. In summary, Dutasteride and Spironolactone exhibit remarkable affinity for PDE5 and possess characteristics that suggest their potential as therapeutic agents for conditions associated with PDE5 dysfunction.

Computational insights into potential marine natural products as selective inhibitors of Mycobacterium tuberculosis InhA: A structure-based virtual screening study

Several factors are associated with the emergence of drug resistance mechanisms, such as impermeable cell walls, gene mutations, and drug efflux systems. Consequently, bacteria acquire resistance, leading to a decrease in drug efficacy. A new and innovative strategy is required to combat drug resistance in tuberculosis (TB) effectively. Therefore, targeting the mycolic acid biosynthesis pathway, which is involved in synthesising mycolic acids (MAs), essential structural components responsible for mycobacterial pathogenicity, has garnered interest in TB research and the concept of drug resistance. In this context, InhA, which plays a crucial role in the fatty acid synthase-II (FAS-II) system of the MA biosynthetic pathway, was selected as a druggable target for screening investigation. To identify potential lead molecules against InhA, diverse marine natural products (MNPs) were collected from the comprehensive marine natural products database (CMNPD). Virtual screening studies aided in selecting potential lead molecules that best fit within the substrate-binding pocket (SBP) of InhA, forming crucial hydrogen bond interaction with the catalytic residue Tyr158. Three MNPs, CMNPD30814, CMNPD1702, and CMNPD27355, were chosen as prospective alternative molecules due to their favorable pharmacokinetic properties and lack of toxicity according to ProTox-II predictions. Additionally, improved reactivity of the MNPs was observed in the results of density functional theory (DFT) studies. Furthermore, comparative molecular dynamics simulation (MDS), principal component (PC)-based free energy landscape (FEL) analysis, and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) were employed to show enhanced structural stability, increased H-bond potential, and high binding affinity toward the target InhA. Moreover, the hot spot residues that contributed to the high binding energy profile and anchored the stability of the complexes were revealed with their individual interaction energy. The computational insights from this study provide potential avenues to combat TB through the multifaceted mode of action of these marine lead molecules, which can be further explored in future experimental investigations.