An atomistic view on human hemoglobin carbon monoxide migration processes

Dynamics of camel and human hemoglobin revealed by molecular simulations

Hemoglobin is one of the most widely studied proteins genetically, biochemically, and structurally. It is an oxygen carrying tetrameric protein that imparts the characteristic red color to blood. Each chain of hemoglobin harbors a heme group embedded in a hydrophobic pocket. Several studies have investigated structural variations present in mammalian hemoglobin and their functional implications. However, camel hemoglobin has not been thoroughly explored, especially from a structural perspective. Importantly, very little is known about how the heme group interacts with hemoglobin under varying conditions of osmolarity and temperature. Several experimental studies have indicated that the tense (T) state is more stable than the relaxed (R) state of hemoglobin under normal physiological conditions. Despite the fact that R state is less stable than the T state, no extensive structural dynamics studies have been performed to investigate global quaternary transitions of R state hemoglobin under normal physiological conditions. To evaluate this, several 500 ns all-atom molecular dynamics simulations were performed to get a deeper understanding of how camel hemoglobin behaves under stress, which it is normally exposed to, when compared to human hemoglobin. Notably, camel hemoglobin was more stable under physiological stress when compared to human hemoglobin. Additionally, when compared to camel hemoglobin, cofactor-binding regions of hemoglobin also exhibited more fluctuations in human hemoglobin under the conditions studied. Several differences were observed between the residues of camel and human hemoglobin that interacted with heme. Importantly, distal residues His58 of α hemoglobin and His63 of β hemoglobin formed more sustained interactions, especially at higher temperatures, in camel hemoglobin. These residues are important for oxygen binding to hemoglobin. Thus, this work provides insights into how camel and human hemoglobin differ in their interactions under stress.

Ordered Motions in the Nitric-Oxide Dioxygenase Mechanism of Flavohemoglobin and Assorted Globins with Tightly Coupled Reductases

Nitric-oxide dioxygenases (NODs) activate and combine O₂ with NO to form nitrate. A variety of oxygen-binding hemoglobins with associated partner reductases or electron donors function as enzymatic NODs. Kinetic and structural investigations of the archetypal two-domain microbial flavohemoglobin-NOD have illuminated an allosteric mechanism that employs selective tunnels for O₂ and NO, gates for NO and nitrate, transient O₂ association with ferric heme, and an O₂ and NO-triggered, ferric heme spin crossover-driven, motion-controlled, and dipole-regulated electron-transfer switch. The proposed mechanism facilitates radical-radical coupling of ferric-superoxide with NO to form nitrate while preventing suicidal ferrous-NO formation. Diverse globins display the structural and functional motifs necessary for a similar allosteric NOD mechanism. In silico docking simulations reveal monomeric erythrocyte hemoglobin alpha-chain and beta-chain intrinsically matched and tightly coupled with NADH-cytochrome b₅ oxidoreductase and NADPH-cytochrome P450 oxidoreductase, respectively, forming membrane-bound flavohemoglobin-like mammalian NODs. The neuroprotective neuroglobin manifests a potential NOD role in a close-fitting ternary complex with membrane-bound NADH-cytochrome b₅ oxidoreductase and cytochrome b₅. Cytoglobin interfaces weakly with cytochrome b₅ for O₂ and NO-regulated electron-transfer and coupled NOD activity. The mechanistic model also provides insight into the evolution of O₂ binding cooperativity in hemoglobin and a basis for the discovery of allosteric NOD inhibitors.

Virtual screening of potential anticancer drugs based on microbial products

The development of microbial products for cancer treatment has been in the spotlight in recent years. In order to accelerate the lengthy and expensive drug development process, in silico screening tools are systematically employed, especially during the initial discovery phase. Moreover, considering the steadily increasing number of molecules approved by authorities for commercial use, there is a demand for faster methods to repurpose such drugs. Here we present a review on virtual screening web tools, such as publicly available databases of molecular targets and libraries of ligands, with the aim to facilitate the discovery of potential anticancer drugs based on microbial products. We provide an entry-level step-by-step description of the workflow for virtual screening of microbial metabolites with known protein targets, as well as two practical examples using freely available web tools. The first case presents a virtual screening study of drugs developed from microbial products using Caver Web, a web tool that performs docking along a tunnel. The second case comprises a comparative analysis between a wild type isocitrate dehydrogenase 1 and a mutant that results in cancer, using the recently developed web tool PredictSNP^Onco. In summary, this review provides the basic and essential background information necessary for virtual screening experiments, which may accelerate the discovery of novel anticancer drugs.

Screening of world approved drugs against highly dynamical spike glycoprotein of SARS-CoV-2 using CaverDock and machine learning

The new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes pathological pulmonary symptoms. Most efforts to develop vaccines and drugs against this virus target the spike glycoprotein, particularly its S1 subunit, which is recognised by angiotensin-converting enzyme 2. Here we use the in-house developed tool CaverDock to perform virtual screening against spike glycoprotein using a cryogenic electron microscopy structure (PDB-ID: 6VXX) and the representative structures of five most populated clusters from a previously published molecular dynamics simulation. The dataset of ligands was obtained from the ZINC database and consists of drugs approved for clinical use worldwide. Trajectories for the passage of individual drugs through the tunnel of the spike glycoprotein homotrimer, their binding energies within the tunnel, and the duration of their contacts with the trimer's three subunits were computed for the full dataset. Multivariate statistical methods were then used to establish structure-activity relationships and select top candidate for movement inhibition. This new protocol for the rapid screening of globally approved drugs (4359 ligands) in a multi-state protein structure (6 states) showed high robustness in the rate of finished calculations. The protocol is universal and can be applied to any target protein with an experimental tertiary structure containing protein tunnels or channels. The protocol will be implemented in the next version of CaverWeb (https://loschmidt.chemi.muni.cz/caverweb/) to make it accessible to the wider scientific community.

Towards understanding non-equivalence of α and β subunits within human hemoglobin in conformational relaxation and molecular oxygen rebinding

Picosecond to millisecond laser time-resolved transient absorption spectroscopy was used to study molecular oxygen (O₂) rebinding and conformational relaxation following O₂ photodissociation in the α and β subunits within human hemoglobin in the quaternary R-like structure. Oxy-cyanomet valency hybrids, α₂(Fe²⁺-O₂)β₂(Fe³⁺-CN) and α₂(Fe³⁺-CN)β₂(Fe²⁺-O₂), were used as models for oxygenated R-state hemoglobin. An extended kinetic model for geminate O₂ rebinding in the ferrous hemoglobin subunits, ligand migration between the primary and secondary docking site(s), and nonexponential tertiary relaxation within the R quaternary structure, was introduced and discussed. Significant functional non-equivalence of the α and β subunits in both the geminate O₂ rebinding and concomitant structural relaxation was revealed. For the β subunits, the rate constant for the geminate O₂ rebinding to the unrelaxed tertiary structure and the tertiary transition rate were found to be greater than the corresponding values for the α subunits. The conformational relaxation following the O₂ photodissociation in the α and β subunits was found to decrease the rate constant for the geminate O₂ rebinding, this effect being more than one order of magnitude greater for the β subunits than for the α subunits. Evidence was provided for the modulation of the O₂ rebinding to the individual α and β subunits within human hemoglobin in the R-state structure by the intrinsic heme reactivity through a change in proximal constraints upon the relaxation of the tertiary structure on a picosecond to microsecond time scale. Our results demonstrate that, for native R-state oxyhemoglobin, O₂ rebinding properties and spectral changes following the O₂ photodissociation can be adequately described as the sum of those for the α and β subunits within the valency hybrids. The isolated β chains (hemoglobin H) show similar behavior to the β subunits within the valency hybrids and can be used as a model for the β subunits within the R-state oxyhemoglobin. At the same time, the isolated α chains behave differently to the α subunits within the valency hybrids.

aquaPELE: A Monte Carlo-Based Algorithm to Sample the Effects of Buried Water Molecules in Proteins

Water is frequently found inside proteins, carrying out important roles in catalytic reactions or molecular recognition tasks. Therefore, computational models that aim to study protein-ligand interactions usually have to include water effects through explicit or implicit approaches to obtain reliable results. While full explicit models might be too computationally daunting for some applications, implicit models are normally faster but omit some of the most important contributions of water. This is the case of our in-house software, called protein energy landscape exploration (PELE), which uses implicit models to speed up conformational explorations as much as possible; the lack of explicit water sampling, however, limits its model. In this work, we confront this problem with the development of aquaPELE. It is a new algorithm that extends the exploration capabilities while keeping efficiency as it employs a mixed implicit/explicit approach to also take into account the effects of buried water molecules. With an additional Monte Carlo (MC) routine, a set of explicit water molecules is perturbed inside protein cavities and their effects are dynamically adjusted to the current state of the system. As a result, this implementation can be used to predict the principal hydration sites or the rearrangement and displacement of conserved water molecules upon the binding of a ligand. We benchmarked this new tool focusing on estimating ligand binding modes and hydration sites in cavities with important interfacial water molecules, according to crystallographic structures. Results suggest that aquaPELE sets a fast and reliable alternative for molecular recognition studies in systems with a strong water-dependency.

Direct observation of ligand migration within human hemoglobin at work

Human hemoglobin is the textbook example of the stereochemistry of an allosteric protein and of the exquisite control that a protein can exert over ligand binding. However, the fundamental basis by which the protein facilitates the ligand movement remains unknown. In this study, we used cryogenic X-ray crystallography and a high-repetition pulsed laser irradiation technique to elucidate the atomic details of ligand migration processes in hemoglobin after photolysis of the bound CO. Our data clarify the distinct CO migration pathways in the individual subunits of hemoglobin and unravel the functional roles of the internal cavities and neighboring amino acid residues in ligand exit and entry. Our results also demonstrate the high gas permeability and porosity of hemoglobin, facilitating O₂ delivery.

Hemoglobin is one of the best-characterized proteins with respect to structure and function, but the internal ligand diffusion pathways remain obscure and controversial. Here we captured the CO migration processes in the tense (T), relaxed (R), and second relaxed (R2) quaternary structures of human hemoglobin by crystallography using a high-repetition pulsed laser technique at cryogenic temperatures. We found that in each quaternary structure, the photodissociated CO molecules migrate along distinct pathways in the α and β subunits by hopping between the internal cavities with correlated side chain motions of large nonpolar residues, such as α14Trp(A12), α105Leu(G12), β15Trp(A12), and β71Phe(E15). We also observe electron density evidence for the distal histidine [α58/β63His(E7)] swing-out motion regardless of the quaternary structure, although less evident in α subunits than in β subunits, suggesting that some CO molecules have escaped directly through the E7 gate. Remarkably, in T-state Fe(II)-Ni(II) hybrid hemoglobins in which either the α or β subunits contain Ni(II) heme that cannot bind CO, the photodissociated CO molecules not only dock at the cavities in the original Fe(II) subunit, but also escape from the protein matrix and enter the cavities in the adjacent Ni(II) subunit even at 95 K, demonstrating the high gas permeability and porosity of the hemoglobin molecule. Our results provide a comprehensive picture of ligand movements in hemoglobin and highlight the relevance of cavities, nonpolar residues, and distal histidines in facilitating the ligand migration.

Fast Screening of Inhibitor Binding/Unbinding Using Novel Software Tool CaverDock

Protein tunnels and channels are attractive targets for drug design. Drug molecules that block the access of substrates or release of products can be efficient modulators of biological activity. Here, we demonstrate the applicability of a newly developed software tool CaverDock for screening databases of drugs against pharmacologically relevant targets. First, we evaluated the effect of rigid and flexible side chains on sets of substrates and inhibitors of seven different proteins. In order to assess the accuracy of our software, we compared the results obtained from CaverDock calculation with experimental data previously collected with heat shock protein 90α. Finally, we tested the virtual screening capabilities of CaverDock with a set of oncological and anti-inflammatory FDA-approved drugs with two molecular targets—cytochrome P450 17A1 and leukotriene A4 hydrolase/aminopeptidase. Calculation of rigid trajectories using four processors took on average 53 min per molecule with 90% successfully calculated cases. The screening identified functional tunnels based on the profile of potential energies of binding and unbinding trajectories. We concluded that CaverDock is a sufficiently fast, robust, and accurate tool for screening binding/unbinding processes of pharmacologically important targets with buried functional sites. The standalone version of CaverDock is available freely at https://loschmidt.chemi.muni.cz/caverdock/ and the web version at https://loschmidt.chemi.muni.cz/caverweb/.

GPathFinder: Identification of Ligand-Binding Pathways by a Multi-Objective Genetic Algorithm

Protein-ligand docking is a widely used method to generate solutions for the binding of a small molecule with its target in a short amount of time. However, these methods provide identification of physically sound protein-ligand complexes without a complete view of the binding process dynamics, which has been recognized to be a major discriminant in binding affinity and ligand selectivity. In this paper, a novel piece of open-source software to approach this problem is presented, called GPathFinder. It is built as an extension of the modular GaudiMM platform and is able to simulate ligand diffusion pathways at atomistic level. The method has been benchmarked on a set of 20 systems whose ligand-binding routes were studied by other computational tools or suggested from experimental "snapshots". In all of this set, GPathFinder identifies those channels that were already reported in the literature. Interestingly, the low-energy pathways in some cases indicate novel possible binding routes. To show the usefulness of GPathFinder, the analysis of three case systems is reported. We believe that GPathFinder is a software solution with a good balance between accuracy and computational cost, and represents a step forward in extending protein-ligand docking capacities, with implications in several fields such as drug or enzyme design.

Dynamics based clustering of globin family members

A methodology to cluster proteins based on their dynamics’ similarity is presented. For each pair of proteins from a dataset, the structures are superimposed, and the Anisotropic Network Model modes of motions are calculated. The twelve slowest modes from each protein are matched using a local mode alignment algorithm based on the local sequence alignment algorithm of Smith–Waterman. The dynamical similarity distance matrix is calculated based on the top scoring matches of each pair and the proteins are clustered using a hierarchical clustering algorithm. The utility of this method is exemplified on a dataset of protein chains from the globin family and a dataset of tetrameric hemoglobins. The results demonstrate the effect of the quaternary structure of globin members on their intrinsic dynamics and show good ability to distinguish between different states of hemoglobin, revealing the dynamical relations between them.

Dynamical comparison between myoglobin and hemoglobin

Myoglobin and hemoglobin are globular hemeproteins, when the former is a monomer and the latter a heterotetramer. Despite the structural similarity of myoglobin to α and β subunits of hemoglobin, there is a functional difference between the two proteins, owing to the quaternary structure of hemoglobin. The effect of the quaternary structure of hemoglobin on the intrinsic dynamics of its subunits is explored by dynamical comparison of the two proteins. Anisotropic Network Model modes of motion were calculated for hemoglobin and myoglobin. Dynamical comparison between the proteins was performed using global and local Anisotropic Network Model mode alignment algorithms based on the algorithms of Smith-Waterman and Needleman-Wunsch for sequence comparison. The results indicate that the quaternary structure of hemoglobin substantially alters the intrinsic dynamics of its subunits, an effect that may contribute to the functional difference between the two proteins. Local dynamics similarity between the proteins is still observed at the major exit route of the ligand.

The ins and outs of vanillyl alcohol oxidase: Identification of ligand migration paths

Vanillyl alcohol oxidase (VAO) is a homo-octameric flavoenzyme belonging to the VAO/PCMH family. Each VAO subunit consists of two domains, the FAD-binding and the cap domain. VAO catalyses, among other reactions, the two-step conversion of p-creosol (2-methoxy-4-methylphenol) to vanillin (4-hydroxy-3-methoxybenzaldehyde). To elucidate how different ligands enter and exit the secluded active site, Monte Carlo based simulations have been performed. One entry/exit path via the subunit interface and two additional exit paths have been identified for phenolic ligands, all leading to the si side of FAD. We argue that the entry/exit path is the most probable route for these ligands. A fourth path leading to the re side of FAD has been found for the co-ligands dioxygen and hydrogen peroxide. Based on binding energies and on the behaviour of ligands in these four paths, we propose a sequence of events for ligand and co-ligand migration during catalysis. We have also identified two residues, His466 and Tyr503, which could act as concierges of the active site for phenolic ligands, as well as two other residues, Tyr51 and Tyr408, which could act as a gateway to the re side of FAD for dioxygen. Most of the residues in the four paths are also present in VAO’s closest relatives, eugenol oxidase and p-cresol methylhydroxylase. Key path residues show movements in our simulations that correspond well to conformations observed in crystal structures of these enzymes. Preservation of other path residues can be linked to the electron acceptor specificity and oligomerisation state of the three enzymes. This study is the first comprehensive overview of ligand and co-ligand migration in a member of the VAO/PCMH family, and provides a proof of concept for the use of an unbiased method to sample this process.

Enzymes are bionanomachines, which speed up chemical reactions in organisms. To understand how they achieve that, we need to study their mechanisms. Computational enzymology can show us what happens in the enzyme’s active site during a reaction. But molecules need first to reach the active site before a reaction can start. The process of substrate entry and product exit to the active site is often neglected when studying enzymes. However, these two events are of fundamental importance to the proper functioning of any enzyme. We are interested in these dynamic processes to complete our understanding of the mode of action of enzymes. In our work, we have studied substrate and product migration in vanillyl alcohol oxidase. This enzyme can produce the flavour vanillin and enantiopure alcohols, but also catalyses other reactions. The named products are of interest to the flavour- and fine-chemical industries.

Tertiary and quaternary structural basis of oxygen affinity in human hemoglobin as revealed by multiscale simulations

Human hemoglobin (Hb) is a benchmark protein of structural biology that shaped our view of allosterism over 60 years ago, with the introduction of the MWC model based on Perutz structures of the oxy(R) and deoxy(T) states and the more recent Tertiary Two-State model that proposed the existence of individual subunit states -“r” and “t”-, whose structure is yet unknown. Cooperative oxygen binding is essential for Hb function, and despite decades of research there are still open questions related to how tertiary and quaternary changes regulate oxygen affinity. In the present work, we have determined the free energy profiles of oxygen migration and for HisE7 gate opening, with QM/MM calculations of the oxygen binding energy in order to address the influence of tertiary differences in the control of oxygen affinity. Our results show that in the α subunit the low to high affinity transition is achieved by a proximal effect that mostly affects oxygen dissociation and is the driving force of the allosteric transition, while in the β subunit the affinity change results from a complex interplay of proximal and distal effects, including an increase in the HE7 gate opening, that as shown by free energy profiles promotes oxygen uptake.

Simulating Substrate Recognition and Oxidation in Laccases: From Description to Design

To meet the very specific requirements demanded by industry, proteins must be appropriately tailored. Engineering laccases, to improve the oxidation of small molecules, with applications in multiple fields, is, however, a difficult task. Most efforts have concentrated on increasing the redox potential of the enzyme, but in recent work, we have pursued an alternate strategy to engineering these biocatalysts. In particular, we have found that redesigning substrate binding at the T1 pocket, guided by in silico methodologies, to be a more consistent option. In this work, we evaluate the robustness of our computational approach to estimate activity, emphasizing the importance of the binding event in laccase reactivity. Strengths and weaknesses of the protocol are discussed along with its potential for scoring large numbers of protein sequences and thus its significance in protein engineering.

Structural and Mechanistic Insights into Hemoglobin-catalyzed Hydrogen Sulfide Oxidation and the Fate of Polysulfide Products

Hydrogen sulfide is a cardioprotective signaling molecule but is toxic at elevated concentrations. Red blood cells can synthesize H₂S but, lacking organelles, cannot dispose of H₂S via the mitochondrial sulfide oxidation pathway. We have recently shown that at high sulfide concentrations, ferric hemoglobin oxidizes H₂S to a mixture of thiosulfate and iron-bound polysulfides in which the latter species predominates. Here, we report the crystal structure of human hemoglobin containing low spin ferric sulfide, the first intermediate in heme-catalyzed sulfide oxidation. The structure provides molecular insights into why sulfide is susceptible to oxidation in human hemoglobin but is stabilized against it in HbI, a specialized sulfide-carrying hemoglobin from a mollusk adapted to life in a sulfide-rich environment. We have also captured a second sulfide bound at a postulated ligand entry/exit site in the α-subunit of hemoglobin, which, to the best of our knowledge, represents the first direct evidence for this site being used to access the heme iron. Hydrodisulfide, a postulated intermediate at the junction between thiosulfate and polysulfide formation, coordinates ferric hemoglobin and, in the presence of air, generated thiosulfate. At low sulfide/heme iron ratios, the product distribution between thiosulfate and iron-bound polysulfides was approximately equal. The iron-bound polysulfides were unstable at physiological glutathione concentrations and were reduced with concomitant formation of glutathione persulfide, glutathione disulfide, and H₂S. Hence, although polysulfides are unlikely to be stable in the reducing intracellular milieu, glutathione persulfide could serve as a persulfide donor for protein persulfidation, a posttranslational modification by which H₂S is postulated to signal.

Structural Plasticity in Globins: Role of Protein Dynamics in Defining Ligand Migration Pathways

Globins are a family of proteins characterized by the presence of the heme prosthetic group and involved in variety of biological functions in the cell. Due to their biological relevance and widespread distribution in all kingdoms of life, intense research efforts have been devoted to disclosing the relationships between structural features, protein dynamics, and function. Particular attention has been paid to the impact of differences in amino acid sequence on the topological features of docking sites and cavities and to the influence of conformational flexibility in facilitating the migration of small ligands through these cavities. Often, tunnels are carved in the interior of globins, and ligand exchange is regulated by gating residues. Understanding the subtle intricacies that relate the differences in sequence with the structural and dynamical features of globins with the ultimate aim of rationalizing the thermodynamics and kinetics of ligand binding continues to be a major challenge in the field. Due to the evolution of computational techniques, significant advances into our understanding of these questions have been made. In this review we focus our attention on the analysis of the ligand migration pathways as well as the function of the structural cavities and tunnels in a series of representative globins, emphasizing the synergy between experimental and theoretical approaches to gain a comprehensive knowledge into the molecular mechanisms of this diverse family of proteins.

Benchmarking Rapid TLES Simulations of Gas Diffusion in Proteins: Mapping O2 Migration and Escape in Myoglobin as a Case Study

Standard molecular dynamics (MD) simulations of gas diffusion consume considerable computational time and resources even for small proteins. To combat this, temperature-controlled locally enhanced sampling (TLES) examines multiple diffusion trajectories per simulation by accommodating multiple noninteracting copies of a gas molecule that diffuse independently, while the protein and water molecules experience an average interaction from all copies. Furthermore, gas migration within a protein matrix can be accelerated without altering protein dynamics by increasing the effective temperature of the TLES copies. These features of TLES enable rapid simulations of gas diffusion within a protein matrix at significantly reduced (∼98%) computational cost. However, the results of TLES and standard MD simulations have not been systematically compared, which limits the adoption of the TLES approach. We address this drawback here by benchmarking TLES against standard MD in the simulation of O2 diffusion in myoglobin (Mb) as a case study since this model system has been extensively characterized. We find that 2 ns TLES and 108 ns standard simulations map the same network of diffusion tunnels in Mb and uncover the same docking sites, barriers, and escape portals. We further discuss the influence of simulation time as well as the number of independent simulations on the O2 population density within the diffusion tunnels and on the sampling of Mb's conformational space as revealed by principal component analysis. Overall, our comprehensive benchmarking reveals that TLES is an appropriate and robust tool for the rapid mapping of gas diffusion in proteins when the kinetic data provided by standard MD are not required. Furthermore, TLES provides explicit ligand diffusion pathways, unlike most rapid methods.

Molecular determinants for selective C25-hydroxylation of vitamins D2and D3by fungal peroxygenases

Regioselective hydroxylation of vitamin D byAgrocybe aegeritaandCoprinopsis cinereaperoxygenases was investigated in an experimental and computational study.

New Monte Carlo Based Technique To Study DNA-Ligand Interactions

We present a new all-atom Monte Carlo technique capable of performing quick and accurate DNA-ligand conformational sampling. In particular, and using the PELE software as a frame, we have introduced an additional force field, an implicit solvent, and an anisotropic network model to effectively map the DNA energy landscape. With these additions, we successfully generated DNA conformations for a test set composed of six DNA fragments of A-DNA and B-DNA. Moreover, trajectories generated for cisplatin and its hydrolysis products identified the best interacting compound and binding site, producing analogous results to microsecond molecular dynamics simulations. Furthermore, a combination of the Monte Carlo trajectories with Markov State Models produced noncovalent binding free energies in good agreement with the published molecular dynamics results, at a significantly lower computational cost. Overall our approach will allow a quick but accurate sampling of DNA-ligand interactions.

Clustering and percolation in protein loop structures

Background

High precision protein loop modelling remains a challenge, both in template based and template independent approaches to protein structure prediction.

Method

We introduce the concepts of protein loop clustering and percolation, to develop a quantitative approach to systematically classify the modular building blocks of loops in crystallographic folded proteins. These fragments are all different parameterisations of a unique kink solution to a generalised discrete nonlinear Schrödinger (DNLS) equation. Accordingly, the fragments are also local energy minima of the ensuing energy function.

Results

We show how the loop fragments cover practically all ultrahigh resolution crystallographic protein structures in Protein Data Bank (PDB), with a 0.2 Ångström root-mean-square (RMS) precision. We find that no more than 12 different loop fragments are needed, to describe around 38 % of ultrahigh resolution loops in PDB. But there is also a large number of loop fragments that are either unique, or very rare, and examples of unique fragments are found even in the structure of a myoglobin.

Conclusions

Protein loops are built in a modular fashion. The loops are composed of fragments that can be modelled by the kink of the DNLS equation. The majority of loop fragments are also common, which are shared by many proteins. These common fragments are probably important for supporting the overall protein conformation. But there are also several fragments that are either unique to a given protein, or very rare. Such fragments are probably related to the function of the protein. Furthermore, we have found that the amino acid sequence does not determine the structure in a unique fashion. There are many examples of loop fragments with an identical amino acid sequence, but with a very different structure.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-015-0049-x) contains supplementary material, which is available to authorized users.

Understanding the kinetics of ligand binding to globins with molecular dynamics simulations: the necessity of multiple state models

Residue conformational changes and internal cavity migration processes play a key role in regulating the kinetics of ligand migration and binding events in globins. Molecular dynamics simulations have demonstrated their value in the study of these processes in different haemoglobins, but derivation of kinetic data demands the use of more complex techniques like enhanced sampling molecular dynamics methods. This review discusses the different methodologies that are currently applied to study the ligand migration process in globins and highlight those specially developed to derive kinetic data.

Steroid hydroxylation by basidiomycete peroxygenases: a combined experimental and computational study

The goal of this study is the selective oxyfunctionalization of steroids under mild and environmentally friendly conditions using fungal enzymes. With this purpose, peroxygenases from three basidiomycete species were tested for the hydroxylation of a variety of steroidal compounds, using H2O2 as the only cosubstrate. Two of them are wild-type enzymes from Agrocybe aegerita and Marasmius rotula, and the third one is a recombinant enzyme from Coprinopsis cinerea. The enzymatic reactions on free and esterified sterols, steroid hydrocarbons, and ketones were monitored by gas chromatography, and the products were identified by mass spectrometry. Hydroxylation at the side chain over the steroidal rings was preferred, with the 25-hydroxyderivatives predominating. Interestingly, antiviral and other biological activities of 25-hydroxycholesterol have been reported recently (M. Blanc et al., Immunity 38:106-118, 2013, http://dx.doi.org/10.1016/j.immuni.2012.11.004). However, hydroxylation in the ring moiety and terminal hydroxylation at the side chain also was observed in some steroids, the former favored by the absence of oxygenated groups at C-3 and by the presence of conjugated double bonds in the rings. To understand the yield and selectivity differences between the different steroids, a computational study was performed using Protein Energy Landscape Exploration (PELE) software for dynamic ligand diffusion. These simulations showed that the active-site geometry and hydrophobicity favors the entrance of the steroid side chain, while the entrance of the ring is energetically penalized. Also, a direct correlation between the conversion rate and the side chain entrance ratio could be established that explains the various reaction yields observed.

Visualizing phosphodiester-bond hydrolysis by an endonuclease

The enzymatic hydrolysis of DNA phosphodiester bonds has been widely studied, but the chemical reaction has not yet been observed. Here we follow the generation of a DNA double-strand break (DSB) by the Desulfurococcus mobilis homing endonuclease I-DmoI, trapping sequential stages of a two-metal-ion cleavage mechanism. We captured intermediates of the different catalytic steps, and this allowed us to watch the reaction by 'freezing' multiple states. We observed the successive entry of two metals involved in the reaction and the arrival of a third cation in a central position of the active site. This third metal ion has a crucial role, triggering the consecutive hydrolysis of the targeted phosphodiester bonds in the DNA strands and leaving its position once the DSB is generated. The multiple structures show the orchestrated conformational changes in the protein residues, nucleotides and metals during catalysis.

Differential control of heme reactivity in alpha and beta subunits of hemoglobin: a combined Raman spectroscopic and computational study

The use of hybrid hemoglobin (Hb), with mesoheme substituted for protoheme, allows separate monitoring of the α or β hemes along the allosteric pathway. Using resonance Raman (rR) spectroscopy in silica gel, which greatly slows protein motions, we have observed that the Fe-histidine stretching frequency, νFeHis, which is a monitor of heme reactivity, evolves between frequencies characteristic of the R and T states, for both α or β chains, prior to the quaternary R-T and T-R shifts. Computation of νFeHis, using QM/MM and the conformational search program PELE, produced remarkable agreement with experiment. Analysis of the PELE structures showed that the νFeHis shifts resulted from heme distortion and, in the α chain, Fe-His bond tilting. These results support the tertiary two-state model of ligand binding (Henry et al., Biophys. Chem. 2002, 98, 149). Experimentally, the νFeHis evolution is faster for β than for α chains, and pump-probe rR spectroscopy in solution reveals an inflection in the νFeHis time course at 3 μs for β but not for α hemes, an interval previously shown to be the first step in the R-T transition. In the α chain νFeHis dropped sharply at 20 μs, the final step in the R-T transition. The time courses are fully consistent with recent computational mapping of the R-T transition via conjugate peak refinement by Karplus and co-workers (Fischer et al., Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 5608). The effector molecule IHP was found to lower νFeHis selectively for α chains within the R state, and a binding site in the α1α2 cleft is suggested.

A model for the flexibility of the distal histidine in dehaloperoxidase-hemoglobin A based on X-ray crystal structures of the carbon monoxide adduct

Dehaloperoxidase hemoglobin A (DHP A) is a multifunctional hemoglobin that appears to have evolved oxidative pathways for the degradation of xenobiotics as a protective function that complements the oxygen transport function. DHP A possesses at least two internal binding sites, one for substrates and one for inhibitors, which include various halogenated phenols and indoles. Herein, we report the X-ray crystallographic structure of the carbonmonoxy complex (DHPCO). Unlike other DHP structures with 6-coordinated heme, the conformation of the distal histidine (H55) in DHPCO is primarily external or solvent exposed, despite the fact that the heme Fe is 6-coordinated. As observed generally in globins, DHP exhibits two distal histidine conformations (one internal and one external). In previous structural studies, we have shown that the distribution of H55 conformations is weighted strongly toward the external position when the DHP heme Fe is 5-coordinated. The large population of the external conformation of the distal histidine observed in DHPCO crystals at pH 6.0 indicates that some structural factor in DHP must account for the difference from other globins, which exhibit a significant external conformation only when pH < 4.5. While the original hypothesis suggested that interaction with a heme-Fe-bound ligand was the determinant of H55 conformation, the current study forces a refinement of that hypothesis. The external or open conformation of H55 is observed to have interactions with two propionate groups in heme, at distances of 3.82 and 2.73 Å, respectively. A relatively weak hydrogen bonding interaction between H55 and CO, combined with strong interactions with heme propionate (position 6), is hypothesized to strengthen the external conformation of H55. Density function theory (DFT) calculations were conducted to test whether there is a weaker hydrogen bond interaction between H55 and heme bonded CO or O₂. Molecular dynamics simulations were conducted to examine how the tautomeric forms of H55 affect the dynamic motions of the distal histidine that govern the switching between open and closed conformations. The calculations support the modified hypothesis suggesting a competition between the strength of interactions with heme ligand and the heme propionates as the factors that determine the conformation of the distal histidine.

Ligand uptake modulation by internal water molecules and hydrophobic cavities in hemoglobins

Internal water molecules play an active role in ligand uptake regulation, since displacement of retained water molecules from protein surfaces or cavities by incoming ligands can promote favorable or disfavorable effects over the global binding process. Detection of these water molecules by X-ray crystallography is difficult given their positional disorder and low occupancy. In this work, we employ a combination of molecular dynamics simulations and ligand rebinding over a broad time range to shed light into the role of water molecules in ligand migration and binding. Computational studies on the unliganded structure of the thermostable truncated hemoglobin from Thermobifida fusca (Tf-trHbO) show that a water molecule is in the vicinity of the iron heme, stabilized by WG8 with the assistance of YCD1, exerting a steric hindrance for binding of an exogenous ligand. Mutation of WG8 to F results in a significantly lower stabilization of this water molecule and in subtle dynamical structural changes that favor ligand binding, as observed experimentally. Water is absent from the fully hydrophobic distal cavity of the triple mutant YB10F-YCD1F-WG8F (3F), due to the lack of residues capable of stabilizing it nearby the heme. In agreement with these effects on the barriers for ligand rebinding, over 97% of the photodissociated ligands are rebound within a few nanoseconds in the 3F mutant case. Our results demonstrate the specific involvement of water molecules in shaping the energetic barriers for ligand migration and binding.

Monte Carlo Free Ligand Diffusion with Markov State Model Analysis and Absolute Binding Free Energy Calculations

Obtaining absolute binding free energies from unbiased ligand diffusion has attracted a significant amount of attention due to its implications in drug design. Several studies have used special purpose computers and software to achieve microsecond molecular dynamics which, combined with a Markov state model analysis, are capable of providing absolute binding free energies. We have recently developed a Monte Carlo based technique, PELE, capable of performing a dynamical exploration of the protein-ligand energy landscape including free ligand diffusion into the active site, at a fraction of the computational cost of molecular dynamics techniques. We demonstrate here the capabilities of our Monte Carlo technique in obtaining absolute binding free energies for a series of benzamidine like inhibitors into trypsin. Our results are in good agreement with experimental data and other molecular dynamics simulations, indicating that PELE can be a useful tool for quick estimates of binding free energies and mechanisms.

Real-time tracking of CO migration and binding in the α and β subunits of human hemoglobin via 150-ps time-resolved Laue crystallography

We have developed the method of picosecond Laue crystallography and used this capability to probe ligand dynamics in tetrameric R-state hemoglobin (Hb). Time-resolved, 2 Å-resolution electron density maps of photolyzed HbCO reveal the time-dependent population of CO in the binding (A) and primary docking (B) sites of both α and β subunits from 100 ps to 10 μs. The proximity of the B site in the β subunit is about 0.25 Å closer to its A binding site, and its k_BA rebinding rate (~300 μs^-1) is six times faster, suggesting distal control of the rebinding dynamics. Geminate rebinding in the β subunit exhibits both prompt and delayed geminate phases. We developed a microscopic model to quantitatively explain the observed kinetics, with three states for the α subunit and four states for the β subunit. This model provides a consistent framework for interpreting rebinding kinetics reported in prior studies of both HbCO and HbO₂.

Ligand migration through hemeprotein cavities: insights from laser flash photolysis and molecular dynamics simulations

The presence of cavities and tunnels in the interior of proteins, in conjunction with the structural plasticity arising from the coupling to the thermal fluctuations of the protein scaffold, has profound consequences on the pathways followed by ligands moving through the protein matrix. In this perspective we discuss how quantitative analysis of experimental rebinding kinetics from laser flash photolysis, trapping of unstable conformational states by embedding proteins within the nanopores of silica gels, and molecular simulations can synergistically converge to gain insight into the migration mechanism of ligands. We show how the evaluation of the free energy landscape for ligand diffusion based on the outcome of computational techniques can assist the definition of sound reaction schemes, leading to a comprehensive understanding of the broad range of chemical events and time scales that encompass the transport of small ligands in hemeproteins.

Comparative analysis of inner cavities and ligand migration in non-symbiotic AHb1 and AHb2

This study reports a comparative analysis of the topological properties of inner cavities and the intrinsic dynamics of non-symbiotic hemoglobins AHb1 and AHb2 from Arabidopsis thaliana. The two proteins belong to the 3/3 globin fold and have a sequence identity of about 60%. However, it is widely assumed that they have distinct physiological roles. In order to investigate the structure-function relationships in these proteins, we have examined the bis-histidyl and ligand-bound hexacoordinated states by atomistic simulations using in silico structural models. The results allow us to identify two main pathways to the distal cavity in the bis-histidyl hexacoordinated proteins. Nevertheless, a larger accessibility to small gaseous molecules is found in AHb2. This effect can be attributed to three factors: the mutation Leu35(AHb1)→Phe32(AHb2), the enhanced flexibility of helix B, and the more favorable energetic profile for ligand migration to the distal cavity. The net effect of these factors would be to facilitate the access of ligands, thus compensating the preference for the fully hexacoordination of AHb2, in contrast to the equilibrium between hexa- and pentacoordinated species in AHb1. On the other hand, binding of the exogenous ligand introduces distinct structural changes in the two proteins. A well-defined tunnel is formed in AHb1, which might be relevant to accomplish the proposed NO detoxification reaction. In contrast, no similar tunnel is found in AHb2, which can be ascribed to the reduced flexibility of helix E imposed by the larger number of salt bridges compared to AHb1. This feature would thus support the storage and transport functions proposed for AHb2. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

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.

Single vs. multiple ligand pathways in globins: a computational view

Diatomic ligand migration in globins has been the subject of numerous studies. Still, a consensus picture for the ligand entrance is not clear, with a growing concern among experimental researchers that computational simulations always show multiple pathways for any globin. Modeling non-biased ligand entrance from conventional molecular dynamics techniques, however, has shown to be difficult (and expensive). Here we use our Monte Carlo methodology, capable of freely mapping ligand diffusion and the description of rare events, to two well-studied systems: myoglobin and the mini-hemoglobin from the sea worm Cerebratulus lacteus. Our results clearly show that the simulations are specific to the system providing a different trend in the entrance pathway, as expected from experiments. While Mb presents multiple entrance pathways, populating the well-known xenon cavities, in CerHb the ligand enters the protein only by one apolar channel. Most of the trajectories (64%) visiting myoglobin's active site though, are gated by the distal histidine. Such detailed information, accessible through the state of the art algorithms in PELE, is computationally inexpensive and available to all non-profit researchers. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.