Exceptionally stable salt bridges in cytochrome P450cam have functional roles

Biophysical chemistry behind sickle cell anemia and the mechanism of voxelotor action

Sickle cell anemia disease has been a great challenge to the world in the present situation. It occurs only due to the polymerization of sickle hemoglobin (HbS) having Pro-Val-Glu typed mutation, while the polymerization does not occur in normal hemoglobin (HbA) having Pro-Glu-Glu peptides. It is also well confirmed that the oxygenated HbS (OHbS) does not participate in the polymerization, while the deoxygenated HbS (dHbS) does, which causes the shape of red blood cells sickled. After polymerization, the blood has a low oxygen affinity. Keeping this fact into consideration, only those drugs are being synthesized that stabilize the OHbS structure so that the polymerization of HbS can be stopped. The literature data showed no systematic description of the changes occurring during the OHbS conversion to dHbS before polymerization. Hence, an innovative reasonable study between HbA and HbS, when they convert into their deoxygenated forms, was done computationally. In this evaluation, physiochemical parameters in HbA/HbS before and after deoxygenation were studied and compared deeply. The computationally collected data was used to understand the abnormal behaviour of dHbS arising due to the replacement of Glu6 with Val6. Consequently, during the presented computational study, the changes occurring in HbS were found opposite/abnormal as compared to HbA after the deoxygenation of both. The mechanism of Voxelotor (GBT-440) action to stop the HbS polymerization was also explained with the help of computationally collected data. Besides, a comparative study between GBT-440 and another suggested drug was also done to know their antisickling strength. Additionally, the effect of pH, CO, CO2, and 2,3-diphosphoglycerate (2,3-DPG) on HbS structure was also studied computationally.

Machine Learning Subtle Conformational Change due to Phosphorylation in Intrinsically Disordered Proteins

Phosphorylation of intrinsically disordered proteins/regions (IDPs/IDRs) has a profound effect in biological functions such as cell signaling, protein folding or unfolding, and long-range allosteric effects. However, here we focus on two IDPs, namely 83-residue IDR transcription factor Ash1 and 92-residue long N-terminal region of CDK inhibitor Sic1 protein, found in Saccharomyces cerevisiae, for which experimental measurements of average conformational properties, namely, radius of gyration and structure factor, indicate negligible changes upon phosphorylation. Here, we show that a judicious dissection of conformational ensemble via combination of unsupervised machine learning and extensive molecular dynamics (MD) trajectories can highlight key differences and similarities among the phosphorylated and wild-type IDP. In particular, we develop Markov state model (MSM) using the latent-space dimensions of an autoencoder, trained using multi-microsecond long MD simulation trajectories. Examination of structural changes among the states, prior to and upon phosphorylation, captured several similarities and differences in their backbone contact maps, secondary structure, and torsion angles. Hydrogen bonding analysis revealed that phosphorylation not only increases the number of hydrogen bonds but also switches the pattern of hydrogen bonding between the backbone and side chain atoms with the phosphorylated residues. We also observe that although phosphorylation introduces salt bridges, there is a loss of the cation-π interaction. Phosphorylation also improved the probability for long-range hydrophobic contacts and also enhanced interaction with water molecules and improved the local structure of water as evident from the geometric order parameters. The observations on these machine-learnt states gave important insights, as it would otherwise be difficult to determine experimentally which is important, if we were to understand the role of phosphorylation of IDPs in their biological functions.

A newly introduced salt bridge cluster improves structural and biophysical properties of de novo TIM barrels

Protein stability can be fine‐tuned by modifying different structural features such as hydrogen‐bond networks, salt bridges, hydrophobic cores, or disulfide bridges. Among these, stabilization by salt bridges is a major challenge in protein design and engineering since their stabilizing effects show a high dependence on the structural environment in the protein, and therefore are difficult to predict and model. In this work, we explore the effects on structure and stability of an introduced salt bridge cluster in the context of three different de novo TIM barrels. The salt bridge variants exhibit similar thermostability in comparison with their parental designs but important differences in the conformational stability at 25°C can be observed such as a highly stabilizing effect for two of the proteins but a destabilizing effect to the third. Analysis of the formed geometries of the salt bridge cluster in the crystal structures show either highly ordered salt bridge clusters or only single salt bridges. Rosetta modeling of the salt bridge clusters results in a good prediction of the tendency on stability changes but not the geometries observed in the three‐dimensional structures. The results show that despite the similarities in protein fold, the salt bridge clusters differently influence the structural and stability properties of the de novo TIM barrel variants depending on the structural background where they are introduced.

PDB Code(s): 7OSU, 7OT7, 7OSV, 7OT8 and 7P12;

Less Unfavorable Salt Bridges on the Enzyme Surface Result in More Organic Cosolvent Resistance

Biocatalysis for the synthesis of fine chemicals is highly attractive but usually requires organic (co-)solvents (OSs). However, native enzymes often have low activity and resistance in OSs and at elevated temperatures. Herein, we report a smart salt bridge design strategy for simultaneously improving OS resistance and thermostability of the model enzyme, Bacillus subtilits Lipase A (BSLA). We combined comprehensive experimental studies of 3450 BSLA variants and molecular dynamics simulations of 36 systems. Iterative recombination of four beneficial substitutions yielded superior resistant variants with up to 7.6-fold (D64K/D144K) improved resistance toward three OSs while exhibiting significant thermostability (thermal resistance up to 137-fold, and half-life up to 3.3-fold). Molecular dynamics simulations revealed that locally refined flexibility and strengthened hydration jointly govern the highly increased resistance in OSs and at 50-100 °C. The salt bridge redesign provides protein engineers with a powerful and likely general approach to design OSs- and/or thermal-resistant lipases and other α/β-hydrolases.

Insight into the Microcosm of the Human Growth Hormone and Its Interactions with Polymers and Copolymers: A Molecular Dynamics Perspective

Therapeutic proteins nowadays have increasingly been applied for their considerable potential in treating a wide variety of diseases. The effectiveness and potency of native therapeutic proteins are limited by various factors (e.g., stability, blood circulation time, specificity). Over the past years, a great deal of effort has been devoted to developing safe and efficient protein delivery systems. Entrapment of protein into polymeric and copolymeric matrices is common among the different types of protein-based drug formulation. However, despite the massive efforts toward developing therapeutic protein delivery in experimental studies and industrial applications, there is relatively little data on the influence of polymers and copolymers on therapeutic proteins at the atomic and molecular levels. Herein, molecular dynamics (MD) simulations are used to study the effects of biocompatible synthetic polymers including methoxy poly(ethylene glycol) (MPEG), poly(lactic acid) (PLA), and poly(lactic acid) copolymers (poly(lactic-co-glycolic acid)) PLGA and MPEG-PLA(PELA)) on the structure and dynamics of the human growth hormone (hGH), and the results are compared with previous experimental findings. Our results indicate that the hGH conformation remains stable both in pure water and in the presence of polymers, and these results are in good agreement with previous experimental data. It is shown that the MPEG chains are self-assembled and folded into a semicrystalline structure; therefore, only a small portion of the protein interacts with the polymer. The other three polymers, however, interact well with the protein and partially cover its surface. Our findings suggest that the use of these polymers for protein encapsulation has the advantage of reducing protein aggregation and thus increasing drug serum half-life. Eventually, we anticipate that the research results will expand the current knowledge about encapsulation mechanisms and the molecular interactions between hGH and the polymers.

Comparative studies on ion-pair energetic, distribution among three domains of life: Archaea, eubacteria, and eukarya

Salt-bridges play a unique role in the structural and functional stability of proteins, especially under harsh environments. How these salt-bridges contribute to the overall thermodynamic stability of protein structure and function across different domains of life is elusive still date. To address the issue, statistical analyses on the energies of salt-bridges, involved in proteins' structure and function, are performed across three domains of life, that is, archaea, eubacteria, and eukarya. Results show that although the majority of salt-bridges are stable and conserved, yet the stability of archaeal proteins (∆∆G_net = -5.06 ± 3.8) is much more than that of eubacteria (∆∆G_net = -3.7 ± 2.9) and eukarya (∆∆G_net = -3.54 ± 3.1). Unlike earlier study with archaea, in eukarya and eubacteria, not all buried salt-bridge in our dataset are stable. Buried salt-bridges play surprising role in protein stability, whose variations are clearly observed among these domains. Greater desolvation penalty of buried salt-bridges is compensated by stable network of salt-bridges apart from equal contribution of bridge and background energy terms. On the basis proteins' secondary structure, topology, and evolution, our observation shows that salt-bridges when present closer to each other in sequence tend to form a greater number. Overall, our comparative study provides insight into the role of specific electrostatic interactions in proteins from different domains of life, which we hope, would be useful for protein engineering and bioinformatics study.

Mapping the Substrate Recognition Pathway in Cytochrome P450

Cytochrome P450s are ubiquitous metalloenzymes involved in the metabolism and detoxification of foreign components via catalysis of the hydroxylation reactions of a vast array of organic substrates. However, the mechanism underlying the pharmaceutically critical process of substrate access to the catalytic center of cytochrome P450 is a long-standing puzzle, further complicated by the crystallographic evidence of a closed catalytic center in both substrate-free and substrate-bound cytochrome P450. Here, we address a crucial question whether the conformational heterogeneity prevalent in cytochrome P450 translates to heterogeneous pathways for substrate access to the catalytic center of these metalloenzymes. By atomistically capturing the full process of spontaneous substrate association from bulk solvent to the occluded catalytic center of an archetypal system P450cam in multi-microsecond-long continuous unbiased molecular dynamics simulations, we here demonstrate that the substrate recognition in P450cam always occurs through a single well-defined dominant pathway. The simulated final bound pose resulting from these unguided simulations is in striking resemblance with the crystallographic bound pose. Each individual binding trajectory reveals that the substrate, initially placed at random locations in bulk solvent, spontaneously lands on a single key channel on the protein-surface of P450cam and resides there for an uncharacteristically long period, before correctly identifying the occluded target-binding cavity. Surprisingly, the passage of substrate to the closed catalytic center is not accompanied by any large-scale opening in protein. Rather, the unbiased simulated trajectories (∼57 μs) and underlying Markov state model, in combination with free-energy analysis, unequivocally show that the substrate recognition process in P450cam needs a substrate-induced side-chain displacement coupled with a complex array of dynamical interconversions of multiple metastable substrate conformations. Further, the work reconciles multiple precedent experimental and theoretical observations on P450cam and establishes a comprehensive view of substrate-recognition in cytochrome P450 that only occurs via substrate-induced structural rearrangements.

Thermodynamics of camphor migration in cytochrome P450cam by atomistic simulations

Understanding the mechanisms of ligand binding to enzymes is of paramount importance for the design of new drugs. Here, we report on the use of a novel biased molecular dynamics (MD) methodology to study the mechanism of camphor binding to cytochrome P450cam. Microsecond-long MD simulations allowed us to observe reaction coordinates characterizing ligand diffusion from the active site of cytochrome P450cam to solvent via three egress routes. These atomistic simulations were used to estimate thermodynamic quantities along the reaction coordinates and indicate diverse binding configurations. The results suggest that the diffusion of camphor along the pathway near the substrate recognition site (SRS) is thermodynamically preferred. In addition, we show that the diffusion near the SRS is triggered by a transition from a heterogeneous collection of closed ligand-bound conformers to the basin comprising the open conformations of cytochrome P450cam. The conformational change accompanying this switch is characterized by the retraction of the F and G helices and the disorder of the B' helix. These results are corroborated by experimental studies and provide detailed insight into ligand binding and conformational behavior of the cytochrome family. The presented methodology is general and can be applied to other ligand-protein systems.

The crystal structure of the versatile cytochrome P450 enzyme CYP109B1 from Bacillus subtilis

The crystal structure of the versatile CYP109B1 enzyme from Bacillus subtilis has been solved at 1.8 Å resolution. This is the first structure of an enzyme from this CYP family, whose members are prevalent across diverse species of bacteria. In the crystal structure the enzyme has an open conformation with an access channel leading from the heme to the surface. The substrate-free structure reveals the location of the key residues in the active site that are responsible for binding the substrate in the correct orientation for regioselective oxidation. Importantly, there are significant differences among these residues in members of the CYP109 and closely related CYP106 families and these likely account for the variations in substrate binding and oxidation profiles observed with these enzymes. A whole-cell oxidation biosystem was developed, which contains CYP109B1 and a phthalate family oxygenase reductase (PFOR), from Pseudomonas putida KT24440, as the electron transfer partner. This electron transfer system is able to support CYP109B1 activity resulting in the regioselective hydroxylation of both α- and β-ionone in vivo and in vitro. The PFOR is therefore a versatile electron transfer partner that is able to support the activity of CYP enzymes from other bacterium. The crystal structure of CYP109B1 has a positively charged proximal face and this explains why it can interact with PFOR and adrenodoxin which are predominantly negatively charged around their [2Fe-2S] clusters.

Salt-bridge energetics in halophilic proteins

Halophilic proteins have greater abundance of acidic over basic and very low bulky hydrophobic residues. Classical electrostatic stabilization was suggested as the key determinant for halophilic adaptation of protein. However, contribution of specific electrostatic interactions (i.e. salt-bridges) to overall stability of halophilic proteins is yet to be understood. To understand this, we use Adaptive-Poison-Boltzmann-Solver Methods along with our home-built automation to workout net as well as associated component energy terms such as desolvation energy, bridge energy and background energy for 275 salt-bridges from 20 extremely halophilic proteins. We then perform extensive statistical analysis on general and energetic attributes on these salt-bridges. On average, 8 salt-bridges per 150 residues protein were observed which is almost twice than earlier report. Overall contributions of salt-bridges are −3.0 kcal mol⁻¹. Majority (78%) of salt-bridges in our dataset are stable and conserved in nature. Although, average contributions of component energy terms are equal, their individual details vary greatly from one another indicating their sensitivity to local micro-environment. Notably, 35% of salt-bridges in our database are buried and stable. Greater desolvation penalty of these buried salt-bridges are counteracted by stable network salt-bridges apart from favorable equal contributions of bridge and background terms. Recruitment of extensive network salt-bridges (46%) with a net contribution of −5.0 kcal mol⁻¹ per salt-bridge, seems to be a halophilic design wherein favorable average contribution of background term (−10 kcal mol⁻¹) exceeds than that of bridge term (−7 kcal mol⁻¹). Interiors of proteins from halophiles are seen to possess relatively higher abundance of charge and polar side chains than that of mesophiles which seems to be satisfied by cooperative network salt-bridges. Overall, our theoretical analyses provide insight into halophilic signature in its specific electrostatic interactions which we hope would help in protein engineering and bioinformatics studies.

Crystal structures and functional characterization of wild-type CYP101D1 and its active site mutants

Although CYP101D1 and P450cam catalyze the same reaction at similar rates and share strikingly similar active site architectures, there are significant functional differences. CYP101D1 thus provides an opportunity to probe what structural and functional features must be shared and what features can differ but maintain the high catalytic efficiency. Crystal structures of the cyanide complex of wild-type CYP101D1 and it active site mutants, D259N and T260A, have been determined. The conformational changes in CYP101D1 upon cyanide binding are very similar to those of P450cam, indicating a similar mechanism for proton delivery during oxygen activation using solvent-assisted proton transfer. The D259N-CN- complex shows a perturbed solvent structure compared to that of the wild type, which is similar to what was observed in the oxy complex of the corresonding D251N mutant in P450cam. As in P450cam, the T260A mutant is highly uncoupled while the D259N mutant gives barely detectable activity. Despite these similarities, CYP101D1 is able to use the P450cam redox partners while P450cam cannot use the CYP101D1 redox partners. Thus, the strict requirement of P450cam for its own redox partner is relaxed in CYP101D1. Differences in the local environment of the essential Asp (Asp259 in CYP101D1) provide a strucutral basis for understanding these functional differences.

Synergistic effects of mutations in cytochrome P450cam designed to mimic CYP101D1

A close orthologue to cytochrome P450cam (CYP101A1) that catalyzes the same hydroxylation of camphor to 5-exo-hydroxycamphor is CYP101D1. There are potentially important differences in and around the active site that could contribute to subtle functional differences. Adjacent to the heme iron ligand, Cys357, is Leu358 in P450cam, whereas this residue is Ala in CYP101D1. Leu358 plays a role in binding of the P450cam redox partner, putidaredoxin (Pdx). On the opposite side of the heme, about 15-20 Å away, Asp251 in P450cam plays a critical role in a proton relay network required for O2 activation but forms strong ion pairs with Arg186 and Lys178. In CYP101D1 Gly replaces Lys178. Thus, the local electrostatic environment and ion pairing are substantially different in CYP101D1. These sites have been systematically mutated in P450cam to the corresponding residues in CYP101D1 and the mutants analyzed by crystallography, kinetics, and UV-vis spectroscopy. Individually, the mutants have little effect on activity or structure, but in combination there is a major drop in enzyme activity. This loss in activity is due to the mutants being locked in the low-spin state, which prevents electron transfer from the P450cam redox partner, Pdx. These studies illustrate the strong synergistic effects on well-separated parts of the structure in controlling the equilibrium between the open (low-spin) and closed (high-spin) conformational states.

The dynamics of camphor in the cytochrome P450 CYP101D2

The recent crystal structures of CYP101D2, a cytochrome P450 protein from the oligotrophic bacterium Novosphingobium aromaticivorans DSM12444 revealed that both the native (substrate-free) and camphor-soaked forms have open conformations. Furthermore, two other potential camphor-binding sites were also identified from electron densities in the camphor-soaked structure, one being located in the access channel and the other in a cavity on the surface near the F-helix side of the F-G loop termed the substrate recognition site. These latter sites may be key intermediate positions on the pathway for substrate access to or product egress from the active site. Here, we show via the use of unbiased atomistic molecular dynamics simulations that despite the open conformation of the native and camphor-bound crystal structures, the underlying dynamics of CYP101D2 appear to be very similar to other CYP proteins. Simulations of the native structure demonstrated that the protein is capable of sampling many different conformational substates. At the same time, simulations with the camphor positioned at various locations within the access channel or recognition site show that movement towards the active site or towards bulk solvent can readily occur on a short timescale, thus confirming many previously reported in silico studies using steered molecular dynamics. The simulations also demonstrate how the fluctuations of an aromatic gate appear to control access to the active site. Finally, comparison of camphor-bound simulations with the native simulations suggests that the fluctuations can be of similar level and thus are more representative of the conformational selection model rather than induced fit.

Coupled flexibility change in cytochrome P450cam substrate binding determined by neutron scattering, NMR, and molecular dynamics simulation

Neutron scattering and nuclear magnetic resonance relaxation experiments are combined with molecular dynamics (MD) simulations in a novel, to our knowledge, approach to investigate the change in internal dynamics on substrate (camphor) binding to a protein (cytochrome P450cam). The MD simulations agree well with both the neutron scattering, which furnishes information on global flexibility, and the nuclear magnetic resonance data, which provides residue-specific order parameters. Decreased fluctuations are seen in the camphor-bound form using all three techniques, dominated by changes in specific regions of the protein. The combined experimental and simulation results permit a detailed description of the dynamical change, which involves modifications in the coupling between the dominant regions and concomitant substrate access channel closing, via specific salt-bridge, hydrogen-bonding, and hydrophobic interactions. The work demonstrates how the combination of complementary experimental spectroscopies with MD simulation can provide an in-depth description of functional dynamical protein changes.

Sequence-dependent pKa shift induced by molecular self-assembly: insights from computer simulation

The control of catalytic activity using molecular self-assembly is of fundamental interest. Recent experiments (Muller et al., Angew. Chem., Int. Ed., 2009, 48, 922-925) have demonstrated that two sequence isomers of β-peptides show remarkably different activity as an amine catalyst for a retro-aldol cleavage reaction, a difference attributed to the ability of one of the sequences to form large aggregates. The self-assembly and catalytic activity of these two isomers are investigated using constant pH molecular dynamics (CPHMD), for an atomistic model of β-peptides in implicit solvent. Simulations show that the globally amphiphilic (GA) isomer, which experimentally has high activity, forms large aggregates, while the non-GA isomer forms aggregates that are at most three or four molecules in size. The pK(a) shift of the βK-residues is significantly higher in the GA isomers that make a large aggregate. Since the decrease in pK(a) of the side-chain ammonium group is the main driving force for amine catalysis, the calculations are consistent with experiment. We find that the buried βK residues become entirely deprotonated, and the pK(a) shift for other titratable βK residues is accompanied mainly by a clustering of solvent exposed βK residues. We conclude that simulations can be used to understand catalytic activity due to self-assembly.

The structure of CYP101D2 unveils a potential path for substrate entry into the active site

The cytochrome P450 CYP101D2 from Novosphingobium aromaticivorans DSM12444 is closely related to CYP101D1 from the same bacterium and to P450cam (CYP101A1) from Pseudomonas putida. All three are capable of oxidizing camphor stereoselectively to 5-exo-hydroxycamphor. The crystal structure of CYP101D2 revealed that the likely ferredoxin-binding site on the proximal face is largely positively charged, similar to that of CYP101D1. However, both the native and camphor-soaked forms of CYP101D2 had open conformations with an access channel. In the active site of the camphor-soaked form, the camphor carbonyl interacted with the haem-iron-bound water. Two other potential camphor-binding sites were also identified from electron densities in the camphor-soaked structure: one located in the access channel, flanked by the B/C and F/G loops and the I helix, and the other in a cavity on the surface of the enzyme near the F helix side of the F/G loop. The observed open structures may be conformers of the CYP101D2 enzyme that enable the substrate to enter the buried active site via a conformational selection mechanism. The second and third binding sites may be intermediate locations of substrate entry and translocation into the active site, and provide insight into a multi-step substrate-binding mechanism.

Enhanced stability of a protein with increasing temperature

The unusual stability of a structured but locally flexible protein, human growth hormone (hGH) at pH 2.7, was investigated using the temperature dependence of the nanosecond-picosecond dynamics of the backbone amide groups obtained from (15)N NMR relaxation data. It is found that the flexibility of the backbone of the helices decreases with temperature in the range from 24 °C to ∼40 °C, corresponding to an increasing stability. A concomitant increase with temperature of the electrostatic interactions between charged residues forming an interhelical network of salt bridges at the center of the four-helix core suggests that these interactions give rise to the decreasing flexibility and increasing stability of the protein. However, numerous hydrophobic interactions in the interior of the four-helix core may also contribute. Above ∼40 °C, where the thermal energy overcomes the electrostatic and hydrophobic interactions, a substantial increase in the flexibility of the helix backbones results in a highly positive contribution from the local conformational heat capacity, C(p, conf), of the helix backbones to the total heat capacity, C(p), of the protein. This reduces the change in heat capacity upon unfolding, ΔC(p), increases the change in the Gibbs free energy, ΔG(unfold), and stabilizes the protein at high temperatures. A similar decrease in flexibility is found near other salt bridges in hGH and in Calmodulin and may be of general importance for the thermostability of proteins and, in particular, of the salt bridge intensive thermophilic proteins.

Electron transfer in the complex of membrane-bound human cytochrome P450 3A4 with the flavin domain of P450BM-3: the effect of oligomerization of the heme protein and intermittent modulation of the spin equilibrium

We studied the kinetics of NADPH-dependent reduction of human CYP3A4 incorporated into Nanodiscs (CYP3A4-ND) and proteoliposomes in order to probe the effect of P450 oligomerization on its reduction. The flavin domain of cytochrome P450-BM3 (BMR) was used as a model electron donor partner. Unlike CYP3A4 oligomers, where only 50% of the enzyme was shown to be reducible by BMR, CYP3A4-ND could be reduced almost completely. High reducibility was also observed in proteoliposomes with a high lipid-to-protein ratio (L/P=910), where the oligomerization equilibrium is displaced towards monomers. In contrast, the reducibililty in proteoliposomes with L/P=76 did not exceed 55+/-6%. The effect of the surface density of CYP3A4 in proteoliposomes on the oligomerization equilibrium was confirmed with a FRET-based assay employing a cysteine-depleted mutant labeled on Cys-468 with BODIPY iodoacetamide. These results confirm a pivotal role of CYP3A4 oligomerization in its functional heterogeneity. Furthermore, the investigation of the initial phase of the kinetics of CYP3A4 reduction showed that the addition of NADPH causes a rapid low-to-high-spin transition in the CYP3A4-BMR complex, which is followed by a partial slower reversal. This observation reveals a mechanism whereby the CYP3A4 spin equilibrium is modulated by the redox state of the bound flavoprotein.

Pathways and mechanisms for product release in the engineered haloalkane dehalogenases explored using classical and random acceleration molecular dynamics simulations

Eight mutants of the DhaA haloalkane dehalogenase carrying mutations at the residues lining two tunnels, previously observed by protein X-ray crystallography, were constructed and biochemically characterized. The mutants showed distinct catalytic efficiencies with the halogenated substrate 1,2,3-trichloropropane. Release pathways for the two dehalogenation products, 2,3-dichloropropane-1-ol and the chloride ion, and exchange pathways for water molecules, were studied using classical and random acceleration molecular dynamics simulations. Five different pathways, denoted p1, p2a, p2b, p2c, and p3, were identified. The individual pathways showed differing selectivity for the products: the chloride ion releases solely through p1, whereas the alcohol releases through all five pathways. Water molecules play a crucial role for release of both products by breakage of their hydrogen-bonding interactions with the active-site residues and shielding the charged chloride ion during its passage through a hydrophobic tunnel. Exchange of the chloride ions, the alcohol product, and the waters between the buried active site and the bulk solvent can be realized by three different mechanisms: (i) passage through a permanent tunnel, (ii) passage through a transient tunnel, and (iii) migration through a protein matrix. We demonstrate that the accessibility of the pathways and the mechanisms of ligand exchange were modified by mutations. Insertion of bulky aromatic residues in the tunnel corresponding to pathway p1 leads to reduced accessibility to the ligands and a change in mechanism of opening from permanent to transient. We propose that engineering the accessibility of tunnels and the mechanisms of ligand exchange is a powerful strategy for modification of the functional properties of enzymes with buried active sites.

Salt bridges in the hyperthermophilic protein Ssh10b are resilient to temperature increases

A double mutant cycle (DMC) approach was employed to estimate the effect of temperature on the contribution of two highly conserved salt bridges to protein stability in the hyperthermophilic protein Ssh10b. The coupling free energy were 2.4 +/- 0.4 kJ/mol at 298 K and 2.2 +/- 0.4 kJ/mol at 353 K for Glu-54/Arg-57, and 6.0 +/- 0.2 kJ/mol at 298 K and 5.9 +/- 0.6 kJ/mol at 353 K for Glu-36/Lys-68. The stability free energy of Ssh10b decrease greatly with increasing temperature, while the direct contribution of these two salt bridges to protein stability remain almost constant, providing evidence supporting the theoretical prediction that salt bridges are extremely resilient to temperature increases and thus are specially suited to improving protein stability at high temperatures. The reason for the difference in coupling free energy between salt bridges Glu-54/Arg-57 and Glu-36/Lys-68 is discussed. Comparing our results with published DMC data for the contribution of salt bridges to stability in other proteins, we found that the energy contribution of a salt bridge formed by two charged residues far apart in the primary sequence is higher than that of those formed between two very close ones. Implications of this finding are useful for engineering proteins with enhanced thermostability.

Quantum and molecular mechanical study of the first proton transfer in the catalytic cycle of cytochrome P450cam and its mutant D251N

In the catalytic cycle of cytochrome P450cam, the hydroperoxo intermediate (Cpd 0) is formed by proton transfer from a reduced oxyheme complex (S5). This process is drastically slowed down when Asp251 is mutated to Asn (D251N). We report quantum mechanical/molecular mechanical (QM/MM) calculations that address this proton delivery in the doublet state through a hydrogen-bond network in the Asp251 channel, both for the wild-type enzyme and the D251N mutant, using four different active-site models. For the wild-type, we find a facile concerted mechanism for proton transfer from protonated Asp251 via Wat901 and Thr252 to the FeOO moiety, with a barrier of about 1 kcal/mol and a high exothermicity of more than 20 kcal/mol. In the D251N mutant with a neutral Asn251 residue, the proton transfer is almost thermoneutral or slightly exothermic in the three models considered. It is still very facile when the Asn251 residue adopts a conformation analogous to Asp251 in the wild-type enzyme, but the barrier increases significantly when the Asn251 side chain flips (as indicated by classical molecular dynamics simulations). This flip disrupts the hydrogen-bond network and hence the proton-transfer pathway, which causes a longer lifetime of S5 in the D251N mutant (consistent with experimental observations). The entry of an additional water molecule into the active site of D251N with flipped Asn251 regenerates the hydrogen-bond network and provides a viable mechanism for proton delivery in the mutant, with a moderate barrier of about 7 kcal/mol.

Structures of prostacyclin synthase and its complexes with substrate analog and inhibitor reveal a ligand-specific heme conformation change

Prostacyclin synthase (PGIS) is a cytochrome P450 (P450) enzyme that catalyzes production of prostacyclin from prostaglandin H(2). PGIS is unusual in that it catalyzes an isomerization rather than a monooxygenation, which is typical of P450 enzymes. To understand the structural basis for prostacyclin biosynthesis in greater detail, we have determined the crystal structures of ligand-free, inhibitor (minoxidil)-bound and substrate analog U51605-bound PGIS. These structures demonstrate a stereo-specific substrate binding and suggest features of the enzyme that facilitate isomerization. Unlike most microsomal P450s, where large substrate-induced conformational changes take place at the distal side of the heme, conformational changes in PGIS are observed at the proximal side and in the heme itself. The conserved and extensive heme propionate-protein interactions seen in all other P450s, which are largely absent in the ligand-free PGIS, are recovered upon U51605 binding accompanied by water exclusion from the active site. In contrast, when minoxidil binds, the propionate-protein interactions are not recovered and water molecules are largely retained. These findings suggest that PGIS represents a divergent evolution of the P450 family, in which a heme barrier has evolved to ensure strict binding specificity for prostaglandin H(2), leading to a radical-mediated isomerization with high product fidelity. The U51605-bound structure also provides a view of the substrate entrance and product exit channels.

Systematic QM/MM investigation of factors that affect the cytochrome P450-catalyzed hydrogen abstraction of camphor

The hydrogen abstraction reaction of camphor in cytochrome P450(cam) has been investigated in the native enzyme environment by combined quantum mechanical/molecular mechanical (QM/MM) calculations and in the gas phase by density functional calculations. This work has been motivated by contradictory published QM/MM results. In an attempt to pinpoint the origin of these discrepancies, we have systematically studied the factors that may affect the computed barriers, including the QM/MM setup, the optimization procedures, and the choice of QM region, basis set, and protonation states. It is found that the ChemShell and QSite programs used in the published QM/MM calculations yield similar results at given geometries, and that the discrepancies mainly arise from two technical issues (optimization protocols and initial system preparation) that need to be well controlled in QM/MM work. In the course of these systematic investigations, new mechanistic insights have been gained. The crystallographic water 903 placed near the oxo atom of Compound I lowers the hydrogen abstraction barrier by ca. 4 kcal/mol, and thus acts as a catalyst for this reaction. Spin density may appear at the A-propionate side chain of the heme if the carboxylate group is not properly screened, which might be expected to happen during protein dynamics, but not in static equilibrium situations. There is no clear correlation between the computed A-propionate spin density and the hydrogen abstraction barrier, and hence, no support for a previously proposed side-chain mediated transition state stabilization mechanism. Standard QM/MM optimizations yield an A-propionate environment close to the X-ray structure only for protonated Asp297, and not for deprotonated Asp297, but the computed barriers are similar in both cases. An X-ray like A-propionate environment can also be obtained when deprotonated Asp297 is included in the QM region and His355 is singly protonated, but this Compound II-type species with a closed-shell porphyrin ring has a higher hydrogen abstraction barrier and should thus not be mechanistically relevant.

Combined quantum mechanical/molecular mechanical study on the pentacoordinated ferric and ferrous cytochrome P450cam complexes

The pentacoordinated ferric and ferrous cytochrome P450(cam) complexes have been investigated by combined quantum mechanical/molecular mechanical (QM/MM) calculations in the presence of a protein/solvent environment and by QM calculations on the isolated QM regions with use of density functional theory. The B3LYP functional has been found more reliable than the BLYP and BHLYP functionals for estimating the relative state energies. The B3LYP/CHARMM calculations with an all-electron basis set for iron give high-spin ground states for the title complexes, in agreement with experiment. The comparison of the B3LYP/CHARMM results of the entire protein system with the B3LYP calculations on the naked QM regions shows that the amount of stabilization by the protein environment is largest for the intermediate-spin states, followed by the high-spin states of the complexes. The calculation of Mössbauer parameters in the presence of the enzyme environment confirms the double occupation of the d(xz) orbital in the quintet spin state of the ferrous complex, consistent with the computed QM/MM energies in the enzyme environment, while the d(x)2(-)(y)2 orbital is doubly occupied in the gas-phase quintet state.

What is the Active Species of Cytochrome P450 during Camphor Hydroxylation? QM/MM Studies of Different Electronic States of Compound I and of Reduced and Oxidized Iron−Oxo Intermediates

We have investigated C-H hydroxylation of camphor by Compound I (Cpd I) of cytochrome P450cam in different electronic states and by its one-electron reduced and oxidized forms, using QM/MM calculations in the native protein/solvent environment. Cpd I species with five unpaired electrons (pentaradicaloids) are ca. 12 kcal/mol higher in energy than the ground state Cpd I species with three unpaired electrons (triradicaloids). The H-abstraction transition states of pentaradicaloids lie ca. 21 (9) kcal/mol above the triradicaloid (pentaradicaloid) reactants. Hydroxylation via pentaradicaloids is thus facile provided that they can react before relaxing to the ground-state triradicaloids. Excited states of Cpd I with an Fe(V)-oxo moiety lie more than 20 kcal/mol above the triradicaloid ground state in single-point gas-phase calculations, but these electronic configurations are not stable upon including the point-charge protein environment which causes SCF convergence to the triradicaloid ground state. One-electron reduced species (Cpd II) show sluggish reactivity compared with Cpd I in agreement with experimental model studies. One-electron oxidized species are more reactive than Cpd I but seem too high in energy to be accessible. The barriers to hydrogen abstraction for the various forms of Cpd I are generally not affected much by the chosen protonation states of the Asp297 and His355 residues near the propionate side chains of the heme or by the appearance of radical character at Asp297, His355, or the propionates.

Understanding antibody–antigen associations by molecular dynamics simulations: Detection of important intra- and inter-molecular salt bridges

1 NSec molecular dynamics (MD) simulation of anti-hen egg white antibody, HyHEL63 (HH63), complexed with HEL reveals important molecular interactions, not revealed in its X-ray crystal structure. These molecular interactions were predicted to be critical for the complex formation, based on structure-function studies of this complex and 3-other anti-HEL antibodies, HH8, HH10 and HH26, HEL complexes. All four antibodies belong to the same structural family, referred to here as HH10 family. Ala scanning results show that they recognize 'coincident epitopes'. 1 NSec explicit, with periodic boundary condition, MD simulation of HH63- HEL reveals the presence of functionally important saltbridges. Around 200 ps in vacuo and an additional 20 ps explicit simulation agree with the observations from 1 Nsec simulation. Intra-molecular salt-bridges predicted to play significant roles in the complex formation, were revealed during MD simulation. A very stabilizing saltbridge network, and another intra-molecular salt-bridge, at the binding site of HEL, revealed during the MD simulation, is proposed to predipose binding site geometry for specific binding. All the revealed saltbridges are present in one or more of the other three complexes and/or involve \"hot-spot\" epitope and paratope residues. Most of these charged epitope residues make large contribution to the binding free energy. The "hot spot" epitope residue Lys97Y, which significantly contributes to the free energy of binding in all the complexes, forms an intermolecular salt-bridge in several MD conformers. Our earlier computations have shown that this inter-molecular salt-bridge plays a significant role in determining specificity and flexibility of binding in the HH8-HEL and HH26-HEL complexes. Using a robust criterion of salt-bridge detection, this intermolecular salt-bridge was detected in the native structures of the HH8-HEL and HH26-HEL complexes, but was not revealed in the crystal structure of HH63-HEL complex. The electrostatic strength of this revealed saltbridge was very strong. During 1 Nsec MD simulation this salt-bridge networks with another inter-molecular salt-bridge to form an inter-molecular salt-bridge triad. Participation of Lys97Y in the formation of inter-molecular triad further validates the functional importance of Lys97Y in HH63-HEL associations. These results demonstrate that many important structural details of biomolecular interactions can be better understood when studied in a dynamic environment, and that MD simulations can complement and expand information obtained from static X-ray structure. This study also highlights "hot-spot" molecular interactions in HyHEL63-HEL complex.

Structural analysis of cytochromes P450 shows differences in flexibility of heme 2- and 4-vinyls

Structural analysis of the orientations of heme vinyl side chains was carried out using the published crystallographic data for different cytochromes P450. Torsional angles (tau, C(alpha)C(beta)-C(a)C(b)) show different distributions for the vinyls in positions 2 and 4. Whereas the orientation of 2-vinyls is rather restricted (tau between -120 degrees and -180 degrees ), the 4-vinyls have a much higher mobility over almost the entire conformational space. On the basis of the empirical correlation recently reported for peroxidases (M.P. Marzocchi, G. Smulevich, Relationship between heme vinyl conformation and the protein matrix in peroxidases, J. Raman Spectrosc. 34 (2003), 725-736), an attempt has been made to compare the observed vinyl orientations with the experimental frequencies of the vinyl stretching vibrational modes. The data for P450 proteins do not exactly match the peroxidase-derived function, although a qualitatively similar relationship is likely to exist. Differences between P450 forms suggest a variability in heme-region flexibility and in communication with the rest of enzyme.

FTIR studies of the redox partner interaction in cytochrome P450: The Pdx–P450cam couple

Recently we have developed a new approach to study protein-protein interactions using Fourier transform infrared spectroscopy in combination with titration experiments and principal component analysis (FTIR-TPCA). In the present paper we review the FTIR-TPCA results obtained for the interaction between cytochrome P450 and the redox partner protein in two P450 systems, the Pseudomonas putida P450cam (CYP101) with putidaredoxin (P450cam-Pdx), and the Bacillus megaterium P450BM-3 (CYP102) heme domain with the FMN domain (P450BMP-FMND). Both P450 systems reveal similarities in the structural changes that occur upon redox partner complex formation. These involve an increase in beta-sheets and alpha-helix content, a decrease in the population of random coil/3(10)-helix structure, a redistribution of turn structures within the interacting proteins and changes in the protonation states or hydrogen-bonding of amino acid carboxylic side chains. We discuss in detail the P450cam-Pdx interaction in comparison with literature data and conclusions drawn from experiments obtained by other spectroscopic techniques. The results are also interpreted in the context of a 3D structural model of the Pdx-P450cam complex.

Common system setup for the entire catalytic cycle of cytochrome P450cam in quantum mechanical/molecular mechanical studies

We describe a system setup that is applicable to all species in the catalytic cycle of cytochrome P450(cam). The chosen procedure starts from the X-ray coordinates of the ferrous dioxygen complex and follows a protocol that includes the careful assignment of protonation states, comparison between different conceivable hydration schemes, and system preparation through a series of classical minimizations and molecular dynamics (MD) simulations. The resulting setup was validated by quantum mechanical/molecular mechanical (QM/MM) calculations on the resting state, the pentacoordinated ferric and ferrous complexes, Compound I, the transition state and hydroxo intermediate of the C--H hydroxylation reaction, and the product complex. The present QM/MM results are generally consistent with those obtained previously with individual setups. Concerning hydration, we find that saturating the protein interior with water is detrimental and leads to higher structural flexibility and catalytically inefficient active-site geometries. The MD simulations favor a low water density around Asp251 that facilitates side chain rotation of protonated Asp251 during the conversion of Compound 0 to Compound I. The QM/MM results for the two preferred hydration schemes (labeled SE-1 and SE-4) are similar, indicating that slight differences in the solvation close to the active site are not critical as long as camphor and the crystallographic water molecules preserve their positions in the experimental X-ray structures.

Mechanisms of reaction in cytochrome P450: Hydroxylation of camphor in P450cam

The fundamental nature of reactivity in cytochrome P450 enzymes is currently controversial. Modelling of bacterial P450cam has suggested an important role for the haem propionates in the catalysis, though this finding has been questioned. Understanding the mechanisms of this enzyme family is important both in terms of basic biochemistry and potentially in the prediction of drug metabolism. We have modelled the hydroxylation of camphor by P450cam, using combined quantum mechanics/molecular mechanics (QM/MM) methods. A set of reaction pathways in the enzyme was determined. We were able to pinpoint the source of the discrepancies in the previous results. We show that when a correct ionization state is assigned to Asp297, no spin density appears on the haem propionates and the protein structure in this region remains preserved. These results indicate that the haem propionates are not involved in catalysis.

The effect of heme environment on the hydrogen abstraction reaction of camphor in P450cam catalysis: a QM/MM study

The discrepancies between the published QM/MM studies (Schöneboom, J. C.; Cohen, S.; Lin, H.; Shaik, S.; Thiel, W. J. Am. Chem. Soc. 2004, 126, 4017; Guallar, V.; Friesner, R. A. J. Am. Chem. Soc. 2004, 126, 8501) on H-abstraction of camphor in P450cam have largely been resolved. The crystallographic water molecule 903 situated near the oxo atom of Compound I acts as a catalyst for H-abstraction, lowering the barrier by about 4 kcal/mol. Spin density at the A-propionate side chain of heme can occur in the case of incomplete screening but has no major effect on the computed barrier.

Contribution of protein-surface ion pairs of a hyperthermophilic protein on thermal and thermodynamic stability

Hyperthermophilic proteins possess many ion pairs on their surface. To reveal the role of the ion pairs, O6-methylguanine-DNA methyltransferase from Thermococcus kodakaraensis KOD1 (Tk-MGMT) was studied as a model protein. The maximum free-energy changes of the protein in 0.1 and 0.5 M NaCl at pH 7.0 were 61.7 kJ mol(-1) at 31.5 degrees C and 77.4 kJ mol(-1) at 39.7 degrees C, respectively. On the other hand, mid points of the thermal unfolding temperatures in 0.1 and 0.5 M NaCl at pH 7.0 were 94.8 degrees C and 90.1 degrees C, respectively. The results suggest that the protein-surface ion pairs contribute to thermal stability (Tm), rather than thermodynamic stability (DeltaG).

Statistical characterization of salt bridges in proteins

The structure and folding mechanism of a given protein are determined by many factors, including the electrostatic interactions between charged residues of protein molecules known in general as salt bridges. In this study, analyses were conducted on 10,370 salt bridges in 2017 proteins and the results compared to previous statistical surveys of 36 protein structures. Although many of the general trends remained consistent with other studies, more detailed information was illuminated by the larger dataset. In particular, it was shown that there is a strong correlation between secondary structure and salt bridge formation, and that salt bridges display preferential formation in an environment of about 30% solvent accessible surface area.

Conformational States of Cytochrome P450cam Revealed by Trapping of Synthetic Molecular Wires

Members of the ubiquitous cytochrome P450 family catalyze a vast range of biologically significant reactions in mammals, plants, fungi, and bacteria. Some P450s display a remarkable promiscuity in substrate recognition, while others are very specific with respect to substrate binding or regio and stereo-selective catalysis. Recent results have suggested that conformational flexibility in the substrate access channel of many P450s may play an important role in controlling these effects. Here, we report the X-ray crystal structures at 1.8A and 1.5A of cytochrome P450cam complexed with two synthetic molecular wires, D-4-Ad and D-8-Ad, consisting of a dansyl fluorophore linked to an adamantyl substrate analog via an alpha,omega-diaminoalkane chain of varying length. Both wires bind with the adamantyl moiety in similar positions at the camphor-binding site. However, each wire induces a distinct conformational response in the protein that differs from the camphor-bound structure. The changes involve significant movements of the F, G, and I helices, allowing the substrate access channel to adapt to the variable length of the probe. Wire-induced opening of the substrate channel also alters the I helix bulge and Thr252 at the active site with binding of water that has been proposed to assist in peroxy bond cleavage. The structures suggest that the coupling of substrate-induced conformational changes to active-site residues may be different in P450cam and recently described mammalian P450 structures. The wire-induced changes may be representative of the conformational intermediates that must exist transiently during substrate entry and product egress, providing a view of how substrates enter the deeply buried active site. They also support observed examples of conformational plasticity that are believed be responsible for the promiscuity of drug metabolizing P450s. Observation of such large changes in P450cam suggests that substrate channel plasticity is a general property inherent to all P450 structures.

Different roles of electrostatics in heat and in cold: adaptation by citrate synthase

Electrostatics plays a major role in heat adaptation by thermophilic proteins. Here we ask whether electrostatics similarly contributes to cold adaptation in psychrophilic proteins. We compare the sequences and structures of citrate synthases from the psychrophile Arthobacter Ds2-3R, from chicken, and from the hyperthermophile Pyrococcus furiosus. The three enzymes share similar packing, burial of nonpolar surface area, and main-chain hydrogen bonding. However, both psychrophilic and hyperthermophilic citrate synthases contain more charged residues, salt bridges, and salt-bridge networks than the mesophile. The electrostatic free-energy contributions toward protein stability by individual charged residues show greater variabilities in the psychrophilic citrate synthase than in the hyperthermophilic enzyme. The charged residues in the active-site regions of the psychrophile are more destabilizing than those in the active-site regions of the hyperthermophile. In the hyperthermophilic enzyme, salt bridges and their networks largely cluster in the active-site regions and at the dimer interface. In contrast, in the psychrophile, they are more dispersed throughout the structure. On average, salt bridges and their networks provide greater electrostatic stabilization to the thermophilic citrate synthase at 100 degrees C than to the psychrophilic enzyme at 0 degrees C. Electrostatics appears to play an important role in both heat and cold adaptation of citrate synthase. However, remarkably, the role may be different in the two types of enzyme: In the hyperthermophile, it may contribute to the integrity of both the protein dimer and the active site by possibly countering conformational disorder at high temperatures. On the other hand, in the psychrophile at low temperatures, electrostatics may contribute to enhance protein solvation and to ensure active-site flexibility.

Theoretical and experimental study of the D2194G mutation in the C2 domain of coagulation factor V

Coagulation factor V (FV) is a large plasma glycoprotein with functions in both the pro- and anticoagulant pathways. In carriers of the so-called R2-FV haplotype, the FV D2194G mutation, in the C2 membrane-binding domain, is associated with low expression levels, suggesting a potential folding/stability problem. To analyze the molecular mechanisms potentially responsible for this in vitro phenotype, we used molecular dynamics (MD) and continuum electrostatic calculations. Implicit solvent simulations were performed on the x-ray structure of the wild-type C2 domain and on a model of the D2194G mutant. Because D2194 is located next to a disulfide bond (S-S bond), MD calculations were also performed on S-S bond depleted structures. D2194 is part of a salt-bridge network and investigations of the stabilizing/destabilizing role of these ionic interactions were carried out. Five mutant FV molecules were created and the expression levels measured with the aim of assessing the tolerance to amino acid changes in this region of molecule. Analysis of the MD trajectories indicated increased flexibility in some areas and energetic comparisons suggested overall destabilization of the structure due to the D2194G mutation. This substitution causes electrostatic destabilization of the domain by approximately 3 kcal/mol. Together these effects likely explain the lowered expression levels in R2-FV carriers.