Modeling of loops in protein structures

Discovery of novel gating checkpoints in the Orai1 calcium channel by systematic analysis of constitutively active mutants of its paralogs and orthologs

In humans, there are three paralogs of the Orai Ca2+ channel that form the core of the store-operated calcium entry (SOCE) machinery. While the STIM-mediated gating mechanism of Orai channels is still under active investigation, several artificial and natural variants are known to cause constitutive activity of the human Orai1 channel. Surprisingly, little is known about the conservation of the gating checkpoints among the different human Orai paralogs and orthologs in other species. In our work, we show that the mutation corresponding to the activating mutation H134A in transmembrane helix 2 (TM2) of human Orai1 also activates Orai2 and Orai3, likely via a similar mechanism. However, this cross-paralog conservation does not apply to the "ANSGA" nexus mutations in TM4 of human Orai1, which is reported to mimic the STIM1-activated state of the channel. In investigating the mechanistic background of these differences, we identified two positions, H171 and F246 in human Orai1, that are not conserved among paralogs and that seem to be crucial for the channel activation triggered by the "ANSGA" mutations in Orai1. However, mutations of the same residues still allow gating of Orai1 by STIM1, suggesting that the ANSGA mutant of Orai1 may not be a surrogate for the STIM1-activated state of the Orai1 channel. Our results shed new light on these important gating checkpoints and show that the gating mechanism of Orai channels is affected by multiple factors that are not necessarily conserved among orai homologs, such as the TM4-TM3 coupling.

The SARS-CoV-2 helicase as a target for antiviral therapy: Identification of potential small molecule inhibitors by in silico modelling

Although vaccines that provide protection against severe illness from coronavirus disease (COVID-19) have been made available, emerging variant strains of severe acute respiratory syndrome 2 coronavirus 2 (SARS-CoV-2) are of concern. A different research direction involves investigation of antiviral therapeutics. In addition to structural proteins, the SARS-CoV-2 non-structural proteins are of interest and this includes the helicase (nsp13). In this study, an initial screen of 300 ligands was performed to identify potential inhibitors of the SARS-CoV-2 nsp13 examining the nucleoside triphosphatase site (NTPase activity) as the target region. The antiviral activity of polyphenols has been previously reported in the literature and as a result, the phenolic compounds and fatty acids from the OliveNet™ library were utilised. Synthetic compounds with antimicrobial and anti-inflammatory properties were also selected. The structures of the SARS-CoV and MERS-CoV helicases, as well as the human RECQ-like DNA helicase, DHX9 helicase, PcrA helicase, hepatitis C NS3 helicase, and mouse Dna2 nuclease-helicase were used for comparison. As expected, sequence and structural homology between the various species was evident. A number of broad-spectrum and well-known inhibitors interacted with the NTPase active site highlighting the need to potentially identify more specific inhibitors for SARS-CoV-2. Acetylcysteine, clavulanic acid and homovanillic acid were identified as potential lead compounds for the SARS-CoV-2 helicase. Molecular dynamics simulations were performed with the leads bound to the SARS-CoV-2 helicase for 200 ns in triplicate, with favourable binding free energies to the NTPase site. Given their availability, further exploration of their potential inhibitory activity could be considered.

Protein-Nucleic Acid Interactions for RNA Polymerase II Elongation Factors by Molecular Dynamics Simulations

RNA polymerase II (Pol II) forms a complex with elongation factors to proceed to the elongation stage of the transcription process. In this work, we studied the elongation factor SPT5 and explored the protein-nucleic acid interactions for the isolated systems of KOW1 and KOW4 domains of SPT5 with DNA and RNA, respectively. We performed molecular dynamics (MD) simulations using three commonly used force fields that are CHARMM c36m, AMBER ff14sb, and ff19sb. Simulations showed strong protein-nucleic acid interactions and low electrostatic binding free energies for all force fields used. RNA was found to be highly dynamic with all force fields, while DNA had relatively more stable conformations with the AMBER force fields compared to that with CHARMM. Furthermore, we performed MD simulations of the complete elongation complex using CHARMM c36m and AMBER ff19sb force fields to compare the dynamics and interactions with the isolated systems. Similarly, strong KOW1 and DNA interactions were observed in the complete elongation complex simulations and DNA was further stabilized by a network of interactions involving SPT5-KOW1, SPT4, and rpb2 of Pol II. Overall, our study showed that the differences between CHARMM and AMBER force fields strongly affect the dynamics of the nucleic acids. CHARMM provides highly flexible DNA, while AMBER largely stabilizes the DNA structure. Although the presence of the entire interaction network stabilized the DNA and decreased the differences in the results from the two force fields, the discrepancies of the force fields for smaller systems may reflect their problems in generating accurate dynamics of nucleic acids.

GalaxyWater-CNN: Prediction of Water Positions on the Protein Structure by a 3D-Convolutional Neural Network

Proteins interact with numerous water molecules to perform their physiological functions in biological organisms. Most water molecules act as solvent media; hence, their roles may be considered implicitly in theoretical treatments of protein structure and function. However, some water molecules interact intimately with proteins and require explicit treatment to understand their effects. Most physics-based computational methods are limited in their ability to accurately locate water molecules on protein surfaces because of inaccurate energy functions. Instead of relying on an energy function, this study attempts to learn the locations of water molecules from structural data. GalaxyWater-convolutional neural network (CNN) predicts water positions on protein chains, protein-protein interfaces, and protein-compound binding sites using a 3D-CNN model that is trained to generate a water score map on a given protein structure. The training data are compiled from high-resolution protein crystal structures resolved together with water molecules. GalaxyWater-CNN shows improved water prediction performance both in the coverage of crystal water molecules and in the accuracy of the predicted water positions when compared with previous energy-based methods. This method shows a superior performance in predicting water molecules that form hydrogen-bond networks precisely. The web service and the source code of this water prediction method are freely available at https://galaxy.seoklab.org/gwcnn and https://github.com/seoklab/GalaxyWater-CNN, respectively.

Prolyl Isomerization-Mediated Conformational Changes Define ATR Subcellular Compartment-Specific Functions

ATR is a PI3K-like kinase protein, regulating checkpoint responses to DNA damage and replication stress. Apart from its checkpoint function in the nucleus, ATR actively engages in an antiapoptotic role at mitochondria following DNA damage. The different functions of ATR in the nucleus and cytoplasm are carried out by two prolyl isomeric forms of ATR: trans- and cis-ATR, respectively. The isomerization occurs at the Pin1 Ser428-Pro429 motif of ATR. Here, we investigated the structural basis of the subcellular location-specific functions of human ATR. Using a mass spectrometry-based footprinting approach, the surface accessibility of ATR lysine residues to sulfo-NHS-LC-biotin modification was monitored and compared between the cis- and the trans-isomers. We have identified two biotin-modified lysine residues, K459 and K469, within the BH3-like domain of cis-ATR that were not accessible in trans-ATR, indicating a conformational change around the BH3 domain between cis- and trans-ATR. The conformational alteration also involved the N-terminal domain and the middle HEAT domain. Moreover, experimental results from an array of complementary assays show that cis-ATR with the accessible BH3 domain was able to bind to tBid while trans-ATR could not. In addition, both cis- and trans-ATR can directly form homodimers via their C-terminal domains without ATRIP, while nuclear (trans-ATR) in the presence of ATRIP forms dimer-dimer complexes involving both N- and C-termini of ATR and ATRIP after UV. Structural characteristics around the Ser428-Pro429 motif and the BH3 domain region are also analyzed by molecular modeling and dynamics simulation. In support, cis conformation was found to be significantly more energetically favorable than trans at the Ser428-Pro429 bond in a 20-aa wild-type ATR peptide. Taken together, our results suggest that the isomerization-induced structural changes of ATR define both its subcellular location and compartment-specific functions and play an essential role in promoting cell survival and DNA damage responses.

PPI-Affinity: A Web Tool for the Prediction and Optimization of Protein-Peptide and Protein-Protein Binding Affinity

Virtual screening of protein–protein and protein–peptide interactions is a challenging task that directly impacts the processes of hit identification and hit-to-lead optimization in drug design projects involving peptide-based pharmaceuticals. Although several screening tools designed to predict the binding affinity of protein–protein complexes have been proposed, methods specifically developed to predict protein–peptide binding affinity are comparatively scarce. Frequently, predictors trained to score the affinity of small molecules are used for peptides indistinctively, despite the larger complexity and heterogeneity of interactions rendered by peptide binders. To address this issue, we introduce PPI-Affinity, a tool that leverages support vector machine (SVM) predictors of binding affinity to screen datasets of protein–protein and protein–peptide complexes, as well as to generate and rank mutants of a given structure. The performance of the SVM models was assessed on four benchmark datasets, which include protein–protein and protein–peptide binding affinity data. In addition, we evaluated our model on a set of mutants of EPI-X4, an endogenous peptide inhibitor of the chemokine receptor CXCR4, and on complexes of the serine proteases HTRA1 and HTRA3 with peptides. PPI-Affinity is freely accessible at https://protdcal.zmb.uni-due.de/PPIAffinity.

Conserved hydrogen-bond motifs of membrane transporters and receptors

Membrane transporters and receptors often rely on conserved hydrogen bonds to assemble transient paths for ion transfer or long-distance conformational couplings. For transporters and receptors that use proton binding and proton transfer for function, inter-helical hydrogen bonds of titratable protein sidechains that could change protonation are of central interest to formulate hypotheses about reaction mechanisms. Knowledge of hydrogen bonds common at sites of potential interest for proton binding could thus inform and guide studies on functional mechanisms of protonation-coupled membrane proteins. Here we apply graph-theory approaches to identify hydrogen-bond motifs of carboxylate and histidine sidechains in a large data set of static membrane protein structures. We find that carboxylate-hydroxyl hydrogen bonds are present in numerous structures of the dataset, and can be part of more extended H-bond clusters that could be relevant to conformational coupling. Carboxylate-carboxyamide and imidazole-imidazole hydrogen bonds are represented in comparably fewer protein structures of the dataset. Atomistic simulations on two membrane transporters in lipid membranes suggest that many of the hydrogen bond motifs present in static protein structures tend to be robust, and can be part of larger hydrogen-bond clusters that recruit additional hydrogen bonds.

Identifying vaccine escape sites via statistical comparisons of short-term molecular dynamics

The identification of viral mutations that confer escape from antibodies is crucial for understanding the interplay between immunity and viral evolution. We describe a molecular dynamics (MD)-based approach that goes beyond contact mapping, scales well to a desktop computer with a modern graphics processor, and enables the user to identify functional protein sites that are prone to vaccine escape in a viral antigen. We first implement our MD pipeline to employ site-wise calculation of Kullback-Leibler divergence in atom fluctuation over replicate sets of short-term MD production runs thus enabling a statistical comparison of the rapid motion of influenza hemagglutinin (HA) in both the presence and absence of three well-known neutralizing antibodies. Using this simple comparative method applied to motions of viral proteins, we successfully identified in silico all previously empirically confirmed sites of escape in influenza HA, predetermined via selection experiments and neutralization assays. Upon the validation of our computational approach, we then surveyed potential hotspot residues in the receptor binding domain of the SARS-CoV-2 virus in the presence of COVOX-222 and S2H97 antibodies. We identified many single sites in the antigen-antibody interface that are similarly prone to potential antibody escape and that match many of the known sites of mutations arising in the SARS-CoV-2 variants of concern. In the Omicron variant, we find only minimal adaptive evolutionary shifts in the functional binding profiles of both antibodies. In summary, we provide an inexpensive and accurate computational method to monitor hotspots of functional evolution in antibody binding footprints.

Subtype-specific responses of hKv7.4 and hKv7.5 channels to polyunsaturated fatty acids reveal an unconventional modulatory site and mechanism

The K_V7.4 and K_V7.5 subtypes of voltage-gated potassium channels play a role in important physiological processes such as sound amplification in the cochlea and adjusting vascular smooth muscle tone. Therefore, the mechanisms that regulate K_V7.4 and K_V7.5 channel function are of interest. Here, we study the effect of polyunsaturated fatty acids (PUFAs) on human K_V7.4 and K_V7.5 channels expressed in Xenopus oocytes. We report that PUFAs facilitate activation of hK_V7.5 by shifting the V₅₀ of the conductance versus voltage (G(V)) curve toward more negative voltages. This response depends on the head group charge, as an uncharged PUFA analogue has no effect and a positively charged PUFA analogue induces positive V₅₀ shifts. In contrast, PUFAs inhibit activation of hK_V7.4 by shifting V₅₀ toward more positive voltages. No effect on V₅₀ of hK_V7.4 is observed by an uncharged or a positively charged PUFA analogue. Thus, the hK_V7.5 channel's response to PUFAs is analogous to the one previously observed in hK_V7.1-7.3 channels, whereas the hK_V7.4 channel response is opposite, revealing subtype-specific responses to PUFAs. We identify a unique inner PUFA interaction site in the voltage-sensing domain of hK_V7.4 underlying the PUFA response, revealing an unconventional mechanism of modulation of hK_V7.4 by PUFAs.

Inhibiting a dynamic viral protease by targeting a non-catalytic cysteine

Viruses are responsible for some of the most deadly human diseases, yet available vaccines and antivirals address only a fraction of the potential viral human pathogens. Here, we provide a methodology for managing human herpesvirus (HHV) infection by covalently inactivating the HHV maturational protease via a conserved, non-catalytic cysteine (C161). Using human cytomegalovirus protease (HCMV Pr) as a model, we screened a library of disulfides to identify molecules that tether to C161 and inhibit proteolysis, then elaborated hits into irreversible HCMV Pr inhibitors that exhibit broad-spectrum inhibition of other HHV Pr homologs. We further developed an optimized tool compound targeted toward HCMV Pr and used an integrative structural biology and biochemical approach to demonstrate inhibitor stabilization of HCMV Pr homodimerization, exploiting a conformational equilibrium to block proteolysis. Irreversible HCMV Pr inhibition disrupts HCMV infectivity in cells, providing proof of principle for targeting proteolysis via a non-catalytic cysteine to manage viral infection.

Seipin transmembrane segments critically function in triglyceride nucleation and lipid droplet budding from the membrane

Lipid droplets (LDs) are organelles formed in the endoplasmic reticulum (ER) to store triacylglycerol (TG) and sterol esters. The ER protein seipin is key for LD biogenesis. Seipin forms a cage-like structure, with each seipin monomer containing a conserved hydrophobic helix and two transmembrane (TM) segments. How the different parts of seipin function in TG nucleation and LD budding is poorly understood. Here, we utilized molecular dynamics simulations of human seipin, along with cell-based experiments, to study seipin’s functions in protein–lipid interactions, lipid diffusion, and LD maturation. An all-atom simulation indicates that seipin TM segment residues and hydrophobic helices residues located in the phospholipid tail region of the bilayer attract TG. Simulating larger, growing LDs with coarse-grained models, we find that the seipin TM segments form a constricted neck structure to facilitate conversion of a flat oil lens into a budding LD. Using cell experiments and simulations, we also show that conserved, positively charged residues at the end of seipin’s TM segments affect LD maturation. We propose a model in which seipin TM segments critically function in TG nucleation and LD growth.

6'-O-Caffeoylarbutin from Que Zui tea ameliorates acetaminophen-induced liver injury via enhancing antioxidant ability and regulating the PI3K signaling pathway

Que Zui tea (QT), a traditional herbal tea in China, has a significant hepatoprotective effect. 6'-O-Caffeoylarbutin (CA) is the most abundant chemical compound in the QT. However, the hepatoprotective effect of CA has not been investigated. This study is aimed to evaluate the protective effect of CA on acetaminophen (APAP) induced hepatotoxicity in vivo and in vitro and its possible underlying mechanism. In APAP-induced HepG-2 cells, CA inhibited intracellular ROS accumulation and cell apoptosis, and improved the expression of antioxidants including SOD, CAT and GSH. In APAP-administrated mice, CA pretreatment remarkably ameliorated the histopathological damage and inflammatory response, and antioxidant enzyme activity in the serum and liver tissues. Moreover, the immunohistochemistry and immunofluorescence assay results revealed that the CA markedly reduced ROS production and apoptosis, and activated antioxidant transcription factor Nrf2 in the liver. Meanwhile, molecular docking results showed that the strong binding force of CA and PI3K was due to the higher number of hydrogen- and π-bonds with active site residues. Notably, CA pretreatment significantly regulated the expression of PI3K, Akt, Nrf2, NQO1, HO-1, Bcl-2, Bax, caspase-3, and caspase-9 proteins in APAP-treated liver tissues. These data demonstrated that CA had a protective effect against APAP-induced hepatotoxicity via regulating the PI3K/Akt and Nrf2 signaling pathway.

Protposer: the web server that readily proposes protein stabilizing mutations with high PPV

Protein stability is a requisite for most biotechnological and medical applications of proteins. As natural proteins tend to suffer from a low conformational stability ex vivo, great efforts have been devoted toward increasing their stability through rational design and engineering of appropriate mutations. Unfortunately, even the best currently used predictors fail to compute the stability of protein variants with sufficient accuracy and their usefulness as tools to guide the rational stabilisation of proteins is limited. We present here Protposer, a protein stabilising tool based on a different approach. Instead of quantifying changes in stability, Protposer uses structure- and sequence-based screening modules to nominate candidate mutations for subsequent evaluation by a logistic regression model, carefully trained to avoid overfitting. Thus, Protposer analyses PDB files in search for stabilization opportunities and provides a ranked list of promising mutations with their estimated success rates (eSR), their probabilities of being stabilising by at least 0.5 kcal/mol. The agreement between eSRs and actual positive predictive values (PPV) on external datasets of mutations is excellent. When Protposer is used with its Optimal kappa selection threshold, its PPV is above 0.7. Even with less stringent thresholds, Protposer largely outperforms FoldX, Rosetta and PoPMusiC. Indicating the PDB file of the protein suffices to obtain a ranked list of mutations, their eSRs and hints on the likely source of the stabilization expected. Protposer is a distinct, straightforward and highly successful tool to design protein stabilising mutations, and it is freely available for academic use at http://webapps.bifi.es/the-protposer.

Can ERAP1 and ERAP2 Form Functional Heterodimers? A Structural Dynamics Investigation

Endoplasmic reticulum aminopeptidases 1 and 2 (ERAP1 and ERAP2) play important roles in the generation of antigenic peptides presented by Major Histocompatibility Class I (MHCI) molecules and indirectly regulate adaptive immune responses. Although the discrete function of these enzymes has been extensively characterized, recent reports have suggested that they can also form heterodimers with functional consequences. However, lack of structural characterization of a putative ERAP1/ERAP2 dimer has limited our understanding of its biological role and significance. To address this, we employed computational molecular dynamics calculations to explore the topology of interactions between these two, based on experimentally determined homo-dimerization interfaces observed in crystal structures of ERAP2 or homologous enzymes. Our analysis of 8 possible dimerization models, suggested that the most likely ERAP1/ERAP2 heterodimerization topology involves the exon 10 loop, a non-conserved loop previously implicated in interactions between ERAP1 and the disulfide-bond shuffling chaperone ERp44. This dimerization topology allows access to the active site of both enzymes and is consistent with a previously reported construct in which ERAP1 and ERAP2 were linked by Fos/Jun zipper tags. The proposed model constitutes a tentative structural template to help understand the physiological role and significance of ERAP1/ERAP2 molecular interactions.

Novel angiotensin-converting enzyme (ACE) inhibitory mechanism of peptides from Macadamia integrifolia antimicrobial protein 2 (MiAMP2)

This work aimed to identify novel angiotensin-converting-enzyme (ACE) inhibitory peptides from Macadamia integrifolia antimicrobial protein 2 (MiAMP2). The MiAMP2 protein was hydrolyzed through in silico digestion, and the generated peptides were screened for ACE inhibitory activity. The in silico enzyme digestion results revealed that 18 unreported peptides were obtained using AHTPDB and BIOPEP-UWM, and none were thought to be toxic based on absorption, distribution, metabolism, and excretion (ADMET) prediction. PGPR, RPLY, MNPQR, and AAPR were predicted to exhibit good biological activity. The molecular docking results revealed that the four peptides tightly bound to the active pocket of ACE via hydrogen bonds and hydrophobic interactions, among which RPLY and MNPQR bound to ACE more strongly. The in vitro assay results confirmed that RPLY and MNPQR peptides inhibited ACE via competitive manner. These results provide theoretical guidance for the development of novel foodborne antihypertensive peptides from Macadamia nut proteins. PRACTICAL APPLICATIONS: This study provides new insight on the inhibitory potential of Macadamia nut peptides against ACE, which may be further applied to the development of antihypertensive peptides in the medical industry.

Aryl Hydrocarbon Receptor-Dependent and -Independent Pathways Mediate Curcumin Anti-Aging Effects

The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor whose activity can be modulated by polyphenols, such as curcumin. AhR and curcumin have evolutionarily conserved effects on aging. Here, we investigated whether and how the AhR mediates the anti-aging effects of curcumin across species. Using a combination of in vivo, in vitro, and in silico analyses, we demonstrated that curcumin has AhR-dependent or -independent effects in a context-specific manner. We found that in Caenorhabditis elegans, AhR mediates curcumin-induced lifespan extension, most likely through a ligand-independent inhibitory mechanism related to its antioxidant activity. Curcumin also showed AhR-independent anti-aging activities, such as protection against aggregation-prone proteins and oxidative stress in C. elegans and promotion of the migratory capacity of human primary endothelial cells. These AhR-independent effects are largely mediated by the Nrf2/SKN-1 pathway.

Bending-torsional elasticity and energetics of the plus-end microtubule tip

SignificanceThe mechanochemical basis of microtubule growth, which is essential for the normal function and division of eukaryotic cells, has remained elusive and controversial, despite extensive work. In particular, recent findings have created the paradox that the microtubule plus-end tips look very similar during both growing and shrinking phases, thereby challenging the traditional textbook picture. Our large-scale atomistic simulations resolve this paradox and explain microtubule growth and shrinkage dynamics as a process governed by energy barriers between protofilament conformations, the heights of which are in turn fine-tuned by different nucleotide states, thus implementing an information-driven Brownian ratchet.

Crystallography-based exploration of non-covalent interactions for the design and synthesis of coumarin for stronger protein binding

Protein molecules are a good target for the inhibition or promotion of biological processes. Different methods like QSAR and molecular docking have been developed to accurately design small binder molecules for target proteins. An alternative model has been developed wherein a statistical method is used to find the propensity of different non-covalent interactions between small molecules and amino acid residues of the protein. The results give hints as to the choice of substituents required at the SM to strongly bind to a protein. In this case, 75 different types of proteins bound with coumarin derivatives have been investigated and the non-covalent interactions observed between the basic coumarin moiety and amino acids have been analyzed. Density functional theory (DFT) calculations were used to identify the electronic features of coumarin to understand the feasibility of the observed non-covalent interactions and to find appropriate groups that can modulate these interactions. The binding affinity towards a protein (β-lactoglobulin (BLG)) and the stability of the protein complex have been investigated through docking and molecular dynamics of 100 ns, respectively. The modeled compounds were synthesized and investigated with regards to their interactions with the model carrier protein. The thermodynamics of the interactions were also investigated and the binding is governed by the Le Chatelier principle.

Structure of the anthrax protective antigen D425A dominant negative mutant reveals a stalled intermediate state of pore maturation

The tripartite protein complex produced by anthrax bacteria (Bacillus anthracis) is a member of the AB family of β-barrel pore-forming toxins. The protective antigen (PA) component forms an oligomeric prepore that assembles on the host cell surface and serves as a scaffold for binding of lethal and edema factors. Following endocytosis, the acidic environment of the late endosome triggers a pH-induced conformational rearrangement to promote maturation of the PA prepore to a functional, membrane spanning pore that facilitates delivery of lethal and edema factors to the cytosol of the infected host. Here, we show that the dominant-negative D425A mutant of PA stalls anthrax pore maturation in an intermediate state at acidic pH. Our 2.7 Å cryo-EM structure of the intermediate state reveals structural rearrangements that involve constriction of the oligomeric pore combined with an intramolecular dissociation of the pore-forming module. In addition to defining the early stages of anthrax pore maturation, the structure identifies asymmetric conformational changes in the oligomeric pore that are influenced by the precise configuration of adjacent protomers.

Binding Affinity Characterization of Four Antennae-Enriched Odorant-Binding Proteins From Harmonia axyridis (Coleoptera: Coccinellidae)

Harmonia axyridis is an important natural enemy that consumes many agricultural and forestry pests. It relies on a sensitive olfactory system to find prey and mates. Odorant-binding proteins (OBPs) as the first-step of recognizing volatiles, transport odors through sensillum lymph to odorant receptors (ORs). However, little is known about the molecular mechanisms of H. axyridis olfaction. In this study, four H. axyridis antenna specific OBP genes, HaxyOBP3, 5, 12, and 15, were bacterially expressed and the binding features of the four recombinant proteins to 40 substances were investigated using fluorescence competitive binding assays. Three-dimensional structure modeling and molecular docking analysis predicted the binding sites between HaxyOBPs and candidate volatiles. Developmental expression analyses showed that the four HaxyOBP genes displayed a variety of expression patterns at different development stages. The expression levels of HaxyOBP3 and HaxyOBP15 were higher in the adult stage than in the other developmental stages, and HaxyOBP15 was significantly transcriptionally enriched in adult stage. Ligand-binding analysis demonstrated that HaxyOBP3 and HaxyOBP12 only combined with two compounds, β-ionone and p-anisaldehyde. HaxyOBP5 protein displayed binding affinities with methyl salicylate, β-ionone, and p-anisaldehyde (Ki = 18.15, 11.71, and 13.45 μM). HaxyOBP15 protein had a broad binding profile with (E)-β-farnesene, β-ionone, α-ionone, geranyl acetate, nonyl aldehyde, dihydro-β-ionone, and linalyl acetate (Ki = 4.33-31.01 μM), and hydrophobic interactions played a key role in the binding of HaxyOBP15 to these substances according to molecular docking. Taken together, HaxyOBP15 exhibited a broader ligand-binding spectrum and a higher expression in adult stage than HaxyOBP3, 5, and 12, indicating HaxyOBP15 may play a greater role in binding volatiles than other three HaxyOBPs. The results will increase our understanding of the molecular mechanism of H. axyridis olfaction and may also result in new management strategies (attractants/repellents) that increase the biological control efficacy of H. axyridis.

An in-frame deletion mutation in the degron tail of auxin coreceptor IAA2 confers resistance to the herbicide 2,4-D in Sisymbrium orientale

Synthetic auxin herbicides intersect basic plant developmental biology and applied weed management. We investigated resistance to 2,4-D in the Australian weed Sisymbrium orientale (Indian hedge mustard). We identified a mechanism involving an in-frame 27-bp deletion in the degron tail of auxin coreceptor IAA2, one member of the gene family of Aux/IAA auxin co-receptors. We show that this deletion in IAA2 is a gain-of-function mutation that confers synthetic auxin resistance. This field-evolved mechanism of resistance to synthetic auxin herbicides confirms previous biochemical studies showing the role of the Aux/IAA degron tail in regulating Aux/IAA protein degradation upon auxin perception. The deletion mutation could be generated in crops using gene-editing approaches for cross-resistance to multiple synthetic auxin herbicides.

The natural auxin indole-3-acetic acid (IAA) is a key regulator of many aspects of plant growth and development. Synthetic auxin herbicides such as 2,4-D mimic the effects of IAA by inducing strong auxinic-signaling responses in plants. To determine the mechanism of 2,4-D resistance in a Sisymbrium orientale (Indian hedge mustard) weed population, we performed a transcriptome analysis of 2,4-D-resistant (R) and -susceptible (S) genotypes that revealed an in-frame 27-nucleotide deletion removing nine amino acids in the degron tail (DT) of the auxin coreceptor Aux/IAA2 (SoIAA2). The deletion allele cosegregated with 2,4-D resistance in recombinant inbred lines. Further, this deletion was also detected in several 2,4-D-resistant field populations of this species. Arabidopsis transgenic lines expressing the SoIAA2 mutant allele were resistant to 2,4-D and dicamba. The IAA2-DT deletion reduced binding to TIR1 in vitro with both natural and synthetic auxins, causing reduced association and increased dissociation rates. This mechanism of synthetic auxin herbicide resistance assigns an in planta function to the DT region of this Aux/IAA coreceptor for its role in synthetic auxin binding kinetics and reveals a potential biotechnological approach to produce synthetic auxin-resistant crops using gene-editing.

Molecular modeling of natural and synthesized inhibitors against SARS-CoV-2 spike glycoprotein

Purpose

Viral diseases increasingly endanger the world public health because of the transient efficacy of antiviral therapies. The novel coronavirus disease 2019 (COVID-19) has been recently identified as caused by a new type of coronaviruses. This type of coronavirus binds to the human receptor through the Spike glycoprotein (S) Receptor Binding Domain (RBD). The spike protein is found in inaccessible (closed) or accessible (open) conformations in which the accessible conformation causes severe infection. Thus, this receptor is a significant target for antiviral drug design.

Methods

An attempt was made to recognize 111 natural and synthesized compounds in order to utilize them against SARS-CoV-2 spike glycoprotein to inhibit Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) using simulation methods, such as molecular docking. The FAF-Drugs3, Pan-Assay Interference Compounds (PAINS), ADME (absorption, distribution, metabolism, excretion) databases along with Lipinski’s rules were used to evaluate the drug-like properties of the identified ligands. In order to analyze and identify the residues critical in the docking process of the spike glycoprotein, the interactions of proposed ligands with both conformations of the spike glycoprotein was simulated.

Results

The results showed that among the available ligands, seven ligands had significant interactions with the binding site of the spike glycoprotein, in which angiotensin-converting enzyme 2 (ACE2) is bounded. Out of seven candidate molecules, six ligands exhibited drug-like characteristics. The results also demonstrated that fluorophenyl and propane groups of ligands had optimal interactions with the binding site of the spike glycoprotein.

Conclusion

According to the results, our findings indicated the ability of six ligands to prevent the binding of the SARS-CoV-2 spike glycoprotein to its cognate receptor, providing novel compounds for the treatment of COVID-19.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42600-020-00122-3.

Lactococcus lactis secreting phage lysins as a potential antimicrobial against multi-drug resistant Staphylococcus aureus

BACKGROUND

Staphylococcus aureus is an opportunistic Gram-positive bacterium that can form biofilm and become resistant to many types of antibiotics. The treatment of multi-drug resistant Staphylococcus aureus (MDRSA) infection is difficult since it possesses multiple antibiotic-resistant mechanisms. Endolysin and virion-associated peptidoglycan hydrolases (VAPGH) enzymes from bacteriophage have been identified as potential alternative antimicrobial agents. This study aimed to assess the ability of Lactococcus lactis NZ9000 secreting endolysin and VAPGH from S. aureus bacteriophage 88 to inhibit the growth of S. aureus PS 88, a MDRSA.

METHOD

Endolysin and VAPGH genes were cloned and expressed in L. lactis NZ9000 after fusion with the SPK1 signal peptide for secretion. The recombinant proteins were expressed and purified, then analyzed for antimicrobial activity using plate assay and turbidity reduction assay. In addition, the spent media of the recombinant lactococcal culture was analyzed for its ability to inhibit the growth of the S. aureus PS 88.

RESULTS

Extracellular recombinant endolysin (Endo88) and VAPGH (VAH88) was successfully expressed and secreted from L. lactis which was able to inhibit S. aureus PS 88, as shown by halozone formation on plate assays as well as inhibition of growth in the turbidity reduction assay. Moreover, it was observed that the spent media from L. lactis NZ9000 expressing Endo88 and VAH88 reduced the viability of PS 88 by up to 3.5-log reduction with Endo88 being more efficacious than VAH88. In addition, Endo88 was able to lyse all MRSA strains tested and Staphylococcus epidermidis but not the other bacteria while VAH88 could only lyse S. aureus PS 88.

CONCLUSION

Recombinant L. lactisNZ9000 expressing phage 88 endolysin may be potentially developed into a new antimicrobial agent for the treatment of MDRSA infection.

Human variation in the protein receptor ACE2 affects its binding affinity to SARS-CoV-2 in a variant-dependent manner

SARS-CoV-2 infection depend on the binding of the viral Spike glycoprotein (S) to the human receptor Angiotensin Converting Enzyme 2 (ACE2) to induce virus-cell membrane fusion. S protein evolved diverse amino acid changes that are possibly linked to more efficient binding to human ACE2, which might explain part of the increase in frequency of SARS-CoV-2 Variants Of Concern (VOCs). In this work, we investigated the role of ACE2 protein variations that are naturally found in human populations and its binding affinity with S protein from SARS-CoV-2 representative genotypes, based on a series of in silico approaches involving molecular modelling, docking and molecular dynamics simulations. Our results indicate that SARS-CoV-2 VOCs bind more efficiently to the human receptor ACE2 than the ancestral Wuhan genotype. Additionally, variations in the ACE2 protein can affect SARS-CoV-2 binding and protein-protein stability, mostly making the interaction weaker and unstable in some cases. We show that some VOCs, such as B.1.1.7 and P.1 are much less sensitive to ACE2 variants, while others like B.1.351 appear to be specifically optimized to bind to the widespread wild-type ACE2 protein.Communicated by Ramaswamy H. Sarma.

Discovery of Dual TAF1-ATR Inhibitors and Ligand-Induced Structural Changes of the TAF1 Tandem Bromodomain

Bromodomains regulate chromatin remodeling and gene transcription through recognition of acetylated lysines on histones and other proteins. Bromodomain-containing protein TAF1, a subunit of general transcription factor TFIID, initiates preinitiation complex formation and cellular transcription. TAF1 serves as a cofactor for certain oncogenic transcription factors and is implicated in regulating the p53 tumor suppressor. Therefore, TAF1 is a potential target to develop small molecule therapeutics for diseases arising from dysregulated transcription, such as cancer. Here, we report the ATR kinase inhibitor AZD6738 (Ceralasertib) and analogues thereof as bona fide inhibitors of TAF1. Crystallographic and small-angle X-ray scattering studies established that newly identified and previously reported inhibitors stabilize distinct structural states of the TAF1 tandem bromodomain through "open-closed" transitions and dimerization. Combined with functional studies on p53 signaling in cancer cell lines, the data provide new insights into the feasibility and challenges of TAF1 inhibitors as chemical probes and therapeutics.

Activation Mechanism of Strigolactone Receptors and Its Impact on Ligand Selectivity between Host and Parasitic Plants

Parasitic weeds such as Striga have led to significant losses in agricultural productivity worldwide. These weeds use the plant hormone strigolactone as a germination stimulant. Strigolactone signaling involves substrate hydrolysis followed by a conformational change of the receptor to a "closed" or "active" state that associates with a signaling partner, MAX2/D3. Crystal structures of active and inactive AtD14 receptors have helped elucidate the structural changes involved in activation. However, the mechanism by which the receptor activates remains unknown. The ligand dependence of AtD14 activation has been disputed by mutagenesis studies showing that enzymatically inactive receptors are able to associate with MAX2 proteins. Furthermore, activation differences between strigolactone receptor in Striga, ShHTL7, and AtD14 could contribute to the high sensitivity to strigolactones exhibited by parasitic plants. Using molecular dynamics simulations, we demonstrate that both AtD14 and ShHTL7 could adopt an active conformation in the absence of ligand. However, ShHTL7 exhibits a higher population in the inactive apo state as compared to the AtD14 receptor. We demonstrate that this difference in inactive state population is caused by sequence differences between their D-loops and interactions with the catalytic histidine that prevent full binding pocket closure in ShHTL7. These results indicate that ligand hydrolysis would enhance the active state population by destabilizing the inactive state in ShHTL7 as compared to AtD14. We also show that the mechanism of activation is more concerted in AtD14 than in ShHTL7 and that the main barrier to activation in ShHTL7 is closing of the binding pocket.

Characterization of the immune resistance of SARS-CoV-2 Mu variant and the robust immunity induced by Mu infection

We have recently revealed that the SARS-CoV-2 Mu variant shows a pronounced resistance to antibodies elicited by natural SARS-CoV-2 infection and vaccination. However, it remains unclear which mutations determine the resistance of SARS-CoV-2 Mu to antiviral sera. Also, it is unclear how SARS-CoV-2 Mu infection induces antiviral immunity. Here we reveal that the two mutations in the SARS-CoV-2 Mu spike protein, YY144-145TSN and E484K, are responsible for the resistance to COVID-19 convalescent sera during early 2020 and vaccine sera. Notably, the convalescent sera of SARS-CoV-2 Mu-infected individuals are broadly antiviral against Mu as well as other SARS-CoV-2 variants of concern/interest.

The Full Model of the pMHC-TCR-CD3 Complex: A Structural and Dynamical Characterization of Bound and Unbound States

The machinery involved in cytotoxic T-cell activation requires three main characters: the major histocompatibility complex class I (MHC I) bound to the peptide (p), the T-cell receptor (TCR), and the CD3 complex, a multidimer interfaced with the intracellular side. The pMHC:TCR interaction has been largely studied by means of both experimental and computational models, giving a contribution in understanding the complexity of the TCR triggering. Nevertheless, a detailed study of the structural and dynamical characterization of the full complex (pMHC:TCR:CD3 complex) is still missing due to a lack of structural information of the CD3-chains arrangement around the TCR. Very recently, the determination of the TCR:CD3 complex structure by means of Cryo-EM technique has given a chance to build the entire system essential in the activation of T-cells, a fundamental mechanism in the adaptive immune response. Here, we present the first complete model of the pMHC interacting with the TCR:CD3 complex, built in a lipid environment. To describe the conformational behavior associated with the unbound and the bound states, all-atom Molecular Dynamics simulations were performed for the TCR:CD3 complex and for two pMHC:TCR:CD3 complex systems, bound to two different peptides. Our data point out that a conformational change affecting the TCR Constant β (Cβ) region occurs after the binding to the pMHC, revealing a key role of this region in the propagation of the signal. Moreover, we found that TCR reduces the flexibility of the MHC I binding groove, confirming our previous results.

In Silico Molecular Characterization of Human TMPRSS2 Protease Polymorphic Variants and Associated SARS-CoV-2 Susceptibility

The 2019 coronavirus disease (COVID-19) pandemic continues to challenge health care systems worldwide. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been responsible for the cause of global pandemic. Type 2 transmembrane serine protease (TMPRSS2) is important in the cell entry and spread of SARS-CoV-2 and plays a crucial role in the proteolytic cleavage of SARS-CoV-2 spike (S) glycoprotein. Here, using reported structural data, we analyzed the molecular complex of TMPRSS2 and the S glycoprotein and further examined intermolecular interactions of natural TMPRSS2 polymorphic variants. We identified several TMPRSS2 variants that could possibly alter host susceptibility to the SARS-CoV-2 infection. Molecular docking analysis revealed that G462D/G462S variants were predicted to be protective variants, whereas Q438E and S339F variants were predicted to increase susceptibility. In addition, we examined intermolecular interactions between TMPRSS2 and its two potential serine protease inhibitors, camostat mesylate and nafamostat. Further, we investigated the effect of TMPRSS2 variants on these interactions. Our structural analysis revealed that G462D, C297S and S460R variants had possibly altered the interactions with the protease inhibitors. Our results identified important TMPRSS2 variations that could be useful to develop high affinity and personalized drugs for treating COVID-19 patients.

Slipknot or Crystallographic Error: A Computational Analysis of the Plasmodium falciparum DHFR Structural Folds

The presence of protein structures with atypical folds in the Protein Data Bank (PDB) is rare and may result from naturally occurring knots or crystallographic errors. Proper characterisation of such folds is imperative to understanding the basis of naturally existing knots and correcting crystallographic errors. If left uncorrected, such errors can frustrate downstream experiments that depend on the structures containing them. An atypical fold has been identified in P. falciparum dihydrofolate reductase (PfDHFR) between residues 20–51 (loop 1) and residues 191–205 (loop 2). This enzyme is key to drug discovery efforts in the parasite, necessitating a thorough characterisation of these folds. Using multiple sequence alignments (MSA), a unique insert was identified in loop 1 that exacerbates the appearance of the atypical fold-giving it a slipknot-like topology. However, PfDHFR has not been deposited in the knotted proteins database, and processing its structure failed to identify any knots within its folds. The application of protein homology modelling and molecular dynamics simulations on the DHFR domain of P. falciparum and those of two other organisms (E. coli and M. tuberculosis) that were used as molecular replacement templates in solving the PfDHFR structure revealed plausible unentangled or open conformations of these loops. These results will serve as guides for crystallographic experiments to provide further insights into the atypical folds identified.

The design of TOPK inhibitors using structure-based pharmacophore modeling and molecular docking based on an MD-refined homology model

The TOPK enzyme (also known as PBK) is a serine-threonine protein kinase that is rarely detected in normal tissues yet is found to be overexpressed and activated in a variety of cancers such as lung, colorectal, breast, and esophageal cancer. Its prevalence in cancerous cells is associated with their poor prognosis and responsiveness to treatment. This enzyme plays a vital role in cell division, specifically in regulating cytokinesis. Unlike drugs targeting early phases in mitosis, inhibition of cytokinesis by targeting biomolecules that are unique to multiplying cells poses no threat to the normal function of non-multiplying cells. Studies have shown that inhibition of cytokinesis is promising in suppressing the growth of proliferating cancerous cells as exemplified by the complete tumor regression seen with the suppression of TOPK. Herein, we report the identification of potent TOPK inhibitors with anticancer potential using a structure-based drug design approach. The only available crystal structure of TOPK corresponds to a double mutant (T9E and T198E) dimer with a distorted N-lobe conformation, thus 3D homology modeling was implemented to rebuild the enzyme's native conformation. The resulting refined model was used to generate 3D pharmacophore models for the virtual screening of small molecules databases. Retrieved hits were filtered, docked into the ATP binding site of the enzyme, rescored, and the binding free energies for the top consensually scoring hits were calculated. Consequently, 45 compounds were selected and their in vitro inhibitory activity against TOPK was tested revealing four potential hits with the most active compound having an IC₅₀ of 3.85 µM. This compound will be chosen as a lead compound to synthesize analogs aiming to enhance potency and drug-like properties and to enrich the SAR data.

The SARS-CoV-2 Lambda variant exhibits enhanced infectivity and immune resistance

SARS-CoV-2 Lambda, a variant of interest, has spread in some South American countries; however, its virological features and evolutionary traits remain unclear. In this study, we use pseudoviruses and reveal that the spike protein of the Lambda variant is more infectious than that of other variants due to the T76I and L452Q mutations. The RSYLTPGD246-253N mutation, a unique 7-amino acid deletion in the N-terminal domain of the Lambda spike protein, is responsible for evasion from neutralizing antibodies and further augments antibody-mediated enhancement of infection. Although this mutation generates a nascent N-linked glycosylation site, the additional N-linked glycan is dispensable for the virological property conferred by this mutation. Since the Lambda variant has dominantly spread according to the increasing frequency of the isolates harboring the RSYLTPGD246-253N mutation, our data suggest that the RSYLTPGD246-253N mutation is closely associated with the substantial spread of the Lambda variant in South America.

PEA-15 engages in allosteric interactions using a common scaffold in a phosphorylation-dependent manner

Phosphoprotein enriched in astrocytes, 15 kDa (PEA-15) is a death-effector domain (DED) containing protein involved in regulating mitogen-activated protein kinase and apoptosis pathways. In this molecular dynamics study, we examined how phosphorylation of the PEA-15 C-terminal tail residues, Ser-104 and Ser-116, allosterically mediates conformational changes of the DED and alters the binding specificity from extracellular-regulated kinase (ERK) to Fas-associated death domain (FADD) protein. We delineated that the binding interfaces between the unphosphorylated PEA-15 and ERK2 and between the doubly phosphorylated PEA-15 and FADD are similarly composed of a scaffold that includes both the DED and the C-terminal tail residues of PEA-15. While the unphosphorylated serine residues do not directly interact with ERK2, the phosphorylated Ser-116 engages in strong electrostatic interactions with arginine residues on FADD DED. Upon PEA-15 binding, FADD repositions its death domain (DD) relative to the DED, an essential conformational change to allow the death-inducing signaling complex (DISC) assembly.

Structure of the ancient TRPY1 channel from Saccharomyces cerevisiae reveals mechanisms of modulation by lipids and calcium

Transient receptor potential (TRP) channels emerged in fungi as mechanosensitive osmoregulators. The Saccharomyces cerevisiae vacuolar TRP yeast 1 (TRPY1) is the most studied TRP channel from fungi, but the structure and details of channel modulation remain elusive. Here, we describe the full-length cryoelectron microscopy structure of TRPY1 at 3.1 Å resolution in a closed state. The structure, despite containing an evolutionarily conserved and archetypical transmembrane domain, reveals distinctive structural folds for the cytosolic N and C termini, compared with other eukaryotic TRP channels. We identify an inhibitory phosphatidylinositol 3-phosphate (PI(3)P) lipid-binding site, along with two Ca2+-binding sites: a cytosolic site, implicated in channel activation and a vacuolar lumen site, implicated in inhibition. These findings, together with data from microsecond-long molecular dynamics simulations and a model of a TRPY1 open state, provide insights into the basis of TRPY1 channel modulation by lipids and Ca2+, and the molecular evolution of TRP channels.

Modeling a unit cell: crystallographic refinement procedure using the biomolecular MD simulation platform Amber

A procedure has been developed for the refinement of crystallographic protein structures based on the biomolecular simulation program Amber. The procedure constructs a model representing a crystal unit cell, which generally contains multiple protein molecules and is fully hydrated with TIP3P water. Periodic boundary conditions are applied to the cell in order to emulate the crystal lattice. The refinement is conducted in the form of a specially designed short molecular-dynamics run controlled by the Amber ff14SB force field and the maximum-likelihood potential that encodes the structure-factor-based restraints. The new Amber-based refinement procedure has been tested on a set of 84 protein structures. In most cases, the new procedure led to appreciably lower R free values compared with those reported in the original PDB depositions or obtained by means of the industry-standard phenix.refine program. In particular, the new method has the edge in refining low-accuracy scrambled models. It has also been successful in refining a number of molecular-replacement models, including one with an r.m.s.d. of 2.15 Å. In addition, Amber-refined structures consistently show superior MolProbity scores. The new approach offers a highly realistic representation of protein-protein interactions in the crystal, as well as of protein-water interactions. It also offers a realistic representation of protein crystal dynamics (akin to ensemble-refinement schemes). Importantly, the method fully utilizes the information from the available diffraction data, while relying on state-of-the-art molecular-dynamics modeling to assist with those elements of the structure that do not diffract well (for example mobile loops or side chains). Finally, it should be noted that the protocol employs no tunable parameters, and the calculations can be conducted in a matter of several hours on desktop computers equipped with graphical processing units or using a designated web service.

In silico study on the effects of disulfide bonds in ORF8 of SARS-CoV-2

The COVID-19 epidemic, caused by virus SARS-CoV-2, has been a pandemic and threatening everyone's health in the past two years. In SARS-CoV-2, ORF8 is one of the most important accessory...

The Therapeutic Antibody Profiler for Computational Developability Assessment

The need to consider an antibody's "developability" (immunogenicity, solubility, specificity, stability, manufacturability, and storability) is now well understood in therapeutic antibody design. Predicting these properties rapidly and inexpensively is critical to industrial workflows, to avoid devoting resources to non-productive candidates. Here, we describe a high-throughput computational developability assessment tool, the Therapeutic Antibody Profiler (TAP), which assesses the physicochemical "druglikeness" of an antibody candidate. Input variable domain sequences are converted to three-dimensional structural models, and then five developability-linked molecular surface descriptors are calculated and compared to advanced-stage clinical therapeutics. Values at the extremes of/outside of the distributions seen in therapeutics imply an increased risk of developability issues. Therefore, TAP, starting only from sequence information, provides a route to rapidly identifying drug candidate antibodies that are likely to have poor developability. Our web application ( opig.stats.ox.ac.uk/webapps/tap ) profiles input antibody sequences against a continually updated reference set of clinical therapeutics.

Structural basis of sphingosine-1-phosphate receptor 1 activation and biased agonism

Sphingosine-1-phosphate receptor 1 (S1PR1) is a master regulator of lymphocyte egress from the lymph node and an established drug target for multiple sclerosis (MS). Mechanistically, therapeutic S1PR1 modulators activate the receptor yet induce sustained internalization through a potent association with β-arrestin. However, a structural basis of biased agonism remains elusive. Here, we report the cryo-electron microscopy (cryo-EM) structures of G_i-bound S1PR1 in complex with S1P, fingolimod-phosphate (FTY720-P) and siponimod (BAF312). In combination with functional assays and molecular dynamics (MD) studies, we reveal that the β-arrestin-biased ligands direct a distinct activation path in S1PR1 through the extensive interplay between the PIF and the NPxxY motifs. Specifically, the intermediate flipping of W269^6.48 and the retained interaction between F265^6.44 and N307^7.49 are the key features of the β-arrestin bias. We further identify ligand-receptor interactions accounting for the S1PR subtype specificity of BAF312. These structural insights provide a rational basis for designing novel signaling-biased S1PR modulators.

Free-Energy Simulations Support a Lipophilic Binding Route for Melatonin Receptors

The effects of the neurohormone melatonin are mediated by the activation of the GPCRs MT₁ and MT₂ in a variety of tissues. Crystal structures suggest ligand access to the orthosteric binding site of MT₁ and MT₂ receptors through a lateral channel between transmembrane (TM) helices IV and V. We investigated the feasibility of this lipophilic entry route for 2-iodomelatonin, a nonselective agonist with a slower dissociation rate from the MT₂ receptor, applying enhanced sampling simulations and free-energy calculations. 2-Iodomelatonin unbinding was investigated with steered molecular dynamics simulations which revealed different trajectories passing through the gap between TM helices IV and V for both receptors. For one of these unbinding trajectories from the MT₁ receptor, an umbrella-sampling protocol with path-collective variables provided a calculated energy barrier consistent with the experimental dissociation rate. The side-chain flexibility of Tyr5.38 was significantly different in the two receptor subtypes, as assessed by metadynamics simulations, and during ligand unbinding it frequently assumes an open conformation in the MT₁ but not in the MT₂ receptor, favoring 2-iodomelatonin egress. Taken together, our simulations are consistent with the possibility that the gap between TM IV and V is a way of connecting the orthosteric binding site and the membrane core for lipophilic melatonin receptor ligands. Our simulations also suggest that the open state of Tyr5.38 generates a small pocket on the surface of MT₁ receptor, which could participate in the recognition of MT₁-selective ligands and may be exploited in the design of new selective compounds.

Modeling the native ensemble of PhuS using enhanced sampling MD and HDX-ensemble reweighting

The cytoplasmic heme binding protein from Pseudomonas aeruginosa, PhuS, plays two essential roles in regulating heme uptake and iron homeostasis. First, PhuS shuttles exogenous heme to heme oxygenase (HemO) for degradation and iron release. Second, PhuS binds DNA and modulates the transcription of the prrF/H small RNAs (sRNAs) involved in the iron-sparing response. Heme binding to PhuS regulates this dual function, as the unliganded form binds DNA, whereas the heme-bound form binds HemO. Crystallographic studies revealed nearly identical structures for apo- and holo-PhuS, and yet numerous solution-based measurements indicate that heme binding is accompanied by large conformational rearrangements. In particular, hydrogen-deuterium exchange mass spectrometry (HDX-MS) of apo- versus holo-PhuS revealed large differences in deuterium uptake, notably in α-helices 6, 7, and 8 (α6,7,8), which contribute to the heme binding pocket. These helices were mostly labile in apo-PhuS but largely protected in holo-PhuS. In contrast, in silico-predicted deuterium uptake levels of α6,7,8 from molecular dynamics (MD) simulations of the apo- and holo-PhuS structures are highly similar, consistent only with the holo-PhuS HDX-MS data. To rationalize this discrepancy between crystal structures, simulations, and observed HDX-MS, we exploit a recently developed computational approach (HDXer) that fits the relative weights of conformational populations within an ensemble of structures to conform to a target set of HDX-MS data. Here, a combination of enhanced sampling MD, HDXer, and dimensionality reduction analysis reveals an apo-PhuS conformational landscape in which α6, 7, and 8 are significantly rearranged compared to the crystal structure, including a loss of secondary structure in α6 and the displacement of α7 toward the HemO binding interface. Circular dichroism analysis confirms the loss of secondary structure, and the extracted ensembles of apo-PhuS and of heme-transfer-impaired H212R mutant, are consistent with known heme binding and transfer properties. The proposed conformational landscape provides structural insights into the modulation by heme of the dual function of PhuS.

Highly synergistic combinations of nanobodies that target SARS-CoV-2 and are resistant to escape

The emergence of SARS-CoV-2 variants threatens current vaccines and therapeutic antibodies and urgently demands powerful new therapeutics that can resist viral escape. We therefore generated a large nanobody repertoire to saturate the distinct and highly conserved available epitope space of SARS-CoV-2 spike, including the S1 receptor binding domain, N-terminal domain, and the S2 subunit, to identify new nanobody binding sites that may reflect novel mechanisms of viral neutralization. Structural mapping and functional assays show that indeed these highly stable monovalent nanobodies potently inhibit SARS-CoV-2 infection, display numerous neutralization mechanisms, are effective against emerging variants of concern, and are resistant to mutational escape. Rational combinations of these nanobodies that bind to distinct sites within and between spike subunits exhibit extraordinary synergy and suggest multiple tailored therapeutic and prophylactic strategies.