Improvement of an antibody neutralizing the anthrax toxin by simultaneous mutagenesis of its six hypervariable loops

The enhancement of antibody affinity by mutagenesis targeting only complementarity determining regions has the advantage of respecting the framework regions, which are important for tolerance if clinical use is envisaged. Here, starting from a Fab (antigen-binding fragment; 35PA₈₃) capable of neutralizing the lethal toxin of anthrax and having an affinity of 3.4 nM for its antigen, a phage-displayed library of variants where all six complementarity determining regions (73 positions) were targeted for mutagenesis was built. This library contained 5 × 10⁸ variants, and each variant carried four mutations on average. The library was first panned with adsorbed antigen and washes of increasing stringency. It was then screened in parallel with either small concentrations of soluble biotinylated antigen or adsorbed antigen and long elution times in the presence of soluble antigen. The stringencies of both selections were pushed as far as possible. Compared with 35PA₈₃, the best selected clone had an affinity enhanced 19-fold, to 180 pM, and its 50% inhibitory concentration was decreased by 40%. The results of the two selection methods were compared, and the generality of these methods was considered.     

Probing the mechanism of GSH activation in Schistosoma haematobium glutathione-S-transferase by site-directed mutagenesis and X-ray crystallography

During turnover, the catalytic tyrosine residue (Tyr10) of the sigma class Schistosoma haematobium wild-type glutathione-S-transferase is expected to switch alternately in and out of the reduced glutathione-binding site (G-site). The Tyrout10 conformer forms a pi-cation interaction with the guanidinium group of Arg21. As in other similar glutathione-S-transferases, the catalytic Tyr has a low pKa of 7.2. In order to investigate the catalytic role of Tyr10, and the structural and functional roles of Arg21, we carried out structural studies on two Arg21 mutants (R21L and R21Q) and a Tyr10 mutant, Y10F. Our crystallographic data for the two Arg21 mutants indicate that only the Tyrout10 conformation is populated, thereby excluding a role of Arg21 in the stabilisation of the out conformation. However, Arg21 was confirmed to be catalytically important and essential for the low pKa of Tyr10. Upon comparison with structural data generated for reduced glutathione-bound and inhibitor-bound wild-type enzymes, it was observed that the orientations of Tyr10 and Arg35 are concerted and that, upon ligand binding, minor rearrangements occur within conserved residues in the active site loop. These rearrangements are coupled to quaternary rigid-body movements at the dimer interface and alterations in the localisation and structural order of the C-terminal domain.     

Probing hydrogen peroxide oxidation kinetics of wild-type Synechocystis catalase-peroxidase (KatG) and selected variants

Catalase-peroxidases (KatGs) are unique bifunctional heme peroxidases that exhibit peroxidase and substantial catalase activities. Nevertheless, the reaction pathway of hydrogen peroxide dismutation, including the electronic structure of the redox intermediate that actually oxidizes H₂O₂, is not clearly defined. Several mutant proteins with diminished overall catalase but wild-type-like peroxidase activity have been described in the last years. However, understanding of decrease in overall catalatic activity needs discrimination between reduction and oxidation reactions of hydrogen peroxide. Here, by using sequential-mixing stopped-flow spectroscopy, we have investigated the kinetics of the transition of KatG compound I (produced by peroxoacetic acid) to its ferric state by trapping the latter as cyanide complex. Apparent bimolecular rate constants (pH 6.5, 20 °C) for wild-type KatG and the variants Trp122Phe (lacks KatG-typical distal adduct), Asp152Ser (controls substrate access to the heme cavity) and Glu253Gln (channel entrance) are reported to be 1.2 × 10⁴ M^− 1 s^− 1, 30 M^− 1 s^− 1, 3.4 × 10³ M^− 1 s^− 1, and 8.6 × 10³ M^− 1 s^− 1, respectively. These findings are discussed with respect to steady-state kinetic data and proposed reaction mechanism(s) for KatG. Assets and drawbacks of the presented method are discussed.     

X-ray structure and mutagenesis of the scorpion depressant toxin LqhIT2 reveals key determinants crucial for activity and anti-insect selectivity

Scorpion depressant beta-toxins show high preference for insect voltage-gated sodium channels (Na(v)s) and modulate their activation. Although their pharmacological and physiological effects were described, their three-dimensional structure and bioactive surface have never been determined. We utilized an efficient system for expression of the depressant toxin LqhIT2 (from Leiurus quinquestriatushebraeus), mutagenized its entire exterior, and determined its X-ray structure at 1.2 A resolution. The toxin molecule is composed of a conserved cysteine-stabilized alpha/beta-core (core-globule), and perpendicular to it an entity constituted from the N and C-terminal regions (NC-globule). The surface topology and overall hydrophobicity of the groove between the core and NC-globules (N-groove) is important for toxin activity and plays a role in selectivity to insect Na(v)s. The N-groove is flanked by Glu24 and Tyr28, which belong to the "pharmacophore" of scorpion beta-toxins, and by the side-chains of Trp53 and Asn58 that are important for receptor site recognition. Substitution of Ala13 by Trp in the N-groove uncoupled activity from binding, suggesting that this region of the molecule is also involved in "voltage-sensor trapping", the mode of action that typifies scorpion beta-toxins. The involvement of the N-groove in recognition of the receptor site, which seems to require a defined topology, as well as in sensor trapping, which involves interaction with a moving channel region, is puzzling. On the basis of the mutagenesis studies we hypothesize that following binding to the receptor site, the toxin undergoes a conformational change at the N-groove region that facilitates the trapping of the voltage-sensor in its activated position.     

Redox control of fast ligand dissociation from Escherichia coli cytochrome bd

Bacterial bd-type quinol oxidases, such as cytochrome bd from Escherichia coli, contain three hemes, but no copper. In contrast to heme-copper oxidases and similarly to globins, single electron-reduced cytochrome bd forms stable complexes with O(2), NO and CO at ferrous heme d. Kinetics of ligand dissociation from heme d(2+) in the single electron- and fully-reduced cytochrome bd from E. coli has been investigated by rapid mixing spectrophotometry at 20 degrees C. Data show that (i) O(2) dissociates at 78 s(-1), (ii) NO and CO dissociation is fast as compared to heme-copper oxidases and (iii) dissociation in the single electron-reduced state is hindered as compared to the fully-reduced enzyme. Presumably, rapid ligand dissociation requires reduced heme b(595). As NO, an inhibitor of respiratory oxidases, is involved in the immune response against microbial infection, the rapid dissociation of NO from cytochrome bd may have important bearings on the patho-physiology of enterobacteria.     

Substrate kringle-mediated catalysis by the streptokinase-plasmin activator complex: Critical contribution of kringle-4 revealed by the mutagenesis approaches

Streptokinase (SK) is a protein co-factor with a potent capability for human plasminogen (HPG) activation. Our previous studies [1] have indicated a major role of long-range protein-protein contacts between the three domains (alpha, beta, and gamma) of SK and the multi-domain HPG substrate (K1-K5CD). To further explore this phenomenon, we prepared truncated derivatives of HPG with progressive removal of kringle domains, like K5CD, K4K5CD, K3-K5CD (K3K4K5CD), K2-K5CD (K2K3K4K5CD) and K1-K5CD (K1K2K3K4K5CD). While urokinase (uPA) cleaved the scissile peptide in the isolated catalytic domain (μPG) with nearly the same rate as with full-length HPG, SK-plasmin showed only 1-2% activity, revealing mutually distinct mechanisms of HPG catalysis between the eukaryotic and prokaryotic activators. Remarkably, with SK.HPN (plasmin), the ‘addition’ of both kringles 4 and 5 onto the catalytic domain showed catalytic rates comparable to full length HPG, thus identifying the dependency of the “long-range” enzyme-substrate interactions onto these two CD-proximal domains. Further, chimeric variants of K5CD were generated by swapping the kringle domains of HPG with those of uPA and TPA (tissue plasminogen activator), separately. Surprisingly, although native-like catalytic turnover rates were retained when either K1, K2 or K4 of HPG was substituted at the K5 position in K5CD, these were invariably lost once substituted with the evolutionarily more distant TPA- and uPA-derived kringles. The present results unveil a novel mechanism of SK.HPN action in which augmented catalysis occurs through enzyme-substrate interactions centered on regions in substrate HPG (kringles 4 and 5) that are spatially distant from the scissile peptide bond.     

Absolute correlation between lag time and growth rate in the spontaneous formation of several amyloid-like aggregates and fibrils

The formation of polypeptide aggregates, including amyloid fibrils and prions, is a biochemical process of considerable interest in the context of its association with ageing and neurodegeneration. Aggregation occurs typically with a lag phase and a growth phase that reflect an underlying nucleation-polymerisation mechanism. While the propensity of nucleation can be estimated from the lag time t(l), the efficiency of growth is represented by the growth rate k(g). Here, I have analysed the absolute k(g) and t(l) values from a total of 298 samples prepared from insulin, glucagon and different sequence variants of the Alzheimer's Abeta(1-40) peptide. Although these samples differ in the conditions of aggregation, systematic comparison reveals an overall similarity in the plot of k(g)versus t(l). The plot fits readily with the simple equation k(g)=alpha/t(l) and by using a proportionality factor alpha of 4.5. In contrast to the individual values of k(g) and t(l) that depend substantially on sequential and environmental parameters, alpha seems much less affected by such factors. These data suggest mechanistic similarities in the nucleation behaviour of different amyloid-like fibrils and aggregates.     

Uncoupling the MgATP-induced inhibition and aggregation of Escherichia coli phosphofructokinase-2 by C-terminal mutations

Binding of MgATP to an allosteric site of Escherichia coli phosphofructokinase-2 (Pfk-2) provoked inhibition and a dimer-tetramer (D-T) conversion of the enzyme. Successive deletions of up to 10 residues and point mutations at the C-terminal end led to mutants with elevated K(Mapp) values for MgATP which failed to show the D-T conversion, but were still inhibited by the nucleotide. Y306 was required for the quaternary packing involved in the D-T conversion and the next residue, L307, was crucial for the ternary packing necessary for the catalytic MgATP-binding site. These results show that the D-T conversion could be uncoupled from the conformational changes that lead to the MgATP-induced allosteric inhibition.     

Enzymatic activity assay of D-hydantoinase by isothermal titration calorimetry. Determination of the thermodynamic activation parameters for the hydrolysis of several substrates

Isothermal titration calorimetry (ITC) has been applied to the determination of the activity of D-hydantoinase (EC 3.5.2.2) with several substrates by monitoring the heat released during the reaction. The method is based on the proportionality between the reaction rate and the thermal power (heat/time) generated. Microcalorimetric assays carried out at different temperatures provided the dependence of the catalytic rate constant on temperature. We show that ITC assay is a nondestructive method that allows the determination of the catalytic rate constant (kcat), Michaelis constant (KM), activation energy and activation Gibbs energy, enthalpy and entropy of this reaction.     

Dual role of FMN in flavodoxin function: Electron transfer cofactor and modulation of the protein-protein interaction surface

Flavodoxin (Fld) replaces Ferredoxin (Fd) as electron carrier from Photosystem I (PSI) to Ferredoxin-NADP⁺ reductase (FNR). A number of Anabaena Fld (AnFld) variants with replacements at the interaction surface with FNR and PSI indicated that neither polar nor hydrophobic residues resulted critical for the interactions, particularly with FNR. This suggests that the solvent exposed benzenoid surface of the Fld FMN cofactor might contribute to it. FMN has been replaced with analogues in which its 7- and/or 8-methyl groups have been replaced by chlorine and/or hydrogen. The oxidised Fld variants accept electrons from reduced FNR more efficiently than Fld, as expected from their less negative midpoint potential. However, processes with PSI (including reduction of Fld semiquinone by PSI, described here for the first time) are impeded at the steps that involve complex re-arrangement and electron transfer (ET). The groups introduced, particularly chlorine, have an electron withdrawal effect on the pyrazine and pyrimidine rings of FMN. These changes are reflected in the magnitude and orientation of the molecular dipole moment of the variants, both factors appearing critical for the re-arrangement of the finely tuned PSI:Fld complex. Processes with FNR are also slightly modulated. Despite the displacements observed, the negative end of the dipole moment points towards the surface that contains the FMN, still allowing formation of complexes competent for efficient ET. This agrees with several alternative binding modes in the FNR:Fld interaction. In conclusion, the FMN in Fld not only contributes to the redox process, but also to attain the competent interaction of Fld with FNR and PSI.     

Effect of different salt ions on the propensity of aggregation and on the structure of Alzheimer's abeta(1-40) amyloid fibrils

The formation of amyloid fibrils and other polypeptide aggregates depends strongly on the physico-chemical environment. One such factor affecting aggregation is the presence and concentration of salt ions. We have examined the effects of salt ions on the aggregation propensity of Alzheimer's Abeta(1-40) peptide and on the structure of the dissolved and of the fibrillar peptide. All salts examined promote aggregation strongly. The most pronounced effect is seen within the cationic series, i.e. for MgCl2. Evaluation of different possible explanations suggests that Abeta(1-40) aggregation depends on direct interaction between ions and Abeta(1-40) peptide, and correlates with ion-induced changes of the surface tension. Salts have profound effects on the fibril structure. In the presence of salts, fibrils are associated with smaller diameters, narrower crossover distances and lower amide I maxima. Since Abeta(1-40) aggregation responds to salts in a manner unlike that for other polypeptides, such as glucagon, beta2-microglobulin or alpha-synuclein; these data argue that there is no fully uniform way in which salts affect aggregation of different polypeptide chains. These observations are important for understanding and predicting aggregation on the basis of simple physico-chemical properties.     

In vitro metabolism of zolmitriptan in rat cytochromes induced with beta-naphthoflavone and the interaction between six drugs and zolmitriptan

Zolmitriptan is a novel and highly selective 5-HT(1B/1D) receptor agonist used as an acute oral treatment for migraine. There are few reports regarding the in vitro metabolism of zolmitriptan. Previous studies indicated zolmitriptan was metabolized via CYP1A2 in human hepatic microsomes. In order to study the enzyme kinetics and drug interaction, the metabolism of zolmitriptan and possible drug-drug interactions were investigated in rat hepatic microsomes induced with different inducers. An active metabolite, N-demethylzolmitriptan, was detected and another minor, inactive metabolite that was reported in human hepatic microsomes was not detected in this study. The enzyme kinetics for the formation of N-demethylzolmitriptan from zolmitriptan in rat liver microsomes pretreated with BNF were 96+/-22 microM (K(m)), 11+/-3 pmol min(-1)mg protein(-1) (V(max)), and 0.12+/-0.02 microl min(-1)mg protein(-1) (CL(int)). Fluvoxamine and diphenytriazol inhibited zolmitriptan N-demethylase activity catalyzed by CYP1A2 (K(i)=3.8+/-0.3 and 3.2+/-0.1 microM, respectively). Diazepam and propranolol elicited a slight inhibitory effect on the metabolism of zolmitriptan (K(i)=70+/-11 and 90+/-18 microM, respectively). Cimetidine and moclobemide produced no significant effect on the metabolism of zolmitriptan. Fluvoxamine yielded a k(inactivation) value of 0.16 min(-1), and K(i) of 57 microM. The results suggest that rat hepatic microsomes are a reasonable model to study the metabolism of zolmitriptan, although there is a difference in the amount of minor, inactive metabolites between human hepatic microsomes and rat liver microsomes. The results of the inhibition experiments provided information for the interactions between zolmitriptan and drugs co-administrated in clinic, and it is helpful to explain the drug-drug interactions of clinical relevance on enzyme level. This study aso demonstrated that fluvoxamine may be a mechanism-based inactivator of CYP1A2.     

The cold denatured state is compact but expands at low temperatures: hydrodynamic properties of the cold denatured state of the C-terminal domain of L9

A point mutation of a small globular protein, the C-terminal domain of L9 destabilizes the protein and leads to observable cold-denaturation at temperatures above zero. The cold denatured state is in slow exchange with the native state on the NMR time scale, and this allows the hydrodynamic properties of the cold unfolded state and the native state to be measured under identical conditions using pulsed-field gradient NMR diffusion measurements. This provides the first experimental measurement of the hydrodynamic properties of a cold unfolded protein and its folded form under identical conditions. Hydrodynamic radii of the cold-induced unfolded states were measured for a set of temperatures ranging from 2 °C to 25 °C at pD 6.6 in the absence of denaturant. The cold unfolded state is compact compared to the urea or acid unfolded state and a trend of increasing radii of hydration is observed as the temperature is lowered. These observations are confirmed by experiments on the same protein at pD 8.0, where it is more stable, in the presence of a modest concentration of urea. The expansion of the cold-denatured state at lower temperatures is consistent with the temperature dependence of hydrophobic interactions.     

Photo-oxidation of P740, the primary electron donor in photosystem I from Acaryochloris marina

Fourier transform infrared spectroscopy (FTIR) difference spectroscopy in combination with deuterium exchange experiments has been used to study the photo-oxidation of P740, the primary electron donor in photosystem I from Acaryochloris marina. Comparison of (P740(+)-P740) and (P700(+)-P700) FTIR difference spectra show that P700 and P740 share many structural similarities. However, there are several distinct differences also: 1), The (P740(+)-P740) FTIR difference spectrum is significantly altered upon proton exchange, considerably more so than the (P700(+)-P700) FTIR difference spectrum. The P740 binding pocket is therefore more accessible than the P700 binding pocket. 2), Broad, "dimer" absorption bands are observed for both P700(+) and P740(+). These bands differ significantly in substructure, however, suggesting differences in the electronic organization of P700(+) and P740(+). 3), Bands are observed at 2727(-) and 2715(-) cm(-1) in the (P740(+)-P740) FTIR difference spectrum, but are absent in the (P700(+)-P700) FTIR difference spectrum. These bands are due to formyl CH modes of chlorophyll d. Therefore, P740 consists of two chlorophyll d molecules. Deuterium-induced modification of the (P740(+)-P740) FTIR difference spectrum indicates that only the highest frequency 13(3) ester carbonyl mode of P740 downshifts, indicating that this ester mode is weakly H-bonded. In contrast, the highest frequency ester carbonyl mode of P700 is free from H-bonding. Deuterium-induced changes in (P740(+)-P740) FTIR difference spectrum could also indicate that one of the chlorophyll d 3(1) carbonyls of P740 is hydrogen bonded.     

Structural consideration of the formation of the activation complex between the staphylokinase-like streptococcal plasminogen activator PadA and bovine plasminogen

The characteristics of a streptococcal plasminogen activator (PA) displaying specificity for ruminant plasminogen (Plg) were defined using molecular approaches. The 16-kDa secreted protein PadA was found to be prevalent in Streptococcus dysgalactiae subspecies dysgalactiae isolated from cases of bovine mastitis and septic arthritis in lambs. PadA was able to activate bovine, ovine and caprine Plg, but not human Plg. Amino acid sequence analysis identified a limited level of homology to other streptococcal PAs, including streptokinase; however, PadA was found to align well with and match in size the staphylococcal PA, staphylokinase. Recombinant PadA was used to investigate interaction with bovine Plg, leading to formation of an activator complex that was capable of recruiting and converting further substrate Plg into plasmin. Individual non-overlapping peptides of PadA or bovine microplasminogen were found to block the interaction between PadA and bovine Plg, preventing the formation of the activation complex. Homology modelling based upon structures of staphylokinase complexed with human microplasminogen supported these findings by placing critical residues in close proximity to the plasmin component of the activation complex.     

Experimental and computational analyses of the energetic basis for dual recognition of immunity proteins by colicin endonucleases

Colicin endonucleases (DNases) are bound and inactivated by immunity (Im) proteins. Im proteins are broadly cross-reactive yet specific inhibitors binding cognate and non-cognate DNases with K_d values that vary between 10^− 4 and 10^− 14 M, characteristics that are explained by a ‘dual-recognition’ mechanism. In this work, we addressed for the first time the energetics of Im protein recognition by colicin DNases through a combination of E9 DNase alanine scanning and double-mutant cycles (DMCs) coupled with kinetic and calorimetric analyses of cognate Im9 and non-cognate Im2 binding, as well as computational analysis of alanine scanning and DMC data. We show that differential ΔΔGs observed for four E9 DNase residues cumulatively distinguish cognate Im9 association from non-cognate Im2 association. E9 DNase Phe86 is the primary specificity hotspot residue in the centre of the interface, which is coordinated by conserved and variable hotspot residues of the cognate Im protein. Experimental DMC analysis reveals that only modest coupling energies to Im9 residues are observed, in agreement with calculated DMCs using the program ROSETTA and consistent with the largely hydrophobic nature of E9 DNase-Im9 specificity contacts. Computed values for the 12 E9 DNase alanine mutants showed reasonable agreement with experimental ΔΔG data, particularly for interactions not mediated by interfacial water molecules. ΔΔG predictions for residues that contact buried water molecules calculated using solvated rotamer models met with mixed success; however, we were able to predict with a high degree of accuracy the location and energetic contribution of one such contact. Our study highlights how colicin DNases are able to utilise both conserved and variable amino acids to distinguish cognate from non-cognate Im proteins, with the energetic contributions of the conserved residues modulated by neighbouring specificity sites.     

Water deficit in field-grown Gossypium hirsutum primarily limits net photosynthesis by decreasing stomatal conductance,increasing photorespiration,and increasing the ratio of dark respiration to gross photosynthesis

Much effort has been expended to improve irrigation efficiency and drought tolerance of agronomic crops; however, a clear understanding of the physiological mechanisms that interact to decrease source strength and drive yield loss has not been attained. To elucidate the underlying mechanisms contributing to inhibition of net carbon assimilation under drought stress, three cultivars of Gossypium hirsutum were grown in the field under contrasting irrigation regimes during the 2012 and 2013 growing season near Camilla, Georgia, USA. Physiological measurements were conducted on three sample dates during each growing season (providing a broad range of plant water status) and included, predawn and midday leaf water potential (Ψ_PD and Ψ_MD), gross and net photosynthesis, dark respiration, photorespiration, and chlorophyll a fluorescence. End-of-season lint yield was also determined. Ψ_PD ranged from −0.31 to −0.95 MPa, and Ψ_MD ranged from −1.02 to −2.67 MPa, depending upon irrigation regime and sample date. G. hirsutum responded to water deficit by decreasing stomatal conductance, increasing photorespiration, and increasing the ratio of dark respiration to gross photosynthesis, thereby limiting P_N and decreasing lint yield (lint yield declines observed during the 2012 growing season only). Conversely, even extreme water deficit, causing a 54% decline in P_N, did not negatively affect actual quantum yield, maximum quantum yield, or photosynthetic electron transport. It is concluded that P_N is primarily limited in drought-stressed G. hirsutum by decreased stomatal conductance, along with increases in respiratory and photorespiratory carbon losses, not inhibition or down-regulation of electron transport through photosystem II. It is further concluded that Ψ_PD is a reliable indicator of drought stress and the need for irrigation in field-grown cotton.     

Substrate water exchange in photosystem II core complexes of the extremophilic red alga Cyanidioschyzon merolae

The binding affinity of the two substrate–water molecules to the water-oxidizing Mn₄CaO₅ catalyst in photosystem II core complexes of the extremophilic red alga Cyanidioschyzon merolae was studied in the S₂ and S₃ states by the exchange of bound ¹⁶O-substrate against ¹⁸O-labeled water. The rate of this exchange was detected via the membrane-inlet mass spectrometric analysis of flash-induced oxygen evolution. For both redox states a fast and slow phase of water-exchange was resolved at the mixed labeled m/z 34 mass peak: k_f = 52 ± 8 s^− 1 and k_s = 1.9 ± 0.3 s^− 1 in the S₂ state, and k_f = 42 ± 2 s^− 1 and k_slow = 1.2 ± 0.3 s^− 1 in S₃, respectively. Overall these exchange rates are similar to those observed previously with preparations of other organisms. The most remarkable finding is a significantly slower exchange at the fast substrate–water site in the S₂ state, which confirms beyond doubt that both substrate–water molecules are already bound in the S₂ state. This leads to a very small change of the affinity for both the fast and the slowly exchanging substrates during the S₂ → S₃ transition. Implications for recent models for water-oxidation are briefly discussed.     

Role of specific residues in coenzyme binding, charge-transfer complex formation, and catalysis in Anabaena ferredoxin NADP-reductase

Two transient charge-transfer complexes (CTC) form prior and upon hydride transfer (HT) in the reversible reaction of the FAD-dependent ferredoxin-NADP⁺ reductase (FNR) with NADP⁺/H, FNR_ox-NADPH (CTC-1), and FNR_rd-NADP⁺ (CTC-2). Spectral properties of both CTCs, as well as the corresponding interconversion HT rates, are here reported for several Anabaena FNR site-directed mutants. The need for an adequate initial interaction between the 2′P-AMP portion of NADP⁺/H and FNR that provides subsequent conformational changes leading to CTC formation is further confirmed. Stronger interactions between the isoalloxazine and nicotinamide rings might relate with faster HT processes, but exceptions are found upon distortion of the active centre. Thus, within the analyzed FNR variants, there is no strict correlation between the stability of the transient CTCs formation and the rate of the subsequent HT. Kinetic isotope effects suggest that, while in the WT, vibrational enhanced modulation of the active site contributes to the tunnel probability of HT; complexes of some of the active site mutants with the coenzyme hardly allow the relative movement of isoalloxazine and nicotinamide rings along the HT reaction. The architecture of the WT FNR active site precisely contributes to reduce the stacking probability between the isoalloxazine and nicotinamide rings in the catalytically competent complex, modulating the angle and distance between the N5 of the FAD isoalloxazine and the C4 of the coenzyme nicotinamide to values that ensure efficient HT processes.