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

Cryo EM structure of the rabies virus ribonucleoprotein complex

Rabies virus is an important zoonotic pathogen. Its bullet shaped particle contains a helical nucleocapsid. We used cryo-electron tomography and subsequent subtomogram averaging to determine the structure of its ribonucleoprotein. The resulting electron density map allowed for confident fitting of the N-protein crystal structure, indicating that interactions between neighbouring N-proteins are only mediated by N- and C-terminal protruding subdomains (aa 1-27 and aa 355-372). Additional connecting densities, likely stabilizing the ribonucleoprotein complex, are present between neighbouring M-protein densities on the same helical turn and between M- and N-protein densities located on neighbouring helical turns, but not between M-proteins of different turns, as is observed for the related Vesicular stomatitis virus (VSV). This insight into the architecture of the rabies virus nucleocapsid highlights the surprising structural divergence of large biological assemblies even if the building blocks - here exemplified by VSV M- and N-protein - are structurally closely related.

Effect of Macromolecular Crowding on the FMN-Heme Intraprotein Electron Transfer in Inducible NO Synthase

Previous biochemical studies of nitric oxide synthase enzymes (NOSs) were conducted in diluted solutions. However, the intracellular milieu where the proteins perform their biological functions is crowded with macromolecules. The effect of crowding on the electron transfer kinetics of multidomain proteins is much less understood. Herein, we investigated the effect of macromolecular crowding on the FMN-heme intraprotein interdomain electron transfer (IET), an obligatory step in NOS catalysis. A noticeable increase in the IET rate in the bidomain oxygenase/FMN (oxyFMN) and the holoprotein of human inducible NOS (iNOS) was observed upon addition of Ficoll 70 in a nonsaturable manner. Additionally, the magnitude of IET enhancement for the holoenzyme is much higher than that that of the oxyFMN construct. The crowding effect is also evident at different ionic strengths. Importantly, the enhancing extent is similar for the iNOS oxyFMN protein with added Ficoll 70 and Dextran 70 that give the same solution viscosity, showing that specific interactions do not exist between the NOS protein and the crowder. Moreover, the population of the docked FMN-heme state is significantly increased upon addition of Ficoll 70 and the fluorescence lifetime values do not correspond to those in the absence of Ficoll 70. The steady-state cytochrome c reduction by the holoenzyme is noticeably enhanced by the crowder, while the ferricyanide reduction is unchanged. The NO production activity of the iNOS holoenzyme is stimulated by Ficoll 70. The effect of macromolecular crowding on the kinetics can be rationalized on the basis of the excluded volume effect, with an entropic origin. The intraprotein electron transfer kinetics, fluorescence lifetime, and steady-state enzymatic activity results indicate that macromolecular crowding modulates the NOS electron transfer through multiple pathways. Such a mechanism should be applicable to electron transfer in other multidomain redox proteins.

Structures of BCL-2 in complex with venetoclax reveal the molecular basis of resistance mutations

Venetoclax is a first-in-class cancer therapy that interacts with the cellular apoptotic machinery promoting apoptosis. Treatment of patients suffering chronic lymphocytic leukaemia with this BCL-2 antagonist has revealed emergence of a drug-selected BCL-2 mutation (G101V) in some patients failing therapy. To understand the molecular basis of this acquired resistance we describe the crystal structures of venetoclax bound to both BCL-2 and the G101V mutant. The pose of venetoclax in its binding site on BCL-2 reveals small but unexpected differences as compared to published structures of complexes with venetoclax analogues. The G101V mutant complex structure and mutant binding assays reveal that resistance is acquired by a knock-on effect of V101 on an adjacent residue, E152, with venetoclax binding restored by a E152A mutation. This provides a framework for considering analogues of venetoclax that might be effective in combating this mutation.

The BCL-2 mutation G101V reduces venetoclax affinity and confers drug resistance in patients with chronic lymphocytic leukaemia. Here, the authors present crystal structures and biochemical analyses of venetoclax bound to BCL-2 and the G101V mutant, revealing the structural basis for venetoclax resistance.

High-Resolution Crystal Structure of Endoplasmic Reticulum Aminopeptidase 1 with Bound Phosphinic Transition-State Analogue Inhibitor

Endoplasmic reticulum aminopeptidase 1 (ERAP1) is an intracellular enzyme that helps generate peptides presented by Major Histocompatibility Complex Class I (MHC class I) molecules and is an emerging target for immunotherapy applications. Despite almost two decades of research on ERAP1, lack of high-resolution crystal structures has hampered drug-development efforts. By optimizing the protein construct, we obtained a high-resolution (1.60 Å) crystal structure of the closed-conformation of ERAP1 with a potent phosphinic pseudopeptide inhibitor bound in its active site. The structure provides key insight on the mechanism of inhibition as well as selectivity toward homologous enzymes and allows detailed mapping of the internal cavity of the enzyme that accommodates peptide-substrates. Bis-tris propane and malic acid molecules, found bound in pockets in the internal cavity, reveal potential druggable secondary binding sites. The ability to obtain high-resolution crystal structures of ERAP1 removes a major bottleneck in the development of compounds that regulate its activity and will greatly accelerate drug-discovery efforts.

A grid-based algorithm in conjunction with a gaussian-based model of atoms for describing molecular geometry

A novel grid-based method is presented, which in conjunction with a smooth Gaussian-based model of atoms, is used to compute molecular volume (MV) and surface area (MSA). The MV and MSA are essential for computing nonpolar component of free energies. The objective of our grid-based approach is to identify solute atom pairs that share overlapping volumes in space. Once completed, this information is used to construct a rooted tree using depth-first method to yield the final volume and SA by using the formulations of the Gaussian model described by Grant and Pickup (J. Phys Chem, 1995, 99, 3503). The method is designed to function uninterruptedly with the grid-based finite-difference method implemented in Delphi, a popular and open-source package used for solving the Poisson-Boltzmann equation (PBE). We demonstrate the time efficacy of the method while also validating its performance in terms of the effect of grid-resolution, positioning of the solute within the grid-map and accuracy in identification of overlapping atom pairs. We also explore and discuss different aspects of the Gaussian model with key emphasis on its physical meaningfulness. This development and its future release with the Delphi package are intended to provide a physically meaningful, fast, robust and comprehensive tool for MM/PBSA based free energy calculations. © 2019 Wiley Periodicals, Inc.

Ice formation and solvent nanoconfinement in protein crystals

Ice formation within protein crystals is a major obstacle to the cryocrystallographic study of protein structure, and has limited studies of how the structural ensemble of a protein evolves with temperature in the biophysically interesting range from ∼260 K to the protein-solvent glass transition near 200 K. Using protein crystals with solvent cavities as large as ∼70 Å, time-resolved X-ray diffraction was used to study the response of protein and internal solvent during rapid cooling. Solvent nanoconfinement suppresses freezing temperatures and ice-nucleation rates so that ice-free, low-mosaicity diffraction data can be reliably collected down to 200 K without the use of cryoprotectants. Hexagonal ice (Ih) forms in external solvent, but internal crystal solvent forms stacking-disordered ice (Isd) with a near-random stacking of cubic and hexagonal planes. Analysis of powder diffraction from internal ice and single-crystal diffraction from the host protein structure shows that the maximum crystallizable solvent fraction decreases with decreasing crystal solvent-cavity size, and that an ∼6 Å thick layer of solvent adjacent to the protein surface cannot crystallize. These results establish protein crystals as excellent model systems for the study of nanoconfined solvent. By combining fast cooling, intense X-ray beams and fast X-ray detectors, complete structural data sets for high-value targets, including membrane proteins and large complexes, may be collected at ∼220-240 K that have much lower mosaicities and comparable B factors, and that may allow more confident identification of ligand binding than in current cryocrystallographic practice.

Structural insight into YcbB-mediated beta-lactam resistance in Escherichia coli

The bacterial cell wall plays a crucial role in viability and is an important drug target. In Escherichia coli, the peptidoglycan crosslinking reaction to form the cell wall is primarily carried out by penicillin-binding proteins that catalyse D,D-transpeptidase activity. However, an alternate crosslinking mechanism involving the L,D-transpeptidase YcbB can lead to bypass of D,D-transpeptidation and beta-lactam resistance. Here, we show that the crystallographic structure of YcbB consists of a conserved L,D-transpeptidase catalytic domain decorated with a subdomain on the dynamic substrate capping loop, peptidoglycan-binding and large scaffolding domains. Meropenem acylation of YcbB gives insight into the mode of inhibition by carbapenems, the singular antibiotic class with significant activity against L,D-transpeptidases. We also report the structure of PBP5-meropenem to compare interactions mediating inhibition. Additionally, we probe the interaction network of this pathway and assay beta-lactam resistance in vivo. Our results provide structural insights into the mechanism of action and the inhibition of L,D-transpeptidation, and into YcbB-mediated antibiotic resistance.

In E. coli, alternate peptidoglycan crosslinking reactions carried out by the L,D-transpeptidase YcbB can lead to beta-lactam resistance. Here, Caveney et al. solve the crystal structure of YcbB and shed light into its mechanism of action and into YcbB-mediated antibiotic resistance.

Engineering enzyme access tunnels

Enzymes are efficient and specific catalysts for many essential reactions in biotechnological and pharmaceutical industries. Many times, the natural enzymes do not display the catalytic efficiency, stability or specificity required for these industrial processes. The current enzyme engineering methods offer solutions to this problem, but they mainly target the buried active site where the chemical reaction takes place. Despite being many times ignored, the tunnels and channels connecting the environment with the active site are equally important for the catalytic properties of enzymes. Changes in the enzymatic tunnels and channels affect enzyme activity, specificity, promiscuity, enantioselectivity and stability. This review provides an overview of the emerging field of enzyme access tunnel engineering with case studies describing design of all the aforementioned properties. The software tools for the analysis of geometry and function of the enzymatic tunnels and channels and for the rational design of tunnel modifications will also be discussed. The combination of new software tools and enzyme engineering strategies will provide enzymes with access tunnels and channels specifically tailored for individual industrial processes.

Acceleration of the autoxidation of nitric oxide by proteins

Lipoproteins and lipid membranes accelerate ^•NO autoxidation by increasing local concentration of ^•NO and O₂. Although the idea that proteins could also accelerate this reaction was presented some time ago, it was largely criticized and dismissed. Herein the effect of proteins on ^•NO autoxidation rates was studied following ^•NO disappearance with a selective electrode. It was found that human serum albumin (HSA) accelerated ^•NO autoxidation by a factor of 9 per g/mL of protein, much less than previously suggested. The acceleration by HSA was sensitive to pH and significantly decreased at pH lower than 4.5 coincident with the acid structure transition of HSA to a partially unfolded and rigid conformation. Other proteins with different surface hydrophobicity also accelerated ^•NO autoxidation and it was found to depend mostly on the protein size and dynamics. Mathematical simulations were performed to assess the physiological importance of this acceleration. It was calculated that in plasma the autoxidation of ^•NO is accelerated 1.38 times by HSA relative to water alone, but this becomes of little relevance when whole blood is simulated because of the rapid rate of ^•NO consumption by red blood cells.

Successful Mesoporous Silica Encapsulation of Optimally Functional EcDOS (E. coli Direct Oxygen Sensor), a Heme-based O2-Sensing Phosphodiesterase

The heme-based O₂ sensor from Escherichia coli, EcDOS, exerts phosphodiesterase activity towards cyclic-di-GMP (c-di-GMP), an important second messenger that regulates biofilm formation, virulence, and other important functions necessary for bacterial survival. EcDOS is a two-domain protein composed of an N-terminal heme-bound O₂-sensing domain and a C-terminal functional domain. O₂ binding to the heme Fe(II) complex in the O₂-sensing domain substantially enhances the catalytic activity of the functional domain, a property with potentially promising medical applications. Mesoporous silica is a useful material with finite-state machine-like features suitable for mediating numerous enzymatic functions. Here, we successfully encapsulated EcDOS into mesoporous silica, and demonstrated that encapsulated EcDOS was substantially activated by CO, an alternative signaling molecule used in place of O₂, exhibiting the same activity as the native enzyme in aqueous solution. Encapsulated EcDOS was sufficiently stable to exert its enzymatic function over several experimental cycles under aerobic conditions at room temperature. Thus, the present study demonstrates the successful encapsulation of the heme-based O₂ sensor EcDOS into mesoporous silica and shows that the native gas-stimulated function of EcDOS is well conserved. As such, this represents the first application of mesoporous silica to an oxygen-sensing-or any gas-sensing-enzyme.

Defective mitochondrial protease LonP1 can cause classical mitochondrial disease

LonP1 is a mitochondrial matrix protease whose selective substrate specificity is essential for maintaining mitochondrial homeostasis. Recessively inherited, pathogenic defects in LonP1 have been previously reported to underlie cerebral, ocular, dental, auricular and skeletal anomalies (CODAS) syndrome, a complex multisystemic and developmental disorder. Intriguingly, although classical mitochondrial disease presentations are well-known to exhibit marked clinical heterogeneity, the skeletal and dental features associated with CODAS syndrome are pathognomonic. We have applied whole exome sequencing to a patient with congenital lactic acidosis, muscle weakness, profound deficiencies in mitochondrial oxidative phosphorylation associated with loss of mtDNA copy number and MRI abnormalities consistent with Leigh syndrome, identifying biallelic variants in the LONP1 (NM_004793.3) gene; c.1693T > C predicting p.(Tyr565His) and c.2197G > A predicting p.(Glu733Lys); no evidence of the classical skeletal or dental defects observed in CODAS syndrome patients were noted in our patient. In vitro experiments confirmed the p.(Tyr565His) LonP1 mutant alone could not bind or degrade a substrate, consistent with the predicted function of Tyr565, whilst a second missense [p.(Glu733Lys)] variant had minimal effect. Mixtures of p.(Tyr565His) mutant and wild-type LonP1 retained partial protease activity but this was severely depleted when the p.(Tyr565His) mutant was mixed with the p.(Glu733Lys) mutant, data consistent with the compound heterozygosity detected in our patient. In summary, we conclude that pathogenic LONP1 variants can lead to a classical mitochondrial disease presentations associated with severe biochemical defects in oxidative phosphorylation in clinically relevant tissues.

Resolving polymorphs and radiation-driven effects in microcrystals using fixed-target serial synchrotron crystallography

The ability to determine high-quality, artefact-free structures is a challenge in micro-crystallography, and the rapid onset of radiation damage and requirement for a high-brilliance X-ray beam mean that a multi-crystal approach is essential. However, the combination of crystal-to-crystal variation and X-ray-induced changes can make the formation of a final complete data set challenging; this is particularly true in the case of metalloproteins, where X-ray-induced changes occur rapidly and at the active site. An approach is described that allows the resolution, separation and structure determination of crystal polymorphs, and the tracking of radiation damage in microcrystals. Within the microcrystal population of copper nitrite reductase, two polymorphs with different unit-cell sizes were successfully separated to determine two independent structures, and an X-ray-driven change between these polymorphs was followed. This was achieved through the determination of multiple serial structures from microcrystals using a high-throughput high-speed fixed-target approach coupled with robust data processing.

Uncovering dehydration in cytochrome c refolding from urea- and guanidine hydrochloride-denatured unfolded state by high pressure spectroscopy

To investigate the dehydration associated with protein folding, the partial molar volume changes for protein unfolding (ΔV_u) in cytochrome c (Cyt c) were determined using high pressure absorption spectroscopy. ΔV_u values for the unfolding to urea- and guanidine hydrochloride (GdnHCl)-denatured Cyt c were estimated to be 56±5 and 29±1 mL mol⁻¹, respectively. Considering that the volume change for hydration of hydrophobic groups is positive and that Cyt c has a covalently bonded heme, a positive ΔV_u reflects the primary contribution of the hydration of heme. Because of the marked tendency of guanidium ions to interact with hydrophobic groups, a smaller number of water molecules were hydrated with hydrophobic groups in GdnHCl-denatured Cyt c than in urea-denatured Cyt c, resulting in the smaller positive ΔV_u. On the other hand, urea is a relatively weak denaturant and urea-denatured Cyt c is not completely hydrated, which retains the partially folded structures. To unfold such partial structures, we introduced a mutation near the heme binding site, His26, to Gln, resulting in a negatively shifted ΔV_u (4±2 mL mol⁻¹) in urea-denatured Cyt c. The formation of the more solvated and less structured state in the urea-denatured mutant enhanced hydration to the hydrophilic groups in the unfolding process. Therefore, we confirmed the hydration of amino acid residues in the protein unfolding of Cyt c by estimating ΔV_u, which allows us to discuss the hydrated structures in the denatured states of proteins.

Toxin ζ Reduces the ATP and Modulates the Uridine Diphosphate-N-acetylglucosamine Pool

Toxin ζ expression triggers a reversible state of dormancy, diminishes the pool of purine nucleotides, promotes (p)ppGpp synthesis, phosphorylates a fraction of the peptidoglycan precursor uridine diphosphate-N-acetylglucosamine (UNAG), leading to unreactive UNAG-P, induces persistence in a reduced subpopulation, and sensitizes cells to different antibiotics. Here, we combined computational analyses with biochemical experiments to examine the mechanism of toxin ζ action. Free ζ toxin showed low affinity for UNAG. Toxin ζ bound to UNAG hydrolyzed ATP·Mg2+, with the accumulation of ADP, Pi, and produced low levels of phosphorylated UNAG (UNAG-P). Toxin ζ, which has a large ATP binding pocket, may temporally favor ATP binding in a position that is distant from UNAG, hindering UNAG phosphorylation upon ATP hydrolysis. The residues D67, E116, R158 and R171, involved in the interaction with metal, ATP, and UNAG, were essential for the toxic and ATPase activities of toxin ζ; whereas the E100 and T128 residues were partially dispensable. The results indicate that ζ bound to UNAG reduces the ATP concentration, which indirectly induces a reversible dormant state, and modulates the pool of UNAG.

Size-Selective Hydroformylation by a Rhodium Catalyst Confined in a Supramolecular Cage

Size-selective hydroformylation of terminal alkenes was attained upon embedding a rhodium bisphosphine complex in a supramolecular metal-organic cage that was formed by subcomponent self-assembly. The catalyst was bound in the cage by a ligand-template approach, in which pyridyl-zinc(II) porphyrin interactions led to high association constants (>105  m-1 ) for the binding of the ligands and the corresponding rhodium complex. DFT calculations confirm that the second coordination sphere forces the encapsulated active species to adopt the ee coordination geometry (i.e., both phosphine ligands in equatorial positions), in line with in situ high-pressure IR studies of the host-guest complex. The window aperture of the cage decreases slightly upon binding the catalyst. As a result, the diffusion of larger substrates into the cage is slower compared to that of smaller substrates. Consequently, the encapsulated rhodium catalyst displays substrate selectivity, converting smaller substrates faster to the corresponding aldehydes. This selectivity bears a resemblance to an effect observed in nature, where enzymes are able to discriminate between substrates based on shape and size by embedding the active site deep inside the hydrophobic pocket of a bulky protein structure.

Reversible multi-electron storage in dual-site redox-active supramolecular cages

Redox-active M₆L₄ cages display multiple reversible redox-events, enabling switching from overall +12 to −4 charged species with reversible storage of 16 electrons.

Alpha-hemolysin nanopore allows discrimination of the microcystins variants

Microcystins (MCs) are a class of cyclic heptapeptides with more than 100 variants produced by cyanobacteria present in surface waters. MCs are potent hepatotoxic agents responsible for fatal poisoning in animals and humans. Several techniques are employed in the detection of MCs, however, there is a shortage of methods capable of discriminating variants of MCs. In this work we demonstrate that the α-hemolysin (αHL) nanopore can detect and discriminate the variants (LR, YR and RR) of MCs in aqueous solution. The discrimination process is based on the analysis of the residence times of each variant of MCs within the unitary nanopore, as well as, on the amplitudes of the blockages in the ionic current flowing through it. Simulations of molecular dynamics and calculation of the electrostatic potential revealed that the variants of MCs present different charge distribution and correlated with the three patterns on the amplitudes of the blockages in the ionic current. Additionally, molecular docking analysis indicates different patterns of interaction of the variants of MCs with two specific regions of the nanopore. We conclude that αHL nanopore can discriminate variants of microcystins by a mechanism based mainly on electrostatic interaction. Finally, we propose the use of nanopore-based technology as a promising method for analyzing microcystins in aqueous solutions.

Microcystins (MCs) are a class of cyclic heptapeptides with more than 100 variants produced by cyanobacteria present in surface waters.

A comprehensive analysis of the protein-ligand interactions in crystal structures of Mycobacterium tuberculosis EthR

The Mycobacterium tuberculosis EthR is a member of the TetR family of repressors, controlling the expression of EthA, a mono-oxygenase responsible for the bioactivation of the prodrug ethionamide. This protein was established as a promising therapeutic target against tuberculosis, allowing, when inhibited by a drug-like molecule, to boost the action of ethionamide. Dozens of EthR crystal structures have been solved in complex with ligands. Herein, we disclose EthR structures in complex with 18 different small molecules and then performed in-depth analysis on the complete set of EthR structures that provides insights on EthR-ligand interactions. The 81 molecules solved in complex with EthR show a large diversity of chemical structures that were split up into several chemical clusters. Two of the most striking common points of EthR-ligand interactions are the quasi-omnipresence of a hydrogen bond bridging compounds with Asn¹⁷⁹ and the high occurrence of π-π interactions involving Phe¹¹⁰. A systematic analysis of the protein-ligand contacts identified eight hot spot residues that defined the basic structural features governing the binding mode of small molecules to EthR. Implications for the design of new potent inhibitors are discussed.

Influence of pH and Salts on Partial Molar Volume of Lysozyme and Bovine Serum Albumin in Aqueous Solutions

The partial molar volume of lysozyme and bovine serum albumin in aqueous solutions at different pH values and in aqueous solutions containing sodium chloride, ammonium chloride, sodium sulfate, or ammonium sulfate at different concentrations at pH 7.0 was investigated experimentally at 298.15 K and 1 bar. It was found that the influence of the pH value and the salts on the partial molar volume of the proteins is small, but trends were measurable. Furthermore, the partial molar volume of lysozyme in pure water at different pH values and in aqueous solutions with different sodium chloride concentrations at pH 7.0 was predicted by molecular simulations. The predictions are in good agreement with the experimental data.

Carboxypeptidase in prolyl oligopeptidase family: Unique enzyme activation and substrate-screening mechanisms

Serine peptidases of the prolyl oligopeptidase (POP) family are of substantial therapeutic importance because of their involvement in diseases such as diabetes, cancer, neurological diseases, and autoimmune disorders. Proper annotation and knowledge of substrate specificity mechanisms in this family are highly valuable. Although endopeptidase, dipeptidyl peptidase, tripeptidyl peptidase, and acylaminoacyl peptidase activities have been reported previously, here we report the first instance of carboxypeptidase activity in a POP family member. We determined the crystal structures of this carboxypeptidase, an S9C subfamily member from Deinococcus radiodurans, in its active and inactive states at 2.3-Å resolution, providing an unprecedented view of assembly and disassembly of the active site mediated by an arginine residue. We observed that this residue is poised to bind substrate in the active structure and disrupts the catalytic triad in the inactive structure. The assembly of the active site is accompanied by the ordering of gating loops, which reduces the effective size of the oligomeric pore. This prevents the entry of larger peptides and constitutes a novel mechanism for substrate screening. Furthermore, we observed structural adaptations that enable its carboxypeptidase activity, with a unique loop and two arginine residues in the active site cavity orienting the peptide substrate for catalysis. Using these structural features, we identified homologs of this enzyme in the POP family and confirmed the presence of carboxypeptidase activity in one of them. In conclusion, we have identified a new type within POP enzymes that exhibits not only unique activity but also a novel substrate-screening mechanism.

Photocatalytic Hydrogen Evolution by a Synthetic [FeFe] Hydrogenase Mimic Encapsulated in a Porphyrin Cage

The design of a biomimetic and fully base metal photocatalytic system for photocatalytic proton reduction in a homogeneous medium is described. A synthetic pyridylphosphole-appended [FeFe] hydrogenase mimic was encapsulated inside a supramolecular zinc porphyrin-based metal-organic cage structure Fe4 (Zn-L)6 . The binding is driven by the selective pyridine-zinc porphyrin interaction and results in the catalyst being bound strongly inside the hydrophobic cavity of the cage. Excitation of the capsule-forming porphyrin ligands with visible light while probing the IR spectrum confirmed that electron transfer takes place from the excited porphyrin cage to the catalyst residing inside the capsule. Light-driven proton reduction was achieved by irradiation of an acidic solution of the caged catalyst with visible light.

Loss of Function in Zeaxanthin Epoxidase of Dunaliella tertiolecta Caused by a Single Amino Acid Mutation within the Substrate-Binding Site

The zea1 mutant of marine microalga Dunaliella tertiolecta accumulates zeaxanthin under normal growth conditions, and its phenotype has been speculated to be related to zeaxanthin epoxidase (ZEP). In this study, we isolated the ZEP gene from both wild-type D. tertiolecta and the mutant. We found that the zea1 mutant has a point mutation of the 1337th nucleotide of the ZEP sequence (a change from guanine to adenine), resulting in a change of glycine to aspartate in a highly conserved region in the catalytic domain. Similar expression levels of ZEP mRNA and protein in both wild-type and zea1 were confirmed by using qRT-PCR and western blot analysis, respectively. Additionally, the enzyme activity analysis of ZEPs in the presence of cofactors showed that the inactivation of ZEP in zea1 was not caused by deficiency in the levels of cofactors. From the predicted three-dimensional ZEP structure of zea1, we observed a conformational change on the substrate-binding site in the ZEP. A comparative analysis of the ZEP structures suggested that the conformational change induced by a single amino acid mutation might impact the interaction between the substrate and substrate-binding site, resulting in loss of zeaxanthin epoxidase function.

Functional analysis of a triplet deletion in the gene encoding the sodium glucose transporter 3, a potential risk factor for ADHD

Sodium-glucose transporters (SGLT) belong to the solute carrier 5 family, which is characterized by sodium dependent transport of sugars and other solutes. In contrast, the human SGLT3 (hSGLT3) isoform, encoded by SLC5A4, acts as a glucose sensor that does not transport sugar but induces membrane depolarization by Na⁺ currents upon ligand binding. Whole-exome sequencing (WES) of several extended pedigrees with high density of attention-deficit/hyperactivity disorder (ADHD) identified a triplet ATG deletion in SLC5A4 leading to a single amino acid loss (ΔM500) in the hSGLT3 protein imperfectly co-segregating with the clinical phenotype of ADHD. Since mutations in homologous domains of hSGLT1 and hSGLT2 were found to affect intestinal and renal function, respectively, we analyzed the functional properties of hSGLT3[wt] and [ΔM500] by voltage clamp and current clamp recordings from cRNA-injected Xenopus laevis oocytes.

The cation conductance of hSGLT3[wt] was activated by application of glucose or the specific agonist 1-desoxynojirimycin (DNJ) as revealed by inward currents in the voltage clamp configuration and cell depolarization in the current clamp mode. Almost no currents and changes in membrane potential were observed when glucose or DNJ were applied to hSGLT3[ΔM500]-injected oocytes, demonstrating a loss of function by this amino acid deletion in hSGLT3. To monitor membrane targeting of wt and mutant hSGLT3, fusion constructs with YFP were generated, heterologously expressed in Xenopus laevis oocytes and analyzed for membrane fluorescence by confocal microscopy. In comparison to hSGLT3[wt] the fluorescent signal of mutant [ΔM500] was reduced by 43% indicating that the mutant phenotype might mainly result from inaccurate membrane targeting. As revealed by homology modeling, residue M500 is located in TM11 suggesting that in addition to the core structure (TM1-TM10) of the transporter, the surrounding TMs are equally crucial for transport/sensor function.

In conclusion, our findings indicate that the deletion [ΔM500] in hSGLT3 inhibits membrane targeting and thus largely disrupts glucose-induced sodium conductance, which may, in interaction with other ADHD risk-related gene variants, influence the risk for ADHD in deletion carriers.

Alcohol reduces muscle fatigue through atomistic interactions with nicotinic receptors

Alcohol consumption affects many organs and tissues, including skeletal muscle. However, the molecular mechanism of ethanol action on skeletal muscle remains unclear. Here, using molecular dynamics simulations and single channel recordings, we show that ethanol interacts with a negatively charged amino acid within an extracellular region of the neuromuscular nicotinic acetylcholine receptor (nAChR), thereby altering its global conformation and reducing the single channel current amplitude. Charge reversal of the negatively charged amino acid abolishes the nAChR-ethanol interaction. Moreover, using transgenic animals harboring the charge-reversal mutation, ex vivo measurements of muscle force production show that ethanol counters fatigue in wild type but not homozygous αE83K mutant animals. In accord, in vivo studies of motor coordination following ethanol administration reveal an approximately twofold improvement for wild type compared to homozygous mutant animals. Together, the converging results from molecular to animal studies suggest that ethanol counters muscle fatigue through its interaction with neuromuscular nAChRs.

Structural insights into oxidation of medium-chain fatty acids and flavanone by myxobacterial cytochrome P450 CYP267B1

Oxidative biocatalytic reactions performed by cytochrome P450 enzymes (P450s) are of high interest for the chemical and pharmaceutical industries. CYP267B1 is a P450 enzyme from myxobacterium Sorangium cellulosum So ce56 displaying a broad substrate scope. In this work, a search for new substrates was performed, combined with product characterization and a structural analysis of substrate-bound complexes using X-ray crystallography and computational docking. The results demonstrate the ability of CYP267B1 to perform in-chain hydroxylations of medium-chain saturated fatty acids (decanoic acid, dodecanoic acid and tetradecanoic acid) and a regioselective hydroxylation of flavanone. The fatty acids are mono-hydroxylated at different in-chain positions, with decanoic acid displaying the highest regioselectivity towards ω-3 hydroxylation. Flavanone is preferably oxidized to 3-hydroxyflavanone. High-resolution crystal structures of CYP267B1 revealed a very spacious active site pocket, similarly to other P450s capable of converting macrocyclic compounds. The pocket becomes more constricted near to the heme and is closed off from solvent by residues of the F and G helices and the B-C loop. The crystal structure of the tetradecanoic acid-bound complex displays the fatty acid bound near to the heme, but in a nonproductive conformation. Molecular docking allowed modeling of the productive binding modes for the four investigated fatty acids and flavanone, as well as of two substrates identified in a previous study (diclofenac and ibuprofen), explaining the observed product profiles. The obtained structures of CYP267B1 thus serve as a valuable prediction tool for substrate hydroxylations by this highly versatile enzyme and will encourage future selectivity changes by rational protein engineering.

Refinement and Analysis of the Mature Zika Virus Cryo-EM Structure at 3.1 Å Resolution

Among the several arthropod-borne human flaviviral diseases, the recent outbreak of Zika virus (ZIKV) has caused devastating birth defects and neurological disorders, challenging the world with another major public health concern. We report here the refined structure of the mature ZIKV at a resolution of 3.1 Å as determined by cryo-electron microscopic single-particle reconstruction. The improvement in the resolution, compared with previous enveloped virus structures, was the result of optimized virus preparation methods and data processing techniques. The glycoprotein interactions and surface properties of ZIKV were compared with other mosquito-borne flavivirus structures. The largest structural differences and sequence variations occur at the glycosylation loop associated with receptor binding. Probable drug binding pockets were identified on the viral surface. These results also provide a structural basis for the design of vaccines against ZIKV.

Nanosecond Tracer Diffusion as a Probe of the Solution Structure and Molecular Mobility of Protein Assemblies: The Case of Ovalbumin

Protein diffusion is not only an important process ensuring biological function but can also be used as a probe to obtain information on structural properties of protein assemblies in liquid solutions. Here, we explore the oligomerization state of ovalbumin at high protein concentrations by means of its short-time self-diffusion. We employ high-resolution incoherent quasielastic neutron scattering to access the self-diffusion on nanosecond timescales, on which interparticle contacts are not altered. Our results indicate that ovalbumin in aqueous (D₂O) solutions occurs in increasingly large assemblies of its monomeric subunits with rising protein concentration. It changes from nearly monomeric toward dimeric and ultimately larger than tetrameric complexes. Simultaneously, we access information on the internal molecular mobility of ovalbumin on the nanometer length scale and compare it with results obtained for bovine serum albumin, immunoglobulin, and β-lactoglobulin.

Differences in substrate specificity of V. cholerae FabH enzymes suggest new approaches for the development of novel antibiotics and biofuels

Vibrio cholerae, the causative pathogen of the life-threatening infection cholera, encodes two copies of β-ketoacyl-acyl carrier protein synthase III (vcFabH1 and vcFabH2). vcFabH1 and vcFabH2 are pathogenic proteins associated with fatty acid synthesis, lipid metabolism, and potential applications in biofuel production. Our biochemical assays characterize vcFabH1 as exhibiting specificity for acetyl-CoA and CoA thioesters with short acyl chains, similar to that observed for FabH homologs found in most gram-negative bacteria. vcFabH2 prefers medium chain-length acyl-CoA thioesters, particularly octanoyl-CoA, which is a pattern of specificity rarely seen in bacteria. Structural characterization of one vcFabH1 and six vcFabH2 structures determined in either apo form or in complex with acetyl-CoA/octanoyl-CoA indicate that the substrate-binding pockets of vcFabH1 and vcFabH2 are of different sizes, accounting for variations in substrate chain-length specificity. An unusual and unique feature of vcFabH2 is its C-terminal fragment that interacts with both the substrate-entrance loop and the dimer interface of the enzyme. Our discovery of the pattern of substrate specificity of both vcFabH1 and vcFabH2 can potentially aid the development of novel antibacterial agents against V. cholerae. Additionally, the distinctive substrate preference of FabH2 in V. cholerae and related facultative anaerobes conceivably make it an attractive component of genetically engineered bacteria used for commercial biofuel production.

Structural and Evolutionary Insights within the Polysaccharide Deacetylase Gene Family of Bacillus anthracis and Bacillus cereus

Functional and folding constraints impose interdependence between interacting sites along the protein chain that are envisaged through protein sequence evolution. Studying the influence of structure in phylogenetic models requires detailed and reliable structural models. Polysaccharide deacetylases (PDAs), members of the carbohydrate esterase family 4, perform mainly metal-dependent deacetylation of O- or N-acetylated polysaccharides such as peptidoglycan, chitin and acetylxylan through a conserved catalytic core termed the NodB homology domain. Genomes of Bacillus anthracis and its relative Bacillus cereus contain multiple genes of putative or known PDAs. A comparison of the functional domains of the recently determined PDAs from B. anthracis and B. cereus and multiple amino acid and nucleotide sequence alignments and phylogenetic analysis performed on these closely related species showed that there were distinct differences in binding site formation, despite the high conservation on the protein sequence, the folding level and the active site assembly. This may indicate that, subject to biochemical verification, the binding site-forming sequence fragments are under functionally driven evolutionary pressure to accommodate and recognize distinct polysaccharide residues according to cell location, use, or environment. Finally, we discuss the suggestion of the paralogous nature of at least two genes of B. anthracis, ba0330 and ba0331, via specific differences in gene sequence, protein structure, selection pressure and available localization patterns. This study may contribute to understanding the mechanisms under which sequences evolve in their structures and how evolutionary processes enable structural variations.

Structural basis for cross-reactivity and conformation fluctuation of the major beech pollen allergen Fag s 1

Fag s 1 is a member of the Pathogen Related protein family 10 (PR-10) and can elicit cross-reaction with IgE antibodies produced against the birch pollen allergen Bet v 1. The Nuclear Magnetic Resonance (NMR) structure of Fag s 1 is presented along with its dynamic properties. It shares 66% identity with Bet v 1 and exhibits the expected three α-helices and seven β-sheets arranged as a semi-beta barrel and exposing the residues mapped as the Bet v 1 IgE epitope. The structural dynamics of Fag s 1 were monitored on the fast and intermediate timescales, using relaxation rates. The complex dynamics of Fag s 1 are closely related to the internal cavity, and they modulate IgE and ligand binding.

Effects of protein-crystal hydration and temperature on side-chain conformational heterogeneity in monoclinic lysozyme crystals

The modulation of main-chain and side-chain conformational heterogeneity and solvent structure in monoclinic lysozyme crystals by dehydration (related to water activity) and temperature is examined. Decreasing the relative humidity (from 99 to 11%) and decreasing the temperature both lead to contraction of the unit cell, to an increased area of crystal contacts and to remodeling of primarily contact and solvent-exposed residues. Both lead to the depopulation of some minor side-chain conformers and to the generation of new conformations. Side-chain modifications and main-chain r.m.s.d.s associated with cooling from 298 to 100 K depend on relative humidity and are minimized at 85% relative humidity (r.h.). Dehydration from 99 to 93% r.h. and cooling from 298 to 100 K result in a comparable number of remodeled residues, with dehydration-induced remodeling somewhat more likely to arise from contact interactions. When scaled to equivalent temperatures based on unit-cell contraction, the evolution of side-chain order parameters with dehydration shows generally similar features to those observed on cooling to T = 100 K. These results illuminate the qualitative and quantitative similarities between structural perturbations induced by modest dehydration, which routinely occurs in samples prepared for 298 and 100 K data collection, and cryocooling. Differences between these perturbations in terms of energy landscapes and occupancies, and implications for variable-temperature crystallography between 180 and 298 K, are discussed. It is also noted that remodeling of a key lysozyme active-site residue by dehydration, which is associated with a radical decrease in the enzymatic activity of lysozyme powder, arises due to a steric clash with the residue of a symmetry mate.

A multi-scale approach to describe electrical impulses propagating along actin filaments in both intracellular and in vitro conditions

An accurate and efficient characterization of the polyelectrolyte properties for cytoskeleton filaments are key to the molecular understanding of electrical signal propagation, bundle and network formation, as well as their potential nanotechnological applications. In this article, we introduce an innovative multi-scale approach able to account for the atomistic details of a protein molecular structure, its biological environment, and their impact on electrical impulses propagating along wild type actin filaments. The formulation includes non-trivial contributions to the ionic electrical conductivity and capacitance coming from the diffuse part of the electrical double layer of G-actins. We utilize this monomer characterization in a non-linear inhomogeneous transmission line prototype model to account for the monomer–monomer interactions, dissipation and damping perturbations along the filament length. A novel, simple, accurate, approximate analytic expression has been obtained for the transmission line model. Our results reveal the propagation of electrical signal impulses in the form of solitons for the range of voltage stimulus and electrolyte solutions typically present for intracellular and in vitro conditions. The approach predicts a lower electrical conductivity with higher linear capacitance and non-linear accumulation of charge for intracellular conditions. Our results show a significant influence of the voltage input on the electrical impulse shape, attenuation and kern propagation velocity. The filament is able to sustain the soliton propagation at almost constant kern velocity for the in vitro condition, whereas the intracellular condition displays a remarkable deceleration. Additionally, the solitons are narrower and travel faster at higher voltage input. As a unique feature, this multi-scale theory is able to account for molecular structure conformation (mutation) and biological environment (protonations/deprotonations) changes often present in pathological conditions. It is also applicable to other highly charged rod-like polyelectrolytes with relevance in biomedicine and biophysics.

An innovative analytic solution accounting for the molecular structure, its biological environment, and their impact on electrical impulses along microfilaments.

Anti-CTLA-4 therapy requires an Fc domain for efficacy

Ipilimumab, an antibody that recognizes cytotoxic T lymphocyte antigen (CTLA)-4, was the first approved “checkpoint”-blocking anticancer therapy. In mice, the response to antibodies against CTLA-4 depends entirely on expression of the Fcγ receptor. We developed H11, an alpaca heavy chain-only antibody fragment against CTLA-4 that lacks an Fc portion and inhibits interactions between CTLA-4 and its ligand. By using H11 to visualize CTLA-4 expression in the whole animal, we found that accessible CTLA-4 is largely confined to the tumor; however, H11 treatment has minimal effects on antitumor responses. Installing the murine IgG2a constant region on H11 greatly enhances antitumor response. We were thus able to dissociate CTLA-4 blockade from CTLA-4–dependent receptor engagement as an explanation for the antitumor effect.

Ipilimumab, a monoclonal antibody that recognizes cytotoxic T lymphocyte antigen (CTLA)-4, was the first approved “checkpoint”-blocking anticancer therapy. In mouse tumor models, the response to antibodies against CTLA-4 depends entirely on expression of the Fcγ receptor (FcγR), which may facilitate antibody-dependent cellular phagocytosis, but the contribution of simple CTLA-4 blockade remains unknown. To understand the role of CTLA-4 blockade in the complete absence of Fc-dependent functions, we developed H11, a high-affinity alpaca heavy chain-only antibody fragment (VHH) against CTLA-4. The VHH H11 lacks an Fc portion, binds monovalently to CTLA-4, and inhibits interactions between CTLA-4 and its ligand by occluding the ligand-binding motif on CTLA-4 as shown crystallographically. We used H11 to visualize CTLA-4 expression in vivo using whole-animal immuno-PET, finding that surface-accessible CTLA-4 is largely confined to the tumor microenvironment. Despite this, H11-mediated CTLA-4 blockade has minimal effects on antitumor responses. Installation of the murine IgG2a constant region on H11 dramatically enhances its antitumor response. Coadministration of the monovalent H11 VHH blocks the efficacy of a full-sized therapeutic antibody. We were thus able to demonstrate that CTLA-4–binding antibodies require an Fc domain for antitumor effect.

Crystal structure of an intramembranal phosphatase central to bacterial cell-wall peptidoglycan biosynthesis and lipid recycling

Undecaprenyl pyrophosphate phosphatase (UppP) is an integral membrane protein that recycles the lipid carrier essential to the ongoing biosynthesis of the bacterial cell wall. Individual building blocks of peptidoglycan are assembled in the cytoplasm on undecaprenyl phosphate (C55-P) before being flipped to the periplasmic face, where they are polymerized and transferred to the existing cell wall sacculus, resulting in the side product undecaprenyl pyrophosphate (C55-PP). Interruption of UppP’s regeneration of C55-P from C55-PP leads to the buildup of cell wall intermediates and cell lysis. We present the crystal structure of UppP from Escherichia coli at 2.0 Å resolution, which reveals the mechanistic basis for intramembranal phosphatase action and substrate specificity using an inverted topology repeat. In addition, the observation of key structural motifs common to a variety of cross membrane transporters hints at a potential flippase function in the specific relocalization of the C55-P product back to the cytosolic space.

Undecaprenyl pyrophosphate phosphatase (UppP) recycles the lipid carrier essential for bacterial cell wall synthesis. Here authors present the crystal structure of UppP from E. coli at 2.0 Å resolution, which sheds light on its phosphatase mechanism and indicates a potential flippase role for UppP.

The 3.5-Å CryoEM Structure of Nanodisc-Reconstituted Yeast Vacuolar ATPase Vo Proton Channel

The molecular mechanism of transmembrane proton translocation in rotary motor ATPases is not fully understood. Here, we report the 3.5-Å resolution cryoEM structure of the lipid nanodisc-reconstituted V_o proton channel of the yeast vacuolar H⁺-ATPase, captured in a physiologically relevant, autoinhibited state. The resulting atomic model provides structural detail for the amino acids that constitute the proton pathway at the interface of the proteolipid ring and subunit a. Based on the structure and previous mutagenesis studies, we propose the chemical basis of transmembrane proton transport. Moreover, we discovered that the C terminus of the assembly factor Voa1 is an integral component of mature V_o. Voa1's C-terminal transmembrane α helix is bound inside the proteolipid ring, where it contributes to the stability of the complex. Our structure rationalizes possible mechanisms by which mutations in human V_o can result in disease phenotypes and may thus provide new avenues for therapeutic interventions.

Protein Fluctuations and Cavity Changes Relationship

Protein cavities and tunnels are critical for function. Ligand recognition and binding, transport, and enzyme catalysis require cavities rearrangements. Therefore, the flexibility of cavities should be guaranteed by protein vibrational dynamics. Molecular dynamics simulations provide a framework to explore conformational plasticity of protein cavities. Herein, we present a novel procedure to characterize the dynamics of protein cavities in terms of their volume gradient vector. For this purpose, we make use of algorithms for calculation of the cavity volume that result robust for numerical differentiations. Volume gradient vector is expressed in terms of principal component analysis obtained from equilibrated molecular dynamics simulations. We analyze contributions of principal component modes to the volume gradient vector according to their frequency and degree of delocalization. In all our test cases, we find that low frequency modes play a critical role together with minor contributions of high frequency modes. These modes involve concerted motions of significant fractions of the total residues lining the cavities. We make use of variations of the potential energy of a protein in the direction of the volume gradient vector as a measure of flexibility of the cavity. We show that proteins whose collective low frequency fluctuations contribute the most to changes of cavity volume exhibit more flexible cavities.

Direct examination of the relevance for folding, binding and electron transfer of a conserved protein folding intermediate

Near the minimum free energy basin of proteins where the native ensemble resides, partly unfolded conformations of slightly higher energy can be significantly populated under native conditions. It has been speculated that they play roles in molecular recognition and catalysis, but they might represent contemporary features of the evolutionary process without functional relevance. Obtaining conclusive evidence on these alternatives is difficult because it requires comparing the performance of a given protein when populating and when not populating one such intermediate, in otherwise identical conditions. Wild type apoflavodoxin populates under native conditions a partly unfolded conformation (10% of molecules) whose unstructured region includes the binding sites for the FMN cofactor and for redox partner proteins. We recently engineered a thermostable variant where the intermediate is no longer detectable. Using the wild type and variant, we assess the relevance of the intermediate comparing folding kinetics, cofactor binding kinetics, cofactor affinity, X-ray structure, intrinsic dynamics, redox potential of the apoflavodoxin-cofactor complex (Fld), its affinity for partner protein FNR, and electron transfer rate within the Fld/FNR physiological complex. Our data strongly suggest the intermediate state, conserved in long-chain apoflavodoxins, is not required for the correct assembly of flavodoxin nor does it contribute to shape its electron transfer properties. This analysis can be applied to evaluate other native basin intermediates.

Anomalous pressure effects on the photoreaction of a light-sensor protein from Synechocystis, PixD (Slr1694), and the compressibility change of its intermediates

SyPixD (Slr1694) is a blue-light receptor that contains a BLUF (blue-light sensor using a flavin chromophore) domain for the function of phototaxis. The key reaction of this protein is a light-induced conformational change and subsequent dissociation reaction from the decamer to the dimer. In this study, anomalous effects of pressure on this reaction were discovered, and changes in the compressibility of its short-lived intermediates were investigated. While the absorption spectra of the dark and light states are not sensitive to pressure, the formation yield of the first intermediate decreases with pressure to about 40% at 150 MPa. Upon blue-light illumination with a sufficiently strong intensity, the transient grating signal, which represents the dissociation of the SyPixD decamer, was observed at 0.1 MPa, and the signal intensity significantly decreased with increasing pressure. This behavior shows that the dissociation of the decamer from the second intermediate state is suppressed by pressure. However, while the decamer undergoes no dissociation upon excitation of one monomer unit at 0.1 MPa, dissociation is gradually induced with increasing pressure. For solving this strange behavior, the compressibility changes of the intermediates were measured as a function of pressure at weak light intensity. Interestingly, the compressibility change was negative at low pressure, but became positive with increasing pressure. Because the compressibility is related to the volume fluctuation, this observation suggests that the driving force for this reaction is fluctuation of the protein. The relationship between the cavities at the interfaces of the monomer units and the reactivity was also discussed.

Modular [FeIII8MII6] n+ (MII = Pd, Co, Ni, Cu) Coordination Cages

The reaction of the simple metalloligand [Fe^IIIL₃] [HL = 1-(4-pyridyl)butane-1,3-dione] with a variety of different M^II salts results in the formation of a family of heterometallic cages of formulae [Fe^III₈Pd^II₆L₂₄]Cl₁₂ (1), [Fe^III₈Cu^II₆L₂₄(H₂O)₄Br₄]Br₈ (2), [Fe^III₈Cu^II₆L₂₄(H₂O)₁₀](NO₃)₁₂ (3), [Fe^III₈Ni^II₆L₂₄(SCN)₁₁Cl] (4), and [Fe^III₈Co^II₆L₂₄(SCN)₁₀(H₂O)₂]Cl₂ (5). The metallic skeleton of each cage describes a cube in which the Fe^III ions occupy the eight vertices and the M^II ions lie at the center of the six faces. Direct-current magnetic susceptibility and magnetization measurements on 3-5 reveal the presence of weak antiferromagnetic exchange between the metal ions in all three cases. Computational techniques known in theoretical nuclear physics as statistical spectroscopy, which exploit the moments of the Hamiltonian to calculate relevant thermodynamic properties, determine J_Fe-Cu = 0.10 cm^-1 for 3 and J_Fe-Ni = 0.025 cm^-1 for 4. Q-band electron paramagnetic resonance spectra of 1 reveal a significantly wider spectral width in comparison to [FeL₃], indicating that the magnitude of the Fe^III zero-field splitting is larger in the heterometallic cage than in the monomer.

Coiled coil protein origami: from modular design principles towards biotechnological applications

This review illustrates the current state in designing coiled-coil-based proteins with an emphasis on coiled coil protein origami structures and their potential.