3V: cavity, channel and cleft volume calculator and extractor

Non-steroidal inhibitors of Drosophila melanogaster steroidogenic glutathione S-transferase Noppera-bo

Insect growth regulators (IGRs) can be developed by elucidating the molecular mechanisms of insect-specific biological events. Because insect molting, and metamorphosis are controlled by ecdysteroids, their biosynthetic pathways can serve as targets for IGR development. The glutathione S-transferase Noppera-bo (Nobo), which is conserved in dipteran and lepidopteran species, plays an essential role in ecdysteroid biosynthesis. Our previous study using 17β-estradiol as a molecular probe revealed that Asp113 of Drosophila melanogaster Nobo (DmNobo) is essential for its biological function. However, to develop IGRs with a greater Nobo inhibitory activity than 17β-estradiol, further structural information is warranted. Here, we report five novel non-steroidal DmNobo inhibitors. Analysis of crystal structures of complexes revealed that DmNobo binds these inhibitors in an Asp113-independent manner. Among amino acid residues at the substrate-recognition site, conformation of conserved Phe39 was dynamically altered upon inhibitor binding. Therefore, these inhibitors can serve as seed compounds for IGR development.

A computational study on the interactions between a layered imine-based COF structure and selected anticancer drugs

The covalent organic frameworks (COFs) are important materials in drug delivery. Herein, the interactions between an imine-based COF with selected commercially available anticancer drugs are studied. Molecular dynamics (MD) simulation studies were used. The studies were carried out in four different temperatures to find out the impact of the temperature on the binding free energies between the drugs and COF structure. It was found that the effect of temperature on binding free energy is ignorable. Between the hydrogen bonding, electrostatic, and van der Waals interactions, the last one is the most important one to keep the drug and COF next to each other. Also, the van der Waals interaction is keeping the layers of COF next to each other to create cavities. The cavities can be loaded with different drugs and the system can be used in drug delivery systems. Based on the obtained results, the drugs that are more lipophilic prefer to adhere more strongly to the COF in comparison with hydrophilic drugs.

Structural Basis of WLS/Evi-Mediated Wnt Transport and Secretion

Wnts are evolutionarily conserved ligands that signal at short range to regulate morphogenesis, cell fate, and stem cell renewal. The first and essential steps in Wnt secretion are their O-palmitoleation and subsequent loading onto the dedicated transporter Wntless/evenness interrupted (WLS/Evi). We report the 3.2 Å resolution cryogenic electron microscopy (cryo-EM) structure of palmitoleated human WNT8A in complex with WLS. This is accompanied by biochemical experiments to probe the physiological implications of the observed association. The WLS membrane domain has close structural homology to G protein-coupled receptors (GPCRs). A Wnt hairpin inserts into a conserved hydrophobic cavity in the GPCR-like domain, and the palmitoleate protrudes between two helices into the bilayer. A conformational switch of highly conserved residues on a separate Wnt hairpin might contribute to its transfer to receiving cells. This work provides molecular-level insights into a central mechanism in animal body plan development and stem cell biology.

Prediction of protein-carbohydrate complex binding affinity using structural features

Protein-carbohydrate interactions play a major role in several cellular and biological processes. Elucidating the factors influencing the binding affinity of protein-carbohydrate complexes and predicting their free energy of binding provide deep insights for understanding the recognition mechanism. In this work, we have collected the experimental binding affinity data for a set of 389 protein-carbohydrate complexes and derived several structure-based features such as contact potentials, interaction energy, number of binding residues and contacts between different types of atoms. Our analysis on the relationship between binding affinity and structural features revealed that the important factors depend on the type of the complex based on number of carbohydrate and protein chains. Specifically, binding site residues, accessible surface area, interactions between various atoms and energy contributions are important to understand the binding affinity. Further, we have developed multiple regression equations for predicting the binding affinity of protein-carbohydrate complexes belonging to six categories of protein-carbohydrate complexes. Our method showed an average correlation and mean absolute error of 0.731 and 1.149 kcal/mol, respectively, between experimental and predicted binding affinities on a jackknife test. We have developed a web server PCA-Pred, Protein-Carbohydrate Affinity Predictor, for predicting the binding affinity of protein-carbohydrate complexes. The web server is freely accessible at https://web.iitm.ac.in/bioinfo2/pcapred/. The web server is implemented using HTML and Python and supports recent versions of major browsers such as Chrome, Firefox, IE10 and Opera.

A Molecular Capsule with Revolving Doors Partitioning Its Inner Space

Covalent capsule 1 was designed to include two molecular baskets linked with three mobile pyridines tucked into its inner space. On the basis of both theory (DFT) and experiments (NMR and X-ray crystallography), we found that the pyridine "doors" split the chamber (380 ų ) of 1 so that two equally sizeable compartments (190 ų ) became joined through a conformationally flexible aromatic barrier. The compartments of such unique host could be populated with CCl₄ (88 ų ; PC=46 %), CBr₄ (106 ų ; 56 %) or their combination CCl₄ /CBr₄ (PC=51 %), with thermodynamic stabilities ΔG° tracking the values of packing coefficients (PC). Halogen (C-X⋅⋅⋅π) and hydrogen bonding (C-H⋅⋅⋅X) contacts held the haloalkane guests in the cavities of 1. The consecutive complexations were found to occur in a negative allosteric manner, which we propose to result from the induced-fit mode of complexation. Newly designed 1 opens a way for probing the effects of inner conformational dynamics on noncovalent interactions, reactivity and intramolecular translation in confined spaces of hollow molecules.

Human trans-editing enzyme displays tRNA acceptor-stem specificity and relaxed amino acid selectivity

Accurate translation of genetic information into proteins is vital for cell sustainability. ProXp-ala prevents proteome-wide Pro-to-Ala mutations by hydrolyzing misacylated Ala-tRNAPro, which is synthesized by prolyl-tRNA synthetase. Bacterial ProXp-ala was previously shown to combine a size-based exclusion mechanism with conformational and chemical selection for the recognition of the alanyl moiety, whereas tRNAPro is selected via recognition of tRNA acceptor-stem elements G72 and A73. The identity of these critical bases changed during evolution with eukaryotic cytosolic tRNAPro possessing a cytosine at the corresponding positions. The mechanism by which eukaryotic ProXp-ala adapted to these changes remains unknown. In this work, recognition of the aminoacyl moiety and tRNA acceptor stem by human (Homo sapiens, or Hs) ProXp-ala was examined. Enzymatic assays revealed that Hs ProXp-ala requires C72 and C73 in the context of Hs cytosolic tRNAPro for efficient deacylation of mischarged Ala-tRNAPro The strong dependence on these bases prevents cross-species deacylation of bacterial Ala-tRNAPro or of Hs mitochondrial Ala-tRNAPro by the human enzyme. Similar to the bacterial enzyme, Hs ProXp-ala showed strong tRNA acceptor-stem recognition but differed in its amino acid specificity profile relative to bacterial ProXp-ala. Changes at conserved residues in both the Hs and bacterial ProXp-ala substrate-binding pockets modulated this specificity. These results illustrate how the mechanism of substrate selection diverged during the evolution of the ProXp-ala family, providing the first example of a trans-editing domain whose specificity evolved to adapt to changes in its tRNA substrate.

Dissecting Dynamic and Hydration Contributions to Sequence-Dependent DNA Minor Groove Recognition

Sequence selectivity is a critical attribute of DNA-binding ligands and underlines the need for detailed molecular descriptions of binding in representative sequence contexts. We investigated the binding and volumetric properties of DB1976, a model bis(benzimidazole)-selenophene diamidine compound with emerging therapeutic potential in acute myeloid leukemia, debilitating fibroses, and obesity-related liver dysfunction. To sample the scope of cognate DB1976 target sites, we evaluated three dodecameric duplexes spanning >10³-fold in binding affinity. The attendant changes in partial molar volumes varied substantially, but not in step with binding affinity, suggesting distinct modes of interactions in these complexes. Specifically, whereas optimal binding was associated with loss of hydration water, low-affinity binding released more hydration water. Explicit-atom molecular dynamics simulations showed that minor groove binding perturbed the conformational dynamics and hydration at the termini and interior of the DNA in a sequence-dependent manner. The impact of these distinct local dynamics on hydration was experimentally validated by domain-specific interrogation of hydration with salt, which probed the charged axial surfaces of oligomeric DNA preferentially over the uncharged termini. Minor groove recognition by DB1976, therefore, generates dynamically distinct domains that can make favorable contributions to hydration release in both high- and low-affinity binding. Because ligand binding at internal sites of DNA oligomers modulates dynamics at the termini, the results suggest both short- and long-range dynamic effects along the DNA target that can influence their effectiveness as low-MW competitors of protein binding.

Druggable Lipid GPCRs: Past, Present, and Prospects

G protein-coupled receptors (GPCRs) have seven transmembrane spanning domains and comprise the largest superfamily with ~800 receptors in humans. GPCRs are attractive targets for drug discovery because they transduce intracellular signaling in response to endogenous ligands via heterotrimeric G proteins or arrestins, resulting in a wide variety of physiological and pathophysiological responses. The endogenous ligands for GPCRs are highly chemically diverse and include ions, biogenic amines, nucleotides, peptides, and lipids. In this review, we follow the KonMari method to better understand druggable lipid GPCRs. First, we have a comprehensive tidying up of lipid GPCRs including receptors for prostanoids, leukotrienes, specialized pro-resolving mediators (SPMs), lysophospholipids, sphingosine 1-phosphate (S1P), cannabinoids, platelet-activating factor (PAF), free fatty acids (FFAs), and sterols. This tidying up consolidates 46 lipid GPCRs and declutters several perplexing lipid GPCRs. Then, we further tidy up the lipid GPCR-directed drugs from the literature and databases, which identified 24 clinical drugs targeting 16 unique lipid GPCRs available in the market and 44 drugs under evaluation in more than 100 clinical trials as of 2019. Finally, we introduce drug designs for GPCRs that spark joy, such as positive or negative allosteric modulators (PAM or NAM), biased agonism, functional antagonism like fingolimod, and monoclonal antibodies (MAbs). These strategic drug designs may increase the efficacy and specificity of drugs and reduce side effects. Technological advances will help to discover more endogenous lipid ligands from the vast number of remaining orphan GPCRs and will also lead to the development novel lipid GPCR drugs to treat various diseases.

Structural elucidation of the cis-prenyltransferase NgBR/DHDDS complex reveals insights in regulation of protein glycosylation

Cis-prenyltransferase (cis-PTase) catalyzes the rate-limiting step in the synthesis of glycosyl carrier lipids required for protein glycosylation in the lumen of endoplasmic reticulum. Here, we report the crystal structure of the human NgBR/DHDDS complex, which represents an atomic resolution structure for any heterodimeric cis-PTase. The crystal structure sheds light on how NgBR stabilizes DHDDS through dimerization, participates in the enzyme's active site through its C-terminal -RXG- motif, and how phospholipids markedly stimulate cis-PTase activity. Comparison of NgBR/DHDDS with homodimeric cis-PTase structures leads to a model where the elongating isoprene chain extends beyond the enzyme's active site tunnel, and an insert within the α3 helix helps to stabilize this energetically unfavorable state to enable long-chain synthesis to occur. These data provide unique insights into how heterodimeric cis-PTases have evolved from their ancestral, homodimeric forms to fulfill their function in long-chain polyprenol synthesis.

Cryo-EM structure of arabinosyltransferase EmbB from Mycobacterium smegmatis

Arabinosyltransferase B (EmbB) belongs to a family of membrane-bound glycosyltransferases that build the lipidated polysaccharides of the mycobacterial cell envelope, and are targets of anti-tuberculosis drug ethambutol. We present the 3.3 Å resolution single-particle cryo-electron microscopy structure of Mycobacterium smegmatis EmbB, providing insights on substrate binding and reaction mechanism. Mutations that confer ethambutol resistance map mostly around the putative active site, suggesting this to be the location of drug binding.

X-ray structure of C-phycocyanin from Galdieria phlegrea: Determinants of thermostability and comparison with a C-phycocyanin in the entire phycobilisome

Galdieria phlegrea is a polyextremophilic red alga belonging to Cyanidiophyceae. Galdieria phlegrea C-phycocyanin (GpPC), an abundant light-harvesting pigment with an important role in energy capture and transfer to photosystems, is the C-phycocyanin (C-PC) with the highest thermal stability described so far. GpPC also presents interesting antioxidant and anticancer activities. The X-ray structure of the protein was here solved. GpPC is a [(αβ)₃]₂ hexamer, with the phycocyanobilin chromophore attached to Cys84α, Cys82β and Cys153β. Details of geometry and interaction with solvent of the chromophores are reported. Comparison with the structure of a C-PC in the entire Porphyridium purpureum phycobilisome system reveals that linker polypeptides have a significant effect on the local structure of the chromophores environment. Comparative analyses with the structures of other purified C-PCs, which were carried out including re-refined models of G. sulphuraria C-PC, reveal that GpPC presents a significantly higher number of inter-trimer salt bridges. Notably, the higher number of salt bridges at the (αβ)₃/(αβ)₃ interface is not due to an increased number of charged residues in this region, but to subtle conformational variations of their side chains, which are the result of mutations of close polar and non-polar residues.

Cryo-EM Structures and Regulation of Arabinofuranosyltransferase AftD from Mycobacteria

Mycobacterium tuberculosis causes tuberculosis, a disease that kills over 1 million people each year. Its cell envelope is a common antibiotic target and has a unique structure due, in part, to two lipidated polysaccharides-arabinogalactan and lipoarabinomannan. Arabinofuranosyltransferase D (AftD) is an essential enzyme involved in assembling these glycolipids. We present the 2.9-Å resolution structure of M. abscessus AftD, determined by single-particle cryo-electron microscopy. AftD has a conserved GT-C glycosyltransferase fold and three carbohydrate-binding modules. Glycan array analysis shows that AftD binds complex arabinose glycans. Additionally, AftD is non-covalently complexed with an acyl carrier protein (ACP). 3.4- and 3.5-Å structures of a mutant with impaired ACP binding reveal a conformational change, suggesting that ACP may regulate AftD function. Mutagenesis experiments using a conditional knockout constructed in M. smegmatis confirm the essentiality of the putative active site and the ACP binding for AftD function.

Caver Web 1.0: identification of tunnels and channels in proteins and analysis of ligand transport

Caver Web 1.0 is a web server for comprehensive analysis of protein tunnels and channels, and study of the ligands’ transport through these transport pathways. Caver Web is the first interactive tool allowing both the analyses within a single graphical user interface. The server is built on top of the abundantly used tunnel detection tool Caver 3.02 and CaverDock 1.0 enabling the study of the ligand transport. The program is easy-to-use as the only required inputs are a protein structure for a tunnel identification and a list of ligands for the transport analysis. The automated guidance procedures assist the users to set up the calculation in a way to obtain biologically relevant results. The identified tunnels, their properties, energy profiles and trajectories for ligands’ passages can be calculated and visualized. The tool is very fast (2–20 min per job) and is applicable even for virtual screening purposes. Its simple setup and comprehensive graphical user interface make the tool accessible for a broad scientific community. The server is freely available at https://loschmidt.chemi.muni.cz/caverweb.

An integrated approach to unravel a crucial structural property required for the function of the insect steroidogenic Halloween protein Noppera-bo

Ecdysteroids are the principal steroid hormones essential for insect development and physiology. In the last 18 years, several enzymes responsible for ecdysteroid biosynthesis encoded by Halloween genes were identified and genetically and biochemically characterized. However, the tertiary structures of these proteins have not yet been characterized. Here, we report the results of an integrated series of in silico, in vitro, and in vivo analyses of the Halloween GST protein Noppera-bo (Nobo). We determined crystal structures of Drosophila melanogaster Nobo (DmNobo) complexed with GSH and 17β-estradiol, a DmNobo inhibitor. 17β-Estradiol almost fully occupied the putative ligand-binding pocket and a prominent hydrogen bond formed between 17β-estradiol and Asp-113 of DmNobo. We found that Asp-113 is essential for 17β-estradiol–mediated inhibition of DmNobo enzymatic activity, as 17β-estradiol did not inhibit and physically interacted less with the D113A DmNobo variant. Asp-113 is highly conserved among Nobo proteins, but not among other GSTs, implying that this residue is important for endogenous Nobo function. Indeed, a homozygous nobo allele with the D113A substitution exhibited embryonic lethality and an undifferentiated cuticle structure, a phenocopy of complete loss-of-function nobo homozygotes. These results suggest that the nobo family of GST proteins has acquired a unique amino acid residue that appears to be essential for binding an endogenous sterol substrate to regulate ecdysteroid biosynthesis. To the best of our knowledge, ours is the first study describing the structural characteristics of insect steroidogenic Halloween proteins. Our findings provide insights relevant for applied entomology to develop insecticides that specifically inhibit ecdysteroid biosynthesis.

2.7 Å cryo-EM structure of rotavirus core protein VP3, a unique capping machine with a helicase activity

In many viruses, including rotavirus (RV), the major pathogen of infantile gastroenteritis, capping of viral messenger RNAs is a pivotal step for efficient translation of the viral genome. In RV, VP3 caps the nascent transcripts synthesized from the genomic dsRNA segments by the RV polymerase VP1 within the particle core. Here, from cryo-electron microscopy, x-ray crystallography, and biochemical analyses, we show that VP3 forms a stable tetrameric assembly with each subunit having a modular domain organization, which uniquely integrates five distinct enzymatic steps required for capping the transcripts. In addition to the previously known guanylyl- and methyltransferase activities, we show that VP3 exhibits hitherto unsuspected RNA triphosphatase activity necessary for initiating transcript capping and RNA helicase activity likely required for separating the RNA duplex formed transiently during endogenous transcription. From our studies, we propose a new mechanism for how VP3 inside the virion core caps the nascent transcripts exiting from the polymerase.

The ABC exporter IrtAB imports and reduces mycobacterial siderophores

Intracellular replication of the deadly pathogen Mycobacterium tuberculosis relies on the production of small organic molecules called siderophores to scavenge iron from host proteins¹. M. tuberculosis produces two classes of siderophores, lipid-bound mycobactin and soluble carboxymycobactin^{2, 3}. Functional studies revealed that iron-loaded carboxymycobactin is imported into the cytoplasm by the ABC transporter IrtAB⁴, which features an additional cytoplasmic siderophore interaction domain (SID)⁵. However, IrtAB’s predicted ABC exporter fold seemingly contradicts its import function. Here, we show that membrane-reconstituted IrtAB is sufficient to import mycobactins, which are then reduced by the SID to facilitate iron release. Structure determination by X-ray crystallography and cryo-EM confirms IrtAB’s ABC exporter fold, but also reveals structural peculiarities at the transmembrane region of IrtAB resulting in a partially collapsed inward-facing substrate binding cavity. The SID is positioned in close proximity to the inner membrane leaflet, which allows the reduction of membrane-inserted mycobactin. Enzymatic ATPase activity and in vivo growth assays show that IrtAB prefers mycobactin over carboxymycobactin as its substrate. Our study provides insights into an unusual ABC exporter that evolved as highly specialized siderophore import machinery in mycobacteria.

Mechanisms of Resistance Associated with the Inhibition of the Dehydration Step of Type II Fatty Acid Synthase in Mycobacterium tuberculosis

Isoxyl (ISO) and thiacetazone (TAC) are two antitubercular prodrugs that abolish mycolic acid biosynthesis and kill Mycobacterium tuberculosis (Mtb) through the inhibition of the essential type II fatty acid synthase (FAS-II) dehydratase HadAB. While mutations preventing ISO and TAC either from being converted to their active form or from covalently modifying their target are the most frequent spontaneous mutations associated with high-level resistance to both drugs, the molecular mechanisms underlying the high-level ISO and TAC resistance of Mtb strains harboring missense mutations in the second, nonessential, FAS-II dehydratase HadBC have remained unexplained. Using a combination of genetic, biochemical, and biophysical approaches and molecular dynamics simulation, we here show that all four reported resistance mutations in the HadC subunit of HadBC alter the stability and/or specific activity of the enzyme, allowing it in two cases (HadBC^V85I and HadBC^K157R) to compensate for a deficiency in HadAB in whole Mtb bacilli. The analysis of the mycolic acid profiles of Mtb strains expressing the mutated forms of HadC further points to alterations in the activity of the mycolic acid biosynthetic complex and suggests an additional contributing resistance mechanism whereby HadC mutations may reduce the accessibility of HadAB to ISO and TAC. Collectively, our results highlight the importance of developing optimized inhibitors of the dehydration step of FAS-II capable of inhibiting both dehydratases simultaneously, a goal that may be achievable given the structural resemblance of the two enzymes and their reliance on the same catalytic subunit HadB.

Discovery, expression, cellular localization, and molecular properties of a novel, alternative spliced HP1γ isoform, lacking the chromoshadow domain

By reading the H3K9Me3 mark through their N-terminal chromodomain (CD), HP1 proteins play a significant role in cancer-associated processes, including cell proliferation, differentiation, chromosomal stability, and DNA repair. Here, we used a combination of bioinformatics-based methodologies, as well as experimentally-derived datasets, that reveal the existence of a novel short HP1γ (CBX3) isoform, named here sHP1γ, generated by alternative splicing of the CBX3 locus. The sHP1γ mRNA encodes a protein composed of 101 residues and lacks the C-terminal chromoshadow domain (CSD) that is required for dimerization and heterodimerization in the previously described 183 a. a HP1γ protein. Fold recognition, order-to-disorder calculations, threading, homology-based molecular modeling, docking, and molecular dynamic simulations show that the sHP1γ is comprised of a CD flanked by intrinsically disordered regions (IDRs) with an IDR-CD-IDR domain organization and likely retains the ability to bind to the H3K9Me3. Both qPCR analyses and mRNA-seq data derived from large-scale studies confirmed that sHP1γ mRNA is expressed in the majority of human tissues at approximately constant ratios with the chromoshadow domain containing isoform. However, sHP1γ mRNA levels appear to be dysregulated in different cancer types. Thus, our data supports the notion that, due to the existence of functionally different isoforms, the regulation of HP1γ-mediated functions is more complex than previously anticipated.

Steered Pull Simulation to Determine Nanomechanical Properties of Cellulose Nanofiber

Cellulose nanofiber (CNF) exhibits excellent mechanical properties, which has been extensively proven through experimental techniques. However, understanding the mechanisms and the inherent structural behavior of cellulose is important in its vastly growing research areas of applications. This study focuses on taking a look into what happens to the atomic molecular interactions of CNF, mainly hydrogen bond, in the presence of external force. This paper investigates the hydrogen bond disparity within CNF structure. To achieve this, molecular dynamics simulations of cellulose Iβ nanofibers are carried out in equilibrated conditions in water using GROMACS software in conjunction with OPLS-AA force field. It is noted that the hydrogen bonds within the CNF are disrupted when a pulling force is applied. The simulated Young’s modulus of CNF is found to be 161 GPa. A simulated shear within the cellulose chains presents a trend with more hydrogen bond disruptions at higher forces.

Structural dynamics and allostery of Rab proteins: strategies for drug discovery and design

Rab proteins represent the largest family of the Rab superfamily guanosine triphosphatase (GTPase). Aberrant human Rab proteins are associated with multiple diseases, including cancers and neurological disorders. Rab subfamily members display subtle conformational variations that render specificity in their physiological functions and can be targeted for subfamily-specific drug design. However, drug discovery efforts have not focused much on targeting Rab allosteric non-nucleotide binding sites which are subjected to less evolutionary pressures to be conserved, hence are likely to offer subfamily specificity and may be less prone to undesirable off-target interactions and side effects. To discover druggable allosteric binding sites, Rab structural dynamics need to be first incorporated using multiple experimentally and computationally obtained structures. The high-dimensional structural data may necessitate feature extraction methods to identify manageable representative structures for subsequent analyses. We have detailed state-of-the-art computational methods to (i) identify binding sites using data on sequence, shape, energy, etc., (ii) determine the allosteric nature of these binding sites based on structural ensembles, residue networks and correlated motions and (iii) identify small molecule binders through structure- and ligand-based virtual screening. To benefit future studies for targeting Rab allosteric sites, we herein detail a refined workflow comprising multiple available computational methods, which have been successfully used alone or in combinations. This workflow is also applicable for drug discovery efforts targeting other medically important proteins. Depending on the structural dynamics of proteins of interest, researchers can select suitable strategies for allosteric drug discovery and design, from the resources of computational methods and tools enlisted in the workflow.

Structural basis of substrate recognition by a polypeptide processing and secretion transporter

The peptidase-containing ATP-binding cassette transporters (PCATs) are unique members of the ABC transporter family that proteolytically process and export peptides and proteins. Each PCAT contains two peptidase domains that cleave off the secretion signal, two transmembrane domains forming a translocation pathway, and two nucleotide-binding domains that hydrolyze ATP. Previously the crystal structures of a PCAT from Clostridium thermocellum (PCAT1) were determined in the absence and presence of ATP, revealing how ATP binding regulates the protease activity and access to the translocation pathway. However, how the substrate CtA, a 90-residue polypeptide, is recognized by PCAT1 remained elusive. To address this question, we determined the structure of the PCAT1-CtA complex by electron cryo-microscopy (cryo-EM) to 3.4 Å resolution. The structure shows that two CtAs are bound via their N-terminal leader peptides, but only one is positioned for cleavage and translocation. Based on these results, we propose a model of how substrate cleavage, ATP hydrolysis, and substrate translocation are coordinated in a transport cycle.

Catalytic mechanism and evolutionary characteristics of thioredoxin from Halobacterium salinarum NRC-1

Thioredoxin (TRX) is an important antioxidant against oxidative stress. TRX from the extremely halophilic archaeon Halobacterium salinarum NRC-1 (HsTRX-A), which has the highest acidic residue content [(Asp + Glu)/(Arg + Lys + His) = 9.0] among known TRXs, was chosen to elucidate the catalytic mechanism and evolutionary characteristics associated with haloadaptation. X-ray crystallographic analysis revealed that the main-chain structure of HsTRX-A is similar to those of homologous TRXs; for example, the root-mean-square deviations on C^α atoms were <2.3 Å for extant archaeal TRXs and <1.5 Å for resurrected Precambrian TRXs. A unique water network was located near the active-site residues (Cys45 and Cys48) in HsTRX-A, which may enhance the proton transfer required for the reduction of substrates under a high-salt environment. The high density of negative charges on the molecular surface (3.6 × 10^-3 e Å^-2) should improve the solubility and haloadaptivity. Moreover, circular-dichroism measurements and enzymatic assays using a mutant HsTRX-A with deletion of the long flexible N-terminal region (Ala2-Pro17) revealed that Ala2-Pro17 improves the structural stability and the enzymatic activity of HsTRX-A under high-salt environments (>2 M NaCl). The elongation of the N-terminal region in HsTRX-A accompanies the increased hydrophilicity and acidic residue content but does not affect the structure of the active site. These observations offer insights into molecular evolution for haloadaptation and potential applications in halophilic protein-related biotechnology.

Metal ion adaptive self-assembly of photoactive lanthanide-based supramolecular hosts

Metal-adaptive self-assembly and post-assembly transmetallation modification of functional lanthanide-porphyrin hosts were presented.

A metal-peptide capsule by multiple ring threading

Cavity creation is a key to the origin of biological functions. Small cavities such as enzyme pockets are created simply through liner peptide folding. Nature can create much larger cavities by threading and entangling large peptide rings, as learned from gigantic virus capsids, where not only chemical structures but the topology of threaded rings must be controlled. Although interlocked molecules are a topic of current interest, they have for decades been explored merely as elements of molecular machines, or as a synthetic challenge. No research has specifically targeted them for, and succesfully achieved, cavity creation. Here we report the emergence of a huge capsular framework via multiple threading of metal-peptide rings. Six equivalent C₄-propeller-shaped rings, each consisting of four oligopeptides and Ag⁺, are threaded by each other a total of twelve times (crossing number: 24) to assemble into a well-defined 4 nm-sized sphere, which acts as a huge molecular capsule.

MOLC. A reversible coarse grained approach using anisotropic beads for the modelling of organic functional materials

We describe the development and implementation of a coarse grained (CG) modelling approach where complex organic molecules, and particularly the π-conjugated ones often employed in organic electronics, are modelled in terms of connected sets of attractive-repulsive biaxial Gay-Berne ellipsoidal beads. The CG model is aimed at reproducing realistically large scale morphologies (e.g. up to 100 nm thick films) for the materials involved, while being able to generate, with a back-mapping procedure, atomistic coordinates suitable, with limited effort, to be applied for charge transport calculations. Detailed methodology and an application to the common hole transporter material α-NPD are provided.

Identifying Polymorphs of Amyloid-β (1-40) Fibrils Using High-Resolution Atomic Force Microscopy

Many amyloid-β fibril preparations are highly polymorphic, and the conditions under which they are formed determine their morphology. This report describes the application of high-resolution atomic force microscopy (HR-AFM), combined with volume-per-length analysis, to define, identify, and quantify the structural components of polymorphic Aβ fibril preparations. Volume-per-length analysis confirms that they are composed of discrete cross-β filaments, and the analysis of HR-AFM images yields the number of striations in each fibril. Compared to mass-per-length analysis by electron microscopy, HR-AFM analysis yields narrower distributions, facilitating rapid and label-free quantitative morphological characterization of Aβ fibril preparations.

Structure and drug resistance of the Plasmodium falciparum transporter PfCRT

The emergence and spread of drug-resistant Plasmodium falciparum impedes global efforts to control and eliminate malaria. For decades, treatment of malaria has relied on chloroquine (CQ), a safe and affordable 4-aminoquinoline that was highly effective against intra-erythrocytic asexual blood-stage parasites, until resistance arose in Southeast Asia and South America and spread worldwide¹. Clinical resistance to the chemically related current first-line combination drug piperaquine (PPQ) has now emerged regionally, reducing its efficacy². Resistance to CQ and PPQ has been associated with distinct sets of point mutations in the P. falciparum CQ-resistance transporter PfCRT, a 49-kDa member of the drug/metabolite transporter superfamily that traverses the membrane of the acidic digestive vacuole of the parasite^3-9. Here we present the structure, at 3.2 Å resolution, of the PfCRT isoform of CQ-resistant, PPQ-sensitive South American 7G8 parasites, using single-particle cryo-electron microscopy and antigen-binding fragment technology. Mutations that contribute to CQ and PPQ resistance localize primarily to moderately conserved sites on distinct helices that line a central negatively charged cavity, indicating that this cavity is the principal site of interaction with the positively charged CQ and PPQ. Binding and transport studies reveal that the 7G8 isoform binds both drugs with comparable affinities, and that these drugs are mutually competitive. The 7G8 isoform transports CQ in a membrane potential- and pH-dependent manner, consistent with an active efflux mechanism that drives CQ resistance⁵, but does not transport PPQ. Functional studies on the newly emerging PfCRT F145I and C350R mutations, associated with decreased PPQ susceptibility in Asia and South America, respectively^6,9, reveal their ability to mediate PPQ transport in 7G8 variant proteins and to confer resistance in gene-edited parasites. Structural, functional and in silico analyses suggest that distinct mechanistic features mediate the resistance to CQ and PPQ in PfCRT variants. These data provide atomic-level insights into the molecular mechanism of this key mediator of antimalarial treatment failures.

Oleoresin defenses in conifers: chemical diversity, terpene synthases and limitations of oleoresin defense under climate change

Conifers have evolved complex oleoresin terpene defenses against herbivores and pathogens. In co-evolved bark beetles, conifer terpenes also serve chemo-ecological functions as pheromone precursors, chemical barcodes for host identification, or nutrients for insect-associated microbiomes. We highlight the genomic, molecular and biochemical underpinnings of the large chemical space of conifer oleoresin terpenes and volatiles. Conifer terpenes are predominantly the products of the conifer terpene synthase (TPS) gene family. Terpene diversity is increased by cytochromes P450 of the CYP720B class. Many conifer TPS are multiproduct enzymes. Multisubstrate CYP720B enzymes catalyse multistep oxidations. We summarise known terpenoid gene functions in various different conifer species with reference to the annotated terpenoid gene space in a spruce genome. Overall, biosynthesis of terpene diversity in conifers is achieved through a system of biochemical radiation and metabolic grids. Expression of TPS and CYP720B genes can be specific to individual cell types of constitutive or traumatic resin duct systems. Induced terpenoid transcriptomes in resin duct cells lead to dynamic changes of terpene composition and quantity to fend off herbivores and pathogens. While terpenoid defenses have contributed much to the evolutionary success of conifers, under new conditions of climate change, these defences may become inconsequential against range-expanding forest pests.

Differences in the path to exit the ribosome across the three domains of life

The ribosome exit tunnel is an important structure involved in the regulation of translation and other essential functions such as protein folding. By comparing 20 recently obtained cryo-EM and X-ray crystallography structures of the ribosome from all three domains of life, we here characterize the key similarities and differences of the tunnel across species. We first show that a hierarchical clustering of tunnel shapes closely reflects the species phylogeny. Then, by analyzing the ribosomal RNAs and proteins, we explain the observed geometric variations and show direct association between the conservations of the geometry, structure and sequence. We find that the tunnel is more conserved in the upper part close to the polypeptide transferase center, while in the lower part, it is substantially narrower in eukaryotes than in bacteria. Furthermore, we provide evidence for the existence of a second constriction site in eukaryotic exit tunnels. Overall, these results have several evolutionary and functional implications, which explain certain differences between eukaryotes and prokaryotes in their translation mechanisms. In particular, they suggest that major co-translational functions of bacterial tunnels were externalized in eukaryotes, while reducing the tunnel size provided some other advantages, such as facilitating the nascent chain elongation and enabling antibiotic resistance.

Elastic continuum stiffness of contractile tail sheaths from molecular dynamics simulations

Contractile tails are key components of the biological nanomachinery involved in cell membrane puncturing, where they provide a means to deliver molecules and ions inside cells. Two intriguing examples of contractile tails are those from bacteriophage T4 and R2-pyocin. Although the two systems are different in terms of biological activity, they share a fascinatingly similar injection mechanism, during which the tail sheaths of both systems contract from a so-called extended state to around half of their length (the contracted state), accompanied by release of elastic energy originally stored in the sheath. Despite the great prevalence and biomedical importance of contractile delivery systems, many fundamental details of their injection machinery and dynamics are still unknown. In this work, we calculate the bending and torsional stiffness constants of a helical tail sheath strand of bacteriophage T4 and R2-pyocin, in both extended and contracted states, using molecular dynamics simulations of about one-sixth of the entire sheath. Differences in stiffness constants between the two systems are rationalized by comparing their all-atom monomer structures, changes in sheath architecture on contraction, and differences in interstrand interactions. The calculated coefficients indicate that the T4 strand is stiffer for both bending and torsion than the corresponding R2-pyocin strands in both extended and contracted conformations. The sheath strands also have greater stiffness in the contracted state for both systems. As the main application of this study, we describe how the stiffness constants can be incorporated in a model to simulate the dynamics of contractile nanoinjection machineries.

Solvent flows, conformation changes and lattice reordering in a cold protein crystal

When protein crystals are abruptly cooled, the unit-cell, protein and solvent-cavity volumes all contract, but the volume of bulk-like internal solvent may expand. Outflow of this solvent from the unit cell and its accumulation in defective interior crystal regions has been suggested as one cause of the large increase in crystal mosaicity on cooling. It is shown that when apoferritin crystals are abruptly cooled to temperatures between 220 and 260 K, the unit cell contracts, solvent is pushed out and the mosaicity grows. On temperature-dependent timescales of 10 to 200 s, the unit-cell and solvent-cavity volume then expand, solvent flows back in, and the mosaicity and B factor both drop. Expansion and reordering at fixed low temperature are associated with small-amplitude but large-scale changes in the conformation and packing of apoferritin. These results demonstrate that increases in mosaicity on cooling arise due to solvent flows out of or into the unit cell and to incomplete, arrested relaxation of protein conformation. They indicate a critical role for time in variable-temperature crystallographic studies, and the feasibility of probing interactions and cooperative conformational changes that underlie cold denaturation in the presence of liquid solvent at temperatures down to ∼200 K.

Sera Antibody Repertoire Analyses Reveal Mechanisms of Broad and Pandemic Strain Neutralizing Responses after Human Norovirus Vaccination

Rapidly evolving RNA viruses, such as the GII.4 strain of human norovirus (HuNoV), and their vaccines elicit complex serological responses associated with previous exposure. Specific correlates of protection, moreover, remain poorly understood. Here, we report the GII.4-serological antibody repertoire—pre- and post-vaccination—and select several antibody clonotypes for epitope and structural analysis. The humoral response was dominated by GII.4-specific antibodies that blocked ancestral strains or by antibodies that bound to divergent genotypes and did not block viral-entry-ligand interactions. However, one antibody, A1431, showed broad blockade toward tested GII.4 strains and neutralized the pandemic GII.P16-GII.4 Sydney strain. Structural mapping revealed conserved epitopes, which were occluded on the virion or partially exposed, allowing for broad blockade with neutralizing activity. Overall, our results provide high-resolution molecular information on humoral immune responses after HuNoV vaccination and demonstrate that infection-derived and vaccine-elicited antibodies can exhibit broad blockade and neutralization against this prevalent human pathogen.

•

Serum vaccine response is dominated by a small number of abundant antibody clonotypes

•

Vaccine-boosted antibodies predominantly target conserved norovirus epitopes

•

Identified cross-genogroup and strain-specific epitopes

•

Discovered a pandemic-genotype neutralizing antibody recognizing a conserved epitope

A Chemosensory GPCR as a Potential Target to Control the Root-Knot Nematode Meloidogyne incognita Parasitism in Plants

Root-knot nematodes (RKN), from the Meloidogyne genus, have a worldwide distribution and cause severe economic damage to many life-sustaining crops. Because of their lack of specificity and danger to the environment, most chemical nematicides have been banned from use. Thus, there is a great need for new and safe compounds to control RKN. Such research involves identifying beforehand the nematode proteins essential to the invasion. Since G protein-coupled receptors GPCRs are the target of a large number of drugs, we have focused our research on the identification of putative nematode GPCRs such as those capable of controlling the movement of the parasite towards (or within) its host. A datamining procedure applied to the genome of Meloidogyne incognita allowed us to identify a GPCR, belonging to the neuropeptide GPCR family that can serve as a target to carry out a virtual screening campaign. We reconstructed a 3D model of this receptor by homology modeling and validated it through extensive molecular dynamics simulations. This model was used for large scale molecular dockings which produced a filtered limited set of putative antagonists for this GPCR. Preliminary experiments using these selected molecules allowed the identification of an active compound, namely C260-2124, from the ChemDiv provider, which can serve as a starting point for further investigations.

CavBench: A benchmark for protein cavity detection methods

Extensive research has been applied to discover new techniques and methods to model protein-ligand interactions. In particular, considerable efforts focused on identifying candidate binding sites, which quite often are active sites that correspond to protein pockets or cavities. Thus, these cavities play an important role in molecular docking. However, there is no established benchmark to assess the accuracy of new cavity detection methods. In practice, each new technique is evaluated using a small set of proteins with known binding sites as ground-truth. However, studies supported by large datasets of known cavities and/or binding sites and statistical classification (i.e., false positives, false negatives, true positives, and true negatives) would yield much stronger and reliable assessments. To this end, we propose CavBench, a generic and extensible benchmark to compare different cavity detection methods relative to diverse ground truth datasets (e.g., PDBsum) using statistical classification methods.

Crystal structure of a bacterial homolog to human lysosomal transporter, spinster

Lysosomes break down various biomolecules and spinster is one of the major efflux carriers removing degradation products from lysosomal lumen to keep it in healthy size and proper function. Although it is well established that a dysfunctional spinster will cause enlarged lysosomes and in turn lead to developmental defects and abnormal behavior in animals, little was known about the transportation mechanism and substrate specificity of spinster. Here, we report a crystal structure of spinster homolog from Hyphomonas neptunium, HnSPNS, in its inward-facing conformation with and without substrate bound. HnSPNS is crystallized in a monomer and a substrate-binding cavity was formed in the center of its transmembrane helices. A blob of electron density corresponding to its substrate was found in the cavity near a conserved residue, R42, which is locked in position by the interactions with conserved residues E129 and R122. Our results suggest that human spinster serves as a transporter translocating negatively-charged lipophilic small molecules and E129 might serve as a switch to control the conformational change via its protonation-deprotonation cycle.

Biochemical and Structural Characterization of XoxG and XoxJ and Their Roles in Lanthanide-Dependent Methanol Dehydrogenase Activity

Lanthanide (Ln)-dependent methanol dehydrogenases (MDHs) have recently been shown to be widespread in methylotrophic bacteria. Along with the core MDH protein, XoxF, these systems contain two other proteins, XoxG (a c-type cytochrome) and XoxJ (a periplasmic binding protein of unknown function), about which little is known. In this work, we have biochemically and structurally characterized these proteins from the methyltroph Methylobacterium extorquens AM1. In contrast to results obtained in an artificial assay system, assays of XoxFs metallated with LaIII , CeIII , and NdIII using their physiological electron acceptor, XoxG, display Ln-independent activities, but the Km for XoxG markedly increases from La to Nd. This result suggests that XoxG's redox properties are tuned specifically for lighter Lns in XoxF, an interpretation supported by the unusually low reduction potential of XoxG (+172 mV). The X-ray crystal structure of XoxG provides a structural basis for this reduction potential and insight into the XoxG-XoxF interaction. Finally, the X-ray crystal structure of XoxJ reveals a large hydrophobic cleft and suggests a role in the activation of XoxF. These studies enrich our understanding of the underlying chemical principles that enable the activity of XoxF with multiple lanthanides in vitro and in vivo.

Predicting Protein-Polymer Block Copolymer Self-Assembly from Protein Properties

Protein–polymer bioconjugate self-assembly has attracted a great deal of attention as a method to fabricate protein nanomaterials in solution and the solid state. To identify protein properties that affect phase behavior in protein–polymer block copolymers, a library of 15 unique protein-b-poly(N-isopropylacrylamide) (PNIPAM) copolymers comprising 11 different proteins was compiled and analyzed. Many attributes of phase behavior are found to be similar among all studied bioconjugates regardless of protein properties, such as formation of micellar phases at high temperature and low concentration, lamellar ordering with increasing temperature, and disordering at high concentration, but several key protein-dependent trends are also observed. In particular, hexagonal phases are only observed for proteins within the molar mass range 20–36 kDa, where ordering quality is also significantly enhanced. While ordering is generally found to improve with increasing molecular weight outside of this range, most large bioconjugates exhibited weaker than predicted assembly, which is attributed to chain entanglement with increasing polymer molecular weight. Additionally, order–disorder transition boundaries are found to be largely uncorrelated to protein size and quality of ordering. However, the primary finding is that bioconjugate ordering can be accurately predicted using only protein molecular weight and percentage of residues contained within β sheets. This model provides a basis for designing protein–PNIPAM bioconjugates that exhibit well-defined self-assembly and a modeling framework that can generalize to other bioconjugate chemistries.

A Java Application to Characterize Biomolecules and Nanomaterials in Electrolyte Aqueous Solutions

The electrostatic, entropic and surface interactions between a macroion (nanoparticle or biomolecule), surrounding ions and water molecules play a fundamental role in the behavior and function of colloidal systems. However, the molecular mechanisms governing these phenomena are still poorly understood. One of the major limitations in procuring this understanding is the lack of appropriate computational tools. Additionally, only experts in the field with an extensive background in programming, who are trained in statistical mechanics, and have access to supercomputers are able to study these systems. To overcome these limitations, in this article, we present a free, multiplatform, portable Java software, which provides experts and non-experts in the field an easy and efficient way to obtain an accurate molecular characterization of electrical and structural properties of aqueous electrolyte mixture solutions around both cylindrical- and spherical-like rigid macroions under multiple conditions. These properties include the normalized ions and water density profile distributions, the mean electrostatic potential, the integrated charge, the zeta potential, the electrostatic potential energy, the particle crowding entropy energy, the ion-ion electrostatic direct correlation energy, and the ionic potential of mean force. The Java software does not require outstanding skills and comes with detailed user-guide documentation. The application is based on the so-called Classical Density Functional Theory Solver (CSDFTS), which was successfully applied to a variety of rod-like biopolymers, rigid-like globular proteins, nanoparticles, and nano-rods. CSDFTS implements several electrolyte and macroion models, uses different levels of approximation and takes advantage of high performance Fortran90 routines and optimized libraries. These features enable the software to run on single processor computers at low-to-moderate computational cost depending on the computer performance, the grid resolution, and the characterization of the macroion and the electrolyte solution, among other factors. As a unique feature, the software comes with a graphical user interface (GUI) that allows users to take advantage of the visually guided setup of the required input data to properly characterize the system and configure the solver. Several examples on nanomaterials and biomolecules are provided to illustrate the use of the GUI and the solver performance.

Construction of Type III-C Rotaxane-Branched Dendrimers and Their Anion-Induced Dimension Modulation Feature

Starting from a novel rotaxane building block with dendrimer growth sites being located at both the wheel and axle component, we realized the successful construction of a new family of rotaxane-branched dendrimers, i.e., Type III-C rotaxane-branched dendrimers, up to fourth generation as a highly branched [46]rotaxane through a controllable divergent approach. In the resultant rotaxane-branched dendrimers, the wheel components of the rotaxane units are located on the branches as well as at the branching points, making them excellent candidates to mimic the amplified collective molecular motions. Thus, taking advantage of the urea moiety inserted into the axle components of the rotaxane units as the binding sites, the addition or removal of acetate anion as stimulus endows the individual rotaxane unit a switchable feature that lead to a collective expansion-contraction motion of the integrated rotaxane-branched dendrimers, thus allowing for the remarkable and reversible size modulation. Such a three-dimensional size switching feature makes Type III-C rotaxane-branched dendrimers a very promising platform toward the fabrication of novel dynamic smart materials.

The effect of high pressure on the NMDA receptor: molecular dynamics simulations

Professional divers exposed to ambient pressures above 11 bar develop the high pressure neurological syndrome (HPNS), manifesting as central nervous system (CNS) hyperexcitability, motor disturbances, sensory impairment, and cognitive deficits. The glutamate-type N-methyl-D-aspartate receptor (NMDAR) has been implicated in the CNS hyperexcitability of HPNS. NMDARs containing different subunits exhibited varying degrees of increased/decreased current at high pressure. The mechanisms underlying this phenomenon remain unclear. We performed 100 ns molecular dynamics (MD) simulations of the NMDAR structure embedded in a dioleoylphosphatidylcholine (DOPC) lipid bilayer solvated in water at 1 bar, hydrostatic 25 bar, and in helium at 25 bar. MD simulations showed that in contrast to hydrostatic pressure, high pressure helium causes substantial distortion of the DOPC membrane due to its accumulation between the two monolayers: reduction of the Sn-1 and Sn-2 DOPC chains and helium-dependent dehydration of the NMDAR pore. Further analysis of important regions of the NMDAR protein such as pore surface (M2 α-helix), Mg²⁺ binding site, and TMD-M4 α-helix revealed significant effects of helium. In contrast with previous models, these and our earlier results suggest that high pressure helium, not hydrostatic pressure per se, alters the receptor tertiary structure via protein-lipid interactions. Helium in divers’ breathing mixtures may partially contribute to HPNS symptoms.