On the accurate estimation of free energies using the jarzynski equality.

The Jarzynski equality is one of the most widely celebrated and scrutinized nonequilibrium work theorems, relating free energy to the external work performed in nonequilibrium transitions. In practice, the required ensemble average of the Boltzmann weights of infinite nonequilibrium transitions is estimated as a finite sample average, resulting in the so-called Jarzynski estimator, ΔF^J . Alternatively, the second-order approximation of the Jarzynski equality, though seldom invoked, is exact for Gaussian distributions and gives rise to the Fluctuation-Dissipation estimator ΔF^FD . Here we derive the parametric maximum-likelihood estimator (MLE) of the free energy ΔF^ML considering unidirectional work distributions belonging to Gaussian or Gamma families, and compare this estimator to ΔF^J . We further consider bidirectional work distributions belonging to the same families, and compare the corresponding bidirectional ΔF^ML∗ to the Bennett acceptance ratio ( ΔF^BAR ) estimator. We show that, for Gaussian unidirectional work distributions, ΔF^FD is in fact the parametric MLE of the free energy, and as such, the most efficient estimator for this statistical family. We observe that ΔF^ML and ΔF^ML∗ perform better than ΔF^J and ΔF^BAR , for unidirectional and bidirectional distributions, respectively. These results illustrate that the characterization of the underlying work distribution permits an optimal use of the Jarzynski equality. © 2018 Wiley Periodicals, Inc.

Mechanism of Sulfide Binding by Ferric Hemeproteins.

The reaction of hydrogen sulfide (H2S) with hemeproteins is a key physiological reaction; still, its mechanism and implications are not completely understood. In this work, we propose a combination of experimental and theoretical tools to shed light on the reaction in model system microperoxidase 11 (MP11-FeIII) and myoglobin (Mb-FeIII), from the estimation of the intrinsic binding constants of the species H2S and hydrosulfide (HS-), and the computational description of the overall binding process. Our results show that H2S and HS- are the main reactive species in Mb-FeIII and MP11-FeIII, respectively, and that the magnitude of their intrinsic binding constants are similar to most of the binding constants reported so far for hemeproteins systems and model compounds. However, while the binding of HS- to Mb-FeIII was negligible, the binding of H2S to MP11-FeIII was significant, providing a frame for a discriminated analysis of both species and revealing differential mechanistic aspects. A joint inspection of the kinetic data and the free energy profiles of the binding processes suggests that a dissociative mechanism with the release of a coordinated water molecule as rate limiting step is operative in the binding of H2S to Mb-FeIII and that the binding of HS- is prevented in the access to the protein matrix. For the MP11-FeIII case, where no access restrictions for the ligands are present, an associative component in the mechanism seems to be operative. Overall, the results suggest that if accessing the active site then both H2S and HS- are capable of binding a ferric heme moiety.

A quantitative model for oxygen uptake and release in a family of hemeproteins.

MOTIVATION: Hemeproteins have many diverse functions that largely depend on the rate at which they uptake or release small ligands, like oxygen. These proteins have been extensively studied using either simulations or experiments, albeit only qualitatively and one or two proteins at a time. RESULTS: We present a physical-chemical model, which uses data obtained exclusively from computer simulations, to describe the uptake and release of oxygen in a family of hemeproteins, called truncated hemoglobins (trHbs). Through a rigorous statistical analysis we demonstrate that our model successfully recaptures all the reported experimental oxygen association and dissociation kinetic rate constants, thus allowing us to establish the key factors that determine the rates at which these hemeproteins uptake and release oxygen. We found that internal tunnels as well as the distal site water molecules control ligand uptake, whereas oxygen stabilization by distal site residues controls ligand release. Because these rates largely determine the functions of these hemeproteins, these approaches will also be important tools in characterizing the trHbs members with unknown functions. CONTACT: lboechi@ic.fcen.uba.ar SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Evolutionary and Functional Relationships in the Truncated Hemoglobin Family.

Predicting function from sequence is an important goal in current biological research, and although, broad functional assignment is possible when a protein is assigned to a family, predicting functional specificity with accuracy is not straightforward. If function is provided by key structural properties and the relevant properties can be computed using the sequence as the starting point, it should in principle be possible to predict function in detail. The truncated hemoglobin family presents an interesting benchmark study due to their ubiquity, sequence diversity in the context of a conserved fold and the number of characterized members. Their functions are tightly related to O2 affinity and reactivity, as determined by the association and dissociation rate constants, both of which can be predicted and analyzed using in-silico based tools. In the present work we have applied a strategy, which combines homology modeling with molecular based energy calculations, to predict and analyze function of all known truncated hemoglobins in an evolutionary context. Our results show that truncated hemoglobins present conserved family features, but that its structure is flexible enough to allow the switch from high to low affinity in a few evolutionary steps. Most proteins display moderate to high oxygen affinities and multiple ligand migration paths, which, besides some minor trends, show heterogeneous distributions throughout the phylogenetic tree, again suggesting fast functional adaptation. Our data not only deepens our comprehension of the structural basis governing ligand affinity, but they also highlight some interesting functional evolutionary trends.

Molecular basis of thermal stability in truncated (2/2) hemoglobins.

BACKGROUND: Understanding the molecular mechanism through which proteins are functional at extreme high and low temperatures is one of the key issues in structural biology. To investigate this phenomenon, we have focused on two instructive truncated hemoglobins from Thermobifida fusca (Tf-trHbO) and Mycobacterium tuberculosis (Mt-trHbO); although the two proteins are structurally nearly identical, only the former is stable at high temperatures. METHODS: We used molecular dynamics simulations at different temperatures as well as thermal melting profile measurements of both wild type proteins and two mutants designed to interchange the amino acid residue, either Pro or Gly, at E3 position. RESULTS: The results show that the presence of a Pro at the E3 position is able to increase (by 8°) or decrease (by 4°) the melting temperature of Mt-trHbO and Tf-trHbO, respectively. We observed that the ProE3 alters the structure of the CD loop, making it more flexible. CONCLUSIONS: This gain in flexibility allows the protein to concentrate its fluctuations in this single loop and avoid unfolding. The alternate conformations of the CD loop also favor the formation of more salt-bridge interactions, together augmenting the protein's thermostability. GENERAL SIGNIFICANCE: These results indicate a clear structural and dynamical role of a key residue for thermal stability in truncated hemoglobins.

Reactivity of inorganic sulfide species toward a heme protein model.

The reactivity of inorganic sulfide species toward heme peptides was explored under biorelevant conditions in order to unravel the molecular details of the reactivity of the endogenous hydrogen sulfide toward heme proteins. Unlike ferric porphyrinates, which are reduced by inorganic sulfide, some heme proteins can form stable Fe(III)-sulfide adducts. To isolate the protein factors ruling the redox chemistry, we used as a system model, the undecapeptide microperoxidase (MP11), a heme peptide derived from cytochrome c proteolysis that retains the proximal histidine bound to the Fe(III) atom. Upon addition of gaseous hydrogen sulfide (H2S) at pH 6.8, the UV-vis spectra of MP11 closely resembled those of the low-spin ferric hydroxo complex (only attained at an alkaline pH) and cysteine or alkylthiol derivatives, suggesting that the Fe(III) reduction was prevented. The low-frequency region of the resonance Raman spectrum revealed the presence of an Fe(III)-S band at 366 cm(-1) and the general features of a low-spin hexacoordinated heme. Anhydrous sodium sulfide (Na2S) was the source of sulfide of choice for the kinetic evaluation of the process. Theoretical calculations showed no distal stabilization mechanisms for bound sulfide species in MP11, highlighting a key role of the proximal histidine for the stabilization of the Fe(III)-S adducts of heme compounds devoid of distal counterparts, which is significant with regard to the biochemical reactivity of endogenous hydrogen sulfide.

Hydrophobic effect drives oxygen uptake in myoglobin via histidine E7.

Since the elucidation of the myoglobin (Mb) structure, a histidine residue on the E helix (His-E7) has been proposed to act as a gate with an open or closed conformation controlling access to the active site. Although it is believed that at low pH, the His-E7 gate is in its open conformation, the full relationship between the His-E7 protonation state, its conformation, and ligand migration in Mb is hotly debated. We used molecular dynamics simulations to first address the effect of His-E7 protonation on its conformation. We observed the expected shift from the closed to the open conformation upon protonation, but more importantly, noted a significant difference between the conformations of the two neutral histidine tautomers. We further computed free energy profiles for oxygen migration in each of the possible His-E7 states as well as in two instructive Mb mutants: Ala-E7 and Trp-E7. Our results show that even in the closed conformation, the His-E7 gate does not create a large barrier to oxygen migration and permits oxygen entry with only a small rotation of the imidazole side chain and movement of the E helix. We identify, instead, a hydrophobic site in the E7 channel that can accommodate an apolar diatomic ligand and enhances ligand uptake particularly in the open His-E7 conformation. This rate enhancement is diminished in the closed conformation. Taken together, our results provide a new conceptual framework for the histidine gate hypothesis.

Small ligand-globin interactions: reviewing lessons derived from computer simulation.

In this work we review the application of classical and quantum-mechanical atomistic computer simulation tools to the investigation of small ligand interaction with globins. In the first part, studies of ligand migration, with its connection to kinetic association rate constants (kon), are presented. In the second part, we review studies for a variety of ligands such as O2, NO, CO, HS(-), F(-), and NO2(-) showing how the heme structure, proximal effects, and the interactions with the distal amino acids can modulate protein ligand binding. The review presents mainly results derived from our previous works on the subject, in the context of other theoretical and experimental studies performed by others. The variety and extent of the presented data yield a clear example of how computer simulation tools have, in the last decade, contributed to our deeper understanding of small ligand interactions with globins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Structure, mechanism, and dynamics of UDP-galactopyranose mutase.

The flavoenzyme UDP-galactopyranose mutase (UGM) is a key enzyme in galactofuranose biosynthesis. The enzyme catalyzes the 6-to-5 ring contraction of UDP-galactopyranose to UDP-galactofuranose. Galactofuranose is absent in humans yet is an essential component of bacterial and fungal cell walls and a cell surface virulence factor in protozoan parasites. Thus, inhibition of galactofuranose biosynthesis is a valid strategy for developing new antimicrobials. UGM is an excellent target in this effort because the product of the UGM reaction represents the first appearance of galactofuranose in the biosynthetic pathway. The UGM reaction is redox neutral, which is atypical for flavoenzymes, motivating intense examination of the chemical mechanism and structural features that tune the flavin for its unique role in catalysis. These studies show that the flavin functions as nucleophile, forming a flavin-sugar adduct that facilitates galactose-ring opening and contraction. The 3-dimensional fold is novel and conserved among all UGMs, however the larger eukaryotic enzymes have additional secondary structure elements that lead to significant differences in quaternary structure, substrate conformation, and conformational flexibility. Here we present a comprehensive review of UGM three-dimensional structure, provide an update on recent developments in understanding the mechanism of the enzyme, and summarize computational studies of active site flexibility.

Protein dynamics and ligand migration interplay as studied by computer simulation.

Since proteins are dynamic systems in living organisms, the employment of methodologies contemplating this crucial characteristic results fundamental to allow revealing several aspects of their function. In this work, we present results obtained using classical mechanical atomistic simulation tools applied to understand the connection between protein dynamics and ligand migration. Firstly, we will present a review of the different sampling schemes used in the last years to obtain both ligand migration pathways and the thermodynamic information associated with the process. Secondly, we will focus on representative examples in which the schemes previously presented are employed, concerning the following: i) ligand migration, tunnels, and cavities in myoglobin and neuroglobin; ii) ligand migration in truncated hemoglobin members; iii) NO escape and conformational changes in nitrophorins; iv) ligand selectivity in catalase and hydrogenase; and v) larger ligand migration: the P450 and haloalkane dehalogenase cases. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Uncovering differential equations from data with hidden variables.

SINDy is a method for learning system of differential equations from data by solving a sparse linear regression optimization problem [Brunton, Proctor, and Kutz, Proc. Natl. Acad. Sci. USA 113, 3932 (2016)PNASA60027-842410.1073/pnas.1517384113]. In this article, we propose an extension of the SINDy method that learns systems of differential equations in cases where some of the variables are not observed. Our extension is based on regressing a higher order time derivative of a target variable onto a dictionary of functions that includes lower order time derivatives of the target variable. We evaluate our method by measuring the prediction accuracy of the learned dynamical systems on synthetic data and on a real data set of temperature time series provided by the Réseau de Transport d'Électricité. Our method provides high quality short-term forecasts and it is orders of magnitude faster than competing methods for learning differential equations with latent variables.

Role of the distal hydrogen-bonding network in regulating oxygen affinity in the truncated hemoglobin III from Campylobacter jejuni.

Oxygen affinity in heme-containing proteins is determined by a number of factors, such as the nature and conformation of the distal residues that stabilize the heme bound-oxygen via hydrogen-bonding interactions. The truncated hemoglobin III from Campylobacter jejuni (Ctb) contains three potential hydrogen-bond donors in the distal site: TyrB10, TrpG8, and HisE7. Previous studies suggested that Ctb exhibits an extremely slow oxygen dissociation rate due to an interlaced hydrogen-bonding network involving the three distal residues. Here we have studied the structural and kinetic properties of the G8(WF) mutant of Ctb and employed state-of-the-art computer simulation methods to investigate the properties of the O(2) adduct of the G8(WF) mutant, with respect to those of the wild-type protein and the previously studied E7(HL) and/or B10(YF) mutants. Our data indicate that the unique oxygen binding properties of Ctb are determined by the interplay of hydrogen-bonding interactions between the heme-bound ligand and the surrounding TyrB10, TrpG8, and HisE7 residues.

Fluoride as a probe for H-bonding interactions in the active site of heme proteins: the case of Thermobifida fusca hemoglobin.

The structural and functional properties of the active site of the bacterial hemoglobin from Thermobifida fusca are largely determined by three polar amino acids: TrpG8, TyrCD1, and TyrB10. We have exploited the availability of a combinatorial set of mutants, in each of which these three amino acids have been singly, doubly, or triply replaced by a Phe residue, to perform a detailed study on H-bonding interactions between the protein and heme-bound fluoride. By appropriate choice of the excitation conditions, ν(Fe-F) stretching bands have been detected in the resonance Raman spectra. In the wild-type protein and one of the mutants, two ν(Fe-F) bands have been observed and assigned to the presence of two protein conformers where fluoride is singly or doubly H-bonded. Furthermore, by plotting the CT1 charge-transfer transition energy vs the ν(Fe-F) wavenumbers, an empirical correlation has been found. The data are well fitted by a straight line with a positive slope. The position along the correlation line can be considered as a novel, general spectroscopic indicator of the extent of H-bonding in the active site of heme proteins. In agreement with the spectroscopic results, we have observed that the rate of ligand dissociation in stopped-flow kinetic measurements progressively increases upon substitution of the H-bonding amino acids. Molecular dynamics simulations have been performed on the fluoride complexes of native and mutated forms, indicating the prevalent interactions at the active site. All the techniques yield evidence that TrpG8 and TyrCD1 can form strong H bonds with fluoride, whereas TyrB10 plays only a minor role in the stabilization of the ligand.

Comparing and combining implicit ligand sampling with multiple steered molecular dynamics to study ligand migration processes in heme proteins.

The ubiquitous heme proteins perform a wide variety of tasks that rely on the subtle regulation of their affinity for small ligands like O2, CO, and NO. Ligand affinity is characterized by kinetic association and dissociation rate constants, that partially depend on ligand migration between the solvent and active site, mediated by the presence of internal cavities or tunnels. Different computational methods have been developed to study these processes which can be roughly divided in two strategies: those costly methods in which the ligand is treated explicitly during the simulations, and the free energy landscape of the process is computed; and those faster methods that use prior computed Molecular Dynamics simulation without the ligand, and incorporate it afterwards, called implicit ligand sampling (ILS) methods. To compare both approaches performance and to provide a combined protocol to study ligand migration in heme proteins, we performed ILS and multiple steered molecular dynamics (MSMD) free energy calculations of the ligand migration process in three representative and well theoretically and experimentally studied cases that cover a wide range of complex situations presenting a challenging benchmark for the aim of the present work. Our results show that ILS provides a good description of the tunnel topology and a reasonable approximation to the free energy landscape, while MSMD provides more accurate and detailed free energy profile description of each tunnel. Based on these results, a combined strategy is presented for the study of internal ligand migration in heme proteins.

Protonation of histidine 55 affects the oxygen access to heme in the alpha chain of the hemoglobin from the Antarctic fish Trematomus bernacchii.

The Root effect describes the drastic drop of oxygen affinity and loss of cooperativity at acidic pH expressed in the hemoglobins (Hb) of certain fish. The comparison between the deoxy structures of the Root effect Hb from the Antarctic fish Trematomus bernacchii (HbTb) at different pHs (pH = 6.2 and pH = 8.4) shows that the most significant differences are localized at the CDα region, where a salt bridge between Asp48 and His55 breaks during the low-to-high pH transition. In order to shed light on the relationship between pH, CDα loop structure and dynamics, and oxygen access to the active site in the alpha chain of HbTb, different computer simulation techniques were performed. Our results highlight the importance of the protonation of His55 in regulating oxygen access, underscoring its pivotal role in the structural and functional properties of HbTb. These data provide further support to the hypothesis that this residue might contribute to the release of Root protons in HbTb and underline the fact that an efficient transport of molecular oxygen in Hbs relies on a subtle balance of tertiary structure and protein conformational flexibility.

Structure and dynamics of Antarctic fish neuroglobin assessed by computer simulations.

Neuroglobin (Ngb) is a heme protein, highly conserved along evolution, predominantly found in the nervous system. It is upregulated by hypoxia and ischemia and may have a neuroprotective role under hypoxic stress. Although many other roles have been proposed, the physiological function is still unclear. Antarctic icefishes lack hemoglobin and some species also lack myoglobin, but all have Ngb and thus may help the elucidation of Ngb function. We present the first theoretically derived structure of fish Ngb and describe its behavior using molecular dynamics simulations. Specifically, we sequenced and analyzed Ngbs from a colorless-blooded Antarctic icefish species Chaenocephalus aceratus and a related red-blooded species (Dissostichus mawsoni). Both fish Ngbs are 6-coordinated but have some peculiarities that differentiate them from mammalian counterparts: they have extensions in the N and C termini that can interact with the EF loop, and a gap in the alignment that changes the CD-region structure/dynamics that has been found to play a key role in human neuroglobin. Our results suggest that a single mutation between both fish Ngbs is responsible for significant difference in the behavior of the proteins. The functional role of these characteristics is discussed.

The peculiar heme pocket of the 2/2 hemoglobin of cold-adapted Pseudoalteromonas haloplanktis TAC125.

The genome of the cold-adapted bacterium Pseudoalteromonas haloplanktis TAC125 contains multiple genes encoding three distinct monomeric hemoglobins exhibiting a 2/2 α-helical fold. In the present work, one of these hemoglobins is studied by resonance Raman, electronic absorption and electronic paramagnetic resonance spectroscopies, kinetic measurements, and different bioinformatic approaches. It is the first cold-adapted bacterial hemoglobin to be characterized. The results indicate that this protein belongs to the 2/2 hemoglobin family, Group II, characterized by the presence of a tryptophanyl residue on the bottom of the heme distal pocket in position G8 and two tyrosyl residues (TyrCD1 and TyrB10). However, unlike other bacterial hemoglobins, the ferric state, in addition to the aquo hexacoordinated high-spin form, shows multiple hexacoordinated low-spin forms, where either TyrCD1 or TyrB10 can likely coordinate the iron. This is the first example in which both TyrCD1 and TyrB10 are proposed to be the residues that are alternatively involved in heme hexacoordination by endogenous ligands.

Exploring the molecular basis of human manganese superoxide dismutase inactivation mediated by tyrosine 34 nitration.

Manganese Superoxide Dismutase (MnSOD) is an essential mitochondrial antioxidant enzyme that protects organisms against oxidative damage, dismutating superoxide radical (O2(.)⁻) into H2O2 and O2. The active site of the protein presents a Mn ion in a distorted trigonal-bipyramidal environment, coordinated by H26, H74, H163, D159 and one ⁻OH ion or H2O molecule. The catalytic cycle of the enzyme is a "ping-pong" mechanism involving Mn³+/Mn²+. It is known that nitration of Y34 is responsible for enzyme inactivation, and that this protein oxidative modification is found in tissues undergoing inflammatory and degenerative processes. However, the molecular basis about MnSOD tyrosine nitration affects the protein catalytic function is mostly unknown. In this work we strongly suggest, using computer simulation tools, that Y34 nitration affects protein function by restricting ligand access to the active site. In particular, deprotonation of 3-nitrotyrosine increases drastically the energetic barrier for ligand entry due to the absence of the proton. Our results for the WT and selected mutant proteins confirm that the phenolic moiety of Y34 plays a key role in assisting superoxide migration.

Hospitalized population diagnosed with Covid-19 in public health centers in the southeastern region of Greater Buenos Aires / [Población hospitalizada con diagnóstico de Covid-19 en los centros de salud públicos de la región sudeste del Gran Buenos Aires]

Introduction: The present work describes the clinical characteristics and interventions to minimize morbidity and mortality in hospitalized patients diagnosed with COVID-19. Methods: It is a prospective cohort investigation of patients who received a response from the Health Centers in the southeast region (RS) of the metropolitan area (AMBA) from April 8 to September 30, 2020. A Situation Room was used epidemiological with two monitoring and follow-up boards, one for bed management and the other for patient management. Results: During the analyzed period, 2,588 patients with confirmed COVID-19 diagnosis were admitted, 1,943 with suspected COVID-19 pathology, and 1,464 subjects with other pathologies. 55% of the patients were men and the mean age was 51 years. There were 82.8% patients with pre-existing diseases, hypertension and diabetes were the most frequent. 14% were hospitalized in the Intensive Care Unit. The mortality of the cohort was 15.05%, mortality was higher for men, with a mean age of 60 years, 92.65% had some pre-existing disease. Conclusion: Our cohort is younger than other published works. Older people, men, and people with comorbidities are at increased risk for COVID-19-related mortality. The public health system was able to respond to the demand without collapsing the hospital institutions.

Introducción: En el presente trabajo se describen las características clínicas y las intervenciones para minimizar la morbimortalidad en pacientes hospitalizados con diagnóstico de COVID-19. Métodos: Es una investigación de cohorte prospectiva de pacientes que recibieron respuesta de los Centros de Salud en la región sudeste (RS) del área metropolitana (AMBA) desde el 8 de abril hasta el 30 de septiembre de 2020. Se utilizó una Sala de Situación epidemiológica con dos tableros de monitoreo y seguimiento, uno de gestión de camas y otro de gestión de pacientes. Resultados: Durante el periodo analizado se internaron2.588pacientes con diagnóstico COVID-19 confirmados, 1.943 con sospecha de patología COVID-19, y 1.464sujetos con otras patologías. El 55% de los pacientes eran hombres y la edad media fue de 51 años. Hubo 82,8% pacientes con enfermedades preexistentes, hipertensión y diabetes fueron las más frecuentes. El 14% fue hospitalizado en la Unidad de Terapia Intensiva. La mortalidad de la cohorte fue del 15,05%, la mortalidad fue mayor para los hombres, con una edad media de 60 años, el 92,65% tenía alguna enfermedad preexistente. Conclusión: Nuestra cohorte es más joven que otros trabajos publicados. Las personas mayores, los hombres y las personas con comorbilidades tienen mayor riesgo de mortalidad relacionada con COVID-19. El sistema de salud público pudo responder a la demanda sin llegar a colapsar las instituciones hospitalarias.

Heme pocket structural properties of a bacterial truncated hemoglobin from Thermobifida fusca.

An acidic surface variant (ASV) of the "truncated" hemoglobin from Thermobifida fusca was designed with the aim of creating a versatile globin scaffold endowed with thermostability and a high level of recombinant expression in its soluble form while keeping the active site unmodified. This engineered protein was obtained by mutating the surface-exposed residues Phe107 and Arg91 to Glu. Molecular dynamics simulations showed that the mutated residues remain solvent-exposed, not affecting the overall protein structure. Thus, the ASV was used in a combinatorial mutagenesis of the distal heme pocket residues in which one, two, or three of the conserved polar residues [TyrB10(54), TyrCD1(67), and TrpG8(119)] were substituted with Phe. Mutants were characterized by infrared and resonance Raman spectroscopy and compared with the wild-type protein. Similar Fe-proximal His stretching frequencies suggest that none of the mutations alters the proximal side of the heme cavity. Two conformers were observed in the spectra of the CO complexes of both wild-type and ASV protein: form 1 with ν(FeC) and ν(CO) at 509 and 1938 cm(-1) and form 2 with ν(FeC) and ν(CO) at 518 and 1920 cm(-1), respectively. Molecular dynamics simulations were performed for the wild-type and ASV forms, as well as for the TyrB10 mutant. The spectroscopic and computational results demonstrate that CO interacts with TrpG8 in form 1 and interacts with both TrpG8 and TyrCD1 in form 2. TyrB10 does not directly interact with the bound CO.