Combination of the CHARMM27 force field with united-atom lipid force fields

Computational insights into allosteric inhibition of focal adhesion kinase: A combined pharmacophore modeling and molecular dynamics approach

Focal adhesion kinase (FAK) is a non-receptor tyrosine kinase that modulates integrin and growth factor signaling pathways and is implicated in cancer cell migration, proliferation, and survival. Over the past decade various, FAK kinase, FERM, and FAT domain inhibitors have been reported and a few kinase domain inhibitors are under clinical consideration. However, few of them were identified as multikinase inhibitors. In kinase drug design selectivity is always a point of concern, to improve selectivity allosteric inhibitor development is the best choice. The current research utilized a pharmacophore modeling (PM) approach to identify novel allosteric inhibitors of FAK. The all-available allosteric inhibitor bound 3D structures with PDB ids 4EBV, 4EBW, and 4I4F were utilized for the pharmacophore modeling. The validated PM models were utilized to map a database of 770,550 compounds prepared from ZINC, EXIMED, SPECS, ASINEX, and InterBioScreen, aiming to identify potential allosteric inhibitors. The obtained compounds from screening step were forwarded to molecular docking (MD) for the prediction of binding orientation inside the allosteric site and the results were evaluated with the known FAK allosteric inhibitor (REF). Finally, 14 FAK-inhibitor complexes were selected from the docking study and were studied under molecular dynamics simulations (MDS) for 500 ns. The complexes were ranked according to binding free energy (BFE) and those demonstrated higher affinity for allosteric site of FAK than REF inhibitors were selected. The selected complexes were further analyzed for intermolecular interactions and finally, three potential allosteric inhibitor candidates for the inhibition of FAK protein were identified. We believe that identified scaffolds may help in drug development against FAK as an anticancer agent.

Synthesis, Characterization, Antimicrobial and In Silico Studies of Isatin Schiff Base Linked 1,2,3-Triazole Hybrids

A new series of isatin-Schiff base linked 1,2,3-triazole hybrids has been synthesized using CuAAC approach from (E)-3-(phenylimino)-1-(prop-2-yn-1-yl)indolin-2-one derivatives in high yield (73-91 %). These synthesized derivatives were characterized using FT-IR, 1H NMR, 13C NMR, 2D-NMR and HRMS spectral techniques. The in vitro antimicrobial activity assay demonstrated that most of the tested hybrids exhibited promising activity. Compound 5 j displayed significant antibacterial efficacy against P. aeruginosa and B. subtilis with MIC value of 0.0062 μmol/mL. While, 5 j also showed better antifungal potency against A. niger with MIC value of 0.0123 μmol/mL. The docking studies of most promising compounds were performed with the well-known antibacterial and antifungal targets i. e. 1KZ1, 5TZ1. Molecular modelling investigations demonstrated that hybrids 5 h and 5 l exhibited good interactions with 1KZN and 5TZ1, with binding energies of -9.6 and -11.0 kcal/mol, respectively. Further, molecular dynamics studies of the compounds showing promising binding interactions were also carried out to study the stability of complexes of these hybrids with both the targets.

In silico study to identify novel NEK7 inhibitors from natural sources by a combination strategy

Cancer poses a significant global health challenge and significantly contributes to mortality. NEK7, related to the NIMA protein kinase family, plays a crucial role in spindle assembly and cell division. The dysregulation of NEK7 is closely linked to the onset and progression of various cancers, especially colon and breast cancer, making it a promising target for cancer therapy. Nevertheless, the shortage of high-quality NEK7 inhibitors highlights the need for new therapeutic strategies. In this study, we utilized a multidisciplinary approach, including virtual screening, molecular docking, pharmacokinetics, molecular dynamics simulations (MDs), and MM/PBSA calculations, to evaluate natural compounds as NEK7 inhibitors comprehensively. Through various docking strategies, we identified three natural compounds: (-)-balanol, digallic acid, and scutellarin. Molecular docking revealed significant interactions at residues such as GLU112 and ALA114, with docking scores of -15.054, -13.059, and -11.547 kcal/mol, respectively, highlighting their potential as NEK7 inhibitors. MDs confirmed the stability of these compounds at the NEK7-binding site. Hydrogen bond analysis during simulations revealed consistent interactions, supporting their strong binding capacity. MM/PBSA analysis identified other crucial amino acids contributing to binding affinity, including ILE20, VAL28, ILE75, LEU93, ALA94, LYS143, PHE148, LEU160, and THR161, crucial for stabilizing the complex. This research demonstrated that these compounds exceeded dabrafenib in binding energy, according to MM/PBSA calculations, underscoring their effectiveness as NEK7 inhibitors. ADME/T predictions showed lower oral toxicity for these compounds, suggesting their potential for further development. This study highlights the promise of these natural compounds as bases for creating more potent derivatives with significant biological activities, paving the way for future experimental validation.

Design of human ACE2 mimic miniprotein binders that interact with RBD of SARS-CoV-2 variants of concerns

The world of medicine demands from the research community solutions to the emerging problem of SARS-CoV-2 variants and other such potential global pandemics. With advantages of specificity over small molecule drugs and designability over antibodies, miniprotein therapeutics offers a unique solution to the threats of rapidly emerging SARS-CoV-2 variants. Unfortunately, most of the promising miniprotein binders are de novo designed and it is not viable to generate molecules for each new variant. Therefore in this study, we demonstrate a method for design of miniprotein mimics from the interaction interphase of human angiotensin converting enzyme 2 (ACE2). ACE2 is the natural interacting partner for the SARS-CoV-2 spike receptor binding domain (RBD) and acts as a recognition molecule for viral entry into the host cells. Starting with ACE2 N-terminal triple helix interaction interphase, we generated more than 70 miniprotein sequences. Employing Rosetta folding and docking scores we selected 10 promising miniprotein candidates amongst which 3 were found to be soluble in lab studies. Further, using molecular mechanics (MM) calculations on molecular dynamics (MD) trajectories we test interaction of miniproteins with RBD from various variants of concern (VOC). Presently, we report two key findings; miniproteins in this study are generated using less than 10 lab testing experiments, yet when tested through in-vitro experiments, they show submicro to nanomolar affinities towards SARS-CoV-2 RBD. Also in simulation studies, when compared with previously developed therapeutics, our miniproteins display remarkable ability to mimic ACE2 interphase; making them an ideal solution to the ever evolving problem of VOCs.Communicated by Ramaswamy H. Sarma.

Comprehensive mapping of mutations in TDP-43 and α-Synuclein that affect stability and binding

Abnormal aggregation and amyloid inclusions of TAR DNA-binding protein 43 (TDP-43) and α-Synuclein (α-Syn) are frequently co-observed in amyotrophic lateral sclerosis, Parkinson's disease, and Alzheimer's disease. Several reports showed TDP-43 C-terminal domain (CTD) and α-Syn interact with each other and the aggregates of these two proteins colocalized together in different cellular and animal models. Molecular dynamics simulation was conducted to elucidate the stability of the TDP-43 and Syn complex structure. The interfacial mutations in protein complexes changes the stability and binding affinity of the protein that may cause diseases. Here, we have utilized the computational saturation mutagenesis approach including structure-based stability and binding energy calculations to compute the systemic effects of missense mutations of TDP-43 CTD and α-Syn on protein stability and binding affinity. Most of the interfacial mutations of CTD and α-Syn were found to destabilize the protein and reduced the protein binding affinity. The results thus shed light on the functional consequences of missense mutations observed in TDP-43 associated proteinopathies and may provide the mechanisms of co-morbidities involving these two proteins.Communicated by Ramaswamy H. Sarma.

Label-Free SERS and MD Analysis of Biomarkers for Rapid Point-of-Care Sensors Detecting Head and Neck Cancer and Infections

For the progress of point-of-care medicine, where individual health status can be easily and quickly monitored using a handheld sensor, saliva serves as one of the best-suited body fluids thanks to its availability and abundance of physiological indicators. Salivary biomarkers, combined with rapid and highly sensitive detection tools, may pave the way to new real-time health monitoring and personalized preventative therapy branches using saliva as a target matrix. Saliva is increasing in importance in liquid biopsy, a non-invasive approach that helps physicians diagnose and characterize specific diseases in patients. Here, we propose a proof-of-concept study combining the unique specificity in biomolecular recognition provided by surface-enhanced Raman spectroscopy (SERS) in combination with molecular dynamics (MD) simulations, which give leave to explore the biomolecular absorption mechanism on nanoparticle surfaces, in order to verify the traceability of two validated salivary indicators, i.e., interleukin-8 (IL-8) and lysozyme (LYZ), implicated in oropharyngeal squamous cell carcinoma (OSCC) and oral infection. This strategy simultaneously assures the detection and interpretation of protein biomarkers in saliva, ultimately opening a new route for the evolution of fast and accurate point-of-care SERS-based sensors of interest in precision medicine diagnostics.

Martini 3 Coarse-Grained Force Field for Cholesterol

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. The proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

Different conformational dynamics of SNARE protein Ykt6 among yeast and mammals

Ykt6 is one of the most conserved SNARE (N-ethylmaleimide-sensitive factor attachment protein receptor) proteins involved in multiple intracellular membrane trafficking processes. The membrane-anchoring function of Ykt6 has been elucidated to result from its conformational transition from a closed state to an open state. Two ways of regulating the conformational transition were proposed: the C-terminal lipidation and the phosphorylation at the SNARE core. Despite many aspects of common properties, Ykt6 displays differential cellular localizations and functional behaviors in different species, such as yeast, mammals, and worms. The structure-function relationship underlying these differences remains elusive. Here, we combined biochemical characterization, single-molecule FRET measurement, and molecular dynamics simulation to compare the conformational dynamics of yeast and rat Ykt6. Compared to rat Ykt6 (rYkt6), yeast Ykt6 (yYkt6) has more open conformations and could not bind dodecylphosphocholine that inhibits rYkt6 in the closed state. A point mutation T46L/Q57A was shown to be able to convert yYkt6 to a more closed and dodecylphosphocholine-bound state, where Leu46 contributes key hydrophobic interactions for the closed state. We also demonstrated that the phospho-mutation S174D could shift the conformation of rYkt6 to a more open state, but the corresponding mutation S176D in yYkt6 leads to a slightly more closed conformation. These observations shed light on the regulatory mechanism underlying the variations of Ykt6 functions across species.

In Silico Analyses of a Promising Drug Candidate for the Treatment of Amyotrophic Lateral Sclerosis Targeting Superoxide Dismutase I Protein

Amyotrophic lateral sclerosis (ALS) is the most prevalent motor neuron disorder in adults, which is associated with a highly disabling condition. To date, ALS remains incurable, and the only drugs approved by the FDA for its treatment confer a limited survival benefit. Recently, SOD1 binding ligand 1 (SBL-1) was shown to inhibit in vitro the oxidation of a critical residue for SOD1 aggregation, which is a central event in ALS-related neurodegeneration. In this work, we investigated the interactions between SOD1 wild-type and its most frequent variants, i.e., A4V (NP_000445.1:p.Ala5Val) and D90A (NP_000445.1:p.Asp91Val), with SBL-1 using molecular dynamics (MD) simulations. The pharmacokinetics and toxicological profile of SBL-1 were also characterized in silico. The MD results suggest that the complex SOD1-SBL-1 remains relatively stable and interacts within a close distance during the simulations. This analysis also suggests that the mechanism of action proposed by SBL-1 and its binding affinity to SOD1 may be preserved upon mutations A4V and D90A. The pharmacokinetics and toxicological assessments suggest that SBL-1 has drug-likeness characteristics with low toxicity. Our findings, therefore, suggested that SBL-1 may be a promising strategy to treat ALS based on an unprecedented mechanism, including for patients with these frequent mutations.

Discovering potential inhibitors against SARS-CoV-2 by targeting Nsp13 Helicase

The rise in the incidence of COVID-19 as a result of SARS-CoV-2 infection has threatened public health globally. Till now, there have been no proper prophylactics available to fight COVID-19, necessitating the advancement and evolution of effective curative against SARS-CoV-2. This study aimed at the nonstructural protein 13 (nsp13) helicase as a promising target for drug development against COVID-19. A unique collection of nucleoside analogs was screened against the SARS-CoV-2 helicase protein, for which a molecular docking experiment was executed to depict the selected ligand's binding affinity with the SARS-CoV-2 helicase proteins. Simultaneously, molecular dynamic simulations were performed to examine the protein's binding site's conformational stability, flexibility, and interaction with the ligands. Key nucleoside ligands were selected for pharmacokinetic analysis based on their docking scores. Selected ligands (cordycepin and pritelivir) showed excellent pharmacokinetics and were well stabilized at the proteins' binding site throughout the MD simulation. We have also performed binding free energy analysis or the binding characteristics of ligands with Nsp13 by using MM-PBSA and MM-GBSA. Free energy calculation by MM-PBSA and MM-GBSA analysis suggests that pritelivir may work as viable therapeutics for efficient drug advancement against SARS-CoV-2 Nsp13 helicase, potentially arresting the SARS-CoV-2 replication.Communicated by Ramaswamy H. Sarma.

The antiviral efficacies of small-molecule inhibitors against respiratory syncytial virus based on the F protein

OBJECTIVES

Respiratory syncytial virus (RSV) infection is one of the three most common causes of death in the infants, pre-schoolers, immunocompromised patients and elderly individuals due to many complications and lack of specific treatment. During RSV infection, the fusion protein (F protein) mediates the fusion of the virus envelope with the host cell membrane. Therefore, the F protein is an effective target for viral inhibition.

METHODS

We identified potential small-molecule inhibitors against RSV-F protein for the treatment of RSV infection using virtual screening and molecular dynamics (MD) simulations. The CCK8 assay was used to determine the cytotoxicity and quantitative RT-PCR and indirect fluorescence assay (IFA) were used to determine the viral replication and RSV-induced inflammation in vitro. An RSV-infected mouse model was established, and viral replication was assayed using real-time quantitative PCR and IFA. Virus-induced complications were also examined using histopathological analysis, airway resistance and the levels of IL-1β, IL-6 and TNF-α.

RESULTS

The top three potential inhibitors against the RSV-F protein were screened from the FDA-approved drug database. Z65, Z85 and Z74 significantly inhibited viral replication and RSV-induced inflammation. They also significantly alleviated RSV infection and RSV-induced complications in vivo. Z65 and Z85 had no cytotoxicity and better anti-RSV effects than Z74.

CONCLUSIONS

Z65 and Z85 may be suitable candidates for the treatment of RSV and serve as the basis for the development of new drugs.

Pharmacophore-Oriented Identification of Potential Leads as CCR5 Inhibitors to Block HIV Cellular Entry

Cysteine-cysteine chemokine receptor 5 (CCR5) has been discovered as a co-receptor for cellular entry of human immunodeficiency virus (HIV). Moreover, the role of CCR5 in a variety of cancers and various inflammatory responses was also discovered. Despite the fact that several CCR5 antagonists have been investigated in clinical trials, only Maraviroc has been licensed for use in the treatment of HIV patients. This indicates that there is a need for novel CCR5 antagonists. Keeping this in mind, the present study was designed. The active CCR5 inhibitors with known IC₅₀ value were selected from the literature and utilized to develop a ligand-based common feature pharmacophore model. The validated pharmacophore model was further used for virtual screening of drug-like databases obtained from the Asinex, Specs, InterBioScreen, and Eximed chemical libraries. Utilizing computational methods such as molecular docking studies, molecular dynamics simulations, and binding free energy calculation, the binding mechanism of selected inhibitors was established. The identified Hits not only showed better binding energy when compared to Maraviroc, but also formed stable interactions with the key residues and showed stable behavior throughout the 100 ns MD simulation. Our findings suggest that Hit1 and Hit2 may be potential candidates for CCR5 inhibition, and, therefore, can be considered for further CCR5 inhibition programs.

An understanding of coronavirus and exploring the molecular dynamics simulations to find promising candidates against the Mpro of nCoV to combat the COVID-19: A systematic review

The first infection case of new coronavirus was reported at the end of 2019 and after then, the cases are reported in all nations across the world in a very short period. Further, the regular news of mutations in the virus has made life restricted with appropriate behavior. To date, a new strain (Omicron and its new subvariant Omicron XE) has brought fear amongst us due to a higher trajectory of increase in the number of cases. The researchers thus started giving attention to this viral infection and discovering drug-like candidates to cure the infections. Finding a drug for any viral infection is not an easy task and takes plenty of time. Therefore, computational chemistry/bioinformatics is followed to get promising molecules against viral infection. Molecular dynamics (MD) simulations are being explored to get drug candidates in a short period. The molecules are screened via molecular docking, which provides preliminary information which can be further verified by molecular dynamics (MD) simulations. To understand the change in structure, MD simulations generated several trajectories such as root mean square deviation (RMSD), root mean square fluctuation (RMSF), hydrogen bonding, and radius of gyration for the main protease (Mpro) of the new coronavirus (nCoV) in the presence of small molecules. Additionally, change in free energy for the formation of complex of Mpro of nCoV with the small molecule can be determined by applying molecular mechanics with generalized born and surface area solvation (MM-GBSA). Thus, the promising molecules can be further explored for clinical trials to combat coronavirus disease-19 (COVID-19).

Identification of potential inhibitors, conformational dynamics, and mechanistic insights into mutant Kirsten rat sarcoma virus (G13D) driven cancers

The mutations at the hotspot region of K-Ras result in the progression of cancer types. Our study aimed to explore the small molecule inhibitors against the G13D mutant K-Ras model with anti-cancerous activity from food and drug administration (FDA)-approved drug compounds. We implemented several computational strategies such as pharmacophore-based virtual screening, molecular docking, absorption, distribution, metabolism and excretion features, and molecular simulation to ensure the identified hit compounds have potential efficacy against G13D K-Ras. We found that the FDA-approved compounds, namely, azelastine, dihydrocodeine, paroxetine, and tramadol, are potential candidates to inhibit the action of G13D mutant K-Ras. All four compounds exhibited similar binding patterns of sotorasib, and a structural binding mechanism with significant hydrophobic contacts. The descriptor features from the QikProp of all four compounds are within allowable limits compared to sotorasib drug. Consequently, a molecular simulation result emphasized that the dihydrocodeine and tramadol exhibited less fluctuation, minimal basin, significant h-bonds, and potent inhibition against G13D K-Ras. As a result, the current research identifies prospective K-Ras inhibitors that could be further improved with biochemical analysis for precision medicine against K-Ras-driven cancers.

Conformational Dynamics of Glucagon-like Peptide-2 with Different Electric Field

Molecular dynamics (MD) simulation was used to study the influence of electric field on Glucagon-like Peptide-2 (GLP-2). Different electric field strengths (0 V/nm ≤ E ≤ 1 V/nm) were mainly carried out on GLP-2. The structural changes in GLP-2 were analyzed by the Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), Radius of Gyration (Rg), Solvent Accessible Surface Area (SASA), Secondary Structure and the number of hydrogen bonds. The stable α—helix structure of GLP-2 was unwound and transformed into an unstable Turn and Coil structure since the stability of the GLP-2 protein structure was reduced under the electric field. Our results show that the degree of unwinding of the GLP-2 structure was not linearly related to the electric field intensity. E = 0.5 V/nm was a special point where the degree of unwinding of the GLP-2 structure reached the maximum at this electric field strength. Under a weak electric field, E < 0.5 V/nm, the secondary structure of GLP-2 becomes loose, and the entropy of the chain increases. When E reaches a certain value (E > 0.5 V/nm), the electric force of the charged residues reaches equilibrium, along the z-direction. Considering the confinement of moving along another direction, the residue is less free. Thus, entropy decreases and enthalpy increases, which enhance the interaction of adjacent residues. It is of benefit to recover hydrogen bonds in the middle region of the protein. These investigations, about the effect of an electric field on the structure of GLP-2, can provide some theoretical basis for the biological function of GLP-2 in vivo.

3D-QSAR-Based Pharmacophore Modeling, Virtual Screening, and Molecular Dynamics Simulations for the Identification of Spleen Tyrosine Kinase Inhibitors

Spleen tyrosine kinase (SYK) is an essential mediator of immune cell signaling and has been anticipated as a therapeutic target for autoimmune diseases, notably rheumatoid arthritis, allergic rhinitis, asthma, and cancers. Significant attempts have been undertaken in recent years to develop SYK inhibitors; however, limited success has been achieved due to poor pharmacokinetics and adverse effects of inhibitors. The primary goal of this research was to identify potential inhibitors having high affinity, selectivity based on key molecular interactions, and good drug-like properties than the available inhibitor, fostamatinib. In this study, a 3D-QSAR model was built for SYK based on known inhibitor IC₅₀ values. The best pharmacophore model was then used as a 3D query to screen a drug-like database to retrieve hits with novel chemical scaffolds. The obtained compounds were subjected to binding affinity prediction using the molecular docking approach, and the results were subsequently validated using molecular dynamics (MD) simulations. The simulated compounds were ranked according to binding free energy (ΔG), and the binding affinity was compared with fostamatinib. The binding mode analysis of selected compounds revealed that the hit compounds form hydrogen bond interactions with hinge region residue Ala451, glycine-rich loop residue Lys375, Ser379, and DFG motif Asp512. Identified hits were also observed to form a desirable interaction with Pro455 and Asn457, the rare feature observed in SYK inhibitors. Therefore, we argue that identified hit compounds ZINC98363745, ZINC98365358, ZINC98364133, and ZINC08789982 may help in drug design against SYK.

Determination of Hybrid TSPO Ligands with Minimal Impact of SNP (rs6971) through Molecular Docking and MD Simulation Study

Background:

In an endeavor to ascertain high-affinity TSPO ligands with minimal single nucleotide polymorphism (SNP), six hybrid molecules have been identified as new leads for future inflammation PET imaging.

Objective:

Genesis for chemical design was encouraged from structural families of well-known ligands FEBMP and PBR28/ DAA1106 that have demonstrated remarkable TSPO binding characteristics.

Methods:

All proposed hybrid ligands 1-6 are subjected to molecular docking and simulation studies with wild and mutant protein to study their interactions, binding, consistency of active conformations and are correlated with well-established TSPO ligands.

Results:

Each hybrid ligand demonstrate better docking score > -11.00 kcal/mol with TSPO with respect to gold standard PK11195, i.e., -11.00 kcal/mol for 4UC3 and -12.94 kcal/mol for 4UC1. On comparison with FEBMP and GE-180 (-12.57, -7.24 kcal/mol for 4UC3 and -14.10, -11.32 kcal/mol for 4UC1), ligand 3 demonstrates maximum docking energy (> -15.50 kcal/mol) with minimum SNP (0.26 kcal/mol).

Discussion:

Presence of strong hydrogen bond Arg148-3.27Å (4UC1) and Trp50-2.43Å, Asp28- 2.57Å (4UC3) apart from short-range interactions, including π–π interactions with the aromatic residues, such as (Trp39, Phe46, Trp135) and (Trp39, Trp108), attributes towards its strong binding.

Conclusion:

Utilizing the results of binding energy, we concluded stable complex formation of these hybrid ligands that could bind to TSPO with the least effect of SNP with similar interactions to known ligands. Overall, ligand 3 stands out as the best ligand having insignificant deviations per residue of protein that can be further explored and assessed in detail for future inflammation PET application after subsequent detailed biological evaluation.

Design of new glycosyl-O-fipronil conjugates with improved hydrolysis efficiency assisted by molecular simulations

BACKGROUND

In a previous study, we showed that two glycosyl-pesticide conjugates with a β-d-glucoside moiety, N-{3-cyano-1-[2,6-dichloro-4-(trifluoromethyl) phenyl]-4-[(trifluoromethyl)-sulfinyl]-1H-pyrazol-5-yl}-2-aminoethyl-β-d-glucopyranoside (GOF) and N-{3-cyano-1-[2,6-dichloro-4-(trifluoromethyl)phenyl]-4-[(trifluoromethyl) sulfinyl]-1H-pyrazol-5-yl}-1-(2-triazolethyl-β-d-glucopyranoside)-1H-1,2,3-triazole-4-methanamine (GOTF), can move in the phloem and be hydrolyzed by β-glucosidase at different rates. Simulations were carried out to investigate differences in the hydrolysis process in GOF, GOTF and p-nitrophenyl β-d-glucopyranoside (pNPG). A new series of glycosyl-O-fipronil conjugates was then designed and synthesized based on the simulation results. The phloem mobilities of the new conjugates were examined using a Ricinus model, and their hydrolysis efficiencies based on β-glucosidase were determined.

RESULTS

New glycosyl-O-fipronil conjugates GOE2-6 were designed and synthesized. To reduce steric hindrance, the conjugating site of the glycone moiety was moved to the 4'-sulfonyl group on the pyrrole ring. As a result, the hydrolysis efficiencies of the new conjugates were significantly improved, with GOE4 having the highest hydrolysis efficiency. All five conjugates could be transported in Ricinus phloem sap, consistent with previously studied glycosyl-O-fipronil conjugates. The insecticidal activities of the conjugates were tested against Plutella xylostella.

CONCLUSION

A strategy for the development of new phloem-mobile pesticides was proposed: linking a glycosyl group to the existing pesticide structure with a linear alkyl connection approximately four carbons in length. The resultant conjugates feature not only good phloem mobility, but also potential high bioactivity due to the efficient release of active pesticide components under the action of glucosidase. © 2022 Society of Chemical Industry.

In silico analyses of acetylcholinesterase (AChE) and its genetic variants in interaction with the anti-Alzheimer drug Rivastigmine

Alzheimer's disease (AD) is the leading cause of dementia worldwide. Despite causing great social and economic impact, there is currently no cure for AD. The most effective therapy to manage AD symptoms is based on acetylcholinesterase inhibitors (AChEi), from which rivastigmine presented numerous benefits. However, mutations in AChE, which affect approximately 5% of the population, can modify protein structure and function, changing the individual response to Alzheimer's treatment. In this study, we performed computer simulations of AChE wild type and variants R34Q, P135A, V333E, and H353N, identified by one or more genome-wide association studies, to evaluate their effects on protein structure and interaction with rivastigmine. The functional effects of AChE variants were predicted using eight machine learning algorithms, while the evolutionary conservation of AChE residues was analyzed using the ConSurf server. Autodock4.2.6 was used to predict the binding modes for the hAChE-rivastigmine complex, which is still unknown. Molecular dynamics (MD) simulations were performed in triplicates for the AChE wild type and mutants using the GROMACS packages. Among the analyzed variants, P135A was classified as deleterious by all the functional prediction algorithms, in addition to occurring at highly conserved positions, which may have harmful consequences on protein function. The molecular docking results suggested that rivastigmine interacts with hAChE at the upper active-site gorge, which was further confirmed by MD simulations. Our MD findings also suggested that the complex hAChE-rivastigmine remains stable over time. The essential dynamics revealed flexibility alterations at the active-site gorge upon mutations P135A, V333E, and H353N, which may lead to strong and nonintuitive consequences to hAChE binding. Nonetheless, similar binding affinities were registered in the MMPBSA analysis for the hAChE wild type and variants when complexed to rivastigmine. Finally, our findings indicated that the rivastigmine binding to hAChE is an energetically favorable process mainly driven by negatively charged amino acids.

Computational Insights of Unfolding of N-Terminal Domain of TDP-43 Reveal the Conformational Heterogeneity in the Unfolding Pathway

TDP-43 proteinopathies is a disease hallmark that characterizes amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD). The N-terminal domain of TDP-43 (NTD) is important to both TDP-43 physiology and TDP-43 proteinopathy. However, its folding and dimerization process is still poorly characterized. In the present study, we have investigated the folding/unfolding of NTD employing all-atom molecular dynamics (MD) simulations in 8 M dimethylsulfoxide (DMSO) at high temperatures. The MD results showed that the unfolding of the NTD at high temperature evolves through the formation of a number of conformational states differing in their stability and free energy. The presence of structurally heterogeneous population of intermediate ensembles was further characterized by the different extents of solvent exposure of Trp80 during unfolding. We suggest that these non-natives unfolded intermediate ensembles may facilitate NTD oligomerization and subsequently TDP-43 oligomerization, which might lead to the formation of irreversible pathological aggregates, characteristics of disease pathogenesis.

Diterpenes/Diterpenoids and Their Derivatives as Potential Bioactive Leads against Dengue Virus: A Computational and Network Pharmacology Study

Dengue fever is a dangerous infectious endemic disease that affects over 100 nations worldwide, from Africa to the Western Pacific, and is caused by the dengue virus, which is transmitted to humans by an insect bite of Aedes aegypti. Millions of citizens have died as a result of dengue fever and dengue hemorrhagic fever across the globe. Envelope (E), serine protease (NS3), RNA-directed RNA polymerase (NS5), and non-structural protein 1 (NS1) are mostly required for cell proliferation and survival. Some of the diterpenoids and their derivatives produced by nature possess anti-dengue viral properties. The goal of the computational study was to scrutinize the effectiveness of diterpenoids and their derivatives against dengue viral proteins through in silico study. Methods: molecular docking was performed to analyze the binding affinity of compounds against four viral proteins: the envelope (E) protein, the NS1 protein, the NS3 protein, and the NS5 protein. Results: among the selected drug candidates, triptolide, stevioside, alepterolic acid, sphaeropsidin A, methyl dodovisate A, andrographolide, caesalacetal, and pyrimethamine have demonstrated moderate to good binding affinities (-8.0 to -9.4 kcal/mol) toward the selected proteins: E protein, NS3, NS5, and NS1 whereas pyrimethamine exerts -7.5, -6.3, -7.8, and -6.6 kcal/mol with viral proteins, respectively. Interestingly, the binding affinities of these lead compounds were better than those of an FDA-approved anti-viral medication (pyrimethamine), which is underused in dengue fever. Conclusion: we can conclude that diterpenoids can be considered as a possible anti-dengue medication option. However, in vivo investigation is recommended to back up the conclusions of this study.

Computational investigation to identify potent inhibitors of the GTPase-Kirsten RAt sarcoma virus (K-Ras) mutants G12C and G12D

K-Ras mutations are frequent in various cancer types, and according to recent research, K-Ras possesses four-drug targeting sites. This increased our interest in finding potential small molecule inhibitors with anticancer activity to treat K-Ras-driven cancers. We utilized integrated bioinformatic strategies, such as XP docking, MM-GBSA, cell-line cytotoxicity prediction, ADMET, and molecular simulation, to discover potential inhibitors of G12C and G12D mutants compared to sotorasib, which is a recent FDA-approved inhibitor of G12C. We identified compounds, such as flupentixol, amlodipine, and fluvoxamine, for the G12C mutant and paroxetine, flupentixol, and zuclopenthixol for the G12D mutant with significant inhibitory functions. All five compounds bound to the H95 cryptic groove of mutant K-Ras with high efficiency and, like sotorasib, retained a novel binding mechanism with additional hydrophobic interactions at the molecular level. Furthermore, the simulation studies suggested that the binding of flupentixol and amlodipine to G12C stabilizes switch I and switch II. In contrast, paroxetine and flupentixol to G12D showed a similar trend compared to sotorasib complexes. Thus, despite the very dynamic functionality of K-Ras switches I and II, the binding of shortlisted compounds is highly stable. Therefore, the reported study provides potential drug candidates for K-Ras inhibition that can be further developed with in vitro and in vivo evidence for targeted therapy.

Identification of CDK7 Inhibitors from Natural Sources Using Pharmacoinformatics and Molecular Dynamics Simulations

The cyclin-dependent kinase 7 (CDK7) plays a crucial role in regulating the cell cycle and RNA polymerase-based transcription. Overexpression of this kinase is linked with various cancers in humans due to its dual involvement in cell development. Furthermore, emerging evidence has revealed that inhibiting CDK7 has anti-cancer effects, driving the development of novel and more cost-effective inhibitors with enhanced selectivity for CDK7 over other CDKs. In the present investigation, a pharmacophore-based approach was utilized to identify potential hit compounds against CDK7. The generated pharmacophore models were validated and used as 3D queries to screen 55,578 natural drug-like compounds. The obtained compounds were then subjected to molecular docking and molecular dynamics simulations to predict their binding mode with CDK7. The molecular dynamics simulation trajectories were subsequently used to calculate binding affinity, revealing four hits-ZINC20392430, SN00112175, SN00004718, and SN00262261-having a better binding affinity towards CDK7 than the reference inhibitors (CT7001 and THZ1). The binding mode analysis displayed hydrogen bond interactions with the hinge region residues Met94 and Glu95, DFG motif residue Asp155, ATP-binding site residues Thr96, Asp97, and Gln141, and quintessential residue outside the kinase domain, Cys312 of CDK7. The in silico selectivity of the hits was further checked by docking with CDK2, the close homolog structure of CDK7. Additionally, the detailed pharmacokinetic properties were predicted, revealing that our hits have better properties than established CDK7 inhibitors CT7001 and THZ1. Hence, we argue that proposed hits may be crucial against CDK7-related malignancies.

Effects of temperature and pH on the structure of a metalloprotease from Lactobacillus fermentum R6 isolated from Harbin dry sausages and molecular docking between protease and meat protein

BACKGROUND

Microbial protease can interact with meat protein in fermented meat products at a certain pH and temperature. To investigate the effects of various pH values and temperatures on the structural characteristics of Lactobacillus fermentum R6 protease, which was isolated from Harbin dry sausages, spectroscopy techniques and molecular dynamics were utilized to evaluate structural changes.

RESULTS

The protease exhibited a stable spatial structure at pH 7 and 40 °C, and the extension of the protease structure was also promoted. Although the structure of the protease could be changed or destroyed by pH 8 and 70 °C, it was mainly determined by the changes of secondary and tertiary structures such as α-helix, β-sheet, β-turn and random coil. In addition, carbonyl vibration, -NH vibration, C-H stretching vibration and disulphide bonds were present in L. fermentum R6 protease under various pH and temperature conditions. Molecular docking showed that the protease can interact with myosin light chain, myosin heavy chain, actin and myoglobin.

CONCLUSION

The protease can maintain stable structure and interact with meat protein, which reflected certain application prospects in the fermentation of Harbin dry sausages. © 2021 Society of Chemical Industry.

Peptide Platform as a Powerful Tool in the Fight against COVID-19

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in a global pandemic causing over 195 million infections and more than 4 million fatalities as of July 2021.To date, it has been demonstrated that a number of mutations in the spike glycoprotein (S protein) of SARS-CoV-2 variants of concern abrogate or reduce the neutralization potency of several therapeutic antibodies and vaccine-elicited antibodies. Therefore, the development of additional vaccine platforms with improved supply and logistic profile remains a pressing need. In this work, we have validated the applicability of a peptide-based strategy focused on a preventive as well as a therapeutic purpose. On the basis of the involvement of the dipeptidyl peptidase 4 (DPP4), in addition to the angiotensin converting enzyme 2 (ACE2) receptor in the mechanism of virus entry, we analyzed peptides bearing DPP4 sequences by protein-protein docking and assessed their ability to block pseudovirus infection in vitro. In parallel, we have selected and synthetized peptide sequences located within the highly conserved receptor-binding domain (RBD) of the S protein, and we found that RBD-based vaccines could better promote elicitation of high titers of neutralizing antibodies specific against the regions of interest, as confirmed by immunoinformatic methodologies and in vivo studies. These findings unveil a key antigenic site targeted by broadly neutralizing antibodies and pave the way to the design of pan-coronavirus vaccines.

A Plasmodium falciparum protein tyrosine phosphatase inhibitor identified from the ChEMBL-NTD database blocks parasite growth

Post‐translational modifications, especially reversible phosphorylation, are among the most common mechanisms that regulate protein function and biological processes in Plasmodium species. Of the Plasmodium phosphatases, phosphatase of regenerating liver (PfPRL) is secreted and is an essential phosphatase. Here, we expressed PfPRL in a heterologous expression system, and then purified and characterized its phosphatase activity. We found that Novartis_003209, a previously identified inhibitor, inhibited the PfPRL phosphatase activity of recombinant PfPRL and blocked in vitro parasite growth in a dose‐dependent manner. Further, in silico docking analysis of Novartis_003209 with all four P. falciparum tyrosine phosphatases (PTP) demonstrated that Novartis_003209 is a Plasmodium PTP inhibitor. Overall, our results identify a scaffold as a potential starting point to design a PTP‐specific inhibitor.

Computational Design of Novel Allosteric Inhibitors for Plasmodium falciparum DegP

The serine protease, DegP exhibits proteolytic and chaperone activities, essential for cellular protein quality control and normal cell development in eukaryotes. The P. falciparum DegP is essential for the parasite survival and required to combat the oscillating thermal stress conditions during the infection, protein quality checks and protein homeostasis in the extra-cytoplasmic compartments, thereby establishing it as a potential target for drug development against malaria. Previous studies have shown that diisopropyl fluorophosphate (DFP) and the peptide SPMFKGV inhibit E. coli DegP protease activity. To identify novel potential inhibitors specific to PfDegP allosteric and the catalytic binding sites, we performed a high throughput in silico screening using Malaria Box, Pathogen Box, Maybridge library, ChEMBL library and the library of FDA approved compounds. The screening helped identify five best binders that showed high affinity to PfDegP allosteric (T0873, T2823, T2801, RJC02337, CD00811) and the catalytic binding site (T0078L, T1524, T2328, BTB11534 and 552691). Further, molecular dynamics simulation analysis revealed RJC02337, BTB11534 as the best hits forming a stable complex. WaterMap and electrostatic complementarity were used to evaluate the novel bio-isosteric chemotypes of RJC02337, that led to the identification of 231 chemotypes that exhibited better binding affinity. Further analysis of the top 5 chemotypes, based on better binding affinity, revealed that the addition of electron donors like nitrogen and sulphur to the side chains of butanoate group are more favoured than the backbone of butanoate group. In a nutshell, the present study helps identify novel, potent and Plasmodium specific inhibitors, using high throughput in silico screening and bio-isosteric replacement, which may be experimentally validated.

Computational Simulations Identified Marine-Derived Natural Bioactive Compounds as Replication Inhibitors of SARS-CoV-2

The rapid spread of COVID-19, caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a worldwide health emergency. Unfortunately, to date, a very small number of remedies have been to be found effective against SARS-CoV-2 infection. Therefore, further research is required to achieve a lasting solution against this deadly disease. Repurposing available drugs and evaluating natural product inhibitors against target proteins of SARS-CoV-2 could be an effective approach to accelerate drug discovery and development. With this strategy in mind, we derived Marine Natural Products (MNP)-based drug-like small molecules and evaluated them against three major target proteins of the SARS-CoV-2 virus replication cycle. A drug-like database from MNP library was generated using Lipinski's rule of five and ADMET descriptors. A total of 2,033 compounds were obtained and were subsequently subjected to molecular docking with 3CLpro, PLpro, and RdRp. The docking analyses revealed that a total of 14 compounds displayed better docking scores than the reference compounds and have significant molecular interactions with the active site residues of SARS-CoV-2 virus targeted proteins. Furthermore, the stability of docking-derived complexes was analyzed using molecular dynamics simulations and binding free energy calculations. The analyses revealed two hit compounds against each targeted protein displaying stable behavior, binding affinity, and molecular interactions. Our investigation identified two hit compounds against each targeted proteins displaying stable behavior, higher binding affinity and key residual molecular interactions, with good in silico pharmacokinetic properties, therefore can be considered for further in vitro studies.

Artificial Intelligence-Based Application to Explore Inhibitors of Neurodegenerative Diseases

Neuroinflammation is a common factor in neurodegenerative diseases, and it has been demonstrated that galectin-3 activates microglia and astrocytes, leading to inflammation. This means that inhibition of galectin-3 may become a new strategy for the treatment of neurodegenerative diseases. Based on this motivation, the objective of this study is to explore an integrated new approach for finding lead compounds that inhibit galectin-3, by combining universal artificial intelligence algorithms with traditional drug screening methods. Based on molecular docking method, potential compounds with high binding affinity were screened out from Chinese medicine database. Manifold artificial intelligence algorithms were performed to validate the docking results and further screen compounds. Among all involved predictive methods, the deep learning-based algorithm made 500 modeling attempts, and the square correlation coefficient of the best trained model on the test sets was 0.9. The XGBoost model reached a square correlation coefficient of 0.97 and a mean square error of only 0.01. We switched to the ZINC database and performed the same experiment, the results showed that the compounds in the former database showed stronger affinity. Finally, we further verified through molecular dynamics simulation that the complex composed of the candidate ligand and the target protein showed stable binding within 100 ns of simulation time. In summary, combined with the application based on artificial intelligence algorithms, we unearthed the active ingredients 1,2-Dimethylbenzene and Typhic acid contained in Crataegus pinnatifida and Typha angustata might be the effective inhibitors of neurodegenerative diseases. The high prediction accuracy of the models shows that it has practical application value on small sample data sets such as drug screening.

Mutational landscape of K-Ras substitutions at 12th position-a systematic molecular dynamics approach

K-Ras is a small GTPase and acts as a molecular switch by recruiting GEFs and GAPs, and alternates between the inert GDP-bound and the dynamic GTP-bound forms. The amino acid at position 12 of K-Ras is a hot spot for oncogenic mutations (G12A, G12C, G12D, G12R, G12S, and G12V), disturbing the active fold of the protein, leading to cancer development. This study aimed to investigate the potential conformational changes induced by these oncogenic mutations at the 12^th position, impairing GAP-mediated GTP hydrolysis. Comprehensive computational tools (iStable, FoldX, SNPeffect, DynaMut, and CUPSAT) were used to evaluate the effect of these six mutations on the stability of wild type K-Ras protein. The docking of GTP with K-Ras was carried out using AutoDock4.2, followed by molecular dynamics simulations. Furthermore, on comparison of binding energies between the wild type K-Ras and the six mutants, we have demonstrated that the G12A and G12V mutants exhibited the strongest binding efficiency compared to the other four mutants. Trajectory analyses of these mutations revealed that G12A encountered the least deviation, fluctuation, intermolecular H-bonds, and compactness compared to the wildtype, which was supported by the lower Gibbs free energy value. Our study investigates the molecular dynamics simulations of the mutant K-Ras forms at the 12^th position, which expects to provide insights about the molecular mechanisms involved in cancer development, and may serve as a platform for targeted therapies against cancer. Communicated by Ramaswamy H. Sarma.

The Relevance of Thimet Oligopeptidase in the Regulation of Energy Metabolism and Diet-Induced Obesity

Thimet oligopeptidase (EC 3.4.24.15; EP24.15; THOP1) is a potential therapeutic target, as it plays key biological functions in processing biologically functional peptides. The structural conformation of THOP1 provides a unique restriction regarding substrate size, in that it only hydrolyzes peptides (optimally, those ranging from eight to 12 amino acids) and not proteins. The proteasome activity of hydrolyzing proteins releases a large number of intracellular peptides, providing THOP1 substrates within cells. The present study aimed to investigate the possible function of THOP1 in the development of diet-induced obesity (DIO) and insulin resistance by utilizing a murine model of hyperlipidic DIO with both C57BL6 wild-type (WT) and THOP1 null (THOP1^−/−) mice. After 24 weeks of being fed a hyperlipidic diet (HD), THOP1^−/− and WT mice ingested similar chow and calories; however, the THOP1^−/− mice gained 75% less body weight and showed neither insulin resistance nor non-alcoholic fatty liver steatosis when compared to WT mice. THOP1^−/− mice had increased adrenergic-stimulated adipose tissue lipolysis as well as a balanced level of expression of genes and microRNAs associated with energy metabolism, adipogenesis, or inflammation. Altogether, these differences converge to a healthy phenotype of THOP1^−/− fed a HD. The molecular mechanism that links THOP1 to energy metabolism is suggested herein to involve intracellular peptides, of which the relative levels were identified to change in the adipose tissue of WT and THOP1^−/− mice. Intracellular peptides were observed by molecular modeling to interact with both pre-miR-143 and pre-miR-222, suggesting a possible novel regulatory mechanism for gene expression. Therefore, we successfully demonstrated the previously anticipated relevance of THOP1 in energy metabolism regulation. It was suggested that intracellular peptides were responsible for mediating the phenotypic differences that are described herein by a yet unknown mechanism of action.

Plasmodium Condensin Core Subunits SMC2/SMC4 Mediate Atypical Mitosis and Are Essential for Parasite Proliferation and Transmission

Condensin is a multi-subunit protein complex regulating chromosome condensation and segregation during cell division. In Plasmodium spp., the causative agent of malaria, cell division is atypical and the role of condensin is unclear. Here we examine the role of SMC2 and SMC4, the core subunits of condensin, during endomitosis in schizogony and endoreduplication in male gametogenesis. During early schizogony, SMC2/SMC4 localize to a distinct focus, identified as the centromeres by NDC80 fluorescence and chromatin immunoprecipitation sequencing (ChIP-seq) analyses, but do not form condensin I or II complexes. In mature schizonts and during male gametogenesis, there is a diffuse SMC2/SMC4 distribution on chromosomes and in the nucleus, and both condensin I and condensin II complexes form at these stages. Knockdown of smc2 and smc4 gene expression reveals essential roles in parasite proliferation and transmission. The condensin core subunits (SMC2/SMC4) form different complexes and may have distinct functions at various stages of the parasite life cycle.

Molecular Docking and Molecular Dynamics (MD) Simulation of Human Anti-Complement Factor H (CFH) Antibody Ab42 and CFH Polypeptide

An understanding of the interaction between the antibody and its targeted antigen and knowing of the epitopes are critical for the development of monoclonal antibody drugs. Complement factor H (CFH) is implied to play a role in tumor growth and metastasis. An autoantibody to CHF is associated with anti-tumor cell activity. The interaction of a human monoclonal antibody Ab42 that was isolated from a cancer patient with CFH polypeptide (pCFH) antigen was analyzed by molecular docking, molecular dynamics (MD) simulation, free energy calculation, and computational alanine scanning (CAS). Experimental alanine scanning (EAS) was then carried out to verify the results of the theoretical calculation. Our results demonstrated that the Ab42 antibody interacts with pCFH by hydrogen bonds through the Tyr315, Ser100, Gly33, and Tyr53 residues on the complementarity-determining regions (CDRs), respectively, with the amino acid residues of Pro441, Ile442, Asp443, Asn444, Ile447, and Thr448 on the pCFH antigen. In conclusion, this study has explored the mechanism of interaction between Ab42 antibody and its targeted antigen by both theoretical and experimental analysis. Our results have important theoretical significance for the design and development of relevant antibody drugs.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Identification of a Novel Anti-cancer Protein, FIP-bbo, from Botryobasidium botryosum and Protein Structure Analysis using Molecular Dynamic Simulation

Fungal immunoregulatory proteins (FIP) are effective small molecule proteins with broad-spectrum immunomodulatory and anti-cancer activities and can be potential agents for the development of clinical drugs and health food additives. In this study, a new member of FIP named FIP-bbo was obtained through Botryobasidium botryosum genome mining. FIP-bbo has the typical characteristics of FIP but is genetically distant from other FIPs. Recombinant FIP-bbo (rFIP-bbo) was produced in an optimized E. coli expression system, and the pure protein was isolated using a Ni-NTA column. Antineoplastic experiments suggested that FIP-bbo is similar to LZ-8 in inhibiting various cancer cells (Hela, Spac-1, and A549) at lower concentrations, but it is not as potent as LZ-8. The molecular mechanism by which FIP-bbo, FIP-fve, and LZ-8 are cytotoxic to cancer cells has been discussed based on molecular dynamics simulation. Point mutations that may improve the thermal stability of FIP-fve and FIP-bbo were predicted. These results not only present a new candidate protein for the development of anticancer adjuvants, but also provide an approach for designing FIPs with high anticancer activity.

Biochemical characterization of Plasmodium complement factors binding protein for its role in immune modulation

Complement system is the first line of human defence against intruding pathogens and is recognized as a potentially useful therapeutic target. Human malaria parasite Plasmodium employs a series of intricate mechanisms that enables it to evade different arms of immune system, including the complement system. Here, we show the expression of a multi-domain Plasmodium Complement Control Protein 1, PfCCp1 at asexual blood stages and its binding affinity with C3b as well as C4b proteins of human complement cascade. Using a biochemical assay, we demonstrate that PfCCp1 binds with complement factors and inhibits complement activation. Active immunization of mice with PfCCp1 followed by challenge with Plasmodium berghei resulted in the loss of biphasic growth of parasites and early death in comparison to the control group. The study also showed a role of PfCCp1 in modulating Toll-like receptor (TLR)-mediated signalling and effector responses on antigen-presenting cells. PfCCp1 binds with dendritic cells that down-regulates the expression of signalling molecules and pro-inflammatory cytokines, thereby dampening the TLR2-mediated signalling; hence acting as a potent immuno-modulator. In summary, PfCCp1 appears to be an important component of malaria parasite directed immuno-modulating strategies that promote the adaptive fitness of pathogens in the host.

A change in structural integrity of c-Kit mutant D816V causes constitutive signaling

Several signaling pathways, ligands, and genes that regulate proliferative and self-renewal properties of the Hematopoietic Stem Cells (HSCs) have been studied meticulously. One of the signaling pathways that play a crucial role in the process of hematopoiesis is the Stem Cell Factor (SCF) mediated c-Kit pathway. The c-Kit is a Receptor Tyrosine Kinase (RTK), which is expressed in the cells including HSCs. It undergoes dimerization upon binding with its cognate ligand SCF. As a result, phosphorylation of the Juxtamembrane (JM) domain of c-Kit takes place at Tyr568 and Tyr570 residues. These phosphorylated residues become the docking sites for protein tyrosine phosphatases (PTPs) namely SHP-1 and SHP-2, which in turn cause dephosphorylation and negative regulation of the downstream signaling responsible for the cell proliferation. Interestingly, it has been reported that the mutation of c-Kit at D816V makes it independent of SCF stimulation and SHP-1/SHP-2 inhibition, thereby, causing its constitutive activation. The present study was commenced to elucidate the structural behavior of this mutation in the JM and A-loop region of c-Kit using Molecular Dynamics (MD) simulations of the wild-type and mutant c-Kit in unphosphorylated and phosphorylated states. The energy difference computed between the wild type and mutant (D816V) c-Kit, and protein-protein docking and complex analysis revealed the impact of this single residue mutation on the integrity dynamics of c-Kit that makes it independent of SHP-1/SHP-2 negative regulation.

High throughput in silico identification and characterization of Plasmodium falciparum PRL phosphatase inhibitors

Kinases and phosphatases are involved in many essential processes in Plasmodium lifecycle. Among the identified 67 Plasmodium falciparum phosphatases, Phosphatase of Regenerating Liver (PRL) family protein homolog, PfPRL, is an essential parasite tyrosine phosphatase. PfPRL is shown to be prenylated, secreted, and involved in the host invasion process. In the present study, a structure-based high throughput in silico screening of PfPRL binders, using ChEMBL-NTD compounds lead to the identification of nine compounds based on binding energy, Lipinski rule of five, and QED score. The most of the shortlisted compounds are known to inhibit parasite growth at a concentration (EC50) ≤2 μm in in vitro P. falciparum culture assays. MD simulations were carried out on the shortlisted nine potential enzyme-inhibitor complexes to analyze specificity, stability, and to calculate the free binding energies of the complexes. The study identifies PfPRL as one of the potential drug targets for selected ChEMBL-NTD compounds that may be exploited as a scaffold to develop novel antimalarials.

Exploring Heme and Hemoglobin Binding Regions of Plasmodium Heme Detoxification Protein for New Antimalarial Discovery

Hemoglobin degradation/hemozoin formation, essential steps in the Plasmodium life cycle, are targets of existing antimalarials. The pathway still offers vast possibilities to be explored for new antimalarial discoveries. Here, we characterize heme detoxification protein, PfHDP, a major protein involved in hemozoin formation, as a novel drug target. Using in silico and biochemical approaches, we identified two heme binding sites and a hemoglobin binding site in PfHDP. Treatment of Plasmodium falciparum 3D7 parasites with peptide corresponding to the hemoglobin binding domain in PfHDP resulted in food vacuole abnormalities similar to that seen with a cysteine protease inhibitor, E-64 (I-1). Screening of compounds that bound the modeled PfHDP structure in the heme/hemoglobin-binding pockets from Maybridge Screening Collection identified a compound, ML-2, that inhibited parasite growth in a dose-dependent manner, thus paving the way for testing its potential as a new drug candidate. These results provide functional insights into the role of PfHDP in Hz formation and further suggest that PfHDP could be an important drug target to combat malaria.

L-Asparaginase of Leishmania donovani: Metabolic target and its role in Amphotericin B resistance

Emergence of Amphotericin B (AmB) resistant Leishmania donovani has posed major therapeutic challenge against the parasite. Consequently, combination therapy aimed at multiple molecular targets, based on proteome wise network analysis has been recommended. In this regard we had earlier identified and proposed L-asparaginase of Leishmania donovani (LdAI) as a crucial metabolic target. Here we report that both LdAI overexpressing axenic amastigote and promastigote forms of L. donovani survives better when challenged with AmB as compared to wild type strain. Conversely, qRT-PCR analysis showed an upregulation of LdAI in both forms upon AmB treatment. Our data demonstrates the importance of LdAI in imparting immediate protective response to the parasite upon AmB treatment. In the absence of structural and functional information, we modeled LdAI and validated its solution structure through small angle X-ray scattering (SAXS) analysis. We identified its specific inhibitors through ligand and structure-based approach and characterized their effects on enzymatic properties (K_m, V_max, K_cat) of LdAI. We show that in presence of two of the inhibitors L1 and L2, the survival of L. donovani is compromised whereas overexpression of LdAI in these cells restores viability. Taken together, our results conclusively prove that LdAI is a crucial metabolic enzyme conferring early counter measure against AmB treatment by Leishmania.

Structural Insights into the Mode of Action of the Peptide Antibiotic Copsin

Defensins make up a class of cysteine-rich antimicrobial peptides, expressed by virtually all eukaryotes as part of their innate immune response. Because of their unique mode of action and rapid killing of pathogenic microbes, defensins are considered promising alternatives to clinically applied antibiotics. Copsin is a defensin-like peptide, previously identified in the mushroom Coprinopsis cinerea. It exerts its activity against a range of Gram-positive bacteria by binding to the peptidoglycan precursor lipid II and prevention of proper cell wall formation. In this study, we present a new workflow for the generation, production, and activity-driven selection of copsin derivatives, based on their expression in Pichia pastoris. One hundred fifty-two single-amino acid mutants and combinations thereof allowed the identification of k-copsin, a peptide variant exhibiting significantly enhanced activity against Bacillus subtilis and Staphylococcus aureus. Furthermore, we performed in silico characterizations of membrane interactions of copsin and k-copsin, in the presence and absence of lipid II. The molecular dynamics data highlighted a high variability in lipid II binding, with a preference for the MurNAc moiety with 47 and 35% of the total contacts for copsin and k-copsin, respectively. Mutated amino acids were located in loop regions of k-copsin and shown to be crucial in the perturbation of the bacterial membrane. These structural studies provide a better understanding of how defensins can be developed toward antibacterial therapies less prone to resistance issues.

Escape from humoral immunity is associated with treatment failure in HIV-1-infected patients receiving long-term antiretroviral therapy

Interindividual heterogeneity in the disease progression of HIV-1-infected patients receiving long-term antiretroviral therapy suggests that some host-related factors may have limited treatment efficacy. To understand the nature of factors contributing to treatment failure, we performed a retrospective cohort study of 45 chronically HIV-1-infected individuals sharing a similar demographics and route of infection, compared the differences between virologically suppressed (VS) and treatment failure (TF) patients with respect to clinical, immunological and virological characteristics. We found that the baseline diversity of HIV-1 env quasispecies was the major difference between VS and TF group, and higher baseline diversity in TF patients. We further predicted TF-related env mutations using a selection pressure-based approach, followed by an analysis of these mutations based on the available three-dimensional structures of gp120/gp41 or their complexes with neutralizing antibodies. Notably, almost all of the identified residues could be mapped to the epitopes of known HIV-1 neutralizing antibodies, especially the epitopes of broadly neutralizing antibodies, and these mutations tended to compromise antibody-antigen interactions. These results indicate that the escape of HIV-1 from host humoral immunity may play a direct role in TF in long-term antiretroviral-experienced patients and that based on env gene sequence of the viruses in the patients.

Molecular dynamics study of di-CF4 based reverse micelles in supercritical CO2

Reverse micelles (RMs) in supercritical CO₂ (scCO₂) are promising alternatives for organic solvents, especially when both polar and non-polar components are involved. Fluorinated surfactants, particularly double-chain fluorocarbon surfactants, are able to form well-structured RMs in scCO₂. The inherent self-assembly mechanisms of surfactants in scCO₂ are still subject to discussion. In this study, molecular dynamics simulations are performed to investigate the self-aggregation behavior of di-CF4 based RMs in scCO₂, and stable and spherical RMs are formed. The dynamics process and the self-assembly structure in the RMs reveal a three-step mechanism to form the RMs, that is, small RMs, rod-like RMs and fusion of the rod-like RMs. Hydrogen-bonds between headgroups and water molecules, and salt bridges linking Na⁺ ions, headgroups and water molecules enhance the interfacial packing efficiency of the surfactant. The results show that di-CF4 molecules have a high surfactant coverage at the RM interface, implying a high CO₂-philicity. This mainly results from bending of the short chain (C-COO-CH₂-(CF2)₃-CF3) due to the flexible carboxyl group. The microscopic insight provided in this study is helpful in understanding surfactant self-assembly phenomena and designing new CO₂-philic surfactants.

Identification of high-affinity anti-CD16A allotype-independent human antibody domains

CD16A (FcγRIIIA) is an activating receptor mostly expressed on natural killer (NK) cells and monocytes/macrophages. It can mediate antibody-dependent cell-mediated cytotoxicity (ADCC) through low-affinity interaction with human immunoglobulin G (IgG) Fc. It can also mediate cell lysis if NK cells are guided by bispecific killer cells engagers (BiKEs). BiKEs showed some success in clinical trials of cancer and are promising candidate therapeutics. However, currently reported BiKEs are based on antibody fragments (scFvs) of relatively large size. The CD16A-specific antibodies are also typically from animal origin. Decreasing the BiKE size could result in enhanced penetration into solid tumor and normal tissues, and using fully human antibodies could decrease the likelihood of immunogenicity. Here we report the identification and characterization of two antibody domains, D6 and E11, isolated from a very large human VH antibody domain library displayed on phage. D6 and E11 bound CD16A with EC₅₀ of 4nM and 8nM, respectively, but not other Fc gamma receptors (FcγRs) such as CD64 (FcγRI), CD32 (FcγRII) and CD16B (FcγRIIIB). They bound to both CD16A allotypes (158F,V) with equal affinity and competed with each other as well as with human IgG1 and the mouse anti-CD16A antibody 3G8. These and other results were used to build a molecular docking model predicting that D6 and E11 may bind to the CD16A membrane proximal D2 domain by interacting with its BC, C'E and EF loops. Importantly, cross-linked (bivalent) D6 and E11 induced secretion of IL-2 after binding to CD16A-expressing Jurkat T cells. The small size of these antibody domains combined with their high-affinity, specific, allotype-independent, activating interactions with CD16A could allow generation of novel highly effective BiKEs and other candidate protein therapeutics.