Effect of inorganic phosphate on FMN binding and loop flexibility in Desulfovibrio desulfuricans apo-flavodoxin

From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

Structural and Functional Characterization of a Short-Chain Flavodoxin Associated with a Noncanonical 1,2-Propanediol Utilization Bacterial Microcompartment

Bacterial microcompartments (BMCs) are proteinaceous organelles that encapsulate enzymes involved in CO₂ fixation (carboxysomes) or carbon catabolism (metabolosomes). Metabolosomes share a common core of enzymes and a distinct signature enzyme for substrate degradation that defines the function of the BMC (e.g., propanediol or ethanolamine utilization BMCs, or glycyl-radical enzyme microcompartments). Loci encoding metabolosomes also typically contain genes for proteins that support organelle function, such as regulation, transport of substrate, and cofactor (e.g., vitamin B₁₂) synthesis and recycling. Flavoproteins are frequently among these ancillary gene products, suggesting that these redox active proteins play an undetermined function in many metabolosomes. Here, we report the first characterization of a BMC-associated flavodoxin (Fld1C), a small flavoprotein, derived from the noncanonical 1,2-propanediol utilization BMC locus (PDU1C) of Lactobacillus reuteri. The 2.0 Å X-ray structure of Fld1C displays the α/β flavodoxin fold, which noncovalently binds a single flavin mononucleotide molecule. Fld1C is a short-chain flavodoxin with redox potentials of -240 ± 3 mV oxidized/semiquinone and -344 ± 1 mV semiquinone/hydroquinone versus the standard hydrogen electrode at pH 7.5. It can participate in an electron transfer reaction with a photoreductant to form a stable semiquinone species. Collectively, our structural and functional results suggest that PDU1C BMCs encapsulate Fld1C to store and transfer electrons for the reactivation and/or recycling of the B₁₂ cofactor utilized by the signature enzyme.

Binding of an RNA pol II Ligand to the WW Domain of Pin1 Using Molecular Dynamics Docking Simulations

A novel docking protocol using a long, all atom molecular dynamics (MD) simulation, in an explicit solvent medium, without using any distance constraints is presented. This MD docking protocol is able to dock ligands, based on the C-terminal domain (CTD) of RNA polymerase II, into the tryptophan-tryptophan (WW) domain of Pin1. In this docking process, a significant loop-bending event occurs in order to encircle the ligand into its solvent exposed binding site, which cannot be simulated using current protocols. The simulations were validated structurally and energetically against an X-ray structure to confirm correct sampling of conformational space. Based on these simulations, and justification of the starting structure as a valid intermediate structure, a potential molecular basis for binding was predicted as well as confirming the key residues involved in the formation of the final strong and stable Pin1 WW domain-ligand complex.

Flavin-Protein Complexes: Aromatic Stacking Assisted by a Hydrogen Bond

Enzyme-catalyzed reactions often rely on a noncovalently bound cofactor whose reactivity is tuned by its immediate environment. Flavin cofactors, the most versatile catalyst encountered in biology, are often maintained within the protein throughout numbers of complex ionic and aromatic interactions. Here, we have investigated the role of π-π stacking and hydrogen bond interactions between a tyrosine and the isoalloxazine moiety of the flavin adenine dinucleotide (FAD) in an FAD-dependent RNA methyltransferase. Combining several static and time-resolved spectroscopies as well as biochemical approaches, we showed that aromatic stacking is assisted by a hydrogen bond between the phenol group and the amide of an adjacent active site loop. A mechanism of recognition and binding of the redox cofactor is proposed.

Structural and functional investigation of flavin binding center of the NqrC subunit of sodium-translocating NADH:quinone oxidoreductase from Vibrio harveyi

Na⁺-translocating NADH:quinone oxidoreductase (NQR) is a redox-driven sodium pump operating in the respiratory chain of various bacteria, including pathogenic species. The enzyme has a unique set of redox active prosthetic groups, which includes two covalently bound flavin mononucleotide (FMN) residues attached to threonine residues in subunits NqrB and NqrC. The reason of FMN covalent bonding in the subunits has not been established yet. In the current work, binding of free FMN to the apo-form of NqrC from Vibrio harveyi was studied showing very low affinity of NqrC to FMN in the absence of its covalent bonding. To study structural aspects of flavin binding in NqrC, its holo-form was crystallized and its 3D structure was solved at 1.56 Å resolution. It was found that the isoalloxazine moiety of the FMN residue is buried in a hydrophobic cavity and that its pyrimidine ring is squeezed between hydrophobic amino acid residues while its benzene ring is extended from the protein surroundings. This structure of the flavin-binding pocket appears to provide flexibility of the benzene ring, which can help the FMN residue to take the bended conformation and thus to stabilize the one-electron reduced form of the prosthetic group. These properties may also lead to relatively weak noncovalent binding of the flavin. This fact along with periplasmic location of the FMN-binding domains in the vast majority of NqrC-like proteins may explain the necessity of the covalent bonding of this prosthetic group to prevent its loss to the external medium.

Mechanism of FMN Binding to the Apoflavodoxin from Helicobacter pylori

Flavodoxins are bacterial electron transport proteins whose redox competence is due to the presence of a tightly but noncovalently bound FMN molecule. While the thermodynamics of the complex are understood, the mechanism of association between the apoflavodoxin and the redox cofactor is not so clear. We investigate here the mechanism of FMN binding to the apoflavodoxin from Helicobacter pylori, an essential protein that is being used as a target to develop antimicrobials. This flavodoxin is structurally peculiar as it lacks the typical bulky residue interacting with the FMN re face but bears instead a small alanine. FMN binding is biphasic, regardless of the presence of phosphate molecules in solution, while riboflavin binding takes place in a single step, the rate constant of which coincides with the fast phase of FMN binding. A mutational study at the isoalloxazine and phosphate subsites for FMN binding clearly indicates that FMN association is always limited by interaction with the isoalloxazine subsite because mutating residues that interact with the phosphate moiety of FMN in the native complex hardly changes the observed rate constants and amplitudes. In contrast, replacing tyr92, which interacts with the isoalloxazine, greatly lowers the rate constants. Our analysis indicates that the two FMN binding phases observed are related neither with alternative or sequential interaction with the two binding subsites nor with the presence of bound phosphate. It is possible that they reflect the intrinsic conformational heterogeneity of the apoflavodoxin ensemble.

Pseudosymmetry, high copy number and twinning complicate the structure determination of Desulfovibrio desulfuricans (ATCC 29577) flavodoxin

The crystal structure of oxidized flavodoxin from Desulfovibrio desulfuricans (ATCC 29577) was determined by molecular replacement in two crystal forms, P3(1)21 and P4(3), at 2.5 and 2.0 A resolution, respectively. Structure determination in space group P3(1)21 was challenging owing to the presence of pseudo-translational symmetry and a high copy number in the asymmetric unit (8). Initial phasing attempts in space group P3(1)21 by molecular replacement using a poor search model (46% identity) and multi-wavelength anomalous dispersion were unsuccessful. It was necessary to solve the structure in a second crystal form, space group P4(3), which was characterized by almost perfect twinning, in order to obtain a suitable search model for molecular replacement. This search model with complementary approaches to molecular replacement utilizing the pseudo-translational symmetry operators determined by analysis of the native Patterson map facilitated the selection and manual placement of molecules to generate an initial solution in the P3(1)21 crystal form. During the early stages of refinement, application of the appropriate twin law, (-h, -k, l), was required to converge to reasonable R-factor values despite the fact that in the final analysis the data were untwinned and the twin law could subsequently be removed. The approaches used in structure determination and refinement may be applicable to other crystal structures characterized by these complicating factors. The refined model shows flexibility of the flavin mononucleotide coordinating loops indicated by the isolation of two loop conformations and provides a starting point for the elucidation of the mechanism used for protein-partner recognition.

Folding, stability and shape of proteins in crowded environments: experimental and computational approaches

How the crowded environment inside cells affects folding, stability and structures of proteins is a vital question, since most proteins are made and function inside cells. Here we describe how crowded conditions can be created in vitro and in silico and how we have used this to probe effects on protein properties. We have found that folded forms of proteins become more compact in the presence of macromolecular crowding agents; if the protein is aspherical, the shape also changes (extent dictated by native-state stability and chemical conditions). It was also discovered that the shape of the macromolecular crowding agent modulates the folding mechanism of a protein; in addition, the extent of asphericity of the protein itself is an important factor in defining its folding speed.

Allophycocyanin trimer stability and functionality are primarily due to polar enhanced hydrophobicity of the phycocyanobilin binding pocket

Allophycocyanin (APC) is the primary pigment-protein component of the cores of the phycobilisome antenna complex. In addition to an extremely high degree of amino acid sequence conservation, the overall structures of APC from both mesophilic and thermophilic species are almost identical at all levels of assembly, yet APC from thermophilic organisms should have structural attributes that prevent thermally induced denaturation. We determined the structure of APC from the thermophilic cyanobacterium Thermosynechococcus vulcanus to 2.9 A, reaffirming the conservation of structural similarity with APC from mesophiles. We provide spectroscopic evidence that T. vulcanus APC is indeed more stable at elevated temperatures in vitro, when compared with the APC from mesophilic species. APC thermal and chemical stability levels are further enhanced when monitored in the presence of high concentrations of buffered phosphate, which increases the strength of hydrophobic interactions, and may mimic the effect of cytosolic crowding. Absorption spectroscopy, size-exclusion HPLC, and native gel electrophoresis also show that the thermally or chemically induced changes in the APC absorption spectra that result in the loss of the prominent 652-nm band in trimeric APC are not a result of physical monomerization. We propose that the bathochromic shift that occurs in APC upon trimerization is due to the coupling of the hydrophobicity of the alpha84 phycocyanobilin cofactor environment created by a deep cleft formed by the beta subunit with highly charged flanking regions. This arrangement also provides the additional stability required by thermophiles at elevated temperatures. The chemical environment that induces the bathochromic shift in APC trimers is different from the source of shifts in the absorption of monomers of the terminal energy acceptors APC(B) and L(CM), as visualized by the building of molecular models.

Redox potential difference between Desulfovibrio vulgaris and Clostridium beijerinckii flavodoxins

The redox potential of the flavin mononucleotide (FMN) hydroquinones for one-electron reduction in the Desulfovibrio vulgaris ( D. vulgaris) flavodoxin ( E sq/hq for FMNH (*)/FMNH (-)) was calculated using the crystal structure of the relevant hydroquinone form and compared to the results of the Clostridium beijerinckii ( C. beijerinckii) flavodoxin. In D. vulgaris and C. beijerinckii flavodoxins, the protein side chain causes significant downshifts of 170 and 240 mV in E sq/hq, respectively. In the C. beijerinckii flavodoxin, the E sq/hq downshift because of the protein side chain is essentially compensated by the counter influence of the protein backbone ( E sq/hq upshift of 260 mV). However, in the D. vulgaris flavodoxin, the corresponding protein backbone influence on E sq/hq is significantly small, i.e., less than half of that in the C. beijerinckii flavodoxin. In particular, there is a significant difference in the influence of the protein backbone of the so-called 60s loop region between the two flavodoxins. The E sq/hq difference can be best explained by the lower compensation of the side chain influence by the backbone influence in the D. vulgaris flavodoxin than in the C. beijerinckii flavodoxin.

Common conformational changes in flavodoxins induced by FMN and anion binding: the structure of Helicobacter pylori apoflavodoxin

Flavodoxins, noncovalent complexes between apoflavodoxins and flavin mononucleotide (FMN), are useful models to investigate the mechanism of protein/flavin recognition. In this respect, the only available crystal structure of an apoflavodoxin (that from Anabaena) showed a closed isoalloxazine pocket and the presence of a bound phosphate ion, which posed many questions on the recognition mechanism and on the potential physiological role exerted by phosphate ions. To address these issues we report here the X-ray structure of the apoflavodoxin from the pathogen Helicobacter pylori. The protein naturally lacks one of the conserved aromatic residues that close the isoalloxazine pocket in Anabaena, and the structure has been determined in a medium lacking phosphate. In spite of these significant differences, the isoallozaxine pocket in H. pylori apoflavodoxin appears also closed and a chloride ion is bound at a native-like FMN phosphate site. It seems thus that it is a general characteristic of apoflavodoxins to display closed, non-native, isoalloxazine binding sites together with native-like, rather promiscuous, phosphate binding sites that can bear other available small anions present in solution. In this respect, both binding energy hot spots of the apoflavodoxin/FMN complex are initially unavailable to FMN binding and the specific spot for FMN recognition may depend on the dynamics of the two candidate regions. Molecular dynamics simulations show that the isoalloxazine binding loops are intrinsically flexible at physiological temperatures, thus facilitating the intercalation of the cofactor, and that their mobility is modulated by the anion bound at the phosphate site.

Molecular crowding enhances native structure and stability of alpha/beta protein flavodoxin

To investigate the consequences of macromolecular crowding on the behavior of a globular protein, we performed a combined experimental and computational study on the 148-residue single-domain alpha/beta protein, Desulfovibrio desulfuricans apoflavodoxin. In vitro thermal unfolding experiments, as well as assessment of native and denatured structures, were probed by using far-UV CD in the presence of various amounts of Ficoll 70, an inert spherical crowding agent. Ficoll 70 has a concentration-dependent effect on the thermal stability of apoflavodoxin (DeltaT(m) of 20 degrees C at 400 mg/ml; pH 7). As judged by CD, addition of Ficoll 70 causes an increase in the amount of secondary structure in the native-state ensemble (pH 7, 20 degrees C) but only minor effects on the denatured state. Theoretical calculations, based on an off-lattice model and hard-sphere particles, are in good agreement with the in vitro data. The simulations demonstrate that, in the presence of 25% volume occupancy of spheres, native flavodoxin is thermally stabilized, and the free energy landscape shifts to favor more compact structures in both native and denatured states. The difference contact map reveals that the native-state compaction originates in stronger interactions between the helices and the central beta-sheet, as well as by less fraying in the terminal helices. This study demonstrates that macromolecular crowding has structural effects on the folded ensemble of polypeptides.

Thermodynamics of ligand binding to P450 2B4 and P450eryF studied by isothermal titration calorimetry

Structural plasticity and cooperativity in ligand recognition are two key aspects of the catalytic diversity of cytochrome P450 enzymes. As more mammalian P450 crystal structures have emerged, computational modeling has become a major tool to predict drug metabolism and interactions. There is a need for real solution thermodynamic data to support modeling and crystallographic observations. Using isothermal titration calorimetry (ITC) we successfully evaluated the conformational flexibility of P450 2B4 in binding imidazole inhibitors of different size and chemistry and dissected the stoichiometry and energetics of ligand binding allostery in P450eryF. Thermodynamic signatures obtained by ITC nicely correlated with structural and modeling results. Thus, ITC is a powerful tool to study structure-function relationships in P450s.

Macromolecular crowding increases structural content of folded proteins

Here we show that increased amount of secondary structure is acquired in the folded states of two structurally-different proteins (alpha-helical VlsE and alpha/beta flavodoxin) in the presence of macromolecular crowding agents. The structural content of flavodoxin and VlsE is enhanced by 33% and 70%, respectively, in 400 mg/ml Ficoll 70 (pH 7, 20 degrees C) and correlates with higher protein-thermal stability. In the same Ficoll range, there are only small effects on the unfolded-state structures of the proteins. This is the first in vitro assessment of crowding effects on the native-state structures at physiological conditions. Our findings imply that for proteins with low intrinsic stability, the functional structures in vivo may differ from those observed in dilute buffers.

Dissecting the Thermodynamics and Cooperativity of Ligand Binding in Cytochrome P450eryF

Conformational flexibility and cooperativity in ligand recognition are two key aspects of the catalytic diversity of cytochrome P450 enzymes. In this study, we dissect the ligand binding stoichiometry and energetics of the soluble bacterial P450eryF by isothermal titration calorimetry (ITC) using three allosteric and two non-allosteric ligands of diverse chemistry. Complementary spectral binding studies and sequential, two-ligand docking simulations were performed to help assign the binding sites. Binding of 4-phenylpyridine (4-PP) or 4-(4-chlorophenyl)imidazole (4-CPI) showed 1:1 stoichiometry in ITC, consistent with the lack of cooperativity observed in spectral binding studies. The larger ligands 9-aminophenanthrene (9-AP), 1-pyrenebutanol (1-PB), and alpha-naphthoflavone (ANF) show cooperative spectral binding and yielded 2:1 stoichiometry. The associated thermodynamic parameters for the sites were calculated using a sequential binding mechanism. The binding constant (KD) for the first site was almost two times lower than that of the second site for all three compounds. Ligand binding at site 1 was entropically favored, whereas binding at site 2 was weak and entropically unfavorable. Simulations showed that two molecules of 9-AP, ANF or 1-PB can be adequately docked to two individual sub-sites within a large binding pocket. The absence of hydrophobic tethering and ligand stacking are consistent with the single low affinity binding site observed for 4-CPI and 4-PP. Competitive binding studies with P450eryF preloaded with either 1-PB or ANF showed a decrease in the affinities for 9-AP at both the sites, with large entropy-enthalpy compensation, indicating the ability of the binding pocket to accommodate two ligands of diverse chemistry and enable cooperativity.

Folding of Desulfovibrio desulfuricans flavodoxin is accelerated by cofactor fly-casting

Folding of cofactor-binding proteins involves ligand binding in addition to polypeptide folding. We here assess the kinetic folding/binding landscape for Desulfovibrio desulfuricans flavodoxin that coordinates an FMN cofactor. The apo-form folds in a two-step process involving a burst-phase intermediate. Studies on Tyr98Ala and Trp60Ala variants reveal that these aromatics-that stack with the FMN in the holo-form-are not participating in the apo-protein folding pathway. However, these residues are essential for FMN interactions with the unfolded protein during refolding of holo-flavodoxin. Unfolding of wild-type holo-flavodoxin is coupled to FMN dissociation whereas for Tyr98Ala and Trp60Ala holo-variants, FMN dissociates before polypeptide unfolding. Both variants refold as apo-proteins before FMN rebinds. In sharp contrast, refolding of unfolded wild-type holo-flavodoxin is over an order of magnitude faster than that of the apo-form, the pathway does not include a burst-phase intermediate, and the speed is independent of FMN excess ratio. These observations demonstrate that FMN binds rapidly to the unfolded polypeptide and guides folding straight to the native state. As this path to functional D. desulfuricans holo-flavodoxin is faster than if the cofactor binds to pre-folded apo-protein, this is one of few examples where molecular recognition via a "fly-casting" mechanism is kinetically favored.

Conformational flexibility of mammalian cytochrome P450 2B4 in binding imidazole inhibitors with different ring chemistry and side chains. Solution thermodynamics and molecular modeling

Recent x-ray structures of cytochrome P450 2B4 (CYP2B4) reveal an open form that undergoes a large-scale structural transition to a closed form upon binding to 4-(4-chlorophenyl)imidazole (4-CPI). Here, we report for the first time a complete solution thermodynamic study using isothermal titration calorimetry supported by spectroscopic studies to elucidate the conformational flexibility of CYP2B4 in binding imidazole inhibitors with different ring chemistry and side chains: 4-CPI, 1-benzylimidazole (1-BI), 1-CPI, 4-phenylimidazole (4-PI), 1-(2-(benzyloxy)ethyl)imidazole (BEI), and 1-PI. Each of the inhibitors induced type II spectral changes, and IC50 values for enzyme inhibition ranged from 0.1 to 2.4 microM, following the order 1-BI < 4-CPI < 1-CPI < 4-PI < BEI < 1-PI. Calorimetric titrations using monomeric enzyme yielded a 1:1 binding stoichiometry, with the associated KD values ranging from 0.3 to 4.8 microM and following the same rank order as the IC50 values. Changes in enthalpy at 25 degrees C ranged from -6.5 to -8.8 kcal mol(-1). The largest difference in binding entropy (+5.9 versus -4.1 cal mol(-1) K(-1)) was observed between 4-CPI and BEI, respectively, with a 2-fold difference in heat capacity changes (-604 versus -331 cal mol(-1) K(-1)), which is inferred to result from the reduction of apolar surface area of the enzyme ensuing from a conformational change upon 4-CPI binding. Accessibility to acrylamide of the only tryptophan (Trp121), which is located in helix C, was greatly decreased only in protein bound to 4-CPI. Steric restrictions hindered the perfect docking of only BEI to the closed conformation of the enzyme. The thermodynamic signature obtained for structurally similar inhibitors suggests remarkable plasticity of CYP2B4.