Conformational selection mechanism governs oxygen ligation to H-NOX proteins

H-NOX proteins are at present the only structural models available for the study of the ligand binding affinity and selectivity of soluble guanylate cyclase, the physiological receptor of nitric oxide. The oxy complex and resting state structures of two bacterial H-NOX proteins of markedly different oxygen affinity, but of quite similar sequence, were studied by molecular dynamics simulations at 300 K and 400 K. Unexpectedly, the different O₂ affinity was found to be reflected in differences of the resting states structures. A conformation containing a pre-formed oxygen-binding cage is the most populated in the resting state equilibrium ensemble of the only successful O₂ binder, Tt H-NOX at 300 K, suggesting that conformational selection governs the interaction.     

Ultrafast dynamics of ligands within heme proteins

Physiological bond formation and bond breaking events between proteins and ligands and their immediate consequences are difficult to synchronize and study in general. However, diatomic ligands can be photodissociated from heme, and thus in heme proteins ligand release and rebinding dynamics and trajectories have been studied on timescales of the internal vibrations of the protein that drive many biochemical reactions, and longer. The rapidly expanding number of characterized heme proteins involved in a large variety of functions allows comparative dynamics-structure-function studies. In this review, an overview is given of recent progress in this field, and in particular on initial sensing processes in signaling proteins, and on ligand and electron transfer dynamics in oxidases and cytochromes.     

A simple method for the determination of reduction potentials in heme proteins

We describe a simple method for the determination of heme protein reduction potentials. We use the method to determine the reduction potentials for the PAS-A domains of the regulatory heme proteins human NPAS2 (E_m = −115 mV ± 2 mV, pH 7.0) and human CLOCK (E_m = −111 mV ± 2 mV, pH 7.0). We suggest that the method can be easily and routinely applied to the determination of reduction potentials across the family of heme proteins.     

The FG loop of a heme-based gas sensor enzyme, Ec DOS, functions in heme binding, autoxidation and catalysis

Ec DOS is a heme-based gas sensor enzyme that catalyzes conversion from cyclic-di-GMP to linear-di-GMP in response to gas molecules, such as oxygen, CO and NO. Ec DOS contains an N-terminal heme-binding PAS domain and C-terminal phosphodiesterase domain. Based on crystal structures of the isolated heme-binding domain, it is suggested that the FG loop is involved in intra-molecular signal transduction to the catalytic domain. We generated nine full-length proteins mutated at ionic and non-ionic polar residues between positions 83 and 96 corresponding to the F-helix and FG loop, and examined the heme binding properties, autoxidation rates, and catalytic activities of mutant proteins. N84A and R85A mutant proteins displayed lower heme binding affinities, consistent with the finding that Asn84 interacts with propionate of protoporphyrin IX, and Arg85 with Asp40 on the heme proximal side. Autoxidation rates (0.058-0.54 min⁻¹) of R91A, S96A and K89A/R91A/E93A mutant proteins were significantly higher than that (0.0053 min⁻¹) of wild-type protein, suggesting that these residues in the FG loop form heme distal architecture conferring stability to the Fe(II)-O₂ complex. Catalytic activities of N84A and R85A mutant proteins with low heme affinity were significantly higher than those of wild-type protein in the absence of gas molecules. Accordingly, we propose that loss of heme binding enhances basal catalysis without the gas molecule, consistent with previous reports on heme inhibition of Ec DOS catalysis.     

NMR characterization of the conformational fluctuations of the human lymphocyte function‐associated antigen‐1 I‐domain

Lymphocyte function‐associated antigen‐1 (LFA‐1) is an integrin protein that transmits information across the plasma membrane through the so‐called inside‐out and outside‐in signaling mechanisms. To investigate these mechanisms, we carried out an NMR analysis of the dynamics of the LFA‐1 I‐domain, which has enabled us to characterize the motions of this domain on a broad range of timescales. We studied first the internal motions on the nanosecond timescale by spin relaxation measurements and model‐free analysis. We then extended this analysis to the millisecond timescale motions by measuring ¹⁵N‐¹H residual dipolar couplings of the backbone amide groups. We analyzed these results in the context of the three major conformational states of the I‐domain using their corresponding X‐ray crystallographic structures. Our results highlight the importance of the low‐frequency motions of the LFA‐1 I‐domain in the inactive apo‐state. We found in particular that α‐helix 7 is in a position in the apo‐closed state that cannot be fully described by any of the existing X‐ray structures, as it appears to be in dynamic exchange between different conformations. This type of motion seems to represent an inherent property of the LFA‐1 I‐domain and might be relevant for controlling the access to the allosteric binding pocket, as well as for the downward displacement of α‐helix 7 that is required for the activation of LFA‐1.     

Coarse-grained molecular dynamics simulations of membrane proteins and peptides

Molecular dynamics (MD) simulations provide a valuable approach to the dynamics, structure, and stability of membrane-protein systems. Coarse-grained (CG) models, in which small groups of atoms are treated as single particles, enable extended (>100 ns) timescales to be addressed. In this study, we explore how CG-MD methods that have been developed for detergents and lipids may be extended to membrane proteins. In particular, CG-MD simulations of a number of membrane peptides and proteins are used to characterize their interactions with lipid bilayers. CG-MD is used to simulate the insertion of synthetic model membrane peptides (WALPs and LS3) into a lipid (PC) bilayer. WALP peptides insert in a transmembrane orientation, whilst the LS3 peptide adopts an interfacial location, both in agreement with experimental biophysical data. This approach is extended to a transmembrane fragment of the Vpu protein from HIV-1, and to the coat protein from fd phage. Again, simulated protein/membrane interactions are in good agreement with solid state NMR data for these proteins. CG-MD has also been applied to an M3-M4 fragment from the CFTR protein. Simulations of CFTR M3-M4 in a detergent micelle reveal formation of an alpha-helical hairpin, consistent with a variety of biophysical data. In an I231D mutant, the M3-M4 hairpin is additionally stabilized via an inter-helix Q207/D231 interaction. Finally, CG-MD simulations are extended to a more complex membrane protein, the bacterial sugar transporter LacY. Comparison of a 200 ns CG-MD simulation of LacY in a DPPC bilayer with a 50 ns atomistic simulation of the same protein in a DMPC bilayer shows that the two methods yield comparable predictions of lipid-protein interactions. Taken together, these results demonstrate the utility of CG-MD simulations for studies of membrane/protein interactions.     

An interaction-based analysis of calcium-induced conformational changes in Ca2+ sensor proteins.

Calcium sensor proteins translate transient increases in intracellular calcium levels into metabolic or mechanical responses, by undergoing dramatic conformational changes upon Ca2+ binding. A detailed analysis of the calcium binding-induced conformational changes in the representative calcium sensors calmodulin (CaM) and troponin C was performed to obtain insights into the underlying molecular basis for their response to the binding of calcium. Distance difference matrices, analysis of interresidue contacts, comparisons of interhelical angles, and inspection of structures using molecular graphics were used to make unbiased comparisons of the various structures. The calcium-induced conformational changes in these proteins are dominated by reorganization of the packing of the four helices within each domain. Comparison of the closed and open conformations confirms that calcium binding causes opening within each of the EF-hands. A secondary analysis of the conformation of the C-terminal domain of CaM (CaM-C) clearly shows that CaM-C occupies a closed conformation in the absence of calcium that is distinct from the semi-open conformation observed in the C-terminal EF-hand domains of myosin light chains. These studies provide insight into the structural basis for these changes and into the differential response to calcium binding of various members of the EF-hand calcium-binding protein family. Factors contributing to the stability of the Ca2+-loaded open conformation are discussed, including a new hypothesis that critical hydrophobic interactions stabilize the open conformation in Ca2+ sensors, but are absent in "non-sensor" proteins that remain closed upon Ca2+ binding. A role for methionine residues in stabilizing the open conformation is also proposed.     

Thermally induced hydrogen exchange processes in small proteins as seen by FTIR spectroscopy

Fourier-transform infrared (FTIR) spectroscopy has been used to study the thermally induced exchange characteristics of those backbone amide protons which persist H-D exchange at ambient conditions in ribonuclease A, in wild type ribonuclease T1 and some of its variants, and in the histone-like protein HBsu. The H-D exchange processes were induced by increasing the thermal energy of the protein solutions in two ways: (i) by linearly increasing the temperature, and (ii) by a temperature jump. To trace the H-D exchange in the proteins, various infrared absorption bands known to be sensitive to H-D exchange were used as specific monitors. Characteristic H-D exchange curves were obtained from which the endpoints (T_H/D) of H-D exchange could be determined. The H-D exchange curves, the T_H/D-values and the phase transition temperatures T_m were used to estimate the structural flexibility and stability of the given proteins. It is suggested that time-resolved FTIR spectroscopy can be used to determine global stability parameters of proteins.     

Oligomerization and pigmentation dependent excitation energy transfer in fucoxanthin-chlorophyll proteins

The ultrafast caroteonid to chlorophyll a energy transfer dynamics of the isolated fucoxanthin-chlorophyll proteins FCPa and FCPb from the diatom Cyclotella meneghiniana was investigated in a comprehensive study using transient absorption in the visible and near infrared spectral region as well as static fluorescence spectroscopy. The altered oligomerization state of both antenna systems results in a more efficient energy transfer for FCPa, which is also reflected in the different chlorophyll a fluorescence quantum yields. We therefore assume an increased quenching in the higher oligomers of FCPb. The influence of the carotenoid composition was investigated using FCPa and FCPb samples grown under different light conditions and excitation wavelengths at the blue (500 nm) and red (550 nm) wings of the carotenoid absorption. The different light conditions yield in altered amounts of the xanthophyll cycle pigments diadinoxanthin and diatoxanthin. Since no significant dynamic changes are observed for high light and low light samples, the contribution of the xanthophyll cycle pigments to the energy transfer is most likely negligible. On the contrary, the observed dynamics change drastically for the different excitation wavelengths. The analyses of the decay associated spectra of FCPb suggest an altered energy transfer pathway. For FCPa even an additional time constant was found after excitation at 500 nm. It is assigned to the intrinsic lifetime of either the xanthophyll cycle carotenoids or more probable the blue absorbing fucoxanthins. Based on our studies we propose a detailed model explaining the different excitation energy transfer pathways in FCPa.     

Heme-binding storage proteins in the Chelicerata

Lipoglycoproteins in the Chelicerata that bind and store heme appear to represent a unique evolutionary strategy to both mitigate the toxicity of heme and utilize the molecule as a prosthetic group. Knowledge of heme-binding storage proteins in these organisms is in its infancy and much of what is known is from studies with vitellogenins (Vg) and more recently the main hemolymph storage protein in ixodid ticks characterized as a hemelipoglyco-carrier protein (CP). Data have also been reported from another arachnid, the black widow spider, Latrodectus mirabilis, and seem to suggest that the heme-binding capability of these large multimeric proteins is not a phenomenon found only in the Acari. CP appears to be most closely related to Vg in ticks in terms of primary structure but post-translational processing is different. Tick CP and L. mirabilis high-density lipoprotein 1 (HDL1) are similar in that they consist of two subunits of approximate molecular masses of 90 and 100 kDa, are found in the hemolymph as the dominant protein, and bind lipids, carbohydrates and cholesterol. CP binds heme which may also be the case for HDL1 since the protein was found to contain a brown pigment when analyzed by native polyacrylamide gel electrophoresis. Vgs in ticks are composed of multiple subunits and are the precursor of the yolk protein, vitellin. The phylogeny of these proteins, regulation of gene expression and putative functions of binding and storing heme throughout reproduction, blood-feeding and development are discussed. Comparisons with non-chelicerate arthropods are made in order to highlight the mechanisms and putative functions of heme-binding storage proteins and their possible critical function in the evolution of hematophagy.     

Outer membrane proteins: comparing X-ray and NMR structures by MD simulations in lipid bilayers

The structures of three bacterial outer membrane proteins (OmpA, OmpX and PagP) have been determined by both X-ray diffraction and NMR. We have used multiple (7 × 15 ns) MD simulations to compare the conformational dynamics resulting from the X-ray versus the NMR structures, each protein being simulated in a lipid (DMPC) bilayer. Conformational drift was assessed via calculation of the root mean square deviation as a function of time. On this basis the ‘quality’ of the starting structure seems mainly to influence the simulation stability of the transmembrane β-barrel domain. Root mean square fluctuations were used to compare simulation mobility as a function of residue number. The resultant residue mobility profiles were qualitatively similar for the corresponding X-ray and NMR structure-based simulations. However, all three proteins were generally more mobile in the NMR-based than in the X-ray simulations. Principal components analysis was used to identify the dominant motions within each simulation. The first two eigenvectors (which account for >50% of the protein motion) reveal that such motions are concentrated in the extracellular loops and, in the case of PagP, in the N-terminal α-helix. Residue profiles of the magnitude of motions corresponding to the first two eigenvectors are similar for the corresponding X-ray and NMR simulations, but the directions of these motions correlate poorly reflecting incomplete sampling on a ∼10 ns timescale.     

Comparison studies of the structural stability of rabbit prion protein with human and mouse prion proteins

Background: Prion diseases are fatal and infectious neurodegenerative diseases affecting humans and animals. Rabbits are one of the few mammalian species reported to be resistant to infection from prion diseases isolated from other species (I. Vorberg et al., Journal of Virology 77 (3) (2003) 2003-2009). Thus the study of rabbit prion protein structure to obtain insight into the immunity of rabbits to prion diseases is very important.Findings: The paper is a straight forward molecular dynamics simulation study of wild-type rabbit prion protein (monomer cellular form) which apparently resists the formation of the scrapie form. The comparison analyses with human and mouse prion proteins done so far show that the rabbit prion protein has a stable structure. The main point is that the enhanced stability of the C-terminal ordered region especially helix 2 through the D177-R163 salt-bridge formation renders the rabbit prion protein stable. The salt bridge D201-R155 linking helixes 3 and 1 also contributes to the structural stability of rabbit prion protein. The hydrogen bond H186-R155 partially contributes to the structural stability of rabbit prion protein.Conclusions: Rabbit prion protein was found to own the structural stability, the salt bridges D177-R163, D201-R155 greatly contribute and the hydrogen bond H186-R155 partially contributes to this structural stability. The comparison of the structural stability of prion proteins from the three species rabbit, human and mouse showed that the human and mouse prion protein structures were not affected by the removing these two salt bridges. Dima et al. (Biophysical Journal 83 (2002) 1268-1280 and Proceedings of the National Academy of Sciences of the United States of America 101 (2004) 15335-15340) also confirmed this point and pointed out that “correlated mutations that reduce the frustration in the second half of helix 2 in mammalian prion proteins could inhibit the formation of PrP^Sc”.     

Dielectric relaxation in proteins: the computational perspective

In photoexcitation and electron transfer, a new dipole or charge is introduced, and the structure is adjusted. This adjustment represents dielectric relaxation, which is the focus of this review. We concentrate on a few selected topics. We discuss linear response theory, as a unifying framework and a tool to describe non-equilibrium states. We review recent, molecular dynamics simulation studies that illustrate the calculation of dynamic and thermodynamic properties, such as Stokes shifts or reorganization free energies. We then turn to the macroscopic, continuum electrostatic view. We recall the physical definition of a dielectric constant and revisit the decomposition of the free energy into a reorganization and a static term. We review some illustrative continuum studies and discuss some difficulties that can arise with the continuum approach. In conclusion, we consider recent developments that will increase the accuracy and broaden the scope of all these methods.     

Role of the tertiary and quaternary structure in the formation of bis-histidyl adducts in cold-adapted hemoglobins

All tetrameric hemoglobins from Antarctic fish, including that from Trematomus bernacchii, HbTb form in the ferric state, promptly and distinctively from all the other tetrameric hemoglobins, a mixture of aquo-met at the α subunits and bis-histidyl adduct (hemichrome) at the β subunits. The role of the tertiary and quaternary structure in the hemichrome formation is unknown. Here we report the cloning, expression, purification, spectroscopic and computational characterization of the β-chain of HbTb (β-HbTb). Similarly to the human β-chains, β-HbTb self-assembles to form the homotetramer β(4)-HbTb; however, the latter quantitatively forms reversible ferric and ferrous bis-histidyl adducts, which are only partially present in the human tetramer (β(4)-HbA). A molecular dynamics study of the isolated β subunit of the two Hbs indicates that the ability to form hemichrome is an intrinsic feature of the chain; moreover, the greater propensity of β-HbTb to form the bis-histidyl adduct is probably linked to the higher flexibility of the CD loop region. On the bases of these experimental and computational results on the isolated chain, the influence of the quaternary structure on the stability of the endogenous ferrous and ferric hexa-coordination is also discussed.     

A comparative molecular dynamics study on thermostability of human and chicken prion proteins

To compare the thermostabilities of human and chicken normal cellular prion proteins (HuPrP(C) and CkPrP(C)), molecular dynamics (MD) simulations were performed for both proteins at an ensemble level (10 parallel simulations at 400 K and 5 parallel simulations at 300 K as a control). It is found that the thermostability of HuPrP(C) is comparable with that of CkPrP(C), which implicates that the non-occurrence of prion diseases in non-mammals cannot be completely attributed to the thermodynamic properties of non-mammalian PrP(C).     

The mechanism of proton translocation in respiratory complex I from molecular dynamics

Abstract

Respiratory complex I, the biggest enzyme of respiratory chain, plays a key role in energy production by the mitochondrial respiratory chain and has been implicated in many human neurodegenerative diseases. Recently, the crystal structure of respiratory complex I is reported. We perform 50?ns molecular dynamics simulations on the membrane domain of respiratory complex I under two hypothetical states (oxidized state and reduced state). We find that the density of water molecules in the trans-membrane domain under reduced state is bigger than that under oxidized state. The connecting elements (helix HL and β-hairpins-helix element) fluctuate stronger under reduced state than that under oxidized state, causing more internal water molecules and facilitating the proton conduction. The conformational changes of helix HL and the crucial charged residue Glu in TM5 play key roles in the mechanism of proton translocation. Our results illustrate the dynamic behavior and the potential mechanism of respiratory complex I, which provides the structural basis for drug design of respiratory complex I.     

Comparative structural analysis of two proteins belonging to quorum sensing system in Vibrio cholerae

Vibrio cholerae uses quorum sensing communication system to interact with other bacteria and for gauzing environmental parameters. This organism dwells equally well in both human host and aquatic environments. Quorum sensing regulates multitude of activities and is one of the lucrative targets presently pursued for drug design in bacteria to encounter virulence. Histidine phosphotransfer protein LuxU and response regulator LuxO of V. cholerae are known to play important roles in biofilms and virulence machinery. In the present study, we used computational methods to model LuxU and LuxO and simulated the interactions of LuxO and LuxU. Since no structural details of the proteins were available, we employed homology modeling to construct the three-dimensional structures and then performed molecular dynamics simulations to study dynamic behavior of the LuxO and LuxU from V. cholerae. The modeled proteins were validated and subjected to molecular docking analyses. This allowed us to predict the binding modes of the proteins to elucidate probable sites of interference.     

Comparison of the dynamics of the membrane-bound form of fd coat protein in micelles and in bilayers by solution and solid-state nitrogen-15 nuclear magnetic resonance spectroscopy

Solid-state and solution 15N nuclear magnetic resonance experiments on uniformly and specifically 15N labeled coat protein in phospholipid bilayers and in detergent micelles are used to describe the dynamics of the membrane-bound form of the protein. The residues in the N- and C-terminal portions of the coat protein in both phospholipid bilayers and in detergent micelles are mobile, while those in the hydrophobic midsection are immobile. There is evidence for a gradient of mobility in the C-terminal region of the coat protein in micelles; at 25 degrees C only the last two residues are mobile on the 10(9)-Hz timescale, while the last six to eight residues appear to be mobile on slower timescales and highly mobile at higher temperatures. Since all of the C-terminal residues are immobile in the virus particles, the mobility of these residues in the membrane-bound form of the protein may be important for the formation of protein-DNA interactions in the assembly process.     

Ultrafast multi-pulse transient absorption spectroscopy of fucoxanthin chlorophyll a protein from Phaeodactylum tricornutum

We have applied femtosecond transient absorption spectroscopy in pump-probe and pump-dump-probe regimes to study energy transfer between fucoxanthin and Chl a in fucoxanthin-Chl a complex from the pennate diatom Phaeodactylum tricornutum. Experiments were carried out at room temperature and 77?K to reveal temperature dependence of energy transfer. At both temperatures, the ultrafast (<100?fs) energy transfer channel from the fucoxanthin S₂ state is active and is complemented by the second pathway via the combined S₁/ICT state. The S₁/ICT-Chl a pathway has two channels, the fast one characterized by sub-picosecond energy transfer, and slow having time constants of 4.5?ps at room temperature and 6.6?ps at 77?K. The overall energy transfer via the S₁/ICT is faster at 77?K, because the fast component gains amplitude upon lowering the temperature. The pump-dump-probe regime, with the dump pulse centered in the spectral region of ICT stimulated emission at 950?nm and applied at 2?ps after excitation, proved that the S₁ and ICT states of fucoxanthin in FCP are individual, yet coupled entities. Analysis of the pump-dump-probe data suggested that the main energy donor in the slow S₁/ICT-Chl a route is the S₁ part of the S₁/ICT potential surface.