Very fast empirical prediction and rationalization of protein pKa values

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

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

Integration of terpesomes loaded Levocetrizine dihydrochloride gel as a repurposed cure for Methicillin-Resistant Staphylococcus aureus (MRSA)-Induced skin infection; D-optimal optimization, ex-vivo, in-silico, and in-vivo studies

The intention of this work is to assess the repurposed antimicrobial impact of Levocetirizine dihydrochloride (LVC), which is a well-known antihistaminic drug, in addition, to augment the antimicrobial effect by using terpene-enriched vesicles (TPs). To investigate how various parameters affect TPs aspects, TPs were made employing the ethanol-injection-method and optimized d-optimal design. The TPs were characterized based on their entrapment efficiency percentage (EE%), particle size (PS), polydispersity index (PDI), and zeta potential (ZP). The optimum TP was submitted to more examinations. The optimum TP (TP12) showed a spherical vesicle having an EE% of 66.39 ± 0.12%, PS of 243.3 ± 4.60 nm, PDI of 0.458 ± 0.003, and ZP of 24.2 ± 0.55 mV. The in-vitro release study results demonstrated that LVC is sustainedly liberated from the optimum TP compared to LVC-solution. The ex-vivo assessment showed that LVC was released in a more sustained manner from TPs-gel related to LVC solution, optimum TP, and LVC gel. Ex-vivo visualization by confocal laser scanning microscopy showed good deposition of the fluorescein-labeled TP. Further, the in-vitro anti-bacterial effect and biofilm inhibition and detachment assessment confirmed the potency of LVC against Methicillin-resistant-Staphylococcus-aureus (MRSA). The in-silico study demonstrated that the LVC has excellent stability with other ingredients combined with it in the TPs, further, it proved that LVC is a potential candidate for treating MRSA. In-vivo assessments revealed a good antimicrobial effect toward MRSA infection. Moreover, the histopathological evaluation confirmed the safety of using TPs-gel topically. In conclusion, MRSA-related skin infections may be treated using the LVC loaded TPs-gel as a promising system.

Identifying Drug Candidates for COVID-19 with Large-Scale Drug Screening

Papain-like protease (PL^pro) is critical to COVID-19 infection. Therefore, it is a significant target protein for drug development. We virtually screened a 26,193 compound library against the PL^pro of SARS-CoV-2 and identified several drug candidates with convincing binding affinities. The three best compounds all had better estimated binding energy than those of the drug candidates proposed in previous studies. By analyzing the docking results for the drug candidates identified in this and previous studies, we demonstrate that the critical interactions between the compounds and PL^pro proposed by the computational approaches are consistent with those proposed by the biological experiments. In addition, the predicted binding energies of the compounds in the dataset showed a similar trend as their IC₅₀ values. The predicted ADME and drug-likeness properties also suggested that these identified compounds can be used for COVID-19 treatment.

Molecular plasticity of CBM3 ancillary domain leads to conformational changes in the cellulose binding interface

Carbohydrate-binding modules (CBMs) constitute independently folded domains typically associated with carbohydrate-active enzymes (CAZymes). These modules are considered to have a rigid structure without notable conformational changes upon ligand binding, exhibiting a complementary topography in relation to the target carbohydrate. Herein, the high-resolution SAD-solved structure of a CBM from family 3 (BsCBM3) that binds to crystalline cellulose is reported in two crystalline forms. This module showed molecular plasticity with structural differences detected between the two crystalline forms and high RMSD values when compared to NMR ensemble of models. Pronounced structural variances were observed in the cellulose binding interface between NMR and XTAL structures, which were corroborated by molecular dynamics simulations. These findings support that family 3 CBMs targeting to cellulose are rather structurally dynamic modules than rigid entities, suggesting a potential role of conformational changes in polysaccharide recognition and modulation of enzyme activity.

The Laccase of Myrothecium roridum VKM F-3565: A New Look at Fungal Laccase Tolerance to Neutral and Alkaline Conditions

Most of the currently known fungal laccases show their maximum activity under acidic environmental conditions. It is known that a decrease in the activity of a typical laccase at neutral or alkaline pH values is the result of an increase in the binding of the hydroxide anion to the T2/T3 copper center, which prevents the transfer of an electron from the T1 Cu to the trinuclear copper center. However, evolutionary pressure has resolved the existing limitations in the catalytic mechanism of laccase, allowing such enzymes to be functionally active under neutral/alkaline pH conditions, thereby giving fungi an advantage for their survival. Combined molecular and biochemical studies, homological modeling, calculation of the electrostatic potential on the Connolly surface at pH 5.0 and 7.0, and structural analysis of the novel alkaliphilic laccase of Myrothecium roridum VKM F-3565 and alkaliphilic and acidophilic fungal laccases with a known structure allowed a new intramolecular channel near the one of the catalytic aspartate residues at T2-copper atom to be found. The amino acid residues of alkaliphilic laccases forming this channel can presumably serve as proton donors for catalytic aspartates under neutral conditions, thus ensuring proper functioning. For the first time for ascomycetous laccases, the production of new trimeric products of phenylpropanoid condensation under neutral conditions has been shown, which could have a potential for use in pharmacology.

Mechanism of Radical Initiation and Transfer in Class Id Ribonucleotide Reductase Based on Density Functional Theory

Class Id ribonucleotide reductase (RNR) is a newly discovered enzyme, which employs the dimanganese cofactor in the superoxidized state (Mn^III/Mn^IV) as the radical initiator. The dimanganese cofactor of class Id RNR in the reduced state (inactive) is clearly based on the crystal structure of the Fj-β subunit. However, the state of the dimanganese cofactor of class Id RNR in the oxidized state (active) is not known. The X-band EPR spectra have shown that the activated Fj-β subunit exists in two distinct complexes, 1 and 2. In this work, quantum mechanical/molecular mechanical calculations were carried out to study class Id RNR. First, we have determined that complex 2 contains a Mn^III-(μ-oxo)₂-Mn^IV cluster, and complex 1 contains a Mn^III-(μ-hydroxo/μ-oxo)-Mn^IV cluster. Then, based on the determined dimanganese cofactors, the mechanism of radical initiation and transfer in class Id RNR is revealed. The Mn^III-(μ-oxo)₂-Mn^IV cluster in complex 2 has not enough reduction potential to initiate radical transfer directly. Instead, it needs to be monoprotonated into Mn^III-(μ-hydroxo/μ-oxo)-Mn^IV (complex 1) before the radical transfer. The protonation state of μ-oxo can be regulated by changing the protein microenvironment, which is induced by the protein aggregation and separation of β subunits with α subunits. The radical transfer between the cluster of Mn^III-(μ-hydroxo/μ-oxo)-Mn^IV and Trp30 in the radical-transfer chain of the Fj-β subunit (Mn^III/Mn^IV ↔ His100 ↔ Asp194 ↔ Trp30 ↔ Arg99) is a water-mediated tri-proton-coupled electron transfer, which transfers proton from the ε-amino group of Lys71 to the carboxyl group of Glu97 via the water molecule Wat551 and the bridging μ-hydroxo ligand through a three-step reaction. This newly discovered proton-coupled electron-transfer mechanism in class Id RNR is different from those reported in the known Ia-Ic RNRs. The ε-amino group of Lys71, which serves as a proton donor, plays an important role in the radical transfer.

Synergistic Antimicrobial Action of Lactoferrin-Derived Peptides and Quorum Quenching Enzymes

Combined use of various antimicrobial peptides (AMPs) with enzymes that hydrolyze the signaling molecules of the resistance mechanism of various microorganisms, quorum sensing (QS), to obtain effective antimicrobials is one of the leading approaches in solving the antimicrobial resistance problem. Our study investigates the lactoferrin-derived AMPs, lactoferricin (Lfcin), lactoferampin and Lf(1-11), as potential partners for combination with enzymes hydrolyzing lactone-containing QS molecules, the hexahistidine-containing organophosphorus hydrolase (His6-OPH) and penicillin acylase, to obtain effective antimicrobial agents with a scope of practical application. The possibility of the effective combination of selected AMPs and enzymes was first investigated in silico using molecular docking method. Based on the computationally obtained results, His6-OPH/Lfcin combination was selected as the most suitable for further research. The study of physical-chemical characteristics of His6-OPH/Lfcin combination revealed the stabilization of enzymatic activity. A notable increase in the catalytic efficiency of action of His6-OPH in combination with Lfcin in the hydrolysis of paraoxon, N-(3-oxo-dodecanoyl)-homoserine lactone and zearalenone used as substrates was established. Antimicrobial efficiency of His6-OPH/Lfcin combination was determined against various microorganisms (bacteria and yeasts) and its improvement was observed as compared to AMP without enzyme. Thus, our findings demonstrate that His6-OPH/Lfcin combination is a promising antimicrobial agent for practical application.

Probing ion binding in the selectivity filter of the Cav1.1 channel with molecular dynamics

Cav1.1 is the voltage-gated calcium channel essential for the contraction of skeletal muscles upon membrane potential changes. Structural determination of the Cav1.1 channel opens the avenue toward understanding of the structure-function relationship of voltage-gated calcium channels. Here, we show that there exist two Ca2+-binding sites, termed S1 and S2, within the selectivity filter of Cav1.1 through extensive molecular dynamics simulations on various initial ion arrangement configurations. The formation of both binding sites is associated with the four Glu residues (Glu292/614/1014/1323) that constitute the so-called EEEE locus. At the S1 site near the extracellular side, the Ca2+ ion is coordinated with the negatively charged carboxylic groups of these Glu residues and of the Asp615 residue either in a direct way or via an intermediate water molecule. At the S2 site, Ca2+ binding shows two distinct states: an upper state involving two out of the four Glu residues in the EEEE locus and a lower state involving only one Glu residue. In addition, there exist two recruitment sites for Ca2+ above the entrance of the filter. These findings promote the understanding of mechanism for ion permeation and selectivity in calcium channels.

Predicting the insecticide-driven mutations in a crop pest insect: Evidence for multiple polymorphisms of acetylcholinesterase gene with potential relevance for resistance to chemicals

The silverleaf whitefly Bemisia tabaci (Gennadius, 1889) (Homoptera: Aleyrodidae) is a serious invasive herbivorous insect pest worldwide. The excessive use of pesticides has progressively selected B. tabaci specimens, reducing the effectiveness of the treatments, and ultimately ending in the selection of pesticide-resistant strains. The management of this crop pest has thus become challenging owing to the level of resistance to all major classes of recommended insecticides. Here, we used in silico techniques for detecting sequence polymorphisms in ace1 gene from naturally occurring B. tabaci variants, and monitor the presence and frequency of the detected putative mutations from 30 populations of the silverleaf whitefly from Egypt and Pakistan. We found several point mutations in ace1-type acetylcholinesterase (ace1) in the studied B. tabaci variants naturally occurring in the field. By comparing ace1 sequence data from an organophosphate-susceptible and an organophosphate-resistant strains of B. tabaci to ace1 sequence data retrieved from GenBank for that species and to nucleotide polymorphisms from other arthropods, we identified novel mutations that could potentially influence insecticide resistance. Homology modeling and molecular docking analyses were performed to determine if the mutation-induced changes in form 1 acetylcholinesterase (AChE1) structure could confer resistance to carbamate and organophosphate insecticides. Mutations had small effects on binding energy (ΔG_b) interactions between mutant AChE1 and insecticides; they altered the conformation of the peripheral anionic site of AChE1, and modified the enzyme surface, and these changes have potential effects on the target-site sensitivity. Altogether, the results from this study provide information on genic variants of B. tabaci ace1 for future monitoring insecticide resistance development and report a potential case of environmentally driven gene variations.

Zaki M, M. Al-Qahtani S, Fayad E, T. Salem E, K. Abdulsahib W, Emam H, M. Hassan H. Integration of Ultrastructural and Computational Approaches Reveals the Protective Effect of Astaxanthin against BPA-Induced Nephrotoxicity

BACKGROUND

Bisphenol A (BPA) is an environmental contaminant that can induce deleterious organ effects. Human Cytochrome P450 CYP2C9 enzyme belongs to the essential xenobiotic-metabolizing enzymes, producing ROS as a byproduct. Astaxanthin (ATX) is a powerful antioxidant that protects organs and tissues from the damaging effects of oxidative stress caused by various diseases.

AIM OF THE STUDY

This study investigated the possible protective impacts of ATX against BPA-induced nephrotoxicity and its underlying mechanism.

MATERIALS AND METHODS

Kidney tissues were isolated and examined microscopically from control, protected, and unprotected groups of rats to examine the potential protective effect of ATX against nephrotoxicity. Moreover, a molecular dynamic (MD) simulation was conducted to predict the performance of ATX upon binding to the active site of P450 CYP2C9 protein receptor as a potential mechanism of ATX protective effect.

RESULTS

Implemented computational methods revealed the possible underlying mechanism of ATX protection; the protective impact of ATX is mediated by inhibiting P450 CYP2C9 through binding to its dimeric state where the RMSF value for apo-protein and ATX-complex system were 5.720.57 and 1.040.41, respectively, implicating the ATX-complex system to have lesser variance in its residues, leading to the prevention of ROS excess production, maintaining the oxidant-antioxidant balance and re-establishing the proper mitochondrial functionality. Furthermore, the experimental methods validated in silico outcomes and revealed that ATX therapy effectively restored the typical histological architecture of pathological kidney tissues.

CONCLUSIONS

ATX prevents BPA-induced nephrotoxicity by controlling oxidative imbalance and reversing mitochondrial dysfunction. These outcomes shed new light on the appropriate use of ATX as a treatment or prophylactic agent for these severe conditions.

Protein-Mineral Composite Particles with Logarithmic Dependence of Anticancer Cytotoxicity on Concentration of Montmorillonite Nanoplates with Adsorbed Cytochrome c

Montmorillonite (MM) colloid nanoplates have high adsorption capacity due to their large size/thickness ratio, which allows them to be used as carriers for drug delivery. Upon adsorption of the mitochondrial protein cytochrome c (cytC) onto MM plates, the composite cytC-MM particles acquire anticancer properties because of the ability of cancer cells to phagocytize submicron particles (in contrast to the normal cells). In this way, exogenous cytC can be introduced into tumor cells, thereby triggering apoptosis-an irreversible cascade of biochemical reactions leading to cell death. In the present study, we investigated the physicochemical properties of cytC-MM particles as a function of the cytC concentration in the suspension, namely, the electrophoretic mobility, the mass increment of MM monoplates upon cytC adsorption, the ratio of the adsorbed to the free cytC in the bulk, the protein density on the MM's surface, the number of cytC globules adsorbed on an MM monoplate, the concentration of cytC-MM composite particles in the suspension, and the dependence of cytotoxicity on the cytC-MM particle concentration. For this purpose, we used microelectrophoresis, static and electric light scattering, and a colon cancer cell culture to test the cytotoxic effects of the cytC-MM suspensions. The results show that the cytotoxicity depends linearly on the logarithm of the particle concentration in the cytC-MM suspension reaching 97%.

The Binding of Different Substrate Molecules at the Docking Site and the Active Site of γ-Secretase Can Trigger Toxic Events in Sporadic and Familial Alzheimer's Disease

Pathogenic changes in γ-secretase activity, along with its response to different drugs, can be affected by changes in the saturation of γ-secretase with its substrate. We analyze the saturation of γ-secretase with its substrate using multiscale molecular dynamics studies. We found that an increase in the saturation of γ-secretase with its substrate could result in the parallel binding of different substrate molecules at the docking site and the active site. The C-terminal domain of the substrate bound at the docking site can interact with the most dynamic presenilin sites at the cytosolic end of the active site tunnel. Such interactions can inhibit the ongoing catalytic activity and increase the production of the longer, more hydrophobic, and more toxic Aβ proteins. Similar disruptions in dynamic presenilin structures can be observed with different drugs and disease-causing mutations. Both, C99-βCTF-APP substrate and its different Aβ products, can support the toxic aggregation. The aggregation depends on the substrate N-terminal domain. Thus, the C99-βCTF-APP substrate and β-secretase path can be more toxic than the C83-αCTF-APP substrate and α-secretase path. Nicastrin can control the toxic aggregation in the closed conformation. The binding of the C99-βCTF-APP substrate to γ-secretase can be controlled by substrate channeling between the nicastrin and β-secretase. We conclude that the presented two-substrate mechanism could explain the pathogenic changes in γ-secretase activity and Aβ metabolism in different sporadic and familial cases of Alzheimer's disease. Future drug-development efforts should target different cellular mechanisms that regulate the optimal balance between γ-secretase activity and amyloid metabolism.

Computational studies on potential small molecule inhibitors of Leishmania pteridine reductase 1

Leishmaniasis is a neglected tropical disease of major public health concern. Challenges with current therapeutics have led to the exploration of plant medicine for potential antileishmanial agents. Despite the promising activity of some antileishmanial natural products, their protein targets have not been explored. The relevance of folate metabolism in the Leishmania parasite's existence presents crucial targets for the development of antileishmanial chemotherapy. Pteridine reductase 1 (PTR1), a crucial enzyme involved in DNA biosynthesis, is a validated target of the Leishmania parasite. Unearthing inhibitors of this enzyme is therefore an active research area. The goal of this work is to unearth small molecule inhibitors of PTR1 using molecular docking and molecular dynamic simulations. Thus, the interactions between selected antileishmanial natural products and PTR1 were examined. The binding affinities obtained from molecular docking ranged from -6.2 to -9.8 kcal/mol. When compared to the natural PTR1 substrate biopterin, compounds such as anonaine, chimanine D, corynantheine, grifolin, licochalcone A, piperogalin and xylopine produced better binding affinities, making interactions catalytic residues - Tyr194, Asp181, Phe113, Arg17 and Ser111. The PTR1- xylopine, -piperogalin, -grifolin, and -licochalcone A complexes exhibited remarkable stability under dynamic conditions during the entire 200 ns simulation period. The overall binding free energy of grifolin, piperogalin, and licochalcone A were observed to be -105.711, -103.567, and -105.646 kJ/mol respectively. The binding of these complexes was observed to be favorable and spontaneous and as such capable of inhibiting Leishmania PTR1. They could therefore be considered as candidates in the development of antileishmanial chemotherapy.Communicated by Ramaswamy H. Sarma.

New Insights into the Role of the Trypanosoma cruzi Aldo-Keto Reductase TcAKR

Chagas disease is a zoonotic infectious disease caused by the protozoan parasite Trypanosoma cruzi. It is distributed worldwide, affecting around 7 million people; there is no effective treatment, and it constitutes a leading cause of disability and premature death in the Americas. Only two drugs are currently approved for the treatment, Benznidazole and Nifurtimox, and both have to be activated by reducing the nitro-group. The T. cruzi aldo-keto reductase (TcAKR) has been related to the metabolism of benznidazole. TcAKR has been extensively studied, being most efforts focused on characterizing its implication in trypanocidal drug metabolism; however, little is known regarding its biological role. Here, we found that TcAKR is confined, throughout the entire life cycle, into the parasite mitochondria providing new insights into its biological function. In particular, in epimastigotes, TcAKR is associated with the kinetoplast, which suggests additional roles of the protein. The upregulation of TcAKR, which does not affect TcOYE expression, was correlated with an increase in PGF₂α, suggesting that this enzyme is related to PGF₂α synthesis in T. cruzi. Structural analysis showed that TcAKR contains a catalytic tetrad conserved in the AKR superfamily. Finally, we found that TcAKR is also involved in Nfx metabolization.

Bioinspired Polyethylene Glycol Coatings for Reduced Nanoparticle-Protein Interactions

Nanoparticles (NPs) and other engineered nanomaterials have great potential as nanodrugs or nanomedical devices for biomedical applications. However, the adsorption of proteins in blood circulation or similar physiological fluids can significantly alter the surface properties and therapeutic response induced by most nanomaterials. For example, interaction with proteins can change the bloodstream circulation time and availability of therapeutic NPs or hinder the accumulation in their desired target organs. Proteins can also trigger or prevent agglomeration. By combining experimental and computational approaches, we have developed NPs carrying polyethylene glycol (PEG) polymeric coatings that mimic the surface charge distribution of proteins typically found in blood, which are known to show low aggregation under normal blood conditions. Here, we show that NPs with coatings based on apoferritin or human serum albumin display better antifouling properties and weaker protein interaction compared to similar NPs carrying conventional PEG polymeric coatings.

Structural models of viral insulin-like peptides and their analogs

The insulin receptor (IR), the insulin-like growth factor-1 receptor (IGF1R), and the insulin/IGF1 hybrid receptors (hybR) are homologous transmembrane receptors. The peptide ligands, insulin and IGF1, exhibit significant structural homology and can bind to each receptor via site-1 and site-2 residues with distinct affinities. The variants of the Iridoviridae virus family show capability in expressing single-chain insulin/IGF1 like proteins, termed viral insulin-like peptides (VILPs), which can stimulate receptors from the insulin family. The sequences of VILPs lacking the central C-domain (dcVILPs) are known, but their structures in unbound and receptor-bound states have not been resolved to date. We report all-atom structural models of three dcVILPs (dcGIV, dcSGIV, and dcLCDV1) and their complexes with the receptors (μIR, μIGF1R, and μhybR), and probed the peptide/receptor interactions in each system using all-atom molecular dynamics (MD) simulations. Based on the nonbonded interaction energies computed between each residue of peptides (insulin and dcVILPs) and the receptors, we provide details on residues establishing significant interactions. The observed site-1 insulin/μIR interactions are consistent with previous experimental studies, and a residue-level comparison of interactions of peptides (insulin and dcVILPs) with the receptors revealed that, due to sequence differences, dcVILPs also establish some interactions distinct from those between insulin and IR. We also designed insulin analogs and report enhanced interactions between some analogs and the receptors.

Molecular dynamics simulations and kinetic measurements provide insights into the structural requirements of substrate size-dependent specificity of oligogalacturonide oxidase 1 (OGOX1)

Oligogalacturonides (OGs) are pectin fragments released from the breakdown of the homogalacturonan during pathogenesis that act as Damage-Associated Molecular Patterns. OG-oxidase 1 (OGOX1) is an Arabidopsis berberine bridge enzyme-like (BBE-l) oligosaccharide oxidase that oxidizes OGs, impairing their elicitor activity and concomitantly releasing H₂O₂. The OG-oxidizing activity of OGOX1 is markedly pH-dependent, with optimum pH around 10, and is higher towards OGs with a degree of polymerization higher than two. Here, the molecular determinants of OGOX1 responsible for the binding of OGs with different lengths have been investigated through molecular dynamics simulations and enzyme kinetics studies. OGOX1 was simulated in complex with OGs with different degree of polymerization such as di-, tri-, tetra- and penta-galacturonide, in water solution at alkaline pH. Our simulations revealed that, among the four OGOX1/OG combinations, the penta-galacturonide (OG5) showed the best conformation in the active site to be efficiently oxidized by OGOX1. The optimal conformation can be stabilized by salt-bridges formed between the carboxyl groups of OG5 and five positively charged amino acids of OGOX1, highly conserved in all OGOX paralogs. Our results suggest that these interactions limit the mobility of OG5 as well as longer OGs, contributing to maintain the terminal monomer of OGs in the optimal orientation in order to be oxidized by the enzyme. In accordance with these results, the enzyme efficiency (K_cat/K_M) of OGOX1 on OG5 (40.04) was found to be significantly higher than that on OG4 (13.05) and OG3 (0.6).

Identification of human phosphoglycerate mutase 1 (PGAM1) inhibitors using hybrid virtual screening approaches

PGAM1 plays a critical role in cancer cell metabolism through glycolysis and different biosynthesis pathways to promote cancer. It is generally known as a crucial target for treating pancreatic ductal adenocarcinoma, the deadliest known malignancy worldwide. In recent years different studies have been reported that strived to find inhibitory agents to target PGAM1, however, no validated inhibitor has been reported so far, and only a small number of different inhibitors have been reported with limited potency at the molecular level. Our in silico studies aimed to identify potential new PGAM1 inhibitors that could bind at the allosteric sites. At first, shape and feature-based models were generated and optimized by performing receiver operating characteristic (ROC) based enrichment studies. The best query model was then employed for performing shape, color, and electrostatics complementarity-based virtual screening of the ChemDiv database. The top two hundred and thirteen hits with greater than 1.2 TanimotoCombo score were selected and then subjected to structure-based molecular docking studies. The hits yielded better docking scores than reported compounds, were selected for subsequent structural similarity-based clustering analysis to select the best hits from each cluster. Molecular dynamics simulations and binding free energy calculations were performed to validate their plausible binding modes and their binding affinities with the PGAM1 enzyme. The results showed that these compounds were binding in the reported allosteric site of the enzyme and can serve as a good starting point to design better active selective scaffolds against PGAM1enzyme.

Formation and stabilization of C4a-hydroperoxy-FAD by the Arg/Asn pair in HadA monooxygenase

HadA monooxygenase catalyses the detoxification of halogenated phenols and nitrophenols via dehalogenation and denitration respectively. C4a-hydroperoxy-FAD is a key reactive intermediate wherein its formation, protonation and stabilization reflect enzyme efficiency. Herein, transient kinetics, site-directed mutagenesis and pH-dependent behaviours of HadA reaction were employed to identify key features stabilizing C4a-adducts in HadA. The formation of C4a-hydroperoxy-FAD is pH independent, whereas its decay and protonation of distal oxygen are associated with pK_a values of 8.5 and 8.4 respectively. These values are correlated with product formation within a pH range of 7.6-9.1, indicating the importance of adduct stabilization to enzymatic efficiency. We identified Arg101 as a key residue for reduced FAD (FADH^- ) binding and C4a-hydroperoxy-FAD formation due to the loss of these abilities as well as enzyme activity in HadA_R101A and HadA_R101Q . Mutations of the neighbouring Asn447 do not affect the rate of C4a-hydroperoxy-FAD formation; however, they impair FADH^- binding. The disruption of Arg101/Asn447 hydrogen bond networking in HadA_N447A increases the pK_a value of C4a-hydroperoxy-FAD decay to 9.5; however, this pK_a was not altered in HadA_N447D (pK_a of 8.5). Thus, Arg101/Asn447 pair should provide important interactions for FADH^- binding and maintain the pK_a associated with H₂ O₂ elimination from C4a-hydroperoxy-FAD in HadA. In the presence of substrate, the formation of C4a-hydroxy-FAD at the hydroxylation step is pH insensitive, and it dehydrates to form the oxidized FAD with pK_a of 7.9. This structural feature might help elucidate how the reactive intermediate was stabilized in other flavin-dependent monooxygenases.

Discovery of allosteric SHP2 inhibitors through ensemble-based consensus molecular docking, endpoint and absolute binding free energy calculations

SHP2 (Src homology-2 domain-containing protein tyrosine phosphatase-2) is a cytoplasmic protein -tyrosine phosphatase encoded by the gene PTPN11. It plays a crucial role in regulating cell growth and differentiation. Specifically, SHP2 is an oncoprotein associated with developmental pathologies and several different cancer types, including gastric, leukemia and breast cancer and is of great therapeutic interest. Given these roles, current research efforts have focused on developing SHP2 inhibitors. Allosteric SHP2 inhibitors have been shown to be more selective and pharmacologically appealing compared to competitive catalytic inhibitors targeting SHP2. Nevertheless, there remains a need for novel allosteric inhibitor scaffolds targeting SHP2 to develop compounds with improved selectivity, cell permeability, and bioavailability. Towards this goal, this study applied various computational tools to screen over 6 million compounds against the allosteric site within SHP2. The top-ranked hits from our in-silico screening were validated using protein thermal shift and biolayer interferometry assays, revealing three potent compounds. Kinetic binding assays were employed to measure the binding affinities of the top-ranked compounds and demonstrated that they all bind to SHP2 with a nanomolar affinity. Hence the compounds and the computational workflow described herein provide an effective approach for identifying and designing a generation of improved allosteric inhibitors of SHP2.

Probing intermolecular interactions and binding stability of antimicrobial peptides with beta-lactamase of Klebsiella aerogenes by comparing FDA approved beta-lactam drugs: a docking and molecular dynamics approach

Hospital pathogens, including Klebsiella aerogenes are becoming increasingly common, with the rise of Beta-lactam-resistant strains, especially in isolates recovered from intensive care rooms. Beta-lactamases participate in both the antibacterial activity and the mediation of the antibiotic resistance of Beta-lactams. The rapid spread of broad-spectrum Beta-lactam antibiotic resistance in pathogenic bacteria has recently become a major global health problem. As a result, new drugs that specifically target Beta-lactamases are urgently needed, and this enzyme has been identified to resolve the problem of bacterial resistance. In previous work, we de-novo developed, synthesized, and studied the in-vitro and in-silico behavior of four novel broad spectrum antimicrobial peptides, namely PEP01 to PEP04. All four peptides had significant antibacterial action against K. aerogenes. The literature evidence strongly suggests that Beta-lactamases are extremely important for bacteria, including K. aerogenes, and hence are therapeutically important and possible targets. Therefore, in this study we incorporated molecular modeling, docking, and simulation studies of the above four AMPs against the Beta-lactamase protein of K. aerogenes. The docking findings were also compared to eight FDA approved Beta-lactam antibiotics. According to our findings, all four peptides have strong binding affinity and interactions with Beta-lactamases and PEP02 has the highest docking score. In MD simulations, the protein-peptide complexes were more stable at 50 ns. We found that the new AMP-PEP02 is the most efficient and suitable drug candidate for inactivating Beta-lactamase protein, and that it is an alternative to or complements existing antibiotics for managing Beta-lactamase related resistance mechanisms based on this computational conclusion.Communicated by Ramaswamy H. Sarma.

Lipid Modulation of a Class B GPCR: Elucidating the Modulatory Role of PI(4,5)P2 Lipids

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂) lipids have been shown to stabilize an active conformation of class A G-protein coupled receptors (GPCRs) through a conserved binding site, not present in class B GPCRs. For class B GPCRs, previous molecular dynamics (MD) simulation studies have shown PI(4,5)P₂ interacting with the Glucagon receptor (GCGR), which constitutes an important target for diabetes and obesity therapeutics. In this work, we applied MD simulations supported by native mass spectrometry (nMS) to study lipid interactions with GCGR. We demonstrate how tail composition plays a role in modulating the binding of PI(4,5)P₂ lipids to GCGR. Specifically, we find the PI(4,5)P₂ lipids to have a higher affinity toward the inactive conformation of GCGR. Interestingly, we find that in contrast to class A GPCRs, PI(4,5)P₂ appear to stabilize the inactive conformation of GCGR through a binding site conserved across class B GPCRs but absent in class A GPCRs. This suggests differences in the regulatory function of PI(4,5)P₂ between class A and class B GPCRs.

Effect of Double Mutation (L452R and E484Q) on the Binding Affinity of Monoclonal Antibodies (mAbs) against the RBD-A Target for Vaccine Development

The COVID-19 pandemic, caused by SARS-CoV-2, emerges as a global health problem, as the viral genome is evolving rapidly to form several variants. Advancement and progress in the development of effective vaccines and neutralizing monoclonal antibodies are promising to combat viral infections. In the current scenario, several lineages containing "co-mutations" in the receptor-binding domain (RBD) region of the spike (S) protein are imposing new challenges. Co-occurrence of some co-mutations includes delta (L452R/T478K), kappa (L452R/E484Q), and a common mutation in both beta and gamma variants (E484K/N501Y). The effect of co-mutants (L452R/E484Q) on human angiotensin-converting enzyme 2 (hACE2) binding has already been elucidated. Here, for the first time, we investigated the role of these RBD co-mutations (L452R/E484Q) on the binding affinity of mAbs by adopting molecular dynamics (MD) simulation and free-energy binding estimation. The results obtained from our study suggest that the structural and dynamic changes introduced by these co-mutations reduce the binding affinity of the viral S protein to monoclonal antibodies (mAbs). The structural changes imposed by L452R create a charged patch near the interfacial surface that alters the affinity towards mAbs. In E484Q mutation, polar negatively charged E484 helps in the formation of electrostatic interaction, while the neutrally charged Q residue affects the interaction by forming repulsive forces. MD simulations along with molecular mechanics-generalized Born surface area (MMGBSA) studies revealed that the REGN 10933, BD-368-2, and S2M11 complexes have reduced binding affinity towards the double-mutant RBD. This indicates that their mutant (MT) structures have a stronger ability to escape from most antibodies than the wild type (WT). However, EY6A Ab showed higher affinity towards the double MT-RBD complex as compared to the WT. However, no significant effect of the per-residue contribution of double-mutated residues was observed, as this mAb does not interact with the region harboring L452 and E484 residues.

Computational insights for predicting the binding and selectivity of peptidomimetic plasmepsin IV inhibitors against cathepsin D

Plasmepsins (Plms) are aspartic proteases involved in the degradation of human hemoglobin by P. falciparum and are essential for the survival and growth of the parasite. Therefore, Plm enzymes are reported as an important antimalarial drug target. Herein, we have applied molecular docking, molecular dynamics (MD) simulations, and binding free energy with the Linear Interaction Energy (LIE) approach to investigate the binding of peptidomimetic PlmIV inhibitors with a particular focus on understanding their selectivity against the human Asp protease cathepsin D (CatD). The residual decomposition analysis results suggest that amino acid differences in the subsite S3 of PlmIV and CatD are responsible for the higher selectivity of the 5a inhibitor. These findings yield excellent agreement with experimental binding data and provide new details regarding van der Waals and electrostatic interactions of subsite residues as well as structural properties of the PlmIV and CatD systems.

ACPred-BMF: bidirectional LSTM with multiple feature representations for explainable anticancer peptide prediction

Cancer has become a major factor threatening human life and health. Under the circumstance that traditional treatment methods such as chemotherapy and radiotherapy are not highly specific and often cause severe side effects and toxicity, new treatment methods are urgently needed. Anticancer peptide drugs have low toxicity, stronger efficacy and specificity, and have emerged as a new type of cancer treatment drugs. However, experimental identification of anticancer peptides is time-consuming and expensive, and difficult to perform in a high-throughput manner. Computational identification of anticancer peptides can make up for the shortcomings of experimental identification. In this study, a deep learning-based predictor named ACPred-BMF is proposed for the prediction of anticancer peptides. This method uses the quantitative and qualitative properties of amino acids, binary profile feature to numerical representation for the peptide sequences. The Bidirectional LSTM network architecture is used in the model, and the attention mechanism is also considered. To alleviate the black-box problem of deep learning model prediction, we visualized the automatically extracted features and used the Shapley additive explanations algorithm to determine the importance of features to further understand the anticancer peptide mechanism. The results show that our method is one of the state-of-the-art anticancer peptide predictors. A web server as the implementation of ACPred-BMF that can be accessed via: http://mialab.ruc.edu.cn/ACPredBMFServer/ .

Proton transfer and conformational changes along the hydrogen bond network in heliorhodopsin

Heliorhodopsin releases a proton from the Schiff base during the L-state to M-state transition but not toward the protein bulk surface. Here we investigate proton transfer and induced structural changes along the H-bond network in heliorhodopsin using a quantum mechanical/molecular mechanical approach and molecular dynamics simulations. Light-induced proton transfer could occur from the Schiff base toward Glu107, reorienting Ser76, followed by subsequent proton transfer toward His80. His80 protonation induces the reorientation of Trp246 on the extracellular surface, originating from the electrostatic interaction that propagates along the transmembrane H-bond network [His80…His23…H₂O_[H23/Q26]…Gln26…Trp246] over a distance of 15 Å. Furthermore, it induces structural fluctuation on the intracellular side in the H-bond network [His80…Asn16…Tyr92…Glu230…Arg104…Glu149], opening the inner cavity at the Tyr92 moiety. These may be a basis of how light-induced proton transfer causes conformational changes during the M-state to O-state transition.

Computational Analysis of Triazole-Based Kojic Acid Analogs as Tyrosinase Inhibitors by Molecular Dynamics and Free Energy Calculations

Molecular docking, molecular dynamics (MD) simulations and the linear interaction energy (LIE) method were used here to predict binding modes and free energy for a set of 1,2,3-triazole-based KA analogs as potent inhibitors of Tyrosinase (TYR), a key metalloenzyme of the melanogenesis process. Initially, molecular docking calculations satisfactorily predicted the binding mode of evaluated KA analogs, where the KA part overlays the crystal conformation of the KA inhibitor into the catalytic site of TYR. The MD simulations were followed by the LIE method, which reproduced the experimental binding free energies for KA analogs with an r² equal to 0.97, suggesting the robustness of our theoretical model. Moreover, the van der Waals contributions performed by some residues such as Phe197, Pro201, Arg209, Met215 and Val218 are responsible for the binding recognition of 1,2,3-triazole-based KA analogs in TYR catalytic site. Finally, our calculations provide suitable validation of the combination of molecular docking, MD, and LIE approaches as a powerful tool in the structure-based drug design of new and potent TYR inhibitors.

1,4-Dideoxy-1,4-imino-D- and L-lyxitol-based inhibitors bind to Golgi α-mannosidase II in different protonation forms

The development of effective inhibitors of Golgi α-mannosidase II (GMII, E.C.3.2.1.114) with minimal off-target effects on phylogenetically-related lysosomal α-mannosidase (LMan, E.C.3.2.1.24) is a complex task due to the complicated structural and chemical properties of their active sites. The pKa values (and also protonation forms in some cases) of several ionizable amino acids, such as Asp, Glu, His or Arg of enzymes, can be changed upon the binding of the inhibitor. Moreover, GMII and LMan work under different pH conditions. The pKa calculations on large enzyme-inhibitor complexes and FMO-PIEDA energy decomposition analysis were performed on the structures of selected inhibitors obtained from docking and hybrid QM/MM calculations. Based on the calculations, the roles of the amino group incorporated in the ring of the imino-D-lyxitol inhibitors and some ionizable amino acids of Golgi-type (Asp270-Asp340-Asp341 of Drosophila melanogaster α-mannosidase dGMII) and lysosomal-type enzymes (His209-Asp267-Asp268 of Canavalia ensiformis α-mannosidase, JBMan) were explained in connection with the observed inhibitory properties. The pyrrolidine ring of the imino-D-lyxitols prefers at the active site of dGMII the neutral form while in JBMan the protonated form, whereas that of imino-L-lyxitols prefers the protonation form in both enzymes. The calculations indicate that the binding mechanism of inhibitors to the active-site of α-mannosidases is dependent on the inhibitor structure and could be used to design new selective inhibitors of GMII. A series of novel synthetic N-substituted imino-D-lyxitols were evaluated with four enzymes from the glycoside hydrolase GH38 family (two of Golgi-type, Drosophila melanogaster GMIIb and Caenorhabditis elegans AMAN-2, and two of lysosomal-type, Drosophila melanogaster LManII and Canavalia ensiformis JBMan, enzymes). The most potent structures [N-9-amidinononyl and N-2-(1-naphthyl)ethyl derivatives] inhibited GMIIb (Ki = 40 nM) and AMAN-2 (Ki = 150 nM) with a weak selectivity index (SI) toward Golgi-type enzymes of IC50(LManII)/IC50(GMIIb) = 35 or IC50(JBMan)/IC50(AMAN-2) = 86. On the other hand, weaker micromolar inhibitors, such as N-2-naphthylmethyl or 4-iodobenzyl derivatives [IC50(GMIIb) = 2.4 μM and IC50 (AMAN-2) = 7.6 μM], showed a significant SI in the range from 111 to 812.

Spontaneous and Electrically Induced Anisotropy of Composite Agarose Gels

Agarose gels containing and not bacteriorhodopsin purple membranes (incorporated before gelling) manifest spontaneous optical anisotropy. The dependencies of the anisotropy on the agarose concentration and time have been studied. The rise in the anisotropy is explained by the predominant orientation of the agarose fibers during the gelling and subsequent deformation of the gel net. In the electric field, additional optical anisotropy rises, which is caused by the orientation of the membranes. A procedure has been developed to separate electrically induced and spontaneous anisotropy in composite gels. The isoelectric points and surface electric potential of bacteriorhodopsin trimer and purple membranes are calculated by the method of protein electrostatics to explain their electric asymmetry, which leads to perpendicular orientation in the direct electric field and longitudinal in the kilohertz sinusoidal field. The results allow for an increase in the separation capability of composite gels of electrophoresis for macromolecules with different sizes by applying an appropriate electric field to modulate the effective pore size.

Crowding modulates the glycation of plasma proteins: In vitro analysis of structural modifications to albumin and transferrin and identification of sites of modification

Protein modification occurs in biological milieus that are characterized by high concentrations of (macro)molecules (i.e. heterogeneous and packed environments). Recent data indicate that crowding can modulate the extent and rate of protein oxidation, however its effect on other post-translational modifications remains to be explored. In this work we hypothesized that crowding would affect the glycation of plasma proteins. Physiologically-relevant concentrations of albumin (35 mg mL^-1) and transferrin (2 mg mL^-1) were incubated with methylglyoxal and glyoxal (5 μM-5 mM), two α-oxoaldehyde metabolites that are elevated in the plasma of people with diabetes. Crowding was induced by adding dextran or ficoll polymers. Electrophoresis, electron microscopy, fluorescence spectroscopy and mass spectrometry were employed to investigate the structural consequences of glycation under crowded conditions. Our data demonstrate that crowding modulates the extent of formation of transferrin cross-links, and also the modification pathways in both albumin and transferrin. Arginine was the most susceptible residue to modification, with lysine and cysteine also affected. Loss of 0.48 and 7.28 arginine residues per protein molecule were determined on incubation with 500 μM methylglyoxal for albumin and transferrin, respectively. Crowding did not influence the extent of loss of arginine and lysine for either protein, but the sites of modification, detected by LC-MS, were different between dilute and crowded conditions. These data confirm the relevance of studying modification processes under conditions that closely mimic biological milieus. These data unveil additional factors that influence the pattern and extent of protein modification, and their structural consequences, in biological systems.

Endowing homodimeric carbamoyltransferase GdmN with iterative functions through structural characterization and mechanistic studies

Iterative enzymes, which catalyze sequential reactions, have the potential to improve the atom economy and diversity of industrial enzymatic processes. Redesigning one-step enzymes to be iterative biocatalysts could further enhance these processes. Carbamoyltransferases (CTases) catalyze carbamoylation, an important modification for the bioactivity of many secondary metabolites with pharmaceutical applications. To generate an iterative CTase, we determine the X-ray structure of GdmN, a one-step CTase involved in ansamycin biosynthesis. GdmN forms a face-to-face homodimer through unusual C-terminal domains, a previously unknown functional form for CTases. Structural determination of GdmN complexed with multiple intermediates elucidates the carbamoylation process and identifies key binding residues within a spacious substrate-binding pocket. Further structural and computational analyses enable multi-site enzyme engineering, resulting in an iterative CTase with the capacity for successive 7-O and 3-O carbamoylations. Our findings reveal a subclade of the CTase family and exemplify the potential of protein engineering for generating iterative enzymes.

Mechanism of the distinct toxicity level of imidacloprid and thiacloprid against honey bees: An in silico study based on cytochrome P450 9Q3

The honey bee, Apis mellifera, shows variation in sensitivity to imidacloprid and thiacloprid, which does not reside at the target site but rather in the rapidly oxidative metabolism mediated by P450s (such as a single P450, CYP9Q3). An in silico study was conducted to investigate the various metabolism of imidacloprid and thiacloprid. The binding potency of thiacloprid was stronger and a stable π-π interaction with Phe121 and the N-H⋯N hydrogen bond with Asn214 are found in the CYP9Q3-thiacloprid system but absent in imidacloprid, which might affect the potential metabolic activity. Moreover, the values of highest occupied molecular orbit (HOMO) energy and the vertical ionization potential (IP) of two compounds demonstrated that thiacloprid is more likely to oxidation. The findings revealed the probable binding modes of imidacloprid and thiacloprid with CYP9Q3 and might facilitate future design of the low bee toxicity neonicotinoid insecticides.

Targeting a Multidrug-Resistant Pathogen: First Generation Antagonists of Burkholderia cenocepacia's BC2L-C Lectin

Multidrug-resistant pathogens such as Burkholderia cenocepacia have become a hazard in the context of healthcare-associated infections, especially for patients admitted with cystic fibrosis or immuno-compromising conditions. Like other opportunistic Gram-negative bacteria, this pathogen establishes virulence and biofilms through lectin-mediated adhesion. In particular, the superlectin BC2L-C is believed to cross-link human epithelial cells to B. cenocepacia during pulmonary infections. We aimed to obtain glycomimetic antagonists able to inhibit the interaction between the N-terminal domain of BC2L-C (BC2L-C-Nt) and its target fucosylated human oligosaccharides. In a previous study, we identified by fragment virtual screening and validated a small set of molecular fragments that bind BC2L-C-Nt in the vicinity of the fucose binding site. Here, we report the rational design and synthesis of bifunctional C- or N-fucosides, generated by connecting these fragments to a fucoside core using a panel of rationally selected linkers. A modular route starting from two key fucoside intermediates was implemented for the synthesis, followed by evaluation of the new compounds as BC2L-C-Nt ligands with a range of techniques (surface plasmon resonance, isothermal titration calorimetry, saturation transfer difference NMR, differential scanning calorimetry, and X-ray crystallography). This study resulted in a hit molecule with an order of magnitude gain over the starting methyl fucoside and in two crystal structures of antagonist/lectin complexes.

Antibody apparent solubility prediction from sequence by transfer learning

Developing therapeutic monoclonal antibodies (mAbs) for the subcutaneous administration requires identifying mAbs with superior solubility that are amenable for high-concentration formulation. However, experimental screening is often material and labor intensive. Here, we present a strategy (named solPredict) that employs the embeddings from pretrained protein language modeling to predict the apparent solubility of mAbs in histidine (pH 6.0) buffer. A dataset of 220 diverse, in-house mAbs were used for model training and hyperparameter tuning through 5-fold cross validation. solPredict achieves high correlation with experimental solubility on an independent test set of 40 mAbs. Importantly, solPredict performs well for both IgG1 and IgG4 subclasses despite the distinct solubility behaviors. This approach eliminates the need of 3D structure modeling of mAbs, descriptor computation, and expert-crafted input features. The minimal computational expense of solPredict enables rapid, large-scale, and high-throughput screening of mAbs using sequence information alone during early antibody discovery.

Novel Calcium-Binding Ablating Mutations Induce Constitutive RET Activity and Drive Tumorigenesis

Distinguishing oncogenic mutations from variants of unknown significance (VUS) is critical for precision cancer medicine. Here, computational modeling of 71,756 RET variants for positive selection together with functional assays of 110 representative variants identified a three-dimensional cluster of VUSs carried by multiple human cancers that cause amino acid substitutions in the calmodulin-like motif (CaLM) of RET. Molecular dynamics simulations indicated that CaLM mutations decrease interactions between Ca2+ and its surrounding residues and induce conformational distortion of the RET cysteine-rich domain containing the CaLM. RET-CaLM mutations caused ligand-independent constitutive activation of RET kinase by homodimerization mediated by illegitimate disulfide bond formation. RET-CaLM mutants possessed oncogenic and tumorigenic activities that could be suppressed by tyrosine kinase inhibitors targeting RET. This study identifies calcium-binding ablating mutations as a novel type of oncogenic mutation of RET and indicates that in silico-driven annotation of VUSs of druggable oncogenes is a promising strategy to identify targetable driver mutations.

SIGNIFICANCE

Comprehensive proteogenomic and in silico analyses of a vast number of VUSs identify a novel set of oncogenic and druggable mutations in the well-characterized RET oncogene.

SARS-CoV-2 Spike Protein (RBD) Subunit Adsorption at Abiotic Surfaces and Corona Formation at Polymer Particles

The adsorption kinetics of the SARS-CoV-2 spike protein subunit with the receptor binding domain at abiotic surfaces was investigated. A combination of sensitive methods was used such as atomic force microscopy yielding a molecular resolution, a quartz microbalance, and optical waveguide lightmode spectroscopy. The two latter methods yielded in situ information about the protein adsorption kinetics under flow conditions. It was established that at pH 3.5-4 the protein adsorbed on mica and silica surfaces in the form of compact quasi-spherical aggregates with an average size of 14 nm. The maximum coverage of the layers was equal to 3 and 1 mg m^-2 at pH 4 and 7.4, respectively. The experimental data were successfully interpreted in terms of theoretical results derived from modeling. The experiments performed for flat substrates were complemented by investigations of the protein corona formation at polymer particles carried out using in situ laser Doppler velocimetry technique. In this way, the zeta potential of the protein layers was acquired as a function of the coverage. Applying the electrokinetic model, these primary data were converted to the dependence of the subunit zeta potential on pH. It was shown that a complete acid-base characteristic of the layer can be acquired only using nanomolar quantities of the protein.

AMOEBA Force Field Trajectories Improve Predictions of Accurate pKa Values of the GFP Fluorophore: The Importance of Polarizability and Water Interactions

Precisely quantifying the magnitude, direction, and biological functions of electric fields in proteins has long been an outstanding challenge in the field. The most widely implemented experimental method to measure such electric fields at a particular residue in a protein has been through changes in pKa of titratable residues. While many computational strategies exist to predict these values, it has been difficult to do this accurately or connect predicted results to key structural or mechanistic features of the molecule. Here, we used experimentally determined pKa values of the fluorophore in superfolder green fluorescent protein (GFP) with amino acid mutations made at position Thr 203 to evaluate the pKa prediction ability of molecular dynamics (MD) simulations using a polarizable force field, AMOEBA. Structure ensembles from AMOEBA were used to calculate pKa values of the GFP fluorophore. The calculated pKa values were then compared to trajectories using a conventional fixed charge force field (Amber03 ff). We found that the position of water molecules included in the pKa calculation had opposite effects on the pKa values between the trajectories from AMOEBA and Amber03 force fields. In AMOEBA trajectories, the inclusion of water molecules within 35 Å of the fluorophore decreased the difference between the predicted and experimental values, resulting in calculated pKa values that were within an average of 0.8 pKa unit from the experimental results. On the other hand, in Amber03 trajectories, including water molecules that were more than 5 Å from the fluorophore increased the differences between the calculated and experimental pKa values. The inaccuracy of pKa predictions determined from Amber03 trajectories was caused by a significant stabilization of the deprotonated chromophore's free energy compared to the result in AMOEBA. We rationalize the cutoffs for explicit water molecules when calculating pKa to better predict the electrostatic environment surrounding the fluorophore buried in GFP. We discuss how the results from this work will assist the prospective prediction of pKa values or other electrostatic effects in a wide variety of folded proteins.

Predicting Experimental B22 Values and the Effects of Histidine Charge States for Monoclonal Antibodies Using Coarse-Grained Molecular Simulations

Static light scattering (SLS) was used to characterize five monoclonal antibodies (MAbs) as a function of total ionic strength (TIS) at pH values between 5.5 and 7.0. Second osmotic virial coefficient (B₂₂) values were determined experimentally for each MAb as a function of TIS using low protein concentration SLS data. Coarse-grained molecular simulations were performed to predict the B₂₂ values for each MAb at a given pH and TIS. To include the effect of charge fluctuations of titratable residues in the B₂₂ calculations, a statistical approach was introduced in the Monte Carlo algorithm based on the protonation probability based on a given pH value and the Henderson-Hasselbalch equation. The charged residues were allowed to fluctuate individually, based on the sampled microstates and the influence of electrostatic interactions on net protein-protein interactions during the simulations. Compared to static charge simulations, the new approach provided improved results compared to experimental B₂₂ values at pH conditions near the pK_a of titratable residues.

Decoding drug resistant mechanism of V32I, I50V and I84V mutations of HIV-1 protease on amprenavir binding by using molecular dynamics simulations and MM-GBSA calculations

Mutations V32I, I50V and I84V in the HIV-1 protease (PR) induce drug resistance towards drug amprenavir (APV). Multiple short molecular dynamics (MSMD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) method were utilized to investigate drug-resistant mechanism of V32I, I50V and I84V towards APV. Dynamic information arising from MSMD simulations suggest that V32I, I50V and I84V highly affect structural flexibility, motion modes and conformational behaviours of two flaps in the PR. Binding free energies calculated by MM-GBSA method suggest that the decrease in binding enthalpy and the increase in binding entropy induced by mutations V32I, I50V and I84V are responsible for drug resistance of the mutated PRs on APV. The energetic contributions of separate residues on binding of APV to the PR show that V32I, I50V and I84V highly disturb the interactions of two flaps with APV and mostly drive the decrease in binding ability of APV to the PR. Thus, the conformational changes of two flaps in the PR caused by V32I, I50V and I84V play key roles in drug resistance of three mutated PR towards APV. This study can provide useful dynamics information for the design of potent inhibitors relieving drug resistance.

Determination of specific and non-specific protein-protein interactions for beta-lactoglobulin by analytical ultracentrifugation and membrane osmometry experiments

Protein-protein interactions are essential for the understanding of biological processes. Specific protein aggregation is an important aspect for many biological systems. In particular, electrostatic interactions play the key role for protein-protein interactions, as many amino acids have pH-dependent charge states. Moreover, protein dissociation is directly related to the solution pH, ionic strength, temperature and protein concentration. The subtle interplay between different specific and non-specific interactions is demonstrated for beta-lactoglobulin (BLG) with a focus on low salt concentrations, thus mimicking technically relevant processing conditions. BLG is a well-characterized model system, proven to attain its monomer-dimer equilibrium strongly dependent upon the pH of the solution. In this manuscript, we present a unique combination of analytical ultracentrifugation and membrane osmometry experiments, which quantifies specific and non-specific interactions, i.e. in terms of the dimer dissociation constants and the second osmotic virial coefficient, at pH 3 and 7 and sodium chloride concentrations of 10 mM and 100 mM. This provides direct insight to protein-protein interactions for a system with a concentration-dependent monomer-dimer equilibrium. Moreover, using a coarse-grained extended DLVO model in combination with molecular dynamics simulations, we quantify non-specific monomer-monomer, monomer-dimer and dimer-dimer interactions as well as the binding free energy of BLG dimerization from theoretical calculations. The experimentally determined interactions are shown to be mainly governed by electrostatic interactions and further agree with free energy calculations. Our experimental protocol aims to determine non-specific and specific interactions for a dynamically interacting system and provides an understanding of protein-protein interactions for BLG at low salt concentrations.

2-(Fluoromethoxy)-4'-(S-methanesulfonimidoyl)-1,1'-biphenyl (UCM-1306), an Orally Bioavailable Positive Allosteric Modulator of the Human Dopamine D1 Receptor for Parkinson's Disease

Tolerance development caused by dopamine replacement with l-DOPA and therapeutic drawbacks upon activation of dopaminergic receptors with orthosteric agonists reveal a significant unmet need for safe and effective treatment of Parkinson’s disease. In search for selective modulators of the D₁ receptor, the screening of a chemical library and subsequent medicinal chemistry program around an identified hit resulted in new synthetic compound 26 [UCM-1306, 2-(fluoromethoxy)-4′-(S-methanesulfonimidoyl)-1,1′-biphenyl] that increases the dopamine maximal effect in a dose-dependent manner in human and mouse D₁ receptors, is inactive in the absence of dopamine, modulates dopamine affinity for the receptor, exhibits subtype selectivity, and displays low binding competition with orthosteric ligands. The new allosteric modulator potentiates cocaine-induced locomotion and enhances l-DOPA recovery of decreased locomotor activity in reserpinized mice after oral administration. The behavior of compound 26 supports the interest of a positive allosteric modulator of the D₁ receptor as a promising therapeutic approach for Parkinson’s disease.