Cholesterol oxidase: ultrahigh-resolution crystal structure and multipolar atom model-based analysis

Peptide bonds revisited

Understanding the structural and chemical properties of peptide bonds within protein secondary structures is vital for elucidating their roles in protein folding, stability and function. This study examines the distinct characteristics of peptide bonds in α-helices and β-strands using a nonredundant data set comprising 1024 high-resolution protein crystal structures from the Protein Data Bank (PDB). The analysis reveals surprising and intriguing insights into bond lengths, angles, dihedral angles, electron-density distributions and hydrogen bonding within α-helices and β-strands. While the respective bond lengths (CN and CO) do not differ much between helices and strands, the bond angles (∠CNCα and ∠OCN) are significantly larger in strands compared with helices. Furthermore, the peptide dihedral angle (ω) in helices clusters around 180° and follows a sharp Gaussian distribution with a standard deviation of 4.1°. In contrast, the distribution of dihedral angles in strands spans a much wider range, with a more flattened Gaussian peak around 180°. This distinct difference in the distribution of dihedral angles reflects the unique structural characteristics of helices and strands, highlighting their respective conformational preferences. Additionally, if the ratio of the electron-density values (2mFo - DFc) at the midpoint of the CO bond and of the CN bond is calculated, a skewed distribution is observed, with the ratio being lower for helices than for strands. Moreover, higher normalized mean atomic displacement parameters (ADPs) for peptide atoms in helices relative to strands suggest increased flexibility or a more dynamic structure within helical regions. Analysis of hydrogen-bond distances between O and N atoms of the main chain reveals larger distances in helices compared with strands, indicative of distinct hydrogen-bonding patterns associated with different secondary structures. All of these observations taken together led us to conclude that peptide bonds in α-helices are different from peptide bonds in β-strands. Overall, α-helical peptide bonds seem to display a more enol-like character. This suggests that peptide oxygen atoms in helices are more likely to be protonated. These findings have several important implications for refining protein structures, particularly in regions susceptible to enol-like transitions or protonation. By recognizing the distinct bond-angle and bond-length variations associated with protonated carbonyl oxygen atoms, current refinement protocols can be adapted to apply more flexible restraints in these regions. This could improve the accuracy of modelling local geometries, where protonation or enol forms lead to subtle structural deviations from the canonical bond parameters typically enforced in refinement strategies.

On the structure refinement of metal complexes against 3D electron diffraction data using multipolar scattering factors

This study examines various methods for modelling the electron density and, thus, the electrostatic potential of an organometallic complex for use in crystal structure refinement against 3D electron diffraction (ED) data. It focuses on modelling the scattering factors of iron(III), considering the electron density distribution specific for coordination with organic linkers. We refined the structural model of the metal-organic complex, iron(III) acetylacetonate (FeAcAc), using both the independent atom model (IAM) and the transferable aspherical atom model (TAAM). TAAM refinement initially employed multipolar parameters from the MATTS databank for acetylacetonate, while iron was modelled with a spherical and neutral approach (TAAM ligand). Later, custom-made TAAM scattering factors for Fe-O coordination were derived from DFT calculations [TAAM-ligand-Fe(III)]. Our findings show that, in this compound, the TAAM scattering factor corresponding to Fe3+ has a lower scattering amplitude than the Fe3+ charged scattering factor described by IAM. When using scattering factors corresponding to the oxidation state of iron, IAM inaccurately represents electrostatic potential maps and overestimates the scattering potential of the iron. In addition, TAAM significantly improved the fitting of the model to the data, shown by improved R1 values, goodness-of-fit (GooF) and reduced noise in the Fourier difference map (based on the residual distribution analysis). For 3D ED, R1 values improved from 19.36% (IAM) to 17.44% (TAAM-ligand) and 17.49% (TAAM-ligand-Fe3+), and for single-crystal X-ray diffraction (SCXRD) from 3.82 to 2.03% and 1.98%, respectively. For 3D ED, the most significant R1 reductions occurred in the low-resolution region (8.65-2.00 Å), dropping from 20.19% (IAM) to 14.67% and 14.89% for TAAM-ligand and TAAM-ligand-Fe(III), respectively, with less improvement in high-resolution ranges (2.00-0.85 Å). This indicates that the major enhancements are due to better scattering modelling in low-resolution zones. Furthermore, when using TAAM instead of IAM, there was a noticeable improvement in the shape of the thermal ellipsoids, which more closely resembled those of an SCXRD-refined model. This study demonstrates the applicability of more sophisticated scattering factors to improve the refinement of metal-organic complexes against 3D ED data, suggesting the need for more accurate modelling methods and highlighting the potential of TAAM in examining the charge distribution of large molecular structures using 3D ED.

Xylitol binding to the M1 site of glucose isomerase induces a conformational change in the substrate binding channel

Glucose isomerase (GI) is extensively used in the food industry for production of high-fructose corn syrup and for the production of biofuels and other renewable chemicals. Structure-based studies on GI inhibitors are important for improving its efficiency in industrial applications. Here, we report the subatomic crystal structure of Streptomyces rubiginosus GI (SruGI) complexed with its inhibitor, xylitol, at 0.99 Å resolution. Electron density map and temperature factor analysis showed partial binding of xylitol to the M1 metal binding site of SruGI, providing two different conformations of the metal binding site and the substrate binding channel. The xylitol molecule induced a conformational change in the M2 metal ion-interacting Asp255 residue, which subsequently led to a conformational change in the side chain of Asp181 residue. This led to the positional shift of Pro25 by 1.71 Å and side chain rotation of Phe26 by 21°, where located on the neighboring protomer in tetrameric SruGI. The conformation change of these two residues affect the size of the substrate-binding channel of GI. Therefore, xylitol binding to M1 site of SruGI induces not only a conformational changes of the metal-binding site, but also conformational change of substrate-binding channel of the tetrameric SruGI. These results expand our knowledge about the mechanism underlying the inhibitory effect of xylitol on GI.

Progressive alignment of crystals: reproducible and efficient assessment of crystal structure similarity

During in silico crystal structure prediction of organic molecules, millions of candidate structures are often generated. These candidates must be compared to remove duplicates prior to further analysis (e.g. optimization with electronic structure methods) and ultimately compared with structures determined experimentally. The agreement of predicted and experimental structures forms the basis of evaluating the results from the Cambridge Crystallographic Data Centre (CCDC) blind assessment of crystal structure prediction, which further motivates the pursuit of rigorous alignments. Evaluating crystal structure packings using coordinate root-mean-square deviation (RMSD) for N molecules (or N asymmetric units) in a reproducible manner requires metrics to describe the shape of the compared molecular clusters to account for alternative approaches used to prioritize selection of molecules. Described here is a flexible algorithm called Progressive Alignment of Crystals (PAC) to evaluate crystal packing similarity using coordinate RMSD and introducing the radius of gyration (R g) as a metric to quantify the shape of the superimposed clusters. It is shown that the absence of metrics to describe cluster shape adds ambiguity to the results of the CCDC blind assessments because it is not possible to determine whether the superposition algorithm has prioritized tightly packed molecular clusters (i.e. to minimize R g) or prioritized reduced RMSD (i.e. via possibly elongated clusters with relatively larger R g). For example, it is shown that when the PAC algorithm described here uses single linkage to prioritize molecules for inclusion in the superimposed clusters, the results are nearly identical to those calculated by the widely used program COMPACK. However, the lower R g values obtained by the use of average linkage are favored for molecule prioritization because the resulting RMSDs more equally reflect the importance of packing along each dimension. It is shown that the PAC algorithm is faster than COMPACK when using a single process and its utility for biomolecular crystals is demonstrated. Finally, parallel scaling up to 64 processes in the open-source code Force Field X is presented.

Electron density based analysis of N–H⋯OC hydrogen bonds and electrostatic interaction energies in high-resolution secondary protein structures: insights from quantum crystallographic approaches

Limiting values of the topological parameters and the electrostatic interaction energies to establish the presence of true N–H⋯OC H-bonds in protein main-chain have been identified using quantitative and qualitative analyses of electron densities.

Expression optimization, purification, and functional characterization of cholesterol oxidase from Chromobacterium sp. DS1

Cholesterol oxidase is a bifunctional bacterial flavoenzyme which catalyzes oxidation and isomerization of cholesterol. This valuable enzyme has attracted a great deal of attention because of its wide application in the clinical laboratory, synthesis of steroid derived drugs, food industries, and its potentially insecticidal activity. Therefore, development of an efficient protocol for overproduction of cholesterol oxidase could be valuable and beneficial in this regard. The present study examined the role of various parameters (host strain, culture media, induction time, isopropyl ß-D-1-thiogalactopyranoside concentration, as well as post-induction incubation time and temperature) on over-expression of cholesterol oxidase from Chromobacterium sp. DS1. Applying the optimized protocol, the yield of recombinant cholesterol oxidase significantly increased from 92 U/L to 2115 U/L. Under the optimized conditions, the enzyme was produced on a large-scale, and overexpressed cholesterol oxidase was purified from cell lysate by column nickel affinity chromatography. K_m and V_max values of the purified enzyme for cholesterol were estimated using Lineweaver-Burk plot. Further, the optimum pH and optimum temperature for the enzyme activity were determined. This study reports a straightforward protocol for cholesterol oxidase production which can be performed in any laboratory.

Neutron scattering in the biological sciences: progress and prospects

The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

Use of an Isotope-Coded Mass Tag (ICMT) Method To Determine the Orientation of Cholesterol Oxidase on Model Membranes

Although the interfacial membrane protein cholesterol oxidase is structurally and kinetically well-characterized, its orientation in and mode of interaction with cholesterol-containing membranes have not been established. Cholesterol oxidase can alter the structure of the cell membrane in pathogenic bacteria and is thus a potential antimicrobial drug target. We recently developed a mass spectrometry-based isotope-coded mass tag (ICMT) labeling method to monitor the real-time solvent-accessible surface of peripheral membrane proteins, such as cholesterol oxidase. The ICMT strategy utilizes maleimide-based isotope tags that covalently react with cysteine residues. In this study, by comparing the ICMT labeling rates of cysteine variants of cholesterol oxidase, we determined which residues of the protein were engaged with the protein-lipid interface. We found that upon addition of cholesterol-containing lipid vesicles, four cysteine residues in a cluster near the substrate entrance channel are labeled more slowly with ICMT probes than in the absence of vesicles, indicating that these four residues were in contact with the membrane surface. From these data, we generated a model of how cholesterol oxidase is oriented when bound to the membrane. In conclusion, this straightforward method, which requires only microgram quantities of protein, offers several advantages over existing methods for the investigation of interfacial membrane proteins and can be applied to a number of different systems.

The electrostatic potential of dynamic charge densities

A procedure to derive the electrostatic potential (ESP) for dynamic charge densities obtained from structure models or maximum-entropy densities is introduced. The ESP essentially is obtained by inverse Fourier transform of the dynamic structure factors of the total charge density corresponding to the independent atom model, the multipole model or maximum-entropy densities, employing dedicated software that will be part of the BayMEM software package. Our approach is also discussed with respect to the Ewald summation method. It is argued that a meaningful ESP can only be obtained if identical thermal smearing is applied to the nuclear (positive) and electronic (negative) parts of the dynamic charge densities. The method is applied to structure models of dl-serine at three different temperatures of 20, 100 and 298 K. The ESP at locations near the atomic nuclei exhibits a drastic reduction with increasing temperature, the largest difference between the ESP from the static charge density and the ESP of the dynamic charge density being at T = 20 K. These features demonstrate that zero-point vibrations are sufficient for changing the spiky nature of the ESP at the nuclei into finite values. On 0.5 e Å-3 isosurfaces of the electron densities (taken as the molecular surface relevant to intermolecular interactions), the dynamic ESP is surprisingly similar at all temperatures, while the static ESP of a single molecule has a slightly larger range and is shifted towards positive potential values.

Validation of X-ray Wavefunction Refinement

In this work, the quality of the electron density in crystals reconstructed by the multipolar model (MM) and by X-ray wavefunction refinement (XWR) is tested on a set of high-resolution X-ray diffraction data sets of four amino acids and six tripeptides. It results in the first thorough validation of XWR. Agreement statistics, figures of merit, residual- and deformation-density maps, as well as atomic displacement parameters are used to measure the quality of the reconstruction relative to the measured structure factors. Topological analysis of the reconstructed density is carried out to obtain atomic and bond-topological properties, which are subsequently compared to the values derived from benchmarking periodic DFT geometry optimizations. XWR is simultaneously in better agreement than the MM with both benchmarking theory and the measured diffraction pattern. In particular, the obvious problems with the description of polar bonds in the MM are significantly reduced by using XWR. Similarly, modeling of electron density in the vicinity of hydrogen atoms with XWR is visibly improved.

A theoretical-electron-density databank using a model of real and virtual spherical atoms

A database describing the electron density of common chemical groups using combinations of real and virtual spherical atoms is proposed, as an alternative to the multipolar atom modelling of the molecular charge density. Theoretical structure factors were computed from periodic density functional theory calculations on 38 crystal structures of small molecules and the charge density was subsequently refined using a density model based on real spherical atoms and additional dummy charges on the covalent bonds and on electron lone-pair sites. The electron-density parameters of real and dummy atoms present in a similar chemical environment were averaged on all the molecules studied to build a database of transferable spherical atoms. Compared with the now-popular databases of transferable multipolar parameters, the spherical charge modelling needs fewer parameters to describe the molecular electron density and can be more easily incorporated in molecular modelling software for the computation of electrostatic properties. The construction method of the database is described. In order to analyse to what extent this modelling method can be used to derive meaningful molecular properties, it has been applied to the urea molecule and to biotin/streptavidin, a protein/ligand complex.

Distribution of valence electrons of the flavin cofactor in NADH-cytochrome b5 reductase

Flavin compounds such as flavin adenine dinucleotide (FAD), flavin mononucleotide and riboflavin make up the active centers in flavoproteins that facilitate various oxidoreductive processes. The fine structural features of the hydrogens and valence electrons of the flavin molecules in the protein environment are critical to the functions of the flavoproteins. However, information on these features cannot be obtained from conventional protein X-ray analyses at ordinary resolution. Here we report the charge density analysis of a flavoenzyme, NADH-cytochrome b₅ reductase (b5R), at an ultra-high resolution of 0.78 Å. Valence electrons on the FAD cofactor as well as the peptide portion, which are clearly visualized even after the conventional refinement, are analyzed by the multipolar atomic model refinement. The topological analysis for the determined electron density reveals the valence electronic structure of the isoalloxazine ring of FAD and hydrogen-bonding interactions with the protein environment. The tetrahedral electronic distribution around the N5 atom of FAD in b5R is stabilized by hydrogen bonding with C_αH of Tyr65 and amide-H of Thr66. The hydrogen bonding network leads to His49 composing the cytochrome b₅-binding site via non-classical hydrogen bonds between N5 of FAD and C_αH of Tyr65 and O of Tyr65 and C_βH of His49.

Transferable aspherical atom model refinement of protein and DNA structures against ultrahigh-resolution X-ray data

In contrast to the independent-atom model (IAM), in which all atoms are assumed to be spherical and neutral, the transferable aspherical atom model (TAAM) takes into account the deformed valence charge density resulting from chemical bond formation and the presence of lone electron pairs. Both models can be used to refine small and large molecules, e.g. proteins and nucleic acids, against ultrahigh-resolution X-ray diffraction data. The University at Buffalo theoretical databank of aspherical pseudo-atoms has been used in the refinement of an oligopeptide, of Z-DNA hexamer and dodecamer duplexes, and of bovine trypsin. The application of the TAAM to these data improves the quality of the electron-density maps and the visibility of H atoms. It also lowers the conventional R factors and improves the atomic displacement parameters and the results of the Hirshfeld rigid-bond test. An additional advantage is that the transferred charge density allows the estimation of Coulombic interaction energy and electrostatic potential.

Charge-density analysis of an iron-sulfur protein at an ultra-high resolution of 0.48 Å

The fine structures of proteins, such as the positions of hydrogen atoms, distributions of valence electrons and orientations of bound waters, are critical factors for determining the dynamic and chemical properties of proteins. Such information cannot be obtained by conventional protein X-ray analyses at 3.0-1.5 Å resolution, in which amino acids are fitted into atomically unresolved electron-density maps and refinement calculations are performed under strong restraints. Therefore, we usually supplement the information on hydrogen atoms and valence electrons in proteins with pre-existing common knowledge obtained by chemistry in small molecules. However, even now, computational calculation of such information with quantum chemistry also tends to be difficult, especially for polynuclear metalloproteins. Here we report a charge-density analysis of the high-potential iron-sulfur protein from the thermophilic purple bacterium Thermochromatium tepidum using X-ray data at an ultra-high resolution of 0.48 Å. Residual electron densities in the conventional refinement are assigned as valence electrons in the multipolar refinement. Iron 3d and sulfur 3p electron densities of the Fe4S4 cluster are visualized around the atoms. Such information provides the most detailed view of the valence electrons of the metal complex in the protein. The asymmetry of the iron-sulfur cluster and the protein environment suggests the structural basis of charge storing on electron transfer. Our charge-density analysis reveals many fine features around the metal complex for the first time, and will enable further theoretical and experimental studies of metalloproteins.

Sub-atomic resolution X-ray crystallography and neutron crystallography: promise, challenges and potential

The International Year of Crystallography saw the number of macromolecular structures deposited in the Protein Data Bank cross the 100000 mark, with more than 90000 of these provided by X-ray crystallography. The number of X-ray structures determined to sub-atomic resolution (i.e. ≤1 Å) has passed 600 and this is likely to continue to grow rapidly with diffraction-limited synchrotron radiation sources such as MAX-IV (Sweden) and Sirius (Brazil) under construction. A dozen X-ray structures have been deposited to ultra-high resolution (i.e. ≤0.7 Å), for which precise electron density can be exploited to obtain charge density and provide information on the bonding character of catalytic or electron transfer sites. Although the development of neutron macromolecular crystallography over the years has been far less pronounced, and its application much less widespread, the availability of new and improved instrumentation, combined with dedicated deuteration facilities, are beginning to transform the field. Of the 83 macromolecular structures deposited with neutron diffraction data, more than half (49/83, 59%) were released since 2010. Sub-mm(3) crystals are now regularly being used for data collection, structures have been determined to atomic resolution for a few small proteins, and much larger unit-cell systems (cell edges >100 Å) are being successfully studied. While some details relating to H-atom positions are tractable with X-ray crystallography at sub-atomic resolution, the mobility of certain H atoms precludes them from being located. In addition, highly polarized H atoms and protons (H(+)) remain invisible with X-rays. Moreover, the majority of X-ray structures are determined from cryo-cooled crystals at 100 K, and, although radiation damage can be strongly controlled, especially since the advent of shutterless fast detectors, and by using limited doses and crystal translation at micro-focus beams, radiation damage can still take place. Neutron crystallography therefore remains the only approach where diffraction data can be collected at room temperature without radiation damage issues and the only approach to locate mobile or highly polarized H atoms and protons. Here a review of the current status of sub-atomic X-ray and neutron macromolecular crystallography is given and future prospects for combined approaches are outlined. New results from two metalloproteins, copper nitrite reductase and cytochrome c', are also included, which illustrate the type of information that can be obtained from sub-atomic-resolution (∼0.8 Å) X-ray structures, while also highlighting the need for complementary neutron studies that can provide details of H atoms not provided by X-ray crystallography.