Ligand migration and cavities within Scapharca Dimeric HbI: studies by time-resolved crystallo-graphy, Xe binding, and computational analysis.

As in many other hemoglobins, no direct route for migration of ligands between solvent and active site is evident from crystal structures of Scapharca inaequivalvis dimeric HbI. Xenon (Xe) and organic halide binding experiments, along with computational analysis presented here, reveal protein cavities as potential ligand migration routes. Time-resolved crystallographic experiments show that photodissociated carbon monoxide (CO) docks within 5 ns at the distal pocket B site and at more remote Xe4 and Xe2 cavities. CO rebinding is not affected by the presence of dichloroethane within the major Xe4 protein cavity, demonstrating that this cavity is not on the major exit pathway. The crystal lattice has a substantial influence on ligand migration, suggesting that significant conformational rearrangements may be required for ligand exit. Taken together, these results are consistent with a distal histidine gate as one important ligand entry and exit route, despite its participation in the dimeric interface.

Computer modeling in biotechnology: a partner in development.

Computational modeling can be a useful partner in biotechnology, in particular, in nanodevice engineering. Such modeling guides development through nanoscale views of biomolecules and devices not available through experimental imaging methods. We illustrate the role of computational modeling, mainly of molecular dynamics, through four case studies: development of silicon bionanodevices for single molecule electrical recording, development of carbon nano-tube-biomolecular systems as in vivo sensors, development of lipoprotein nanodiscs for assays of single membrane proteins, and engineering of oxygen tolerance into the enzyme hydrogenase for photosynthetic hydrogen gas production. The four case studies show how molecular dynamics approaches were adapted to the specific technical uses through (i) multi-scale extensions, (ii) fast quantum chemical force field evaluation, (iii) coarse graining, and (iv) novel sampling methods. The adapted molecular dynamics simulations provided key information on device behavior and revealed development opportunities, arguing that the "computational microscope" is an indispensable nanoengineering tool.

Finding gas migration pathways in proteins using implicit ligand sampling.

Implicit ligand sampling is a practical, efficient, and accurate method for finding the gas migration pathways for small hydrophobic gas molecules, such as oxygen, inside proteins. The method infers the gas migration pathways by calculating the potential of mean force for the gas molecule everywhere inside the protein by means of a molecular dynamics simulation of the protein in the absence of the gas molecule. Pathways can be constructed by connecting the areas of the protein that are favorable to the presence of gas. This method has the advantage of providing a comprehensive overview of all possible gas migration pathways and barriers in a given protein from a single simulation run. Implicit ligand sampling has been applied to a large number of hemoproteins. The example of the truncated hemoglobin from Paramecium caudatum is given to illustrate the method.

O2 migration pathways are not conserved across proteins of a similar fold.

Recent advances in computational biology have made it possible to map the complete network and energy profile of gas migration pathways inside proteins. Although networks of O(2) pathways have already been characterized for a small number of proteins, the general properties and locations of these pathways have not been previously compared between proteins. In this study, maps of the O(2) pathways inside 12 monomeric globins were computed. It is found that, despite the conserved tertiary structure fold of the studied globins, the shape and topology of O(2) pathway networks exhibit a large variability between different globins, except when two globins are nearly identical. The locations of the O(2) pathways are, however, found to be correlated with the location of large hydrophobic residues, and a similar correlation is observed in two unrelated protein families: monomeric globins and copper-containing amine oxidases. The results have implications for the evolution of gas pathways in proteins and for protein engineering applications involving modifications of these pathways.

Exploring molecular oxygen pathways in Hansenula polymorpha copper-containing amine oxidase.

The accessibility of large substrates to buried enzymatic active sites is dependent upon the utilization of proteinaceous channels. The necessity of these channels in the case of small substrates is questionable because diffusion through the protein matrix is often assumed. Copper amine oxidases contain a buried protein-derived quinone cofactor and a mononuclear copper center that catalyze the conversion of two substrates, primary amines and molecular oxygen, to aldehydes and hydrogen peroxide, respectively. The nature of molecular oxygen migration to the active site in the enzyme from Hansenula polymorpha is explored using a combination of kinetic, x-ray crystallographic, and computational approaches. A crystal structure of H. polymorpha amine oxidase in complex with xenon gas, which serves as an experimental probe for molecular oxygen binding sites, reveals buried regions of the enzyme suitable for transient molecular oxygen occupation. Calculated O(2) free energy maps using copper amine oxidase crystal structures in the absence of xenon correspond well with later experimentally observed xenon sites in these systems, and allow the visualization of O(2) migration routes of differing probabilities within the protein matrix. Site-directed mutagenesis designed to block individual routes has little effect on overall k(cat)/K(m) (O(2)), supporting multiple dynamic pathways for molecular oxygen to reach the active site.

Exploring gas permeability of cellular membranes and membrane channels with molecular dynamics.

Aquaporins are a family of membrane proteins specialized in rapid water conduction across biological membranes. Whether these channels also conduct gas molecules and the physiological significance of this potential function have not been well understood. Here we report 140 ns of molecular dynamics simulations of membrane-embedded AQP1 and of a pure POPE bilayer addressing these questions. The permeability of AQP1 to two types of gas molecules, O2 and CO2, was investigated using two complementary methods, namely, explicit gas diffusion simulation and implicit ligand sampling. The simulations show that the central (tetrameric) pore of AQP1 can be readily used by either gas molecule to permeate the channel. The two approaches produced similar free energy profiles associated with gas permeation through the central pore: a -0.4 to -1.7 kcal/mol energy well in the middle, and a 3.6-4.6 kcal/mol energy barrier in the periplasmic vestibule. The barrier appears to be mainly due to a dense cluster of water molecules anchored in the periplasmic mouth of the central pore by four aspartate residues. Water pores show a very low permeability to O2, but may contribute to the overall permeation of CO2 due to its more hydrophilic nature. Although the central pore of AQP1 is found to be gas permeable, the pure POPE bilayer provides a much larger cross-sectional area, thus exhibiting a much lower free energy barrier for CO2 and O2 permeation. As such, gas conduction through AQP1 may only be physiologically relevant either in membranes of low gas permeability, or in cells where a major fraction of the cellular membrane is occupied by AQPs.

Imaging the migration pathways for O2, CO, NO, and Xe inside myoglobin.

Myoglobin (Mb) is perhaps the most studied protein, experimentally and theoretically. Despite the wealth of known details regarding the gas migration processes inside Mb, there exists no fully conclusive picture of these pathways. We address this deficiency by presenting a complete map of all the gas migration pathways inside Mb for small gas ligands (O2, NO, CO, and Xe). To accomplish this, we introduce a computational approach for studying gas migration, which we call implicit ligand sampling. Rather than simulating actual gas migration events, we infer the location of gas migration pathways based on a free-energy perturbation approach applied to simulations of Mb's dynamical fluctuations at equilibrium in the absence of ligand. The method provides complete three-dimensional maps of the potential of mean force of gas ligand placement anywhere inside a protein-solvent system. From such free-energy maps we identify each gas docking site, the pathways between these sites, to the heme and to the external solution. Our maps match previously known features of these pathways in Mb, but also point to the existence of additional exits from the protein matrix in regions that are not easily probed by experiment. We also compare the pathway maps of Mb for different gas ligands and for different animal species.

Finding gas diffusion pathways in proteins: application to O2 and H2 transport in CpI [FeFe]-hydrogenase and the role of packing defects.

We report on a computational investigation of the passive transport of H2 and O2 between the external solution and the hydrogen-producing active site of CpI [FeFe]-hydrogenase from Clostridium pasteurianum. Two distinct methodologies for studying gas access are discussed and applied: (1) temperature-controlled locally enhanced sampling, and (2) volumetric solvent accessibility maps, providing consistent results. Both methodologies confirm the existence and function of a previously hypothesized pathway and reveal a second major pathway that had not been detected by previous analyses of CpI's static crystal structure. Our results suggest that small hydrophobic molecules, such as H2 and O2, diffusing inside CpI, take advantage of well-defined preexisting packing defects that are not always apparent from the protein's static structure, but that can be predicted from the protein's dynamical motion. Finally, we describe two contrasting modes of intraprotein transport for H2 and O2, which in our model are differentiated only by their size.

Mechanism of anionic conduction across ClC.

ClC chloride channels are voltage-gated transmembrane proteins that have been associated with a wide range of regulatory roles in vertebrates. To accomplish their function, they allow small inorganic anions to efficiently pass through, while blocking the passage of all other particles. Understanding the conduction mechanism of ClC has been the subject of many experimental investigations, but until now, the detailed dynamic mechanism was not known despite the availability of crystallographic structures. We investigate Cl(-) conduction by means of an all-atom molecular dynamics simulation of the ClC channel in a membrane environment. Based on our simulation results, we propose a king-of-the-hill mechanism for permeation, in which a lone ion bound to the center of the ClC pore is pushed out by a second ion that enters the pore and takes its place. Although the energy required to extract the single central ion from the pore is enormous, by resorting to this two-ion process, the largest free energy barrier for conduction is reduced to 4 kcal/mol. At the narrowest part of the pore, residues Tyr-445 and Ser-107 stabilize the central ion. There, the bound ion blocks the pore, disrupting the formation of a continuous water file that could leak protons, possibly preventing the passage of uncharged solutes.