PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD Simulations

Integrating Dynamic Network Analysis with AI for Enhanced Epitope Prediction in PD-L1:Affibody Interactions

Understanding binding epitopes involved in protein-protein interactions and accurately determining their structure are long-standing goals with broad applicability in industry and biomedicine. Although various experimental methods for binding epitope determination exist, these approaches are typically low throughput and cost-intensive. Computational methods have potential to accelerate epitope predictions; however, recently developed artificial intelligence (AI)-based methods frequently fail to predict epitopes of synthetic binding domains with few natural homologues. Here we have developed an integrated method employing generalized-correlation-based dynamic network analysis on multiple molecular dynamics (MD) trajectories, initiated from AlphaFold2Multimer structures, to unravel the structure and binding epitope of the therapeutic PD-L1:Affibody complex. Both AlphaFold2 and conventional molecular dynamics trajectory analysis were ineffective in distinguishing between two proposed binding models, parallel and perpendicular. However, our integrated approach, utilizing dynamic network analysis, demonstrated that the perpendicular mode was significantly more stable. These predictions were validated using a suite of experimental epitope mapping protocols, including cross-linking mass spectrometry and next-generation sequencing-based deep mutational scanning. Conversely, AlphaFold3 failed to predict a structure bound in the perpendicular pose, highlighting the necessity for exploratory research in the search for binding epitopes and challenging the notion that AI-generated protein structures can be accepted without scrutiny. Our research underscores the potential of employing dynamic network analysis to enhance AI-based structure predictions for more accurate identification of protein-protein interaction interfaces.

Simulation-driven design of stabilized SARS-CoV-2 spike S2 immunogens

The full-length prefusion-stabilized SARS-CoV-2 spike (S) is the principal antigen of COVID-19 vaccines. Vaccine efficacy has been impacted by emerging variants of concern that accumulate most of the sequence modifications in the immunodominant S1 subunit. S2, in contrast, is the most evolutionarily conserved region of the spike and can elicit broadly neutralizing and protective antibodies. Yet, S2's usage as an alternative vaccine strategy is hampered by its general instability. Here, we use a simulation-driven approach to design S2-only immunogens stabilized in a closed prefusion conformation. Molecular simulations provide a mechanistic characterization of the S2 trimer's opening, informing the design of tryptophan substitutions that impart kinetic and thermodynamic stabilization. Structural characterization via cryo-EM shows the molecular basis of S2 stabilization in the closed prefusion conformation. Informed by molecular simulations and corroborated by experiments, we report an engineered S2 immunogen that exhibits increased protein expression, superior thermostability, and preserved immunogenicity against sarbecoviruses.

Enhanced sensitivity in electrochemical detection of ochratoxin A within food samples using ferrocene- and aptamer-tethered gold nanoparticles on disposable electrodes

Ensuring food security is crucial for public health, and the presence of mycotoxins, produced by fungi in improperly stored processed or unprocessed food, poses a significant threat. This research introduces a novel approach - a disposable aptasensing platform designed for the detection of ochratoxin A (OTA). The platform employs gold-nanostructured screen-printed carbon electrodes functionalized with a ferrocene derivative, serving as an integrated faradaic transducing system, and an anti-OTA aptamer as a bioreceptor site. Detection relies on the ferrocene electrochemical signal changes induced by the aptamer folding in the presence of the target molecule. Remarkably sensitive, the platform detects OTA within the range of 0.5 to 70 ng mL-1 and a detection limit of 11 pg mL-1. This limit is approximately 200 times below the levels stipulated by the European Commission for agricultural commodities. Notably, the sensing device exhibits efficacy in detecting OTA in complex media, such as roasted coffee beans and wine, without the need for sample pretreatment, yielding accurate recoveries. Furthermore, while label-free electrochemical aptasensors have proliferated, this study addresses a gap in understanding the binding mechanisms of some aptasensors. To enhance the experimental findings, a theoretical study was conducted to underscore the specificity of the anti-OTA aptamer as a donor for OTA detection. The molecular docking technique was employed to unveil the key binding region of the aptamer, providing valuable insights into the aptasensor specificity.

Comparative analysis of structural dynamics and allosteric mechanisms of RecA/Rad51 family proteins: Integrated atomistic MD simulation and network-based analysis

Homologous recombination plays a key role in double-strand break repair, stalled replication fork repair, and meiosis. The RecA/Rad51 family recombinases catalyze the DNA strand invasion reaction that occurs during homologous recombination. However, the high sequence differences between homologous groups have hindered the thoroughly studies of this ancient protein family. The dynamic mechanisms of the family, particularly at the residual level, remain poorly understood. In this work, five representative RecA/Rad51 recombinase family members from all major kingdoms of living organisms: prokaryotes, eukaryotes, archaea, and viruses, were selected to explore the molecular mechanisms behind their conserved biological significance. A variety of techniques, including all-atom molecular dynamics simulation, perturbation response scanning, and protein structure network analysis, were used to examine the flexibility and correlation of protein domains, distribution of sensors and effectors and conserved hub residues. Furthermore, the potential communication routes between the ATP-binding region and the DNA-binding region of each recombinase were identified. Our results demonstrate the conserved molecular dynamics of these recombinases in the early stage of homologous recombination, including cooperative motions between regions, conserved sensing and effecting functional residue distribution, and conserved hub residues. Meanwhile, the unique ATP-DNA communication routes of each recombinase was also revealed. These results provide new insights into the mechanism of RecA/Rad51 family proteins, and provide new theoretical guidance for the development of allosteric inhibitors and the application of RecA/Rad51 family proteins.

Key interaction networks: Identifying evolutionarily conserved non-covalent interaction networks across protein families

Protein structure (and thus function) is dictated by non-covalent interaction networks. These can be highly evolutionarily conserved across protein families, the members of which can diverge in sequence and evolutionary history. Here we present KIN, a tool to identify and analyze conserved non-covalent interaction networks across evolutionarily related groups of proteins. KIN is available for download under a GNU General Public License, version 2, from https://www.github.com/kamerlinlab/KIN. KIN can operate on experimentally determined structures, predicted structures, or molecular dynamics trajectories, providing insight into both conserved and missing interactions across evolutionarily related proteins. This provides useful insight both into protein evolution, as well as a tool that can be exploited for protein engineering efforts. As a showcase system, we demonstrate applications of this tool to understanding the evolutionary-relevant conserved interaction networks across the class A β-lactamases.

Binding order of substrate and cofactor in sulfonamide monooxygenase during sulfa drug degradation: in silico studies

For decades, sulfonamide antibiotics have been used across industries such as agriculture and animal husbandry. However, the use and inadvertent misuse of these antibiotics have resulted in the advent of sulfonamide-drug-resistant strains due to antibiotic pollution. Enzymatic bioremediation of antibiotics remains a potential emerging solution to combat antibiotic pollution. Here, we propose an enzymatic model for the degradation of sulfonamides by Microbacterium sp. We have employed a multi-pronged computational strategy involving - protein structure modelling, ligand docking and molecular dynamics simulations to decipher a plausible binding order for the enzymatic degradation of sulfonamides by the bacterial sulfonamide monooxygenase, SulX. Our results enable us to predict that this degradation is achieved through the sequential binding of the antibiotic sulfonamide followed by the reduced flavin cofactor FMNH2, thereby laying the computational foundation for further advancements in enzyme-mediated degradation of the antibiotic. We also provide a list of experiments which may be performed to verify and follow-up on our in-silico studies.Communicated by Ramaswamy H. Sarma.

Discovery of a new dihydroeugenol-chalcone hybrid with cytotoxic and anti-migratory potential: A dual-action hit for cancer therapeutics

Cancer still represents a serious public health problem and one of the main problems related to the worsening of this disease is the ability of some tumors to develop metastasis. In this work, we synthesized a new series of chalcones and isoxazoles derived from eugenol and analogues as molecular hybrids and these compounds were evaluated against different tumor cell lines. This structural pattern was designed considering the cytotoxic potential already known for eugenol, chalcones and isoxazoles. Notably, chalcones 7, 9, 10, and 11 displayed significant activity (4.2-14.5 µM) against two cancer cell lines, surpassing the potency of the control drug doxorubicin. The reaction of chalcones with hydroxylamine hydrochloride provided the corresponding isoxazoles that were inactive against these cancer cells. The dihydroeugenol chalcone 7 showed the most promising results, demonstrating higher potency against HepG2 (CC50: 4.2 µM) and TOV-21G (CC50: 7.2 µM). Chalcone 7 was also three times less toxic than doxorubicin considering HepG2 cells, with a selectivity index greater than 11. Further investigations including clonogenic survival, cell cycle progression and cell migration assays confirmed the compelling antitumoral potential of chalcone 7, as it reduced long-term survival due to DNA fragmentation, inducing cell death and inhibiting HepG2 cells migration. Moreover, in silico studies involving docking and molecular dynamics revealed a consistent binding mode of chalcone 7 with metalloproteinases, particularly MMP-9, shedding light on its potential mechanism of action related to anti-migratory effects. These significant findings suggest the inclusion of compound 7 as a promising candidate for future studies in the field of cancer therapeutics.

Structural and thermodynamic insights into antibody light chain tetramer formation through 3D domain swapping

Overexpression of antibody light chains in small plasma cell clones can lead to misfolding and aggregation. On the other hand, the formation of amyloid fibrils from antibody light chains is related to amyloidosis. Although aggregation of antibody light chain is an important issue, atomic-level structural examinations of antibody light chain aggregates are sparse. In this study, we present an antibody light chain that maintains an equilibrium between its monomeric and tetrameric states. According to data from X-ray crystallography, thermodynamic and kinetic measurements, as well as theoretical studies, this antibody light chain engages in 3D domain swapping within its variable region. Here, a pair of domain-swapped dimers creates a tetramer through hydrophobic interactions, facilitating the revelation of the domain-swapped structure. The negative cotton effect linked to the β-sheet structure, observed around 215 nm in the circular dichroism (CD) spectrum of the tetrameric variable region, is more pronounced than that of the monomer. This suggests that the monomer contains less β-sheet structures and exhibits greater flexibility than the tetramer in solution. These findings not only clarify the domain-swapped structure of the antibody light chain but also contribute to controlling antibody quality and advancing the development of future molecular recognition agents and drugs.

A cystatin C similar protein from Musa acuminata that inhibits cathepsin B involved in rheumatoid arthritis using in silico approach and in vitro cathepsin B inhibition by protein extract

Rheumatoid arthritis (RA) is an auto-immune disease that affects the synovial lining of the joints, causes synovitis and culminates to joint destruction. Cathepsin B is responsible for digesting unwanted proteins in extracellular matrix but its hyper expression could implicate in pathological diseases like RA. Available treatments for RA are classified into non-steroidal anti-inflammatory drugs (NSAIDs), disease-modifying anti-rheumatic drugs (DMARDs), and steroids, but the severe side effects associated with these drugs is one of concerns and cannot be ignored. Thus, any alternative therapy with minimum or no side effects would be a cornerstone. In our in silico studies a cystatin C similar protein (CCSP) has been identified from Musa acuminata that could effectively inhibit the cathepsin B activity. In silico and molecular dynamics studies showed that the identified CCSP and cathepsin B complex has binding energy -66.89 kcal/mol as compared to cystatin C - cathepsin B complex with binding energy of -23.38 kcal/mol. These results indicate that CCSP from Musa acuminata has better affinity towards cathepsin B as compared to its natural inhibitor cystatin C. Hence, CCSP may be suggested as an alternative therapeutic in combating RA by inhibiting its one of the key proteases cathepsin B. Further, in vitro experiments with fractionated protein extracts from Musa sp. peel inhibited cathepsin B to 98.30% at 300 µg protein concentration and its IC50 was found to be 45.92 µg indicating the presence of cathepsin B inhibitor(s) in protein extract of peel which was further confirmed by reverse zymography.Communicated by Ramaswamy H. Sarma.

Molecular dynamics of JUNO-IZUMO1 complexation suggests biologically relevant mechanisms in fertilization

JUNO-IZUMO1 binding is the first known physical link created between the sperm and egg membranes in fertilization, however, how this initiates sperm-egg fusion remains elusive. As advanced structural insights will help to combat the infertility crisis, or advance fertility control, we employed all-atom Molecular Dynamics (MD) to derive dynamic structural insights that are difficult to obtain experimentally. We found that the hydrated JUNO-IZUMO1 interface is composed of a large set of short-lived non-covalent interactions. The contact interface is destabilized by strategically located point mutations, as well as by Zn2+ ions, which shift IZUMO1 into the non-binding "boomerang" conformation. We hypothesize that the latter might explain how the transient zinc spark, as released after sperm entry into the oocyte, might contribute to block polyspermy. To address a second mystery, we performed another set of simulations, as it was previously suggested that JUNO in solution is unable to bind to folate despite it belonging to the folate receptor family. MD now suggests that JUNO complexation with IZUMO1 opens up the binding pocket thereby enabling folate insertion. Our MD simulations thus provide crucial new hypotheses how the dynamics of the JUNO-IZUMO1 complex upon solvation might regulate fertility.

Molecular interaction analysis of anti-IL-8 scFv-10F8-6His against IL-8 monomer through molecular docking and molecular dynamic simulations

Elevated interleukin 8 (IL-8) expression has been linked to unfavorable outcomes in a range of inflammatory conditions, such as rheumatoid arthritis, psoriasis, and cancer. The human monoclonal antibody (HuMab) 10F8 and the hybridoma 35B11-B bind to an epitope on human IL-8, respectively. 10F8 inhibited interaction between IL-8 and neutrophils in eczema and pustulosis palmoplantaris patients while 35B11-B decreased size lesion in rat model. The binding interaction of monoclonal antibodies and IL-8, especially how complementarity-determining region (CDR) loops could bind the N-terminal of IL-8, has not been fully deliberated at molecular-level. Here, we used a combination of molecular docking, heated and long coarse-grained molecular dynamics simulations to identify key residues of established interaction. Based on heated MD simulation, docked pose of complexes generated by ClusPro showed good binding stability throughout of 70 ns simulation. Based on long molecular dynamic simulations, key residues for the binding were identified throughout of 1000 ns simulation. TYR-53, ASP-99, and ARG-100 of heavy chain CDR together with TYR-33 of light chain CDR are among the highest contributing energy residues within the binding interaction. Meanwhile, LYS11 and TYR13 of IL-8 are important for the determination of overall binding energy. Furthermore, the result of decomposition residues analysis is in good agreement with the interaction analysis data. Current study provides a list of important interacting residues and further scrutiny on these residues is essential for future development and design of a new and stable recombinant antibody against IL-8.Communicated by Ramaswamy H. Sarma.

Insights into substrate coordination and glycosyl transfer of poplar cellulose synthase-8

Cellulose is an abundant cell wall component of land plants. It is synthesized from UDP-activated glucose molecules by cellulose synthase, a membrane-integrated processive glycosyltransferase. Cellulose synthase couples the elongation of the cellulose polymer with its translocation across the plasma membrane. Here, we present substrate- and product-bound cryogenic electron microscopy structures of the homotrimeric cellulose synthase isoform-8 (CesA8) from hybrid aspen (poplar). UDP-glucose binds to a conserved catalytic pocket adjacent to the entrance to a transmembrane channel. The substrate's glucosyl unit is coordinated by conserved residues of the glycosyltransferase domain and amphipathic interface helices. Site-directed mutagenesis of a conserved gating loop capping the active site reveals its critical function for catalytic activity. Molecular dynamics simulations reveal prolonged interactions of the gating loop with the substrate molecule, particularly across its central conserved region. These transient interactions likely facilitate the proper positioning of the substrate molecule for glycosyl transfer and cellulose translocation.

Molecular dynamics simulations support a multistep pathway for activation of branched actin filament nucleation by Arp2/3 complex

Actin-related protein 2/3 complex (Arp2/3 complex) catalyzes the nucleation of branched actin filaments that push against membranes in processes like cellular motility and endocytosis. During activation by WASP proteins, the complex must bind WASP and engage the side of a pre-existing (mother) filament before a branched filament is nucleated. Recent high-resolution structures of activated Arp2/3 complex revealed two major sets of activating conformational changes. How these activating conformational changes are triggered by interactions of Arp2/3 complex with actin filaments and WASP remains unclear. Here we use a recent high-resolution structure of Arp2/3 complex at a branch junction to design all-atom molecular dynamics simulations that elucidate the pathway between the active and inactive states. We ran a total of ∼4.6 microseconds of both unbiased and steered all-atom molecular dynamics simulations starting from three different binding states, including Arp2/3 complex within a branch junction, bound only to a mother filament, and alone in solution. These simulations indicate that the contacts with the mother filament are mostly insensitive to the massive rigid body motion that moves Arp2 and Arp3 into a short pitch helical (filament-like) arrangement, suggesting actin filaments alone do not stimulate the short pitch conformational change. In contrast, contacts with the mother filament stabilize subunit flattening in Arp3, an intrasubunit change that converts Arp3 from a conformation that mimics an actin monomer to one that mimics a filamentous actin subunit. Our results support a multistep activation pathway that has important implications for understanding how WASP-mediated activation allows Arp2/3 complex to assemble force-producing actin networks.

Usefulness of Docking and Molecular Dynamics in Selecting Tumor Neoantigens to Design Personalized Cancer Vaccines: A Proof of Concept

Personalized cancer vaccines based on neoantigens are a new and promising treatment for cancer; however, there are still multiple unresolved challenges to using this type of immunotherapy. Among these, the effective identification of immunogenic neoantigens stands out, since the in silico tools used generate a significant portion of false positives. Inclusion of molecular simulation techniques can refine the results these tools produce. In this work, we explored docking and molecular dynamics to study the association between the stability of peptide-HLA complexes and their immunogenicity, using as a proof of concept two HLA-A2-restricted neoantigens that were already evaluated in vitro. The results obtained were in accordance with the in vitro immunogenicity, since the immunogenic neoantigen ASTN1 remained bound at both ends to the HLA-A2 molecule. Additionally, molecular dynamic simulation suggests that position 1 of the peptide has a more relevant role in stabilizing the N-terminus than previously proposed. Likewise, the mutations may have a "delocalized" effect on the peptide-HLA interaction, which means that the mutated amino acid influences the intensity of the interactions of distant amino acids of the peptide with the HLA. These findings allow us to propose the inclusion of molecular simulation techniques to improve the identification of neoantigens for cancer vaccines.

Crystal structures of glycoprotein D of equine alphaherpesviruses reveal potential binding sites to the entry receptor MHC-I

Cell entry of most alphaherpesviruses is mediated by the binding of glycoprotein D (gD) to different cell surface receptors. Equine herpesvirus type 1 (EHV-1) and EHV-4 gDs interact with equine major histocompatibility complex I (MHC-I) to initiate entry into equine cells. We have characterized the gD-MHC-I interaction by solving the crystal structures of EHV-1 and EHV-4 gDs (gD1, gD4), performing protein-protein docking simulations, surface plasmon resonance (SPR) analysis, and biological assays. The structures of gD1 and gD4 revealed the existence of a common V-set immunoglobulin-like (IgV-like) core comparable to those of other gD homologs. Molecular modeling yielded plausible binding hypotheses and identified key residues (F213 and D261) that are important for virus binding. Altering the key residues resulted in impaired virus growth in cells, which highlights the important role of these residues in the gD-MHC-I interaction. Taken together, our results add to our understanding of the initial herpesvirus-cell interactions and will contribute to the targeted design of antiviral drugs and vaccine development.

Rare deleterious mutations in Bruton's tyrosine kinase as biomarkers for ibrutinib-based therapy: an in silico insight

CONTEXT

Squamous cell carcinoma (SCC) is the second most common type of skin cancer caused by malignant keratinocytes. Multiple studies have shown that protein mutations have a significant impact on the development and progression of cancer, including SCC. We attempted to decode the effect of single amino acid mutations in the Bruton's tyrosine kinase (BTK) protein in this study. Molecular dynamic (MD) simulations were performed on selected deleterious mutations of the BTK protein, revealing that the variants adversely affect the protein, indicating that they may contribute to the prognosis of SCC by making the protein unstable. Then, we investigated the interaction between the protein and its mutants with ibrutinib, a drug designed to treat SCC. Even though the mutations have deleterious effects on protein structure, they bind to ibrutinib similarly to their wild type counterpart. This study demonstrates that the effect of detected missense mutations is unfavorable and can result in function loss, which is severe for SCC, but that ibrutinib-based therapy can still be effective on them, and the mutations can be used as biomarkers for Ibrutinib-based treatment.

METHODS

Seven different computational techniques were used to compute the effect of SAVs in accordance with the experimental requirements of this study. To understand the differences in protein and mutant dynamics, MD simulation and trajectory analysis, including RMSD, RMSF, PCA, and contact analysis, were performed. The free binding energy and its decomposition for each protein-drug complex were determined using docking, MM-GBSA, MM-PBSA, and interaction analysis (wild and mutants).

Insights into substrate coordination and glycosyl transfer of poplar cellulose synthase-8

Cellulose is an abundant cell wall component of land plants. It is synthesized from UDP-activated glucose molecules by cellulose synthase, a membrane-integrated processive glycosyltransferase. Cellulose synthase couples the elongation of the cellulose polymer with its translocation across the plasma membrane. Here, we present substrate and product-bound cryogenic electron microscopy structures of the homotrimeric cellulose synthase isoform-8 (CesA8) from hybrid aspen (poplar). UDP-glucose binds to a conserved catalytic pocket adjacent to the entrance to a transmembrane channel. The substrate's glucosyl unit is coordinated by conserved residues of the glycosyltransferase domain and amphipathic interface helices. Site-directed mutagenesis of a conserved gating loop capping the active site reveals its critical function for catalytic activity. Molecular dynamics simulations reveal prolonged interactions of the gating loop with the substrate molecule, particularly across its central conserved region. These transient interactions likely facilitate the proper positioning of the substrate molecule for glycosyl transfer and cellulose translocation.

Highlights

Cryo-EM structures of substrate and product bound poplar cellulose synthase provide insights into substrate selectivitySite directed mutagenesis signifies a critical function of the gating loop for catalysisMolecular dynamics simulations support persistent gating loop - substrate interactionsGating loop helps in positioning the substrate molecule to facilitate cellulose elongationConserved cellulose synthesis substrate binding mechanism across the kingdoms.

Structural modelling and dynamics of full-length of TLR10 sheds light on possible modes of dimerization, ligand binding and mechanism of action

Toll like receptors (TLRs) play a pivotal role in innate and adaptive immunity. There are 10 TLRs in the human genome, of which TLR10 is the least characterized. Genetic polymorphism of TLR10 has been shown to be associated with multiple diseases including tuberculosis and rheumatoid arthritis. TLR10 consists of an extracellular domain (ECD), a single-pass transmembrane (TM) helix and intracellular TIR (Toll/Interleukin-1 receptor) domain. ECD is employed for ligand recognition and the intracellular domain interacts with other TIR domain-containing adapter proteins for signal transduction. Experimental structure of ECD or TM domain is not available for TLR10. In this study, we have modelled multiple forms of TLR10-ECD dimers, such as closed and open forms, starting from available structures of homologues. Subsequently, multiple full-length TLR10 homodimer models were generated by utilizing homology modelling and protein-protein docking. The dynamics of these models in membrane-aqueous environment revealed the global motion of ECD and TIR domain towards membrane bilayer. The TIR domain residues exhibited high root mean square fluctuation compared to ECD. The 'closed form' model was observed to be energetically more favorable than 'open form' model. The evaluation of persistent interchain interactions, along with their conservation score, unveiled critical residues for each model. Further, the binding of dsRNA to TLR10 was modelled by defined and blind docking approaches. Differential binding of dsRNA to the protomers of TLR10 was observed upon simulation that could provide clues on ligand disassociation. Dynamic network analysis revealed that the 'open form' model can be the functional form while 'closed form' model can be the apo form of TLR10.

Two-pore channel blockade by phosphoinositide kinase inhibitors YM201636 and PI-103 determined by a histidine residue near pore-entrance

Human two-pore channels (TPCs) are endolysosomal cation channels and play an important role in NAADP-evoked Ca²⁺ release and endomembrane dynamics. We found that YM201636, a PIKfyve inhibitor, potently inhibits PI(3,5)P₂-activated human TPC2 with an IC₅₀ of 0.16 μM. YM201636 also effectively inhibits NAADP-activated TPC2 and a constitutively-open TPC2 L690A/L694A mutant channel; whereas it exerts little effect when applied in the channel’s closed state. PI-103, a YM201636 analog and an inhibitor of PI3K and mTOR, also inhibits human TPC2 with an IC₅₀ of 0.64 μM. With mutational, virtual docking, and molecular dynamic simulation analyses, we found that YM201636 and PI-103 directly block the TPC2’s open-state channel pore at the bundle-cross pore-gate region where a nearby H699 residue is a key determinant for channel’s sensitivity to the inhibitors. H699 likely interacts with the blockers around the pore entrance and facilitates their access to the pore. Substitution of a Phe for H699 largely accounts for the TPC1 channel’s insensitivity to YM201636. These findings identify two potent TPC2 channel blockers, reveal a channel pore entrance blockade mechanism, and provide an ion channel target in interpreting the pharmacological effects of two commonly used phosphoinositide kinase inhibitors.

YM201636 and PI-103 are potent inhibitors of human two-pore channel 2 that act through a channel pore entrance blockade mechanism.

The human RNA polymerase I structure reveals an HMG-like docking domain specific to metazoans

We characterize the human RNA polymerase I by evolutionary biochemistry and cryo-EM revealing a built-in structural domain that apparently serves as transcription factor–binding platform in metazoans.

Transcription of the ribosomal RNA precursor by RNA polymerase (Pol) I is a major determinant of cellular growth, and dysregulation is observed in many cancer types. Here, we present the purification of human Pol I from cells carrying a genomic GFP fusion on the largest subunit allowing the structural and functional analysis of the enzyme across species. In contrast to yeast, human Pol I carries a single-subunit stalk, and in vitro transcription indicates a reduced proofreading activity. Determination of the human Pol I cryo-EM reconstruction in a close-to-native state rationalizes the effects of disease-associated mutations and uncovers an additional domain that is built into the sequence of Pol I subunit RPA1. This “dock II” domain resembles a truncated HMG box incapable of DNA binding which may serve as a downstream transcription factor–binding platform in metazoans. Biochemical analysis, in situ modelling, and ChIP data indicate that Topoisomerase 2a can be recruited to Pol I via the domain and cooperates with the HMG box domain–containing factor UBF. These adaptations of the metazoan Pol I transcription system may allow efficient release of positive DNA supercoils accumulating downstream of the transcription bubble.

Targeting COVID-19 pandemic: in silico evaluation of 2-hydroxy-1, 2-diphenylethanone N(4)-methyl-N(4)-phenylthiosemicarbazone as a potential inhibitor of SARS-CoV-2

The global spread of the COVID-19 pandemic caused by the etiological agent, severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), triggered researchers to identify and develop novel antiviral therapeutics. Herein, we report a new molecule 2-hydroxy-1,2-diphenylethanone N(4)-methyl-N(4)-phenyl thiosemicarbazone (BMPTSC), as a potential inhibitor of SARS-CoV-2. BMPTSC was synthesized, characterized by IR and NMR studies, and the structural parameters were analyzed computationally by B3LYP/cc-pVDZ method. Molecular docking studies were performed to get insights into the energetics and compatibility of BMPTSC against various SARS-CoV-2 drug targets. The best docking poses of target protein-BMPTSC complex structures were further subjected to molecular dynamics (MD) simulations. Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations on the binding of BMPTSC with the target proteins viz. spike glycoprotein and ACE-2 protein showed energy values of -179.87 and -145.61 kJ/mol, respectively. Moreover, BMPTSC obeys Lipinski's rule, and further in silico assessment of oral bioavailability, bioactivity scores, ADME, drug-likeness, and medicinal chemistry friendliness suggests that this molecule is a promising candidate for the COVID-19 drug discovery process.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11224-022-02033-8.

On the Tracks of the Aggregation Mechanism of the PHF6 Peptide from Tau Protein: Molecular Dynamics, Energy, and Interaction Network Investigations

The formation of neurofibrillary tangles (NFTs), composed of tau protein aggregates, is a hallmark of some neurodegenerative diseases called tauopathies. NFTs are composed of paired helical filaments (PHFs) of tau protein with a dominant β-sheet secondary structuration. The NFT formation mechanism is not known yet. This study focuses on PHF6, a crucial hexapeptide responsible for tau aggregation. A 2 μs molecular dynamics simulation was launched to determine the keys of the PHF6 aggregation mechanism. Hydrogen bonding, van der Waals, and other non-covalent interactions as π-stacking were investigated. Parallel aggregation was slightly preferred due to its adaptability, but antiparallel aggregation remained widely present during the PHF6 aggregation. The analysis highlighted the leading role of hydrogen bonds identified at the atomic level for each aggregation process. The aggregation study emphasized the importance of Tyr310 during the β-sheets' complexation through π-stacking.

PITRM1 interaction studies with amyloidogenic nonapeptide mutants of familial Alzheimer's disease

Amyloid β-protein (ABP) is found to be the major cause for the development of neurodegeneration which leads to Alzheimer's. The Aβ nonapeptide segment, QKLVFFAED (amino acids 15-23) is the highly amyloidogenic central region of Aβ. Familial mutation in Aβ increases the aggregation property of the peptide compared to the Native (Wild) amyloid-beta (Aβ) and these mutations fall on the Aβ nonapeptide segment. The catalytic activity of pitrilysin metallopeptidase 1(PITRM1) with familial mutant Aβ (Flemish, Arctic, Dutch, Italian and Iowa) during interaction is examined using molecular dynamic simulation. The molecular dynamics simulation of PITRM1 and the Aβ nonapeptide segment showed similar RMSD with respect to stability. The active site amino acid (AA) H108, hydrophobic pocket AA residues L111, F123, F124, and L127 and the basic pocket AA residues R888 and H896 showed similar interactions with both wild and familial Aβ. The molecular level interaction between amyloid beta and PITRM1 were similar in the wild and familial mutants except for the Arctic mutant. The hydrophobic interaction was commonly observed between the S1 hydrophobic pocket and the LVFF region, the Arctic mutant showed less hydrogen bond formation consistently when compared to other complexes. This molecular information on catalytic activity suggests that modulating inactive PITRM1 or an increase in expression of PITRM1 can help in eliminating different kinds of familial mutant Aβ in neurodegenerative cells.Communicated by Ramaswamy H. Sarma.

Computational discovery of small drug-like compounds as potential inhibitors of PD-1/PD-L1 interactions

The programmed cell death ligand protein 1 (PD-L1) is a strong immunosuppressive molecule that inactivates tumor-specific T cells by binding to the programmed cell death- 1 protein (PD-1). Cancer immunotherapy based on the monoclonal antibodies targeting the PD-1/PD-L1 pathway has demonstrated therapeutic responses without precedent over a wide range of cancers. However, the antibody-based immunotherapies have several limitations such as high production cost or the induction of severe immune-related adverse effects. Small-molecule inhibitors of the PD-1/PD-L1 pathway are a promising alternative or complementary therapeutic to antibodies. Currently, the field of developing anti-PD-1/PD-L1 small-molecule inhibitors is intensively explored. In the present study a pharmacophore model was generated based on previously developed compounds and their atomistic structures with the PD-L1 dimer. Structure-based affinity-based virtual screening of small-molecule inhibitors of the PD-1/PD-L1 pathway according to the pharmacophore model followed by a screening in terms of drug-likeness resulted in ten hit compounds of high affinity towards the PD-L1 dimer and the satisfaction to all of the drug-likeness rules. Molecular dynamics (MD) simulations showed that nine of ten compounds formed stable complexes with the PD-L1 dimer as evidenced by the analysis of MD trajectories. Molecular mechanics Poisson- Boltzmann surface area (MM-PBSA) calculation revealed very low binding energies (<-46 kcal/mol) for the interactions of these ligands with the PD-L1 dimer, suggesting that identified compounds may serve as good scaffolds for the design of novel agents of antitumor immunotherapy able to target the PD-1/PD-L1 interactionCommunicated by Ramaswamy H. Sarma.

Identification of some dietary flavonoids as potential inhibitors of TMPRSS2 through protein-ligand interaction studies and binding free energy calculations

The continuing threat of COVID-19 and deaths need an urgent cost-effective pharmacological approach. Here, we examine the inhibitory activity of a group of dietary bioactive flavonoids against the human protease TMPRSS2, which plays a major role in SARS CoV-2 viral entry. After the molecular docking studies of a large number of flavonoids, four compounds with high binding scores were selected and studied in detail. The binding affinities of these four ligands, Amentoflavone, Narirutin, Eriocitrin, and Naringin, at the active site of the TMPRSS2 target, were investigated using MD simulations followed by MM-PBSA binding energy calculations. From the studies, a number of significant hydrophobic and hydrogen bonding interactions between the ligands and binding site amino residues of TMPRSS2 are identified which showcase their excellent inhibitory activity against TMPRSS2. Among these ligands, Amentoflavone and Narirutin showed MM-PBSA binding energy values of -155.57 and -139.71 kJ/mol, respectively. Our previous studies of the inhibitory activity of these compounds against the main protease of SARS-COV2 and the present study on TMPRSS2 strongly highlighted that Amentoflavone and Naringin can exhibit promising multi-target activity against SARS-CoV-2. Moreover, due to their wide availability, no side effects, and low cost, these compounds could be recommended as dietary supplements for COVID patients or for the development of SARS-CoV-2 treatments.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11224-022-01955-7.

Mechanism behind Polysorbates' Inhibitory Effect on P-Glycoprotein

Much effort has been invested in the search for modulators of membrane transport proteins such as P-glycoprotein (P-gp) to improve drug bioavailability and reverse multidrug resistance in cancer. Nonionic surfactants, a class of pharmaceutical excipients, are known to inhibit such proteins, but knowledge about the exact mechanism of this inhibition is scarce. Here, we perform multiscale molecular dynamics simulations of one of these surfactants, polysorbate 20 (PS20), to reveal the behavior of such compounds on the molecular level and thereby discover the molecular mechanism of the P-gp inhibition. We show that the amphiphilic headgroup of PS20 is too hydrophobic to partition in the water phase, which drives the binding of PS20 to the amphiphilic drug-binding domain of P-gp and thereby causes the inhibition of the protein. Based on our findings, we conclude that PS20 primarily inhibits P-gp through direct binding to the drug-binding domain (DBD) from the extracellular leaflet.

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

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

SPARC: Structural properties associated with residue constraints

SPARC facilitates the generation of plausible hypotheses regarding underlying biochemical mechanisms by structurally characterizing protein sequence constraints. Such constraints appear as residues co-conserved in functionally related subgroups, as subtle pairwise correlations (i.e., direct couplings), and as correlations among these sequence features or with structural features. SPARC performs three types of analyses. First, based on pairwise sequence correlations, it estimates the biological relevance of alternative conformations and of homomeric contacts, as illustrated here for death domains. Second, it estimates the statistical significance of the correspondence between directly coupled residue pairs and interactions at heterodimeric interfaces. Third, given molecular dynamics simulated structures, it characterizes interactions among constrained residues or between such residues and ligands that: (a) are stably maintained during the simulation; (b) undergo correlated formation and/or disruption of interactions with other constrained residues; or (c) switch between alternative interactions. We illustrate this for two homohexameric complexes: the bacterial enhancer binding protein (bEBP) NtrC1, which activates transcription by remodeling RNA polymerase (RNAP) containing σ⁵⁴, and for DnaB helicase, which opens DNA at the bacterial replication fork. Based on the NtrC1 analysis, we hypothesize possible mechanisms for inhibiting ATP hydrolysis until ADP is released from an adjacent subunit and for coupling ATP hydrolysis to restructuring of σ⁵⁴ binding loops. Based on the DnaB analysis, we hypothesize that DnaB 'grabs' ssDNA by flipping every fourth base and inserting it into cavities between subunits and that flipping of a DnaB-specific glutamine residue triggers ATP hydrolysis.

SINAPs: A Software Tool for Analysis and Visualization of Interaction Networks of Molecular Dynamics Simulations

As long as the structural study of molecular mechanisms requires multiple molecular dynamics reflecting contrasted bioactive states, the subsequent analysis of molecular interaction networks remains a bottleneck to be fairly treated and requires a user-friendly 3D view of key interactions. Structural Interaction Network Analysis Protocols (SINAPs) is a proprietary python tool developed to (i) quickly solve key interactions able to distinguish two protein states, either from two sets of molecular dynamics simulations or from two crystallographic structures, and (ii) render a user-friendly 3D view of these key interactions through a plugin of UCSF Chimera, one of the most popular open-source viewing software for biomolecular systems. Through two case studies, glucose transporter-1 (GLUT-1) and A2A adenosine receptor (A2AR), SINAPs easily pinpointed key interactions observed experimentally and relevant for their bioactivities. This very effective tool was thus applied to identify the amino acids involved in the molecular enzymatic mechanisms ruling the activation of an immunomodulator drug candidate, P28 glutathione-S-transferase (P28GST). SINAPs is freely available at https://github.com/ParImmune/SINAPs.

Successes and challenges in multiscale modelling of artificial metalloenzymes: the case study of POP-Rh2 cyclopropanase

Molecular modelling applications in metalloenzyme design are still scarce due to a series of challenges. On top of that, the simulations of metal-mediated binding and the identification of catalytic competent geometries require both large conformational exploration and simulation of fine electronic properties. Here, we demonstrate how the incorporation of new tools in multiscale strategies, namely substrate diffusion exploration, allows taking a step further. As a showcase, the enantioselective profiles of the most outstanding variants of an artificial Rh₂-based cyclopropanase (GSH, HFF and RFY) developed by Lewis and co-workers (Nat. Commun., 2015, 6, 7789 and Nat. Chem., 2018, 10, 318-324) have been rationalized. DFT calculations on the free-cofactor-mediated process identify the carbene insertion and the cyclopropanoid formation as crucial events, the latter being the enantiodetermining step, which displays up to 8 competitive orientations easily altered by the protein environment. The key intermediates of the reaction were docked into the protein scaffold showing that some mutated residues have direct interaction with the cofactor and/or the co-substrate. These interactions take the form of a direct coordination of Rh in GSH and HFF and a strong hydrophobic patch with the carbene moiety in RFY. Posterior molecular dynamics sustain that the cofactor induces global re-arrangements of the protein. Finally, massive exploration of substrate diffusion, based on the GPathFinder approach, defines this event as the origin of the enantioselectivity in GSH and RFY. For HFF, fine molecular dockings suggest that it is likely related to local interactions upon diffusion. This work shows how modelling of long-range mutations on the catalytic profiles of metalloenzymes may be unavoidable and software simulating substrate diffusion should be applied.

Allosteric control of Ubp6 and the proteasome via a bidirectional switch

The proteasome recognizes ubiquitinated proteins and can also edit ubiquitin marks, allowing substrates to be rejected based on ubiquitin chain topology. In yeast, editing is mediated by deubiquitinating enzyme Ubp6. The proteasome activates Ubp6, whereas Ubp6 inhibits the proteasome through deubiquitination and a noncatalytic effect. Here, we report cryo-EM structures of the proteasome bound to Ubp6, based on which we identify mutants in Ubp6 and proteasome subunit Rpt1 that abrogate Ubp6 activation. The Ubp6 mutations define a conserved region that we term the ILR element. The ILR is found within the BL1 loop, which obstructs the catalytic groove in free Ubp6. Rpt1-ILR interaction opens the groove by rearranging not only BL1 but also a previously undescribed network of three interconnected active-site-blocking loops. Ubp6 activation and noncatalytic proteasome inhibition are linked in that they are eliminated by the same mutations. Ubp6 and ubiquitin together drive proteasomes into a unique conformation associated with proteasome inhibition. Thus, a multicomponent allosteric switch exerts simultaneous control over both Ubp6 and the proteasome.

Optimization of the Urea Removal in a Wearable Dialysis Device Using Nitrogen-Doped and Phosphorus-Doped Graphene

Dialysis has been recognized as an essential treatment for end-stage renal disease (ESRD). This therapy, however, suffers from several limitations leading to numerous complications in the patients. As dialysis cannot completely substitute healthy kidney functions, the health condition of an ESRD patient is ultimately affected. Wearable artificial kidney (WAK) can resolve the restrictions of blood purification by the dialysis method. However, absorbing large amounts of urea produced in the body is one of the main challenges of these WAK and overcoming this is necessary to improve both functionality and footprint of the device. This study investigates the adsorption capabilities of N- and P-doped graphene nanosorbents for the first time by using molecular dynamic simulation. Urea removal on carbon nanosheets was simulated with different percentages of phosphorus and nitrogen dopants along with the pristine graphene. Specifically, the effects of interaction energy, adsorption percentage, gyration radius, hydrogen bonding, and other molecular dynamic analyses on urea removal were also investigated. The results from this study match well with the existing research, demonstrating the accuracy of the model. The results further suggest that graphene nanosheets doped by 10% nitrogen are likely the most effective in removing urea given that it is associated with the maximum radial distribution function (RDF), the maximum reduction in gyration radius, a high number of hydrogen bonds, and the most negative adsorption energy. This molecular study offers attractive suggestions for the novel adsorbents of artificial kidney devices and paves the way for the development of novel and enhanced urea adsorbents.

OUP accepted manuscript

Correct identification and effective visualization of interactions in biomolecular structures facilitate understanding of their functions and molecular design. In response to the practical needs of structure-based analysis, we have created a Mapiya web server. The Mapiya integrates four main functionalities: (i) generation of contact maps - intramolecular and intermolecular-for proteins, nucleic acids, and their complexes; (ii) characterization of the interactions physicochemical nature, (iii) interactive visualization of biomolecular conformations with automatic zoom on selected contacts using Molstar and (iv) additional sequence- and structure-based analyses performed with third-party software and in-house algorithms combined into an easy-to-use interface. Thus, Mapiya offers a highly customized analysis of the molecular interactions' in various biological systems. The web server is available at: http://mapiya.lcbio.pl/.

Drug Repurposing: Deferasirox Inhibits the Anti-Apoptotic Activity of Mcl-1

Introduction

With the aim of repositioning commercially available drugs for the inhibition of the anti-apoptotic myeloid cell leukemia protein, Mcl-1, implied in various cancers, five molecules, highlighted from a published theoretical screening, were selected to experimentally validate their affinity toward Mcl-1.

Results

A detailed NMR study revealed that only two of the five tested drugs, Torsemide and Deferasirox, interacted with Mcl-1. NMR data analysis allowed the complete characterization of the binding mode of both drugs to Mcl-1, including the estimation of their affinity for Mcl-1. Biological assays evidenced that the biological activity of Torsemide was lower as compared to the Deferasirox, which was able to efficiently and selectively inhibit the anti-apoptotic activity of Mcl-1. Finally, docking and molecular dynamics led to a 3D model for the Deferasirox:Mcl-1 complex and revealed the positioning of the drug in the Mcl-1 P2/P3 pockets as well as almost all synthetic Mcl-1 inhibitors. Interestingly, contrary to known synthetic Mcl-1 inhibitors which interact through Arg263, Deferasirox, establishes a salt bridge with Lys234.

Conclusion

Deferasirox could be a potential candidate for drug repositioning as Mcl-1 inhibitor.

PIP2-induced membrane binding of the vinculin tail competes with its other binding partners

Vinculin plays a key role during the first phase of focal adhesion formation and interacts with the plasma membrane through specific binding of its tail domain to the lipid phosphatidylinositol 4,5-bisphosphate (PIP₂). Our understanding of the PIP₂-vinculin interaction has been hampered by contradictory biochemical and structural data. Here, we used a multiscale molecular dynamics simulation approach, in which unbiased coarse-grained molecular dynamics were used to generate starting structures for subsequent microsecond-long all-atom simulations. This allowed us to map the interaction of the vinculin tail with PIP₂-enriched membranes in atomistic detail. In agreement with experimental data, we have shown that membrane binding is sterically incompatible with the intramolecular interaction between vinculin's head and tail domain. Our simulations further confirmed biochemical and structural results, which identified two positively charged surfaces, the basic collar and the basic ladder, as the main PIP₂ interaction sites. By introducing a valency-disaggregated binding network analysis, we were able to map the protein-lipid interactions in unprecedented detail. In contrast to the basic collar, in which PIP₂ is specifically recognized by an up to hexavalent binding pocket, the basic ladder forms a series of low-valency binding sites. Importantly, many of these PIP₂ binding residues are also involved in maintaining vinculin in a closed, autoinhibited conformation. These findings led us to propose a molecular mechanism for the coupling between vinculin activation and membrane binding. Finally, our refined binding site suggests an allosteric relationship between PIP₂ and F-actin binding that disfavors simultaneous interaction with both ligands, despite nonoverlapping binding sites.

CDR1 Composition Can Affect Nanobody Recombinant Expression Yields

The isolation of nanobodies from pre-immune libraries by means of biopanning is a straightforward process. Nevertheless, the recovered candidates often require optimization to improve some of their biophysical characteristics. In principle, CDRs are not mutated because they are likely to be part of the antibody paratope, but in this work, we describe a mutagenesis strategy that specifically addresses CDR1. Its sequence was identified as an instability hot spot by the PROSS program, and the available structural information indicated that four CDR1 residues bound directly to the antigen. We therefore modified the loop flexibility with the addition of an extra glycine rather than by mutating single amino acids. This approach significantly increased the nanobody yields but traded-off with moderate affinity loss. Accurate modeling coupled with atomistic molecular dynamics simulations enabled the modifications induced by the glycine insertion and the rationale behind the engineering design to be described in detail.

Polymorphisms at site 469 of B-RAF protein associated with skin melanoma may be correlated with dabrafenib resistance: An in silico study

Melanoma is a type of skin cancer. Numerous genes and their proteins are strongly associated with melanoma susceptibility. This study aims to use an in silico method to identify genetic variants in the melanoma susceptibility gene. The COSMIC database was queried for genes and cross-referenced with three environment-gene interaction databases (EGP, SeattleSNPs and CTD) to identify shared genes. The majority of approved skin melanoma drugs were found to act on the protein serine/threonine-protein kinase (B-RAF) encoded by the BRAF gene, which was also present in all three referenced databases. Comprehensive computational analysis was performed to predict deleterious genetic variants associated with skin melanoma, and the nsSNPs G469V and G469E were prioritized based on their predicted deleterious effects. Molecular dynamic simulation analysis of the B-RAF protein mutants G469V and G469E reveals that variations in the amino acid conformation at the drug binding site result in inconsistency in drug interaction. Additionally, this analysis showed that the G469V and G469E mutants have lower binding energy for dabrafenib than the wild type. The population with the highest frequency of each deleterious and pathogenic variant has been determined. The study's findings would support the development of more effective treatment strategies for skin melanoma. Communicated by Ramaswamy H. Sarma.

VoroContacts: a tool for the analysis of interatomic contacts in macromolecular structures

SUMMARY

VoroContacts is a versatile tool for computing and analyzing contact surface areas (CSAs) and solvent accessible surface areas (SASAs) for 3 D structures of proteins, nucleic acids and their complexes at the atomic resolution. CSAs and SASAs are derived using Voronoi tessellation of 3 D structure, represented as a collection of atomic balls. VoroContacts web server features a highly configurable query interface, which enables on-the-fly analysis of contacts for selected set of atoms and allows filtering interatomic contacts by their type, surface areas, distance between contacting atoms and sequence separation between contacting residues. The VoroContacts functionality is also implemented as part of the standalone Voronota package, enabling batch processing.

AVAILABILITY AND IMPLEMENTATION

https://bioinformatics.lt/wtsam/vorocontacts.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Structural Design and Analysis of the RHOA-ARHGEF1 Binding Mode: Challenges and Applications for Protein-Protein Interface Prediction

The interaction between two proteins may involve local movements, such as small side-chains re-positioning or more global allosteric movements, such as domain rearrangement. We studied how one can build a precise and detailed protein-protein interface using existing protein-protein docking methods, and how it can be possible to enhance the initial structures using molecular dynamics simulations and data-driven human inspection. We present how this strategy was applied to the modeling of RHOA-ARHGEF1 interaction using similar complexes of RHOA bound to other members of the Rho guanine nucleotide exchange factor family for comparative assessment. In parallel, a more crude approach based on structural superimposition and molecular replacement was also assessed. Both models were then successfully refined using molecular dynamics simulations leading to protein structures where the major data from scientific literature could be recovered. We expect that the detailed strategy used in this work will prove useful for other protein-protein interface design. The RHOA-ARHGEF1 interface modeled here will be extremely useful for the design of inhibitors targeting this protein-protein interaction (PPI).

Time-resolved spectroscopic and electrophysiological data reveal insights in the gating mechanism of anion channelrhodopsin

Channelrhodopsins are widely used in optogenetic applications. High photocurrents and low current inactivation levels are desirable. Two parallel photocycles evoked by different retinal conformations cause cation-conducting channelrhodopsin-2 (CrChR2) inactivation: one with efficient conductivity; one with low conductivity. Given the longer half-life of the low conducting photocycle intermediates, which accumulate under continuous illumination, resulting in a largely reduced photocurrent. Here, we demonstrate that for channelrhodopsin-1 of the cryptophyte Guillardia theta (GtACR1), the highly conducting C = N-anti-photocycle was the sole operating cycle using time-resolved step-scan FTIR spectroscopy. The correlation between our spectroscopic measurements and previously reported electrophysiological data provides insights into molecular gating mechanisms and their role in the characteristic high photocurrents. The mechanistic importance of the central constriction site amino acid Glu-68 is also shown. We propose that canceling out the poorly conducting photocycle avoids the inactivation observed in CrChR2, and anticipate that this discovery will advance the development of optimized optogenetic tools.

Molecular mechanisms of ion conduction and ion selectivity in TMEM16 lipid scramblases

TMEM16 lipid scramblases transport lipids and also operate as ion channels with highly variable ion selectivities and various physiological functions. However, their molecular mechanisms of ion conduction and selectivity remain largely unknown. Using computational electrophysiology simulations at atomistic resolution, we identified the main ion-conductive state of TMEM16 lipid scramblases, in which an ion permeation pathway is lined by lipid headgroups that directly interact with permeating ions in a voltage polarity-dependent manner. We found that lipid headgroups modulate the ion-permeability state and regulate ion selectivity to varying degrees in different scramblase isoforms, depending on the amino-acid composition of the pores. Our work has defined the structural basis of ion conduction and selectivity in TMEM16 lipid scramblases and uncovered the mechanisms responsible for the direct effects of membrane lipids on the conduction properties of ion channels.

Seeking the exclusive binding region of phenylalkylamine derivatives on human T-type calcium channels via homology modeling and molecular dynamics simulation approach

Pharmaceutical features of phenylalkylamine derivatives (PAAs) binding to calcium channels have been studied extensively in the past decades. Only a few PAAs have the binding specificity on calcium channels, for example, NNC 55‐0396. Here, we created the homology models of human Ca_v3.2, Ca_v3.3 and use them as a receptor on the rigid docking tests. The nonspecific calcium channel blocker mibefradil showed inconsistent docking preference across four domains; however, NNC 55‐0396 had a unique binding pattern on domain II specifically. The subsequent molecular dynamics (MD) simulations identified that Ca_v3.1, Ca_v3.2, and Ca_v3.3 share domain II when Ca²⁺ appearing in the neighbor region of selective filters (SFs). Moreover, free‐energy perturbation analysis suggests single mutation of lysine at P‐loop domain III, or threonine at the P‐loop domain II largely reduced the total amount of hydration‐free energy in the system. All these findings suggest that P‐loop and segment six domain II in the T‐type calcium channels (TCCs) are crucial for attracting the PAAs with specificity as the antagonist.

The Ras dimer structure

Oncogenic mutated Ras is a key player in cancer, but despite intense and expensive approaches its catalytic center seems undruggable. The Ras dimer interface is a possible alternative drug target. Dimerization at the membrane affects cell growth signal transduction. In vivo studies indicate that preventing dimerization of oncogenic mutated Ras inhibits uncontrolled cell growth. Conventional computational drug-screening approaches require a precise atomic dimer model as input to successfully access drug candidates. However, the proposed dimer structural models are controversial. Here, we provide a clear-cut experimentally validated N-Ras dimer structural model. We incorporated unnatural amino acids into Ras to enable the binding of labels at multiple positions via click chemistry. This labeling allowed the determination of multiple distances of the membrane-bound Ras-dimer measured by fluorescence and electron paramagnetic resonance spectroscopy. In combination with protein-protein docking and biomolecular simulations, we identified key residues for dimerization. Site-directed mutations of these residues prevent dimer formation in our experiments, proving our dimer model to be correct. The presented dimer structure enables computational drug-screening studies exploiting the Ras dimer interface as an alternative drug target.

g_elpot: A Tool for Quantifying Biomolecular Electrostatics from Molecular Dynamics Trajectories

Electrostatic forces drive a wide variety of biomolecular processes by defining the energetics of the interaction between biomolecules and charged substances. Molecular dynamics (MD) simulations provide trajectories that contain ensembles of structural configurations sampled by biomolecules and their environment. Although this information can be used for high-resolution characterization of biomolecular electrostatics, it has not yet been possible to calculate electrostatic potentials from MD trajectories in a way allowing for quantitative connection to energetics. Here, we present g_elpot, a GROMACS-based tool that utilizes the smooth particle mesh Ewald method to quantify the electrostatics of biomolecules by calculating potential within water molecules that are explicitly present in biomolecular MD simulations. g_elpot can extract the global distribution of the electrostatic potential from MD trajectories and measure its time course in functionally important regions of a biomolecule. To demonstrate that g_elpot can be used to gain biophysical insights into various biomolecular processes, we applied the tool to MD trajectories of the P2X3 receptor, TMEM16 lipid scramblases, the secondary-active transporter Glt_Ph, and DNA complexed with cationic polymers. Our results indicate that g_elpot is well suited for quantifying electrostatics in biomolecular systems to provide a deeper understanding of its role in biomolecular processes.

An anti-HER2 nanobody binds to its antigen HER2 via two independent paratopes

High-resolution structural data of complexes between antibodies and membrane receptors still represent a demanding task. In this study, we used complementary sets of experimental data to obtain a structural model of the complex formed by the human epidermal growth factor receptor 2 (HER2) and its specific nanobody A10. First we identified by NMR the residues that bind or rearrange as a consequence of the complex formation. In parallel, the complex was cross-linked, digested and the resulting peptides were characterized by mass-spectrometry to define maximal distance restraints between HER2 and A10 amino acids in their complex. These independent datasets guided a docking process, refined by molecular dynamics simulations, to develop a model of the complex and estimate per-residue free-energy contributions. Such a model explains the experimental data and identifies a second, non-canonical paratope, located in the region opposite to the conventional nanobody paratope, formed by the hypervariable loop regions LH1 and LH3. Both paratopes contributed substantially to the overall affinity by binding to independent HER2 epitopes. Nanobody mutants with substitution of key interaction residues, as indicated by the model, possess significantly lower affinity for HER2. This is the first described case of a "natural" biparatopic nanobody, directly selected by in-vitro panning.

In silico investigation of critical binding pattern in SARS-CoV-2 spike protein with angiotensin-converting enzyme 2

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is a newly-discovered coronavirus and responsible for the spread of coronavirus disease 2019 (COVID-19). SARS-CoV-2 infected millions of people in the world and immediately became a pandemic in March 2020. SARS-CoV-2 belongs to the beta-coronavirus genus of the large family of Coronaviridae. It is now known that its surface spike glycoprotein binds to the angiotensin-converting enzyme-2 (ACE2), which is expressed on the lung epithelial cells, mediates the fusion of the cellular and viral membranes, and facilitates the entry of viral genome to the host cell. Therefore, blocking the virus-cell interaction could be a potential target for the prevention of viral infection. The binding of SARS-CoV-2 to ACE2 is a protein-protein interaction, and so, analyzing the structure of the spike glycoprotein of SARS-CoV-2 and its underlying mechanism to bind the host cell receptor would be useful for the management and treatment of COVID-19. In this study, we performed comparative in silico studies to deeply understand the structural and functional details of the interaction between the spike glycoprotein of SARS-CoV-2 and its cognate cellular receptor ACE2. According to our results, the affinity of the ACE2 receptor for SARS-CoV-2 was higher than SARS-CoV. According to the free energy decomposition of the spike glycoprotein-ACE2 complex, we found critical points in three areas which are responsible for the increased binding affinity of SARS-CoV-2 compared with SARS-CoV. These mutations occurred at the receptor-binding domain of the spike glycoprotein that play an essential role in the increasing the affinity of coronavirus to ACE2. For instance, mutations Pro462Ala and Leu472Phe resulted in the altered binding energy from - 2 kcal mol-1 in SARS-COV to - 6 kcal mol-1 in SARS-COV-2. The results demonstrated that some mutations in the receptor-binding motif could be considered as a hot-point for designing potential drugs to inhibit the interaction between the spike glycoprotein and ACE2.

Rounding Out the Understanding of ACD Toxicity with the Discovery of Cyclic Forms of Actin Oligomers

Actin is an essential element of both innate and adaptive immune systems and can aid in motility and translocation of bacterial pathogens, making it an attractive target for bacterial toxins. Pathogenic Vibrio and Aeromonas genera deliver actin cross-linking domain (ACD) toxin into the cytoplasm of the host cell to poison actin regulation and promptly induce cell rounding. At early stages of toxicity, ACD covalently cross-links actin monomers into oligomers (AOs) that bind through multivalent interactions and potently inhibit several families of actin assembly proteins. At advanced toxicity stages, we found that the terminal protomers of linear AOs can get linked together by ACD to produce cyclic AOs. When tested against formins and Ena/VASP, linear and cyclic AOs exhibit similar inhibitory potential, which for the cyclic AOs is reduced in the presence of profilin. In coarse-grained molecular dynamics simulations, profilin and WH2-motif binding sites on actin subunits remain exposed in modeled AOs of both geometries. We speculate, therefore, that the reduced toxicity of cyclic AOs is due to their reduced configurational entropy. A characteristic feature of cyclic AOs is that, in contrast to the linear forms, they cannot be straightened to form filaments (e.g., through stabilization by cofilin), which makes them less susceptible to neutralization by the host cell.

Computational and network pharmacology analysis of bioflavonoids as possible natural antiviral compounds in COVID-19

Bioflavonoids are the largest group of plant-derived polyphenolic compounds with diverse biological potential and have also been proven efficacious in the treatment of Severe Acute Respiratory Syndrome (SARS) and Middle East Respiratory Syndrome (MERS). The present investigation validates molecular docking, simulation, and MM-PBSA studies of fifteen bioactive bioflavonoids derived from plants as a plausible potential antiviral in the treatment of COVID-19. Molecular docking studies for 15 flavonoids on the three SARS CoV-2 proteins, non-structural protein-15 Endoribonuclease (NSP15), the receptor-binding domain of spike protein (RBD of S protein), and main protease (Mpro/3CLpro) were performed and selected protein-ligand complexes were subjected to Molecular Dynamics simulations. The molecular dynamics trajectories were subjected to free energy calculation by the MM-PBSA method. All flavonoids were further assessed for their effectiveness as adjuvant therapy by network pharmacology analysis on the target proteins. The network pharmacology analysis suggests the involvement of selected bioflavonoids in the modulation of multiple signaling pathways like p53, FoxO, MAPK, Wnt, Rap1, TNF, adipocytokine, and leukocyte transendothelial migration which plays a significant role in immunomodulation, minimizing the oxidative stress and inflammation. Molecular docking and molecular dynamics simulation studies illustrated the potential of glycyrrhizic acid, amentoflavone, and mulberroside in inhibiting key SARS-CoV-2 proteins and these results could be exploited further in designing future ligands from natural sources.

A structural model of the PriB-DnaT complex in Escherichia coli replication restart

In Escherichia coli, DNA replication is restarted following DNA repair by the PriA-dependent pathway, in which the binding and dissociation of proteins such as PriA, PriB, and DnaT on ssDNA lead to the formation of a protein-DNA complex for recruiting the DnaB-DnaC replication protein complex. However, the structure of the PriB-DnaT complex, which is an essential step in the PriA-dependent pathway, remains elusive. In this study, the importance of His26 in PriB for replication restart was reconfirmed using plasmid complementation. Furthermore, we used NMR to examine the DnaT interaction sites on PriB. We also evaluated the PriB-DnaT peptide complex model, which was prepared by in silico docking, using molecular dynamic simulation. From these data, we propose a structural model that provides insight into the PriB-DnaT interaction.