Accurate bond and angle parameters for X-ray protein structure refinement

Sultam based Carbonic Anhydrase VII inhibitors for the management of neuropathic pain

We report a series of compounds 1-17 derived from the antiepileptic drug Sulthiame (SLT) from which both the benzenesulfonamide and the sultam moiety were retained. All compounds were tested in vitro for their inhibition activity against the human (h) Carbonic Anhydrase (CA; EC 4.2.1.1) I, II, VII, IX and XII isoforms. Among the series, derivatives 1 and 11 showed great enhancement of both inhibition potency and selectivity towards the hCA VII isoform, when compared to the reference SLT drug. The binding mode of 11 within the hCA VII active site was deciphered by means of X-ray crystallography and revealed the sultam moiety being exposed to the rim of the active site. In vivo experiments on a model of neuropathic pain induced by oxaliplatin clearly showed 11 being an effective pain relieving agent and therefore worth of further exploitation towards the validation of the hCA VII as new target for the management of neuropathies.

Towards Consistency in Geometry Restraints for Carbohydrates in the Pyranose form: Modern Dictionary Generators Reviewed

Macromolecular restrained refinement is nowadays the most used method for improving the agreement between an atomic structural model and experimental data. Restraint dictionaries, a key tool behind the success of the method, allow fine-tuning geometric properties such as distances and angles between atoms beyond simplistic expectations. Dictionary generators can provide restraint target estimates derived from different sources, from fully theoretical to experimental and any combination in between. Carbohydrates are stereochemically complex biomolecules and, in their pyranose form, have clear conformational preferences. As such, they pose unique problems to dictionary generators and in the course of this study, require special attention from software developers. Functional differences between restraint generators will be discussed, as well as the process of achieving consistent results with different software designs. The study will conclude a set of practical considerations, as well as recommendations for the generation of new restraint dictionaries, using the improved software alternatives discussed.

Machine learning to estimate the local quality of protein crystal structures

Low-resolution electron density maps can pose a major obstacle in the determination and use of protein structures. Herein, we describe a novel method, called quality assessment based on an electron density map (QAEmap), which evaluates local protein structures determined by X-ray crystallography and could be applied to correct structural errors using low-resolution maps. QAEmap uses a three-dimensional deep convolutional neural network with electron density maps and their corresponding coordinates as input and predicts the correlation between the local structure and putative high-resolution experimental electron density map. This correlation could be used as a metric to modify the structure. Further, we propose that this method may be applied to evaluate ligand binding, which can be difficult to determine at low resolution.

Solvent Effects in the Crystalline Sponge Method: Importance of Co-solvents for Ordering Absorbed Guests

In the crystalline sponge method, the crucial step for ordering the absorbed guest is soaking of the guest into the pores of the crystalline sponge. Here, we find that the choice of solvent is particularly important for smooth guest soaking and ordering. Moderately polar solvents, such as ketones and esters, which we have previously avoided for the guest-soaking process, efficiently promote diffusion and guest ordering by filling the gaps in the pores through co-crystallization with the guests. Using this modified protocol, we successfully demonstrate the structural analysis of various steroids.

Ten things I `hate' about refinement

Macromolecular refinement is an optimization process that aims to produce the most likely macromolecular structural model in the light of experimental data. As such, macromolecular refinement is one of the most complex optimization problems in wide use. Macromolecular refinement programs have to deal with the complex relationship between the parameters of the atomic model and the experimental data, as well as a large number of types of prior knowledge about chemical structure. This paper draws attention to areas of unfinished business in the field of macromolecular refinement. In it, we describe ten refinement topics that we think deserve attention and discuss directions leading to macromolecular refinement software that would make the best use of modern computer resources to meet the needs of structural biologists of the twenty-first century.

Structural basis for the Pr-Pfr long-range signaling mechanism of a full-length bacterial phytochrome at the atomic level

Phytochromes constitute a widespread photoreceptor family that typically interconverts between two photostates called Pr (red light–absorbing) and Pfr (far-red light–absorbing). The lack of full-length structures solved at the (near-)atomic level in both pure Pr and Pfr states leaves gaps in the structural mechanisms involved in the signal transmission pathways during the photoconversion. Here, we present the crystallographic structures of three versions from the plant pathogen Xanthomonas campestris virulence regulator XccBphP bacteriophytochrome, including two full-length proteins, in the Pr and Pfr states. The structures show a reorganization of the interaction networks within and around the chromophore-binding pocket, an α-helix/β-sheet tongue transition, and specific domain reorientations, along with interchanging kinks and breaks at the helical spine as a result of the photoswitching, which subsequently affect the quaternary assembly. These structural findings, combined with multidisciplinary studies, allow us to describe the signaling mechanism of a full-length bacterial phytochrome at the atomic level.

Insights on the structure-function relationship of human multidrug resistance protein 7 (MRP7/ABCC10) from molecular dynamics simulations and docking studies

ATP-binding cassette (ABC) transporters superfamily mediates multidrug resistance in cancer by extruding structurally distinct chemotherapeutic agents, causing failure in chemotherapy. Among the 49 ABC transporters, multidrug resistance protein 7 (MRP7 or ABCC10) is relatively new and has been identified as the efflux pump of multiple anticancer agents including Vinca alkaloids and taxanes. Herein, we construct and validate a homology model for human MRP7 based on the cryo-EM structures of MRP1. Structure-function relationship of MRP7 was obtained from molecular dynamics simulations and docking studies and was in accordance with previous studies of ABC transporters. The motion patterns correlated with efflux mechanism were discussed. Additionally, predicted substrate- and modulator-binding sites of MRP7 were described for the first time, which provided rational insights in understanding the drug binding and functional regulation in MRP7. Our findings will benefit the high-throughput virtual screening and development of MRP7 modulators in the future.

Protein conformational transitions explored by a morphing approach based on normal mode analysis in internal coordinates

Large-scale conformational changes are essential for proteins to function properly. Given that these transition events rarely occur, however, it is challenging to comprehend their underlying mechanisms through experimental and theoretical approaches. In this study, we propose a new computational methodology called internal coordinate normal mode-guided elastic network interpolation (ICONGENI) to predict conformational transition pathways in proteins. Its basic approach is to sample intermediate conformations by interpolating the interatomic distance between two end-point conformations with the degrees of freedom constrained by the low-frequency dynamics afforded by normal mode analysis in internal coordinates. For validation of ICONGENI, it is applied to proteins that undergo open-closed transitions, and the simulation results (i.e., simulated transition pathways) are compared with those of another technique, to demonstrate that ICONGENI can explore highly reliable pathways in terms of thermal and chemical stability. Furthermore, we generate an ensemble of transition pathways through ICONGENI and investigate the possibility of using this method to reveal the transition mechanisms even when there are unknown metastable states on rough energy landscapes.

Molecular insights into the Y-domain of hepatitis E virus using computational analyses

Background

Hepatitis E virus (HEV) of the family Hepeviridae is a major causative agent of acute hepatitis in developing countries. The Y-domain is derived from multi-domain non-structural polyprotein encoded by open reading frame 1 (ORF1). Previous studies have demonstrated the essentiality of Y-domain sequences in HEV life cycle; however, its function remains completely unexplored. The following study was thus conceptualized to examine the detailed computational investigation for the putative Y-domain to estimate its phylogenetic assessment, physiochemical properties, structural and functional characteristics using in silico analyses.

Results

The phylogenetic assessment of Y-domain with a vast range of hosts indicated that the protein was very well conserved throughout the course of evolution. The Y-domain was found to be unstable, hydrophilic and basic in nature with high thermostability value. Structural analysis of Y-domain revealed mixed α/β structural fold of the protein having higher percentage of alpha-helices. The three-dimensional (3D) protein model generated through homology modelling revealed the presence of clefts, tunnels and pore. Gene ontology analysis predicted Y-domain protein’s involvement in several binding and catalytic activities as well as significant biological processes. Mutations in the conserved amino acids of the Y-domain suggested that it may stabilize or de-stabilize the protein structure that might affect its structure–function relationship.

Conclusions

This theoretical study will facilitate towards deciphering the role of unexplored Y-domain, thereby providing better understanding towards the pathogenesis of HEV infection.

DNA crosslinking and recombination-activating genes 1/2 (RAG1/2) are required for oncogenic splicing in acute lymphoblastic leukemia

Background

Abnormal alternative splicing is frequently associated with carcinogenesis. In B‐cell acute lymphoblastic leukemia (B‐ALL), double homeobox 4 fused with immunoglobulin heavy chain (DUX4/IGH) can lead to the aberrant production of E‐26 transformation‐specific family related gene abnormal transcript (ERG_alt) and other splicing variants. However, the molecular mechanism underpinning this process remains elusive. Here, we aimed to know how DUX4/IGH triggers abnormal splicing in leukemia.

Methods

The differential intron retention analysis was conducted to identify novel DUX4/IGH‐driven splicing in B‐ALL patients. X‐ray crystallography, small angle X‐ray scattering (SAXS), and analytical ultracentrifugation were used to investigate how DUX4/IGH recognize double DUX4 responsive element (DRE)‐DRE sites. The ERG_alt biogenesis and B‐cell differentiation assays were performed to characterize the DUX4/IGH crosslinking activity. To check whether recombination‐activating gene 1/2 (RAG1/2) was required for DUX4/IGH‐driven splicing, the proximity ligation assay, co‐immunoprecipitation, mammalian two hybrid characterizations, in vitro RAG1/2 cleavage, and shRNA knock‐down assays were performed.

Results

We reported previously unrecognized intron retention events in C‐type lectin domain family 12, member A abnormal transcript (CLEC12A_alt) and chromosome 6 open reading frame 89 abnormal transcript (C6orf89_alt), where also harbored repetitive DRE‐DRE sites. Supportively, X‐ray crystallography and SAXS characterization revealed that DUX4 homeobox domain (HD)1‐HD2 might dimerize into a dumbbell‐shape trans configuration to crosslink two adjacent DRE sites. Impaired DUX4/IGH‐mediated crosslinking abolishes ERG_alt, CLEC12A_alt, and C6orf89_alt biogenesis, resulting in marked alleviation of its inhibitory effect on B‐cell differentiation. Furthermore, we also observed a rare RAG1/2‐mediated recombination signal sequence‐like DNA edition in DUX4/IGH target genes. Supportively, shRNA knock‐down of RAG1/2 in leukemic Reh cells consistently impaired the biogenesis of ERG_alt, CLEC12A_alt, and C6orf89_alt.

Conclusions

All these results suggest that DUX4/IGH‐driven DNA crosslinking is required for RAG1/2 recruitment onto the double tandem DRE‐DRE sites, catalyzing V(D)J‐like recombination and oncogenic splicing in acute lymphoblastic leukemia.

Loops around the Heme Pocket Have a Critical Role in the Function and Stability of BsDyP from Bacillus subtilis

Bacillus subtilis BsDyP belongs to class I of the dye-decolorizing peroxidase (DyP) family of enzymes and is an interesting biocatalyst due to its high redox potential, broad substrate spectrum and thermostability. This work reports the optimization of BsDyP using directed evolution for improved oxidation of 2,6-dimethoxyphenol, a model lignin-derived phenolic. After three rounds of evolution, one variant was identified displaying 7-fold higher catalytic rates and higher production yields as compared to the wild-type enzyme. The analysis of X-ray structures of the wild type and the evolved variant showed that the heme pocket is delimited by three long conserved loop regions and a small α helix where, incidentally, the mutations were inserted in the course of evolution. One loop in the proximal side of the heme pocket becomes more flexible in the evolved variant and the size of the active site cavity is increased, as well as the width of its mouth, resulting in an enhanced exposure of the heme to solvent. These conformational changes have a positive functional role in facilitating electron transfer from the substrate to the enzyme. However, they concomitantly resulted in decreasing the enzyme’s overall stability by 2 kcal mol⁻¹, indicating a trade-off between functionality and stability. Furthermore, the evolved variant exhibited slightly reduced thermal stability compared to the wild type. The obtained data indicate that understanding the role of loops close to the heme pocket in the catalysis and stability of DyPs is critical for the development of new and more powerful biocatalysts: loops can be modulated for tuning important DyP properties such as activity, specificity and stability.

Tandem domain structure determination based on a systematic enumeration of conformations

Protein structure determination is undergoing a change of perspective due to the larger importance taken in biology by the disordered regions of biomolecules. In such cases, the convergence criterion is more difficult to set up and the size of the conformational space is a obstacle to exhaustive exploration. A pipeline is proposed here to exhaustively sample protein conformations using backbone angle limits obtained by nuclear magnetic resonance (NMR), and then to determine the populations of conformations. The pipeline is applied to a tandem domain of the protein whirlin. An original approach, derived from a reformulation of the Distance Geometry Problem is used to enumerate the conformations of the linker connecting the two domains. Specifically designed procedure then permit to assemble the domains to the linker conformations and to optimize the tandem domain conformations with respect to two sets of NMR measurements: residual dipolar couplings and paramagnetic resonance enhancements. The relative populations of optimized conformations are finally determined by fitting small angle X-ray scattering (SAXS) data. The most populated conformation of the tandem domain is a semi-closed one, fully closed and more extended conformations being in minority, in agreement with previous observations. The SAXS and NMR data show different influences on the determination of populations.

Structural Analysis of Strigolactone-Related Gene Products

Structural knowledge of biological macromolecules is essential for understanding their function and for modifying that function by engineering. Protein crystallography is a powerful method for elucidating molecular structures of proteins, but it is essential that the investigator has a basic knowledge of good practices and of the major pitfalls in the technique. Here we describe issues specific for the case of structural studies of strigolactone (SL) receptor structure and function, and in particular the difficulties associated with capturing complexes of SL receptors with the SL hormone ligand in the crystal.

Exploring ligand dynamics in protein crystal structures with ensemble refinement

Understanding the dynamics of ligands bound to proteins is an important task in medicinal chemistry and drug design. However, the dominant technique for determining protein-ligand structures, X-ray crystallography, does not fully account for dynamics and cannot accurately describe the movements of ligands in protein binding sites. In this article, an alternative method, ensemble refinement, is used on six protein-ligand complexes with the aim of understanding the conformational diversity of ligands in protein crystal structures. The results show that ensemble refinement sometimes indicates that the flexibility of parts of the ligand and some protein side chains is larger than that which can be described by a single conformation and atomic displacement parameters. However, since the electron-density maps are comparable and Rfree values are slightly increased, the original crystal structure is still a better model from a statistical point of view. On the other hand, it is shown that molecular-dynamics simulations and automatic generation of alternative conformations in crystallographic refinement confirm that the flexibility of these groups is larger than is observed in standard refinement. Moreover, the flexible groups in ensemble refinement coincide with groups that give high atomic displacement parameters or non-unity occupancy if optimized in standard refinement. Therefore, the conformational diversity indicated by ensemble refinement seems to be qualitatively correct, indicating that ensemble refinement can be an important complement to standard crystallographic refinement as a tool to discover which parts of crystal structures may show extensive flexibility and therefore are poorly described by a single conformation. However, the diversity of the ensembles is often exaggerated (probably partly owing to the rather poor force field employed) and the ensembles should not be trusted in detail.

The missing link: covalent linkages in structural models

Covalent linkages between constituent blocks of macromolecules and ligands have been subject to inconsistent treatment during the model-building, refinement and deposition process. This may stem from a number of sources, including difficulties with initially detecting the covalent linkage, identifying the correct chemistry, obtaining an appropriate restraint dictionary and ensuring its correct application. The analysis presented herein assesses the extent of problems involving covalent linkages in the Protein Data Bank (PDB). Not only will this facilitate the remediation of existing models, but also, more importantly, it will inform and thus improve the quality of future linkages. By considering linkages of known type in the CCP4 Monomer Library (CCP4-ML), failure to model a covalent linkage is identified to result in inaccurate (systematically longer) interatomic distances. Scanning the PDB for proximal atom pairs that do not have a corresponding type in the CCP4-ML reveals a large number of commonly occurring types of unannotated potential linkages; in general, these may or may not be covalently linked. Manual consideration of the most commonly occurring cases identifies a number of genuine classes of covalent linkages. The recent expansion of the CCP4-ML is discussed, which has involved the addition of over 16 000 and the replacement of over 11 000 component dictionaries using AceDRG. As part of this effort, the CCP4-ML has also been extended using AceDRG link dictionaries for the aforementioned linkage types identified in this analysis. This will facilitate the identification of such linkage types in future modelling efforts, whilst concurrently easing the process involved in their application. The need for a universal standard for maintaining link records corresponding to covalent linkages, and references to the associated dictionaries used during modelling and refinement, following deposition to the PDB is emphasized. The importance of correctly modelling covalent linkages is demonstrated using a case study, which involves the covalent linkage of an inhibitor to the main protease in various viral species, including SARS-CoV-2. This example demonstrates the importance of properly modelling covalent linkages using a comprehensive restraint dictionary, as opposed to just using a single interatomic distance restraint or failing to model the covalent linkage at all.

Modelling covalent linkages in CCP4

In this contribution, the current protocols for modelling covalent linkages within the CCP4 suite are considered. The mechanism used for modelling covalent linkages is reviewed: the use of dictionaries for describing changes to stereochemistry as a result of the covalent linkage and the application of link-annotation records to structural models to ensure the correct treatment of individual instances of covalent linkages. Previously, linkage descriptions were lacking in quality compared with those of contemporary component dictionaries. Consequently, AceDRG has been adapted for the generation of link dictionaries of the same quality as for individual components. The approach adopted by AceDRG for the generation of link dictionaries is outlined, which includes associated modifications to the linked components. A number of tools to facilitate the practical modelling of covalent linkages available within the CCP4 suite are described, including a new restraint-dictionary accumulator, the Make Covalent Link tool and AceDRG interface in Coot, the 3D graphical editor JLigand and the mechanisms for dealing with covalent linkages in the CCP4i2 and CCP4 Cloud environments. These integrated solutions streamline and ease the covalent-linkage modelling workflow, seamlessly transferring relevant information between programs. Current recommended practice is elucidated by means of instructive practical examples. By summarizing the different approaches to modelling linkages that are available within the CCP4 suite, limitations and potential pitfalls that may be encountered are highlighted in order to raise awareness, with the intention of improving the quality of future modelled covalent linkages in macromolecular complexes.

Critical evaluation of a crystal structure of nitrogenase with bound N2 ligands

Recently, a 1.83 Å crystallographic structure of nitrogenase was suggested to show N2-derived ligands at three sites in the catalytic FeMo cluster, replacing the three [Formula: see text] bridging sulfide ligands (two in one subunit and the third in the other subunit) (Kang et al. in Science 368: 1381-1385, 2020). Naturally, such a structure is sensational, having strong bearings on the reaction mechanism of the enzyme. Therefore, it is highly important to ensure that the interpretation of the structure is correct. Here, we use standard crystallographic refinement and quantum refinement to evaluate the structure. We show that the original crystallographic raw data are strongly anisotropic, with a much lower resolution in certain directions than others. This, together with the questionable use of anisotropic B factors, give atoms an elongated shape, which may look like diatomic atoms. In terms of standard electron-density maps and real-space Z scores, a resting-state structure with no dissociated sulfide ligands fits the raw data better than the interpretation suggested by the crystallographers. The anomalous electron density at 7100 eV is weaker for the putative N2 ligands, but not lower than for several of the [Formula: see text] bridging sulfide ions and not lower than what can be expected from a statistical analysis of the densities. Therefore, we find no convincing evidence for any N2 binding to the FeMo cluster. Instead, a standard resting state without any dissociated ligands seems to be the most likely interpretation of the structure. Likewise, we find no support that the homocitrate ligand should show monodentate binding.

NMR solution structures of Runella slithyformis RNA 2'-phosphotransferase Tpt1 provide insights into NAD+ binding and specificity

Tpt1, an essential component of the fungal and plant tRNA splicing machinery, catalyzes transfer of an internal RNA 2′-PO₄ to NAD⁺ yielding RNA 2′-OH and ADP-ribose-1′,2′-cyclic phosphate products. Here, we report NMR structures of the Tpt1 ortholog from the bacterium Runella slithyformis (RslTpt1), as apoenzyme and bound to NAD⁺. RslTpt1 consists of N- and C-terminal lobes with substantial inter-lobe dynamics in the free and NAD⁺-bound states. ITC measurements of RslTpt1 binding to NAD⁺ (K_D ∼31 μM), ADP-ribose (∼96 μM) and ADP (∼123 μM) indicate that substrate affinity is determined primarily by the ADP moiety; no binding of NMN or nicotinamide is observed by ITC. NAD⁺-induced chemical shift perturbations (CSPs) localize exclusively to the RslTpt1 C-lobe. NADP⁺, which contains an adenylate 2′-PO₄ (mimicking the substrate RNA 2′-PO₄), binds with lower affinity (K_D ∼1 mM) and elicits only N-lobe CSPs. The RslTpt1·NAD⁺ binary complex reveals C-lobe contacts to adenosine ribose hydroxyls (His99, Thr101), the adenine nucleobase (Asn105, Asp112, Gly113, Met117) and the nicotinamide riboside (Ser125, Gln126, Asn163, Val165), several of which are essential for RslTpt1 activity in vivo. Proximity of the NAD⁺ β-phosphate to ribose-C1″ suggests that it may stabilize an oxocarbenium transition-state during the first step of the Tpt1-catalyzed reaction.

Probing the role of Val228 on the catalytic activity of Scytalidium catalase

Scytalidium catalase is a homotetramer including heme d in each subunit. Its primary function is the dismutation of H₂O₂ to water and oxygen, but it is also able to oxidase various small organic compounds including catechol and phenol. The crystal structure of Scytalidium catalase reveals the presence of three linked channels providing access to the exterior like other catalases reported so far. The function of these channels has been extensively studied, revealing the possible routes for substrate flow and product release. In this report, we have focussed on the semi-conserved residue Val228, located near to the vinyl groups of the heme at the opening of the lateral channel. Its replacement with Ala, Ser, Gly, Cys, Phe and Ile were tested. We observed a significant decrease in catalytic efficiency in all mutants with the exception of a remarkable increase in oxidase activity when Val228 was mutated to either Ala, Gly or Ser. The reduced catalytic efficiencies are characterized in terms of the restriction of hydrogen peroxide as electron acceptor in the active centre resulting from the opening of lateral channel inlet by introducing the smaller side chain residues. On the other hand, the increased oxidase activity is explained by allowing the suitable electron donor to approach more closely to the heme. The crystal structures of V228C and V228I were determined at 1.41 and 1.47 Å resolution, respectively. The lateral channels of the V228C and V228I presented a broadly identical chain of arranged waters to that observed for wild-type enzyme.

Pr-favoured variants of the bacteriophytochrome from the plant pathogen Xanthomonas campestris hint on light regulation of virulence-associated mechanisms

Red/far-red light-sensing bacteriophytochrome photoreceptor (BphP) pathways play key roles in bacterial physiology and ecology. These bilin-binding proteins photoswitch between two states, Pr (red absorbing) and Pfr (far-red absorbing). The isomerization of the chromophore and the downstream structural changes result in the light signal transduction. The agricultural pathogen Xanthomonas campestris pv. campestris (Xcc) code for a single bathy-like type BphP (XccBphP), previously shown to negatively regulate several light-mediated biological processes involved in virulence. Here, we generated three different full-length variants with single amino acid changes within its GAF domain that affect the XccBphP photocycle favouring its Pr state: L193Q, L193N and D199A. While D199A recombinant protein locks XccBphP in a Pr-like state, L193Q and L193N exhibit a significant enrichment of the Pr form in thermal equilibrium. The X-ray crystal structures of the three variants were solved, resembling the wild-type protein in the Pr state. Finally, we studied the effects of altering the XccBphP photocycle on the exopolysaccharide xanthan production and stomatal aperture assays as readouts of its bacterial signalling pathway. Null-mutant complementation assays show that the photoactive Pr-favoured XccBphP variants L193Q and L193N tend to negatively regulate xanthan production in vivo. In addition, our results indicate that strains expressing these variants also promote stomatal apertures in challenged plant epidermal peels, compared to wild-type Xcc. The findings presented in this work provide new evidence on the Pr state of XccBphP as a negative regulator of the virulence-associated mechanisms by light in Xcc.

Macromolecular refinement of X-ray and cryoelectron microscopy structures with Phenix/OPLS3e for improved structure and ligand quality

With the advent of the resolution revolution in cryoelectron microscopy (cryo-EM), low-resolution refinement is common, and likewise increases the need for a reliable force field. Here, we report on the incorporation of the OPLS3e force field with the VSGB2.1 solvation model in the structure determination package Phenix. Our results show significantly improved structure quality and reduced ligand strain at lower resolution for X-ray refinement. For refinement of cryo-EM-based structures, we find comparable quality structures, goodness-of-fit, and reduced ligand strain. We also show how structure quality and ligand strain are related to the map-model cross-correlation as a function of data weight, and how that can detect overfitting. Signs of overfitting are found in over half of our cryo-EM dataset, which can be remedied by a re-refinement at a lower data weight. Finally, a start-to-end script for refining structures with Phenix/OPLS3e is available in the Schrödinger 2020-3 distribution.

Crystal structure of the GTP-binding protein-like domain of AGAP1

AGAP1 is often considered to regulate membrane trafficking, protein transport and actin cytoskeleton dynamics. Recent studies have shown that aberrant expression of AGAP1 is associated with many diseases, including neurodevelopmental disorders and acute lymphoblastic leukemia. It has been proposed that the GTP-binding protein-like domain (GLD) is involved in the binding of cofactors and thus regulates the catalytic activity of AGAP1. To obtain a better understanding of the pathogenic mechanism underpinning AGAP1-related diseases, it is essential to obtain structural information. Here, the GLD (residues 70-235) of AGAP1 was overexpressed in Escherichia coli BL21 (DE3) cells. Affinity and gel-filtration chromatography were used to obtain AGAP1GLD with high purity for crystallization. Using the hanging-drop vapor-diffusion method with the protein at a final concentration of 20 mg ml-1, AGAP1GLD protein crystals of suitable size were obtained. The crystals were found to diffract to 3.0 Å resolution and belonged to space group I4, with unit-cell parameters a = 100.39, b = 100.39, c = 48.08 Å. The structure of AGAP1GLD exhibits the highly conserved functional G1-G5 loops and is generally similar to other characterized ADP-ribosylation factor (Arf) GTPase-activating proteins (GAPs), implying an analogous function to Arf GAPs. Additionally, this study indicates that AGAP1 could be classified as a type of NTPase, the activity of which might be regulated by protein partners or by its other domains. Taken together, these results provide insight into the regulatory mechanisms of AGAP1 in cell signaling.

Flat, Cα,β -Didehydroalanine Foldamers with Ferrocene Pendants: Assessing the Role of α-Peptide Dipolar Moments

The foldamer field is continuously expanding as it allows to produce molecules endowed with 3D-structures and functions never observed in nature. We synthesized flat foldamers based on the natural, but non-coded, Cα,β -didehydroalanine α-amino acid, and covalently linked to them two ferrocene (Fc) moieties, as redox probes. These conjugates retain the flat and extended conformation of the 2.05 -helix, both in solution and in the crystal state (X-ray diffraction). Cyclic voltammetry measurements agree with the adoption of the 2.05 -helix, characterized by a negligible dipole moment. Thus, elongated α-peptide stretches of this type are insulators rather than charge conductors, the latter being constituted by peptide α-helices. Also, our homo-tetrapeptide has a N-to-C length of about 18.2 Å, almost double than that (9.7 Å) of an α-helical α-tetrapeptide.

Protein Structure Prediction Using Population-Based Algorithm Guided by Information Entropy

Ab initio protein structure prediction is one of the most challenging problems in computational biology. Multistage algorithms are widely used in ab initio protein structure prediction. The different computational costs of a multistage algorithm for different proteins are important to be considered. In this study, a population-based algorithm guided by information entropy (PAIE), which includes exploration and exploitation stages, is proposed for protein structure prediction. In PAIE, an entropy-based stage switch strategy is designed to switch from the exploration stage to the exploitation stage. Torsion angle statistical information is also deduced from the first stage and employed to enhance the exploitation in the second stage. Results indicate that an improvement in the performance of protein structure prediction in a benchmark of 30 proteins and 17 other free modeling targets in CASP.

Structural interpretation of cryo-EM image reconstructions

The productivity of single-particle cryo-EM as a structure determination method has rapidly increased as many novel biological structures are being elucidated. The ultimate result of the cryo-EM experiment is an atomic model that should faithfully represent the computed image reconstruction. Although the principal approach of atomic model building and refinement from maps resembles that of the X-ray crystallographic methods, there are important differences due to the unique properties resulting from the 3D image reconstructions. In this review, we discuss the practiced work-flow from the cryo-EM image reconstruction to the atomic model. We give an overview of (i) resolution determination methods in cryo-EM including local and directional resolution variation, (ii) cryo-EM map contrast optimization including complementary map types that can help in identifying ambiguous density features, (iii) atomic model building and (iv) refinement in various resolution regimes including (v) their validation and (vi) discuss differences between X-ray and cryo-EM maps. Based on the methods originally developed for X-ray crystallography, the path from 3D image reconstruction to atomic coordinates has become an integral and important part of the cryo-EM structure determination work-flow that routinely delivers atomic models.

Benzyl alcohol inhibits carbonic anhydrases by anchoring to the zinc coordinated water molecule

Up to date alcohols have been scarcely investigated as carbonic anhydrase (CA) inhibitors. To get more insights into the CA inhibition properties of this class of molecules, in this paper, by means of inhibition assays and X-ray crystallographic studies we report a detailed characterization of the CA inhibition properties and the binding mode to human CA II of benzyl alcohol. Results show that, although possessing a very simple scaffold, this molecule acts as a micromolar CA II inhibitor, which anchors to the enzyme active site by means of an H-bond interaction with the zinc bound solvent molecule. Taken together our results clearly indicate primary alcohols as a class of CA inhibitors that deserve to be more investigated.

Iron uptake and transport by the carboxymycobactin-mycobactin siderophore machinery of Mycobacterium tuberculosis is dependent on the iron-regulated protein HupB

Iron-starved Mycobacterium tuberculosis utilises the carboxymycobactin-mycobactin siderophore machinery to acquire iron. These two siderophores have high affinity for ferric iron and can withdraw the metal ion from insoluble iron hydroxides and iron-binding proteins. We first reported HupB, a multi-functional mycobacterial protein to be associated with iron acquisition in M. tuberculosis. This 28 kDa cell wall protein, up regulated upon iron limitation functions as a transcriptional activator of mycobactin biosynthesis and is essential for the pathogen to survive inside macrophages. The focus of this study is to understand the role of HupB in iron uptake and transport by the carboxmycobactin-mycobactin siderophore machinery in M. tuberculosis. Experimental approaches included radiolabelled iron uptake studies by viable organisms and protein-ligand binding studies using the purified HupB and the two siderophores. Uptake of ⁵⁵Fe-carboxymycobactin by wild type M. tuberculosis (WT M.tb.H37Rv) and not by the hupB KO mutant (M.tb.ΔhupB) showed that HupB is necessary for the uptake of ferri-carboxymycobactin. Additionally, the radiolabel recovery was high in HupB-incorporated liposomes upon addition of the labelled siderophore. Bioinformatic and experimental studies using spectrofluorimetry, CD analysis and surface plasmon resonance not only confirmed the binding of HupB with ferri-carboxymycobactin and ferri-mycobactin but also with free iron. In conclusion, HupB is established as a ferri- carboxymycobactin receptor and by virtue of its property to bind ferric iron, functions as a transporter of the ferric iron from the extracellular siderophore to mycobactin within the cell envelope.

High-resolution structural profile of hylaseptin-4: Aggregation, membrane topology and pH dependence of overall membrane binding process

Hylaseptin-4 (HSP-4, GIGDILKNLAKAAGKAALHAVGESL-NH₂) is an antimicrobial peptide originally isolated from Hypsiboas punctatus tree frog. The peptide has been chemically synthetized for structural investigations by CD and NMR spectroscopies. CD experiments reveal the high helical content of HSP-4 in biomimetic media. Interestingly, the aggregation process seems to occur at high peptide concentrations either in aqueous solution or in presence of biomimetic membranes, indicating an increase in the propensity of the peptide for adopting a helical conformation. High-resolution NMR structures determined in presence of DPC-d₃₈ micelles show a highly ordered α-helix from amino acid residues I2 to S24 and a smooth bend near G14. A large separation between hydrophobic and hydrophilic residues occurs up to the A16 residue, from which a shift in the amphipathicity is noticed. Oriented solid-state NMR spectroscopy show a roughly parallel orientation of the helical structure along the POPC lipid bilayer surface, with an insertion of the hydrophobic N-terminus into the bilayer core. Moreover, a noticeable pH dependence of the aggregation process in both aqueous and in biomimetic membrane environments is attributed to a single histidine residue (H19). The protonation degree of the imidazole side-chain might help in modulating the peptide-peptide or peptide-lipid interactions. Finally, molecular dynamics simulations confirm the orientation and preferential helical conformation and in addition, show that HSP-4 tends to self-aggregate in order to stabilize its active conformation in aqueous or phospholipid bilayer environments.

Precise control of synthetic hydrogel network structure via linear, independent synthesis-swelling relationships

Hydrogel physical properties are tuned by altering synthesis conditions such as initial polymer concentration and polymer-cross-linker stoichiometric ratios. Traditionally, differences in hydrogel synthesis schemes, such as end-linked poly(ethylene glycol) diacrylate hydrogels and cross-linked poly(vinyl alcohol) hydrogels, limit structural comparison between hydrogels. In this study, we use generalized synthesis variables for hydrogels that emphasize how changes in formulation affect the resulting network structure. We identify two independent linear correlations between these synthesis variables and swelling behavior. Analysis through recently updated swollen polymer network models suggests that synthesis-swelling correlations can be used to make a priori predictions of the stiffness and solute diffusivity characteristics of synthetic hydrogels. The same experiments and analyses performed on methacrylamide-modified gelatin hydrogels demonstrate that complex biopolymer structures disrupt the linear synthesis-swelling correlations. These studies provide insight into the control of hydrogel physical properties through structural design and can be used to implement and optimize biomedically relevant hydrogels.

A complete Fourier-synthesis-based backbone-conformation-dependent library for proteins

While broadening the applicability of (φ/ψ)-dependent target values for the bond angles in the peptide backbone, sequence/conformation categories with too few residues to analyze via previous methods were encountered. Here, a method of describing a conformation-dependent library (CDL) using two-dimensional Fourier coefficients is reported where the number of coefficients for individual categories is determined via complete cross-validation. Sample sizes are increased further by selective blending of categories with similar patterns of conformational dependence. An additional advantage of the Fourier-synthesis-based CDL is that it uses continuous functions and has no artifactual steps near the edges of populated regions of φ/ψ space. A set of libraries for the seven main-chain bond angles, along with the ω and ζ angles, was created based on a set of Fourier analyses of 48 368 residues selected from high-resolution models in the wwPDB. This new library encompasses both trans- and cis-peptide bonds and outperforms currently used discrete CDLs.

Transferred-Rotational-Echo Double Resonance

Internuclear distance determination is the foundation for NMR-based structure calculation. However, high-precision distance measurement is a laborious process requiring lengthy data acquisitions due to the large set of multidimensional spectra needed at different mixing times. This prevents application to large or challenging molecular systems. Here, we present a new approach, transferred-rotational-echo double resonance (TREDOR), a heteronuclear transfer method in which we simultaneously detect both starting and transferred signals in a single spectrum. This co-acquisition is used to compensate for coherence decay, resulting in accurate and precise distance determination by a single parameter fit using a single spectrum recorded at an ideal mixing time. We showcase TREDOR with the microcrystalline SH3 protein using 3D spectra to resolve resonances. By combining the measured N-C and H-C distances, we calculate the structure of SH3, which converges to the correct fold, with a root-mean-square deviation of 2.1 Å compared to a reference X-ray structure. The TREDOR data used in the structure calculation were acquired in only 4 days on a 600 MHz instrument. This is achieved due to the more than 2-fold time saving afforded by co-acquisition of additional information and demonstrates TREDOR as a fast and straightforward method for determining structures via magic-angle spinning NMR.

Heme Binding to HupZ with a C-Terminal Tag from Group A Streptococcus

HupZ is an expected heme degrading enzyme in the heme acquisition and utilization pathway in Group A Streptococcus. The isolated HupZ protein containing a C-terminal V5-His6 tag exhibits a weak heme degradation activity. Here, we revisited and characterized the HupZ-V5-His6 protein via biochemical, mutagenesis, protein quaternary structure, UV-vis, EPR, and resonance Raman spectroscopies. The results show that the ferric heme-protein complex did not display an expected ferric EPR signal and that heme binding to HupZ triggered the formation of higher oligomeric states. We found that heme binding to HupZ was an O2-dependent process. The single histidine residue in the HupZ sequence, His111, did not bind to the ferric heme, nor was it involved with the weak heme-degradation activity. Our results do not favor the heme oxygenase assignment because of the slow binding of heme and the newly discovered association of the weak heme degradation activity with the His6-tag. Altogether, the data suggest that the protein binds heme by its His6-tag, resulting in a heme-induced higher-order oligomeric structure and heme stacking. This work emphasizes the importance of considering exogenous tags when interpreting experimental observations during the study of heme utilization proteins.

ProMod3-A versatile homology modelling toolbox

Computational methods for protein structure modelling are routinely used to complement experimental structure determination, thus they help to address a broad spectrum of scientific questions in biomedical research. The most accurate methods today are based on homology modelling, i.e. detecting a homologue to the desired target sequence that can be used as a template for modelling. Here we present a versatile open source homology modelling toolbox as foundation for flexible and computationally efficient modelling workflows. ProMod3 is a fully scriptable software platform that can perform all steps required to generate a protein model by homology. Its modular design aims at fast prototyping of novel algorithms and implementing flexible modelling pipelines. Common modelling tasks, such as loop modelling, sidechain modelling or generating a full protein model by homology, are provided as production ready pipelines, forming the starting point for own developments and enhancements. ProMod3 is the central software component of the widely used SWISS-MODEL web-server.

Covid-19.bioreproducibility.org: A web resource for SARS-CoV-2-related structural models

The COVID-19 pandemic has triggered numerous scientific activities aimed at understanding the SARS-CoV-2 virus and ultimately developing treatments. Structural biologists have already determined hundreds of experimental X-ray, cryo-EM, and NMR structures of proteins and nucleic acids related to this coronavirus, and this number is still growing. To help biomedical researchers, who may not necessarily be experts in structural biology, navigate through the flood of structural models, we have created an online resource, covid19.bioreproducibility.org, that aggregates expert-verified information about SARS-CoV-2-related macromolecular models. In this article, we describe this web resource along with the suite of tools and methodologies used for assessing the structures presented therein.

Diflunisal Derivatives as Modulators of ACMS Decarboxylase Targeting the Tryptophan-Kynurenine Pathway

In the kynurenine pathway for tryptophan degradation, an unstable metabolic intermediate, α-amino-β-carboxymuconate-ε-semialdehyde (ACMS), can nonenzymatically cyclize to form quinolinic acid, the precursor for de novo biosynthesis of nicotinamide adenine dinucleotide (NAD⁺). In a competing reaction, ACMS is decarboxylated by ACMS decarboxylase (ACMSD) for further metabolism and energy production. Therefore, the inhibition of ACMSD increases NAD⁺ levels. In this study, an Food and Drug Administration (FDA)-approved drug, diflunisal, was found to competitively inhibit ACMSD. The complex structure of ACMSD with diflunisal revealed a previously unknown ligand-binding mode and was consistent with the results of inhibition assays, as well as a structure-activity relationship (SAR) study. Moreover, two synthesized diflunisal derivatives showed half-maximal inhibitory concentration (IC₅₀) values 1 order of magnitude better than diflunisal at 1.32 ± 0.07 μM (22) and 3.10 ± 0.11 μM (20), respectively. The results suggest that diflunisal derivatives have the potential to modulate NAD⁺ levels. The ligand-binding mode revealed here provides a new direction for developing inhibitors of ACMSD.

Structural modelling and thermostability of a serine protease inhibitor belonging to the Kunitz-BPTI family from the Rhipicephalus microplus tick

rBmTI-A is a recombinant serine protease inhibitor that belongs to the Kunitz-BPTI family and that was cloned from Rhipicephalus microplus tick. rBmTI-A has inhibitory activities on bovine trypsin, human plasma kallikrein, human neutrophil elastase and plasmin with dissociation constants in nM range. It is characterized by two inhibitory domains and each domain presents six cysteines that form three disulfide bonds, which contribute to the high stability of its structure. Previous studies suggest that serine protease inhibitor rBmTI-A has a protective potential against pulmonary emphysema in mice and anti-inflammatory potential. Besides that, rBmTI-A presented a potent inhibitory activity against in vitro vessel formation. In this study, the tertiary structure of rBmTI-A was modeled. The structure stabilization was evaluated by molecular dynamics analysis. Circular dichroism spectroscopy data corroborated the secondary structure found by the homology modelling. Also, in circular dichroism data it was shown a thermostability of rBmTI-A until approximately 70 °C, corroborated by inhibitory assays toward trypsin.

FASPR: an open-source tool for fast and accurate protein side-chain packing

MOTIVATION

Protein structure and function are essentially determined by how the side-chain atoms interact with each other. Thus, accurate protein side-chain packing (PSCP) is a critical step toward protein structure prediction and protein design. Despite the importance of the problem, however, the accuracy and speed of current PSCP programs are still not satisfactory.

RESULTS

We present FASPR for fast and accurate PSCP by using an optimized scoring function in combination with a deterministic searching algorithm. The performance of FASPR was compared with four state-of-the-art PSCP methods (CISRR, RASP, SCATD and SCWRL4) on both native and non-native protein backbones. For the assessment on native backbones, FASPR achieved a good performance by correctly predicting 69.1% of all the side-chain dihedral angles using a stringent tolerance criterion of 20°, compared favorably with SCWRL4, CISRR, RASP and SCATD which successfully predicted 68.8%, 68.6%, 67.8% and 61.7%, respectively. Additionally, FASPR achieved the highest speed for packing the 379 test protein structures in only 34.3 s, which was significantly faster than the control methods. For the assessment on non-native backbones, FASPR showed an equivalent or better performance on I-TASSER predicted backbones and the backbones perturbed from experimental structures. Detailed analyses showed that the major advantage of FASPR lies in the optimal combination of the dead-end elimination and tree decomposition with a well optimized scoring function, which makes FASPR of practical use for both protein structure modeling and protein design studies.

AVAILABILITY AND IMPLEMENTATION

The web server, source code and datasets are freely available at https://zhanglab.ccmb.med.umich.edu/FASPR and https://github.com/tommyhuangthu/FASPR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The structure of unliganded sterol carrier protein 2 from Yarrowia lipolytica unveils a mechanism for binding site occlusion

Isolated or as a part of multidomain proteins, Sterol Carrier Protein 2 (SCP2) exhibits high affinity and broad specificity for different lipidic and hydrophobic compounds. A wealth of structural information on SCP2 domains in all forms of life is currently available; however, many aspects of its ligand binding activity are poorly understood. ylSCP2 is a well-characterized single domain SCP2 from the yeast Yarrowia lipolytica. Herein, we report the X-ray structure of unliganded ylSCP2 refined to 2.0 Å resolution. Comparison with the previously solved liganded ylSCP2 structure unveiled a novel mechanism for binding site occlusion. The liganded ylSCP2 binding site is a large cavity with a volume of more than 800 Å³. In unliganded ylSCP2 the binding site is reduced to about 140 Å³. The obliteration is caused by a swing movement of the C-terminal α helix 5 and a subtle compaction of helices 2-4. Previous pairwise comparisons were between homologous SCP2 domains with a uncertain binding status. The reported unliganded ylSCP2 structure allows for the first time a fully controlled comparative analysis of the conformational effects of ligand occupation dispelling several doubts regarding the architecture of SCP2 binding site.

Arginine off-kilter: guanidinium is not as planar as restraints denote

Crystallographic refinement of macromolecular structures relies on stereochemical restraints to mitigate the typically poor data-to-parameter ratio. For proteins, each amino acid has a unique set of geometry restraints which represent stereochemical information such as bond lengths, valence angles, torsion angles, dihedrals and planes. It has been shown that the geometry in refined structures can differ significantly from that present in libraries; for example, it was recently reported that the guanidinium moiety in arginine is not symmetric. In this work, the asymmetry of the N^ϵ-C^ζ-N^η1 and N^ϵ-C^ζ-N^η2 valence angles in the guanidinium moiety is confirmed. In addition, it was found that the C^δ atom can deviate significantly (more than 20°) from the guanidinium plane. This requires the relaxation of the planar restraint for the C^δ atom, as it otherwise causes the other atoms in the group to compensate by distorting the guanidinium core plane. A new set of restraints for the arginine side chain have therefore been formulated, and are available in the software package Phenix, that take into account the asymmetry of the group and the planar deviation of the C^δ atom. This is an example of the need to regularly revisit the geometric restraint libraries used in macromolecular refinement so that they reflect the best knowledge of the structural chemistry of their components available at the time.

Real-space quantum-based refinement for cryo-EM: Q|R#3

Electron cryo-microscopy (cryo-EM) is rapidly becoming a major competitor to X-ray crystallography, especially for large structures that are difficult or impossible to crystallize. While recent spectacular technological improvements have led to significantly higher resolution three-dimensional reconstructions, the average quality of cryo-EM maps is still at the low-resolution end of the range compared with crystallography. A long-standing challenge for atomic model refinement has been the production of stereochemically meaningful models for this resolution regime. Here, it is demonstrated that including accurate model geometry restraints derived from ab initio quantum-chemical calculations (HF-D3/6-31G) can improve the refinement of an example structure (chain A of PDB entry 3j63). The robustness of the procedure is tested for additional structures with up to 7000 atoms (PDB entry 3a5x and chain C of PDB entry 5fn5) using the less expensive semi-empirical (GFN1-xTB) model. The necessary algorithms enabling real-space quantum refinement have been implemented in the latest version of qr.refine and are described here.

Phylogenetics-based identification and characterization of a superior 2,3-butanediol dehydrogenase for Zymomonas mobilis expression

BACKGROUND

Zymomonas mobilis has recently been shown to be capable of producing the valuable platform biochemical, 2,3-butanediol (2,3-BDO). Despite this capability, the production of high titers of 2,3-BDO is restricted by several physiological parameters. One such bottleneck involves the conversion of acetoin to 2,3-BDO, a step catalyzed by 2,3-butanediol dehydrogenase (Bdh). Several Bdh enzymes have been successfully expressed in Z. mobilis, although a highly active enzyme is yet to be identified for expression in this host. Here, we report the application of a phylogenetic approach to identify and characterize a superior Bdh, followed by validation of its structural attributes using a mutagenesis approach.

RESULTS

Of the 11 distinct bdh genes that were expressed in Z. mobilis, crude extracts expressing Serratia marcescens Bdh (SmBdh) were found to have the highest activity (8.89 µmol/min/mg), when compared to other Bdh enzymes (0.34-2.87 µmol/min/mg). The SmBdh crystal structure was determined through crystallization with cofactor (NAD⁺) and substrate (acetoin) molecules bound in the active site. Active SmBdh was shown to be a tetramer with the active site populated by a Gln247 residue contributed by the diagonally opposite subunit. SmBdh showed a more extensive supporting hydrogen-bond network in comparison to the other well-studied Bdh enzymes, which enables improved substrate positioning and substrate specificity. This protein also contains a short α6 helix, which provides more efficient entry and exit of molecules from the active site, thereby contributing to enhanced substrate turnover. Extending the α6 helix to mimic the lower activity Enterobacter cloacae (EcBdh) enzyme resulted in reduction of SmBdh function to nearly 3% of the total activity. In great contrast, reduction of the corresponding α6 helix of the EcBdh to mimic the SmBdh structure resulted in ~ 70% increase in its activity.

CONCLUSIONS

This study has demonstrated that SmBdh is superior to other Bdhs for expression in Z. mobilis for 2,3-BDO production. SmBdh possesses unique structural features that confer biochemical advantage to this protein. While coordinated active site formation is a unique structural characteristic of this tetrameric complex, the smaller α6 helix and extended hydrogen network contribute towards improved activity and substrate promiscuity of the enzyme.

Quantum refinement with multiple conformations: application to the P-cluster in nitrogenase

X-ray crystallography is the main source of atomistic information on the structure of proteins. Normal crystal structures are obtained as a compromise between the X-ray scattering data and a set of empirical restraints that ensure chemically reasonable bond lengths and angles. However, such restraints are not always available or accurate for nonstandard parts of the structure, for example substrates, inhibitors and metal sites. The method of quantum refinement, in which these empirical restraints are replaced by quantum-mechanical (QM) calculations, has previously been suggested for small but interesting parts of the protein. Here, this approach is extended to allow for multiple conformations in the QM region by performing separate QM calculations for each conformation. This approach is shown to work properly and leads to improved structures in terms of electron-density maps and real-space difference density Z-scores. It is also shown that the quality of the structures can be gauged using QM strain energies. The approach, called ComQumX-2QM, is applied to the P-cluster in two different crystal structures of the enzyme nitrogenase, i.e. an Fe8S7Cys6 cluster, used for electron transfer. One structure is at a very high resolution (1.0 Å) and shows a mixture of two different oxidation states, the fully reduced PN state (Fe82+, 20%) and the doubly oxidized P2+ state (80%). In the original crystal structure the coordinates differed for only two iron ions, but here it is shown that the two states also show differences in other atoms of up to 0.7 Å. The second structure is at a more modest resolution, 2.1 Å, and was originally suggested to show only the one-electron oxidized state, P1+. Here, it is shown that it is rather a 50/50% mixture of the P1+ and P2+ states and that many of the Fe-Fe and Fe-S distances in the original structure were quite inaccurate (by up to 0.8 Å). This shows that the new ComQumX-2QM approach can be used to sort out what is actually seen in crystal structures with dual conformations and to give locally improved coordinates.

Crystal structures of peanut lectin in the presence of synthetic β-N- and β-S-galactosides disclose evidence for the recognition of different glycomimetic ligands

Carbohydrate-lectin interactions are involved in important cellular recognition processes, including viral and bacterial infections, inflammation and tumor metastasis. Hence, structural studies of lectin-synthetic glycan complexes are essential for understanding lectin-recognition processes and for the further design of promising chemotherapeutics that interfere with sugar-lectin interactions. Plant lectins are excellent models for the study of the molecular-recognition process. Among them, peanut lectin (PNA) is highly relevant in the field of glycobiology because of its specificity for β-galactosides, showing high affinity towards the Thomsen-Friedenreich antigen, a well known tumor-associated carbohydrate antigen. Given this specificity, PNA is one of the most frequently used molecular probes for the recognition of tumor cell-surface O-glycans. Thus, it has been extensively used in glycobiology for inhibition studies with a variety of β-galactoside and β-lactoside ligands. Here, crystal structures of PNA are reported in complex with six novel synthetic hydrolytically stable β-N- and β-S-galactosides. These complexes disclosed key molecular-binding interactions of the different sugars with PNA at the atomic level, revealing the roles of specific water molecules in protein-ligand recognition. Furthermore, binding-affinity studies by isothermal titration calorimetry showed dissociation-constant values in the micromolar range, as well as a positive multivalency effect in terms of affinity in the case of the divalent compounds. Taken together, this work provides a qualitative structural rationale for the upcoming synthesis of optimized glycoclusters designed for the study of lectin-mediated biological processes. The understanding of the recognition of β-N- and β-S-galactosides by PNA represents a benchmark in protein-carbohydrate interactions since they are novel synthetic ligands that do not belong to the family of O-linked glycosides.

In silico improvement of the cyanobacterial lectin microvirin and mannose interaction

Lectins that bind to HIV envelope glycoprotein can inhibit virus-cell fusion and be used for rational drug design. This paper presents the results of an in silico approach to improve affinity interaction between the cyanobacterial lectin microvirin and its ligand Manα(1-2)Man. Comparative modeling and molecular dynamics tools were used. Additionally, the alanine scanning webserver was used to study the importance of protein residues in the binding site and to guide mutant production. The model obtained presented two homologous domains designated as domains A and B, each consisting of a single strand with triple and antiparallel β-sheets of (β1-β3 and β6-β8). Disulfide bonds between the cysteines (Cys60-Cys80, Cys63-Cys78 and Cys8-Cys24) were also found. The highly conserved binding site, including residues Asn44, Ile45, Asp46, Gln54, Asn55, Glu58, Thr59, Gln81, Thr82 and Met83. The RMSD values of the di-mannose and the interaction site were very stable during the molecular dynamics. Calculations of the occupation time of the hydrogen bonds were made for the residues that showed interaction in the complex lectin and ligand. The residue that contributed most to the interaction with Manα(1-2)Man was Asn55. After validation, the model generated remained stable during the entire simulation. Despite its structural similarity with the template we used, our mutant (Thr82Arg) showed a higher affinity interaction with Manα(1-2)Man. Communicated by Ramaswamy H. Sarma.

Systematic design of biomolecular force fields

Force fields for the study of biomolecules have been developed in a predominantly organic manner by regular updates over half a century. Together with better algorithms and advances in computer technology, force fields have improved to yield more robust predictions. However, there are also indications to suggest that intramolecular energy functions have not become better and that there still is room for improvement. In this review, systematic efforts toward development of novel force fields from scratch are described. This includes an estimate of the complexity of the problem and the prerequisites in the form of data and algorithms. It is suggested that in order to make progress, an effort is needed to standardize reference data and force field validation benchmarks.

A reinforcement-learning-based approach to enhance exhaustive protein loop sampling

MOTIVATION

Loop portions in proteins are involved in many molecular interaction processes. They often exhibit a high degree of flexibility, which can be essential for their function. However, molecular modeling approaches usually represent loops using a single conformation. Although this conformation may correspond to a (meta-)stable state, it does not always provide a realistic representation.

RESULTS

In this paper, we propose a method to exhaustively sample the conformational space of protein loops. It exploits structural information encoded in a large library of three-residue fragments, and enforces loop-closure using a closed-form inverse kinematics solver. A novel reinforcement-learning-based approach is applied to accelerate sampling while preserving diversity. The performance of our method is showcased on benchmark datasets involving 9-, 12- and 15-residue loops. In addition, more detailed results presented for streptavidin illustrate the ability of the method to exhaustively sample the conformational space of loops presenting several meta-stable conformations.

AVAILABILITY AND IMPLEMENTATION

We are developing a software package called MoMA (for Molecular Motion Algorithms), which includes modeling tools and algorithms to sample conformations and transition paths of biomolecules, including the application described in this work. The binaries can be provided upon request and a web application will also be implemented in the short future.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cofactor Complexes of DesD, a Model Enzyme in the Virulence-related NIS Synthetase Family

The understudied nonribosomal-peptide-synthetase-independent siderophore (NIS) synthetase family has been increasingly associated with virulence in bacterial species due to its key role in the synthesis of hydroxamate and carboxylate "stealth" siderophores. We have identified a model family member, DesD, from Streptomyces coelicolor, to structurally characterize using a combination of a wild-type and a Arg306Gln variant in apo, cofactor product AMP-bound, and cofactor reactant ATP-bound complexes. The kinetics in the family has been limited by solubility and reporter assays, so we have developed a label-free kinetics assay utilizing a single-injection isothermal-titration-calorimetry-based method. We report second-order rate constants that are 50 times higher than the previous estimations for DesD. Our Arg306Gln DesD variant was also tested under identical buffer and substrate conditions, and its undetectable activity was confirmed. These are the first reported structures for DesD, and they describe the critical cofactor coordination. This is also the first label-free assay to unambiguously determine the kinetics for an NIS synthetase.