The Hunt for the Closed Conformation of the Fruit-Ripening Enzyme 1-Aminocyclopropane-1-carboxylic Oxidase: A Combined Electron Paramagnetic Resonance and Molecular Dynamics Study.

1-Aminocyclopropane-1-carboxylic oxidase (ACCO) is a non-heme iron(II)-containing enzyme involved in the biosynthesis of the phytohormone ethylene, which regulates fruit ripening and flowering in plants. The active conformation of ACCO, and in particular that of the C-terminal part, remains unclear and open and closed conformations have been proposed. In this work, a combined experimental and computational study to understand the conformation and dynamics of the C-terminal part is reported. Site-directed spin-labeling coupled to electron paramagnetic resonance (SDSL-EPR) spectroscopy was used. Mutagenesis experiments were performed to generate active enzymes bearing two paramagnetic labels (nitroxide radicals) anchored on cysteine residues, one in the main core and one in the C-terminal part. Inter-spin distance distributions were measured by pulsed EPR spectroscopy and compared with the results of molecular dynamics simulations. The results reveal the existence of a flexibility of the C-terminal part. This flexibility generates several conformations of the C-terminal part of ACCO that correspond neither to the existing crystal structures nor to the modelled structures. This highly dynamic region of ACCO raises questions on its exact function during enzymatic activity.

Prediction of Reduction Potentials of Copper Proteins with Continuum Electrostatics and Density Functional Theory.

Blue copper proteins, such as azurin, show dramatic changes in Cu2+ /Cu+ reduction potential upon mutation over the full physiological range. Hence, they have important functions in electron transfer and oxidation chemistry and have applications in industrial biotechnology. The details of what determines these reduction potential changes upon mutation are still unclear. Moreover, it has been difficult to model and predict the reduction potential of azurin mutants and currently no unique procedure or workflow pattern exists. Furthermore, high-level computational methods can be accurate but are too time consuming for practical use. In this work, a novel approach for calculating reduction potentials of azurin mutants is shown, based on a combination of continuum electrostatics, density functional theory and empirical hydrophobicity factors. Our method accurately reproduces experimental reduction potential changes of 30 mutants with respect to wildtype within experimental error and highlights the factors contributing to the reduction potential change. Finally, reduction potentials are predicted for a series of 124 new mutants that have not yet been investigated experimentally. Several mutants are identified that are located well over 10 Šfrom the copper center that change the reduction potential by more than 85 mV. The work shows that secondary coordination sphere mutations mostly lead to long-range electrostatic changes and hence can be modeled accurately with continuum electrostatics.

Interactions between the protein barnase and co-solutes studied by NMR.

Protein solubility and stability depend on the co-solutes present. There is little theoretical basis for selection of suitable co-solutes. Some guidance is provided by the Hofmeister series, an empirical ordering of anions according to their effect on solubility and stability; and by osmolytes, which are small organic molecules produced by cells to allow them to function in stressful environments. Here, NMR titrations of the protein barnase with Hofmeister anions and osmolytes are used to measure and locate binding, and thus to separate binding and bulk solvent effects. We describe a rationalisation of Hofmeister (and inverse Hofmeister) effects, which is similar to the traditional chaotrope/kosmotrope idea but based on solvent fluctuation rather than water withdrawal, and characterise how co-solutes affect protein stability and solubility, based on solvent fluctuations. This provides a coherent explanation for solute effects, and points towards a more rational basis for choice of excipients.

Improved methodology for protein NMR structure calculation using hydrogen bond restraints and ANSURR validation: The SH2 domain of SH2B1.

Protein structures calculated using NMR data are less accurate and less well-defined than they could be. Here we use the program ANSURR to show that this deficiency is at least in part due to a lack of hydrogen bond restraints. We describe a protocol to introduce hydrogen bond restraints into the structure calculation of the SH2 domain from SH2B1 in a systematic and transparent way and show that the structures generated are more accurate and better defined as a result. We also show that ANSURR can be used as a guide to know when the structure calculation is good enough to stop.

The accuracy of protein structures in solution determined by AlphaFold and NMR.

In the recent Critical Assessment of Structure Prediction (CASP) competition, AlphaFold2 performed outstandingly. Its worst predictions were for nuclear magnetic resonance (NMR) structures, which has two alternative explanations: either the NMR structures were poor, implying that Alpha-Fold may be more accurate than NMR, or there is a genuine difference between crystal and solution structures. Here, we use the program Accuracy of NMR Structures Using RCI and Rigidity (ANSURR), which measures the accuracy of solution structures, and show that one of the NMR structures was indeed poor. We then compare Alpha-Fold predictions to NMR structures and show that Alpha-Fold tends to be more accurate than NMR ensembles. There are, however, some cases where the NMR ensembles are more accurate. These tend to be dynamic structures, where Alpha-Fold had low confidence. We suggest that Alpha-Fold could be used as the model for NMR-structure refinements and that Alpha-Fold structures validated by ANSURR may require no further refinement.

The accuracy of NMR protein structures in the Protein Data Bank.

The program ANSURR measures the accuracy of NMR structures by comparing rigidity obtained from experimental backbone chemical shifts and from structures. We report on ANSURR analysis of 7,000 PDB NMR ensembles within the Protein Data Bank, which can be found at ansurr.com. The accuracy of NMR structures progressively improved up until 2005, but since then, it has plateaued. Most structures have accurate secondary structure, but are generally too floppy, particularly in loops. Thus, there is a need for more experimental restraints in loops. Currently, the best predictors of accuracy are Ramachandran distribution and the number of NOE restraints per residue. The precision of structures within the ensemble correlates well with accuracy, as does the number of hydrogen bond restraints per residue. Structure accuracy is improved when other components (such as additional polypeptide chains or ligands) are included.

A method for validating the accuracy of NMR protein structures.

We present a method that measures the accuracy of NMR protein structures. It compares random coil index [RCI] against local rigidity predicted by mathematical rigidity theory, calculated from NMR structures [FIRST], using a correlation score (which assesses secondary structure), and an RMSD score (which measures overall rigidity). We test its performance using: structures refined in explicit solvent, which are much better than unrefined structures; decoy structures generated for 89 NMR structures; and conventional predictors of accuracy such as number of restraints per residue, restraint violations, energy of structure, ensemble RMSD, Ramachandran distribution, and clashscore. Restraint violations and RMSD are poor measures of accuracy. Comparisons of NMR to crystal structures show that secondary structure is equally accurate, but crystal structures are typically too rigid in loops, whereas NMR structures are typically too floppy overall. We show that the method is a useful addition to existing measures of accuracy.

Conformational exchange in the potassium channel blocker ShK.

ShK is a 35-residue disulfide-linked polypeptide produced by the sea anemone Stichodactyla helianthus, which blocks the potassium channels Kv1.1 and Kv1.3 with pM affinity. An analogue of ShK has been developed that blocks Kv1.3 > 100 times more potently than Kv1.1, and has completed Phase 1b clinical trials for the treatment of autoimmune diseases such as psoriasis and rheumatoid arthritis. Previous studies have indicated that ShK undergoes a conformational exchange that is critical to its function, but this has proved difficult to characterise. Here, we have used high hydrostatic pressure as a tool to increase the population of the alternative state, which is likely to resemble the active form that binds to the Kv1.3 channel. By following changes in chemical shift with pressure, we have derived the chemical shift values of the low- and high-pressure states, and thus characterised the locations of structural changes. The main difference is in the conformation of the Cys17-Cys32 disulfide, which is likely to affect the positions of the critical Lys22-Tyr23 pair by twisting the 21-24 helix and increasing the solvent exposure of the Lys22 sidechain, as indicated by molecular dynamics simulations.

Does Substrate Positioning Affect the Selectivity and Reactivity in the Hectochlorin Biosynthesis Halogenase?

In this work we present the first computational study on the hectochlorin biosynthesis enzyme HctB, which is a unique three-domain halogenase that activates non-amino acid moieties tethered to an acyl-carrier, and as such may have biotechnological relevance beyond other halogenases. We use a combination of small cluster models and full enzyme structures calculated with quantum mechanics/molecular mechanics methods. Our work reveals that the reaction is initiated with a rate-determining hydrogen atom abstraction from substrate by an iron (IV)-oxo species, which creates an iron (III)-hydroxo intermediate. In a subsequent step the reaction can bifurcate to either halogenation or hydroxylation of substrate, but substrate binding and positioning drives the reaction to optimal substrate halogenation. Furthermore, several key residues in the protein have been identified for their involvement in charge-dipole interactions and induced electric field effects. In particular, two charged second coordination sphere amino acid residues (Glu223 and Arg245) appear to influence the charge density on the Cl ligand and push the mechanism toward halogenation. Our studies, therefore, conclude that nonheme iron halogenases have a chemical structure that induces an electric field on the active site that affects the halide and iron charge distributions and enable efficient halogenation. As such, HctB is intricately designed for a substrate halogenation and operates distinctly different from other nonheme iron halogenases.

Features of reactive cysteines discovered through computation: from kinase inhibition to enrichment around protein degrons.

Large-scale characterisation of cysteine modification is enabling study of the physicochemical determinants of reactivity. We find that location of cysteine at the amino terminus of an α-helix, associated with activity in thioredoxins, is under-represented in human protein structures, perhaps indicative of selection against background reactivity. An amino-terminal helix location underpins the covalent linkage for one class of kinase inhibitors. Cysteine targets for S-palmitoylation, S-glutathionylation, and S-nitrosylation show little correlation with pKa values predicted from structures, although flanking sequences of S-palmitoylated sites are enriched in positively-charged amino acids, which could facilitate palmitoyl group transfer to substrate cysteine. A surprisingly large fraction of modified sites, across the three modifications, would be buried in native protein structure. Furthermore, modified cysteines are (on average) closer to lysine ubiquitinations than are unmodified cysteines, indicating that cysteine redox biology could be associated with protein degradation and degron recognition.