Structural characterization of the tunnels of Mycobacterium tuberculosis truncated hemoglobin N from molecular dynamics simulations

The structure of oxygenated trHbN from Mycobacterium tuberculosis shows an extended heme distal hydrogen‐bond network that includes Tyr33(B10), Gln58(E11), and the bound O₂. In addition, trHbN structure shows a network of hydrophobic cavities organized in two orthogonal branches. In the present work, the structure and the dynamics of oxygenated and deoxygenated trHbN in explicit water was investigated from 100 ns molecular dynamics (MD) simulations. Results show that, depending on the presence or the absence of a coordinated O₂, the Tyr33(B10) and Gln58(E11) side chains adopt two different configurations in concert with hydrogen bond network rearrangement. In addition, our data indicate that Tyr33(B10) and Gln58(E11) control the dynamics of Phe62(E15). In deoxy‐trHbN, Phe62(E15) is restricted to one conformation. Upon O₂ binding, the conformation of Gln58(E11) changes and residue Phe62(E15) fluctuates between two conformations. We also conducted a systematic study of trHbN tunnels by analyzing thousands of MD snapshots with CAVER. The results show that tunnel formation is the result of the dynamic reshaping of short‐lived hydrophobic cavities. The analyses indicate that the presence of these cavities is likely linked to the rigid structure of trHbN and also reveal two tunnels, EH and BE, that link the protein surface to the buried distal heme pocket and not present in the crystallographic structure. The cavities are sufficiently large to accomodate and store ligands. Tunnel dynamics in trHbN was found to be controlled by the side‐chain conformation of the Tyr33(B10), Gln58(E11), and Phe62(E15) residues. Importantly, in contrast to recently published works, our extensive systematic studies show that the presence or absence of a coordinated dioxygen does not control the opening of the long tunnel but rather the opening of the EH tunnel. In addition, the data lead to new and distinctly different conclusion on the impact of the Phe62(E15) residue on trHbN tunnels. We propose that the EH and the long tunnels are used for apolar ligands storage. The trajectories bring important new structural insights related to trHbN function and to ligand diffusion in proteins. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Conformation of the Ras-binding domain of Raf studied by molecular dynamics and free energy simulations

Recognition of Ras by its downstream target Raf is mediated by a Ras-recognition region in the Ras-binding domain (RBD) of Raf. Residues 78–89 in this region occupy two different conformations in the ensemble of NMR solution structures of the RBD: a fully α-helical one, and one where 87–90 form a type IV β-turn. Molecular dynamics simulations of the RBD in solution were performed to explore the stability of these and other possible conformations of both the wild-type RBD and the R89K mutant, which does not bind Ras. The simulations sample a fully helical conformation for residues 78–89 similar to the NMR helical structures, a conformation where 85–89 form a 3₁₀-helical turn, and a conformation where 87–90 form a type I |iB-turn, whose free energies are all within 0.3 kcal/mol of each other. NOE patterns and H^α chemical shifts from the simulations are in reasonable agreement with experiment. The NMR turn structure is calculated to be 3 kcal/mol higher than the three above conformations. In a simulation with the same implicit solvent model used in the NMR structure generation, the turn conformation relaxes into the fully helical conformation, illustrating possible structural artifacts introduced by the implicit solvent model. With the Raf R89K mutant, simulations sample a fully helical and a turn conformation, the turn being 0.9 kcal/mol more stable. Thus, the mutation affects the population of RBD conformations, and this is expected to affect Ras binding. For example, if the fully helical conformation of residues 78–89 is required for binding, its free energy increase in R89K will increase the binding free energy by about 0.6 kcal/mol. Proteins 31:186–200, 1998. © 1998 Wiley-Liss, Inc.     

Theoretical Investigations of Nitric Oxide Channeling in Mycobacterium tuberculosis Truncated Hemoglobin N

Mycobacterium tuberculosis group I truncated hemoglobin trHbN catalyzes the oxidation of nitric oxide (•NO) to nitrate with a second-order rate constant k ≈ 745 μM⁻¹ s⁻¹ at 23°C (nitric oxide dioxygenase reaction). It was proposed that this high efficiency is associated with the presence of hydrophobic tunnels inside trHbN structure that allow substrate diffusion to the distal heme pocket. In this work, we investigated the mechanisms of •NO diffusion within trHbN tunnels in the context of the nitric oxide dioxygenase reaction using two independent approaches. Molecular dynamics simulations of trHbN were performed in the presence of explicit •NO molecules. Successful •NO diffusion from the bulk solvent to the distal heme pocket was observed in all simulations performed. The simulations revealed that •NO interacts with trHbN at specific surface sites, composed of hydrophobic residues located at tunnel entrances. The entry and the internal diffusion of •NO inside trHbN were performed using the Long, Short, and EH tunnels reported earlier. The Short tunnel was preferentially used by •NO to reach the distal heme pocket. This preference is ascribed to its hydrophobic funnel-shape entrance, covering a large area extending far from the tunnel entrance. This funnel-shape entrance triggers the frequent formation of solvent-excluded cavities capable of hosting up to three •NO molecules, thereby accelerating •NO capture and entry. The importance of hydrophobicity of entrances for •NO capture is highlighted by a comparison with a polar mutant for which residues at entrances were mutated with polar residues. A complete map of •NO diffusion pathways inside trHbN matrix was calculated, and •NO molecules were found to diffuse from Xe cavity to Xe cavity. This scheme was in perfect agreement with the three-dimensional free-energy distribution calculated using implicit ligand sampling. The trajectories showed that •NO significantly alters the dynamics of the key amino acids of Phe⁶²(E15), a residue proposed to act as a gate controlling ligand traffic inside the Long tunnel, and also of Ile¹¹⁹(H11), at the entrance of the Short tunnel. It is noteworthy that •NO diffusion inside trHbN tunnels is much faster than that reported previously for myoglobin. The results presented in this work shed light on the diffusion mechanism of apolar gaseous substrates inside protein matrix.     

Structure, dynamics, and Hck interaction of full-length HIV-1 Nef

Nef is an HIV accessory protein that plays an important role in the progression of disease after viral infection. It interferes with numerous signaling pathways, one of which involves serine/threonine kinases. Here, we report the results of an NMR structural investigation on full-length Nef and its interaction with the entire regulatory domain of Hck (residues 72-256; Hck32L). A helical conformation was found at the N-terminus for residues 14-22, preceding the folded core domain. In contrast to the previously studied truncated Nef (Nef Δ1-39), the full-length Nef did not show any interactions of Trp57/Leu58 with the hydrophobic patch formed by helices α1 and α2. Upon Hck32L binding, the N-terminal anchor domain as well as the well-known SH3-binding site of Nef exhibited significant chemical shift changes. Upon Nef binding, resonance changes in the Hck spectrum were confined mostly to the SH3 domain, with additional effects seen for the connector between SH3 and SH2, the N-terminal region of SH2 and the linker region that contains the regulatory polyproline motif. The binding data suggest that in full-length Nef more than the core domain partakes in the interaction. The solution conformation of Hck32L was modeled using RDC data and compared with the crystal structure of the equivalent region in the inactivated, full-length Hck, revealing a notable difference in the relative orientations of the SH3 and SH2 domains. The RDC-based model combined with (15)N backbone dynamics data suggest that Hck32L adopts an open conformation without binding of the polyproline motif in the linker to the SH3 domain.     

Quantitative analysis of ligand migration from transition networks

In this work we use transition network analysis for the first time to investigate ligand migration in truncated hemoglobin (trHbN) and obtain kinetic information about the docking-site dynamics in the protein. A comparison with explicit water molecular dynamics simulations (100 ns in total) shows that the rate constants derived from the network analysis are realistic. The transition network analysis provides 1) The time-resolved connectivity network in the protein; 2) The half-lives of the docking sites; 3) The transition timescales between two given docking sites; and 4) The extent of population transfer among different docking sites of the protein as a function of lag time. We investigate the role of the Tyr³³ and Gln⁵⁸ residues in ligand migration by studying ligand migration in four mutants of trHbN. The mutation study suggests that residues Tyr³³ and Gln⁵⁸ stabilize the NO ligand in the Xe2 docking site of trHbN, thus facilitating the efficiency of the NO detoxification reaction.     

Ligand-induced dynamical regulation of NO conversion in Mycobacterium tuberculosis truncated hemoglobin-N

Mycobacterium tuberculosis, the causative agent of human tuberculosis, is forced into latency by nitric oxide produced by macrophages during infection. In response to nitrosative stress M. tuberculosis has evolved a defense mechanism that relies on the oxygenated form of "truncated hemoglobin" N (trHbN), formally acting as NO-dioxygenase, yielding the harmless nitrate ion. X-ray crystal structures have shown that trHbN hosts a two-branched protein matrix tunnel system, proposed to control diatomic ligand migration to the heme, as the rate-limiting step in NO conversion to nitrate. Extended molecular dynamics simulations (0.1 micros), employed here to characterize the factors controlling diatomic ligand diffusion through the apolar tunnel system, suggest that O2 migration in deoxy-trHbN is restricted to a short branch of the tunnel, and that O2 binding to the heme drives conformational and dynamical fluctuations promoting NO migration through the long tunnel branch. The simulation results suggest that trHbN has evolved a dual-path mechanism for migration of O2 and NO to the heme, to achieve the most efficient NO detoxification.     

NMR structures and localization of the potential fusion peptides and the pre-transmembrane region of SARS-CoV: Implications in membrane fusion

Severe acute respiratory syndrome-associated coronavirus (SARS-CoV) poses a serious public health hazard. The S2 subunit of the S glycoprotein of SARS-CoV carries out fusion between the virus and the host cells. However, the exact mechanism of the cell fusion process is not well understood. Current model suggests that a conformational transition, upon receptor recognition, of the two heptad core regions of S2 may expose the hydrophobic fusogenic peptide or fusion peptide for membrane insertion. Three regions of the S2 subunit have been proposed to be involved in cell–cell fusion. The N-terminal fusion peptide (FP, residues 770–788), an internal fusion peptide (IFP, residues 873–888) and the pre-transmembrane region (PTM, residues 1185–1202) demonstrated interactions with model lipid membranes and potentially involved in the fusion process. Here, we have determined atomic resolution structures of these three peptides in DPC detergent micelles by solution NMR. FP assumes α-helical conformation with significant distortion at the central Gly residues; enabling a close packing among sidechains of aromatic residues including W, Y and F. The 3-D structure of PMT is characterized by a helix–loop–helix with extensive aromatic interactions within the helices. IFP adopts a rather straight α-helical conformation defined by packing among sidechains of aromatic and aliphatic residues. Paramagnetic spin labeled NMR has demonstrated surface localization of PMT whereas FP and IFP inserted into the micelles. Collectively, data presented in this study will aid in understanding fusion mechanism of SARS-CoV.     

Multiple Tight Phospholipid-Binding Modes of α-Synuclein Revealed by Solution NMR Spectroscopy

'In dopaminergic neurons, α-synuclein (αS) partitions between a disordered cytosolic state and a lipid-bound state. Binding of αS to membrane phospholipids is implicated in its functional role in synaptic regulation, but also impacts fibril formation associated with Parkinson's disease. We describe here a solution NMR study in which αS is added to small unilamellar vesicles of a composition mimicking synaptic vesicles; the results provide evidence for multiple distinct phospholipid-binding modes of αS. Exchange between the free state and the lipid-bound αS state, and between different bound states is slow on the NMR timescale, being in the range of 1-10 s^− 1. Partitioning of the binding modes is dependent on lipid/αS stoichiometry, and tight binding with slow-exchange kinetics is observed at stoichiometries as low as 2:1. In all lipid-bound states, a segment of residues starting at the N-terminus of αS adopts an α-helical conformation, while succeeding residues retain the characteristics of a random coil. The 40 C-terminal residues remain dynamically disordered, even at high-lipid concentrations, but can also bind to lipids to an extent that appears to be determined by the fraction of cis X-Pro peptide bonds in this region. While lipid-bound αS exhibits dynamic properties that preclude its direct observation by NMR, its exchange with the NMR-visible free form allows for its indirect characterization. Rapid amide-amide nuclear Overhauser enhancement buildup points to a large α-helical conformation, and a distinct increase in fluorescence anisotropy attributed to Tyr39 indicates an ordered environment for this “dark state.” Titration of αS with increasing amounts of lipids suggests that the binding mode under high-lipid conditions remains qualitatively similar to that in the low-lipid case. The NMR data appear incompatible with the commonly assumed model where αS lies in an α-helical conformation on the membrane surface and instead suggest that considerable remodeling of the vesicles is induced by αS.     

Conformational study of cyclolinopeptide A. A distance geometry and molecular dynamics approach

The conformation of cyclolinopeptide A, c(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val), a naturally occurring peptide with remarkable cytoprotective activity, has been investigated by means of distance geometry calculations and molecular dynamics simulations. The starting points for all the calculations were an X-ray structure and other structures obtained from distance geometry calculations based on NMR data. Restrained and unrestrained molecular dynamics simulations are reported in vacuo and in CCl4. Structural and dynamic properties are investigated and compared with those experimentally determined. The conformation obtained from the MD simulations which best reproduces the NMR parameters is at the same time one of the most stable ones and is also fairly similar to the crystal structure. An explanation for the occurrence of multiple conformations in solution at room temperature is given.     

Structure of crambin in solution, crystal and in the trajectories of molecular dynamics simulations

The mechanisms of the three-dimensional crambin structure alterations in the crystalline environments and in the trajectories of the molecular dynamics simulations in the vacuum and crystal surroundings have been analyzed. In the crystalline state and in the solution the partial regrouping of remote intramolecular packing contacts, involved in the formation and stabilization of the tertiary structure of the crambin molecule, occurs in NMR structures. In the crystalline state it is initiated by the formation of the intermolecular contacts, the conformational influence of its appearance is distributed over the structure. The changes of the conformations and positions of the residues of the loop segments, where the intermolecular contacts of the crystal surroundings are preferably concentrated, are most observable. Under the influence of these contacts the principal change of the regular secondary structure of crambin is taking place: extension of the two-strand β structure to the three-strand structure with the participation of the single last residue N46 of the C-terminal loop. In comparison with the C-terminal loop the more profound changes are observed in the conformation and the atomic positions of the backbone atoms and in the solvent accessibility of the residues of the interhelical loop. In the solution of the ensemble of the 8 NMR structures relative accessibility to the solvent differs more noticeably also in the region of the loop segments and rather markedly in the interhelical loop. In the crambin cryogenic crystal structures the positions of the atoms of the backbone and/or side chain of 14–18 of 46 residues are discretely disordered. The disorganizations of at least 8 of 14 residues occur directly in the regions of the intermolecular contacts and another 5 residues are disordered indirectly through the intramolecular contacts with the residues of the intermolecular contacts. Upon the molecular dynamics simulation in the vacuum surrounding as in the solution of the crystalline structure of crambin the essential changes of the backbone conformation are caused by the intermolecular contacts absence, but partly masked by the structure changes owing to the nonpolar H atoms absence on the simulated structure. The intermolecular contact absence is partly manifested upon the molecular dynamics simulation of the crambin crystal with one protein molecule. Compared to the crystal structure the lengths of the interpeptide hydrogen bonds and other interresidue contacts in an average solution NMR structure are somewhat shorter and accordingly the energy of the interpeptide hydrogen bonds is better. This length shortening can occur at the stage of the refinement of the NMR structures of the crambin and other proteins by its energy minimizations in the vacuum surroundings and not exist in the solution protein structures.     

Restriction in the conformational flexibility of apoproteins in the presence of organic cosolvents: A consequence of the formation of “native-like conformation”

The influence of n-propanol on the overall α-helical conformation of β-globin, apocytochrome C, and the functional domain of streptococcal M49 protein (pepM49) and its consequence on the proteolysis of the respective proteins has been investigated. A significant amount of α-helical conformation is induced into these proteins atpH 6.0 and 4°C in the presence of relatively low concentrations of n-propanol. The induction of α-helical conformation into the proteins increased as a function of the propanol concentration, the maximum induction occurring around 30% n-propanol. In the case of α-globin, the fluorescence of its tryptophyl residues also increased as a function of n-propanol concentration, the midpoint of this transition being around 20% n-propanol. Furthermore, concomitant with the induction of helical conformation into these proteins, the proteolysis of their polypeptide chain by V8 protease also gets restricted. The α-helical conformation induced into α- and β-globin by n-propanol decreased as the temperature is raised from 4 to 24°C. In contrast, the α-helical conformation of both α- and β-chain (i.e., globin with noncovalently bound heme) did not exhibit such a sensitivity to this change in temperature. However, distinct differences exist between the n-propanol induced “α-helical conformation” of globins and the “α-helical conformation” of α- and β-chains. A cross-correlation of the n-propanol induced increase in the fluorescence of β-globin with the corresponding increase in the α-helical conformation of the polypeptide chain suggested that the fluorescence increase represents a structural change of the protein that is secondary to the induction of the α-helical conformation into the protein (i.e., an integration of the helical conformation induced to the segments of the polypeptide chain to influence the microenvironment of the tryptophyl residues). Presumably, the fluorescence increase is a consequence of the packing of the helical segments of globin to generate a “native-like structure.” The induction of α-helical conformation into these proteins in the presence of n-propanol and the consequent generation of “native-like conformation” is not unique to n-propanol. Trifluoroethanol, another helix-inducing organic solvent, also behaves in the same fashion as n-propanol. However, in contrast to the proteins described above, n-propanol could neither induce an α-helical conformation into performic acid oxidized RNAse-A nor restrict its proteolysis by proteases. Thus, the high sensitivity of apoproteins and the protein domains to assume α-helical conformation in the presence of low concentration of n-propanol with a concomitant restriction of the proteolytic susceptibility of their polypeptide chain appears to be unique to those proteins that exhibit high α-helical propensities. Apparently, this phenomenon of helix induction and the restriction of proteolysis reflects the formation of rudimentary tertiary interaction of the native protein and is unique to apoproteins or structural domains of α-helical proteins. Consistent with this concept, the induction of α-helical conformation into shorter polypeptide fragments of 30 residues, (e.g., α_1-30, which exists in an α-helical conformation in hemoglobin) is very low. Besides, this peptide exhibited neither the high sensitivity to the low concentrations of n-propanol seen with the apoproteins/protein domains nor the resistance toward proteolysis. The results suggest that the organic cosolvent induced decrease in the conformational flexibility of the apoprotein, and the consequent restriction of their proteolytic cleavage provides an opportunity to develop new strategies for protease catalyzed segment condensation reactions.     

Ligand binding to truncated hemoglobin N from Mycobacterium tuberculosis is strongly modulated by the interplay between the distal heme pocket residues and internal water

The survival of Mycobacterium tuberculosis requires detoxification of host *NO. Oxygenated Mycobacterium tuberculosis truncated hemoglobin N catalyzes the rapid oxidation of nitric oxide to innocuous nitrate with a second-order rate constant (k'(NOD) approximately 745 x 10(6) m(-1) x s(-1)), which is approximately 15-fold faster than the reaction of horse heart myoglobin. We ask what aspects of structure and/or dynamics give rise to this enhanced reactivity. A first step is to expose what controls ligand/substrate binding to the heme. We present evidence that the main barrier to ligand binding to deoxy-truncated hemoglobin N (deoxy-trHbN) is the displacement of a distal cavity water molecule, which is mainly stabilized by residue Tyr(B10) but not coordinated to the heme iron. As observed in the Tyr(B10)/Gln(E11) apolar mutants, once this kinetic barrier is lowered, CO and O(2) binding is very rapid with rates approaching 1-2 x 10(9) m(-1) x s(-1). These large values almost certainly represent the upper limit for ligand binding to a heme protein and also indicate that the iron atom in trHbN is highly reactive. Kinetic measurements on the photoproduct of the *NO derivative of met-trHbN, where both the *NO and water can be directly followed, revealed that water rebinding is quite fast (approximately 1.49 x 10(8) s(-1)) and is responsible for the low geminate yield in trHbN. Molecular dynamics simulations, performed with trHbN and its distal mutants, indicated that in the absence of a distal water molecule, ligand access to the heme iron is not hindered. They also showed that a water molecule is stabilized next to the heme iron through hydrogen-bonding with Tyr(B10) and Gln(E11).     

Lactam constraints provide insights into the receptor-bound conformation of secretin and stabilize a receptor antagonist

The natural ligands for family B G protein-coupled receptors are moderate-length linear peptides having diffuse pharmacophores. The amino-terminal regions of these ligands are critical for biological activity, with their amino-terminal truncation leading to production of orthosteric antagonists. The carboxyl-terminal regions of these peptides are thought to occupy a ligand-binding cleft within the disulfide-bonded amino-terminal domains of these receptors, with the peptides in amphipathic helical conformations. In this work, we have characterized the binding and activity of a series of 11 truncated and lactam-constrained secretin(5-27) analogues at the prototypic member of this family, the secretin receptor. One peptide in this series with lactam connecting residues 16 and 20 [c[E(16),K(20)][Y(10)]sec(5-27)] improved the binding affinity of its unconstrained parental peptide 22-fold while retaining the absence of endogenous biological activity and competitive antagonist characteristics. Homology modeling with molecular mechanics and molecular dynamics simulations established that this constrained peptide occupies the ligand-binding cleft in an orientation similar to that of natural full-length secretin and provided insights into why this peptide was more effective than other truncated conformationally constrained peptides in the series. This lactam bridge is believed to stabilize an extended α-helical conformation of this peptide while in solution and not to interfere with critical residue-residue approximations while docked to the receptor.     

Right-handed 14-helix in beta 3-peptides from L-aspartic acid monomers

beta-Peptides made from L-aspartic acid monomers form a new class of beta 3-peptides. Here we report the first three-dimensional NMR solution structure of a beta 3-hexapeptide (1) from L-aspartic acid monomers in 2,2,2-trifluoroethanol (TFE). We show that 1 forms a right-handed 14-helical structure in TFE. alpha-peptides from naturally occurring L-amino acids adopt a right-handed alpha-helix whereas beta 3-peptides formed from beta 3-amino acids derived from naturally occurring L-amino acids form left-handed 14-helices. The right-handed 14-helical conformation of 1 is a better mimic of alpha-peptide conformations. Using the NMR structure of 1 in TFE, we further study the conformation of 1 in water, as well as two similar beta 3-peptides (2 and 3) in water and TFE by molecular dynamics (MD) simulations. NMR and MD results suggest loss of secondary structure of 1 in water and show that it forms a fully extended structure. 2 and 3 contain residues with oppositely charged side chains that engage in salt-bridge interactions and dramatically stabilize the 14-helical conformation in aqueous media.     

Proline-induced constraints in alpha-helices

The disrupting effect of a prolyl residue on an α-helix has been analyzed by means of conformational energy computations. In the preferred, nearly α-helical conformations of Ac-Ala₄-Pro-NHMe and of Ac-Ala₇-Pro-Ala₇-NHMe, only the residue preceding Pro is not α-helical, while all other residues can occur in the α-helical A conformation; i.e., it is sufficient to introduce a conformational change of only one residue in order to accommodate proline in a distorted α-helix. Other low-energy conformations exist in which the conformational state of three residues preceding proline is altered considerably; on the other hand, another conformation in which these three residues retain the near-α-helical A-conformational state (with up to 26° changes of their dihedral angles ? and ψ, and a 48° change in one ω from those of the ideal α-helix) has a considerably higher energy. These conclusions are not altered by the substitution of other residues in the place of the Ala preceding Pro. The conformations of the peptide chain next to prolyl residues in or near an α-helix have been analyzed in 58 proteins of known structure, based on published atomic coordinates. Of 331 α-helices, 61 have a Pro at or next to their N-terminus, 21 have a Pro next to their C-terminus, and 30 contain a Pro inside the helix. Of the latter, 16 correspond to a break in the helix, 9 are located inside distorted first turns of the helix, and 5 are parts of irregular helices. Thus, the reported occurrence of prolyl residues next to or inside observed α-helices in proteins is consistent with the computed steric and energetic requirements of prolyl peptides.     

The structure and dynamics of rat apo-cellular retinol-binding protein II in solution: comparison with the X-ray structure

The structure and dynamics of rat apo-cellular retinol binding protein II (apo-CRBP II) in solution has been determined by multidimensional NMR analysis of uniformly enriched recombinant rat 13C, 15N-apo-CRBP II and 15N-apo-CRBP II. The final ensemble of 24 NMR structures has been calculated from 3274 conformational restraints or 24.4 restraints/residue. The average root-mean-square deviation of the backbone atoms for the final 24 structures relative to their mean structure is 1.06 A. Although the average solution structure is very similar to the crystal structure, it differs at the putative entrance to the binding cavity, which is formed by the helix-turn-helix motif, the betaC-betaD turn and the betaE-betaF turn. The mean coordinates of the main-chain atoms of amino acid residues 28-38 are displaced in the solution structure relative to the crystal structure. The side-chain of F58, located on the betaC-betaD turn, is reoriented such that it interacts with L37 and no longer blocks entry into the ligand-binding pocket. Residues 28-35, which form the second helix of the helix-turn-helix motif in the crystal structure, do not exhibit a helical conformation in the solution structure. The solution structure of apo-CRBP II exhibits discrete regions of backbone disorder which are most pronounced at residues 28-32, 37-38 and 73-76 in the betaE-betaF turn as evaluated by the consensus chemical shift index, the root-mean-square deviation, amide 1H exchange rates and 15N relaxation studies. These studies indicate that fluctuations in protein conformation occur on the microseconds to ms time-scale in these regions of the protein. Some of these exchange processes can be directly observed in the three-dimensional 15N-resolved NOESY spectrum. These results suggest that in solution, apo-CRBP II undergoes conformational changes on the microseconds to ms time-scale which result in increased access to the binding cavity.     

The Dynameomics rotamer library: amino acid side chain conformations and dynamics from comprehensive molecular dynamics simulations in water

We have recently completed systematic molecular dynamics simulations of 807 different proteins representing 95% of the known autonomous protein folds in an effort we refer to as Dynameomics. Here we focus on the analysis of side chain conformations and dynamics to create a dynamic rotamer library. Overall this library is derived from 31,000 occurrences of each of 86,217 different residues, or 2.7 × 10(9) rotamers. This dynamic library has 74% overlap of rotamer distributions with rotamer libraries derived from static high-resolution crystal structures. Seventy-five percent of the residues had an assignable primary conformation, and 68% of the residues had at least one significant alternate conformation. The average correlation time for switching between rotamers ranged from 22 ps for Met to over 8 ns for Cys; this time decreased 20-fold on the surface of the protein and modestly for dihedral angles further from the main chain. Side chain S(2) axis order parameters were calculated and they correlated well with those derived from NMR relaxation experiments (R = 0.9). Relationships relating the S(2) axis order parameters to rotamer occupancy were derived. Overall the Dynameomics rotamer library offers a comprehensive depiction of side chain rotamer preferences and dynamics in solution, and more realistic distributions for dynamic proteins in solution at ambient temperature than libraries derived from crystal structures, in particular charged surface residues are better represented. Details of the rotamer library are presented here and the library itself can be downloaded at http://www.dynameomics.org.     

Probing the influence of sequence-dependent interactions upon α-helix stability in alanine-based linear peptides

The observation that short, linear alanine-based polypeptides form stable α-helices in aqueous solution has allowed the development of well-defined experimental systems with which to study the influence of amino acid sequence upon the stability of secondary structure. We have performed detailed conformational searches upon six alanine-based peptides in order to rationalize the observed variation in the α-helical stability in terms of side-chain-backbone and side-chain-side-chain interactions. Although a simple, gas-phase, potential model was used to obtain the conformational energies for these peptides, good agreement was obtained with experiment regarding their relative α-helical stabilities. Our calculations clearly indicate that valine, isoleucine, and phenylalanine residues should destabilize the α-helical conformation when included within alanine-based peptides because of energetically unfavorable side-chain-backbone interactions, which tend to result in the formation of regions of 3₁₀-helix. In the case of valine, the destabilization most probably arises from entropic effects as the isopropyl side chain can assume more orientations in the 3₁₀-helical form of the peptide. A detailed examination of very short-range interactions in these peptides has also indicated that an interaction, involving fewer than five consecutive residues, whose stabilizing effect reinforces that of the (i, i + 4) hydrogen bond may be the basis of the requirement for increased nucleation (σ) and propagation parameters (s) required by Zimm–Bragg theory to predict the α-helical content for compounds in this class of short peptides. Our calculations complement recent work using modified Zimm–Bragg and Lifson–Roig theories of the helix–coil transition, and are consistent with molecular dynamics simulations upon linear peptides in aqueous solution. © 1993 John Wiley & Sons, Inc.     

Ligand migration in the truncated hemoglobin of Mycobacterium tuberculosis

The truncated hemoglobin of Mycobacterium tuberculosis (Mt-trHbO) is a small heme protein belonging to the hemoglobin superfamily. Truncated hemoglobins (trHbs) are believed to have functional roles such as terminal oxidases and oxygen sensors involved in the response to oxidative and nitrosative stress, nitric oxide (NO) detoxification, O?/NO chemistry, O? delivery under hypoxic conditions, and long-term ligand storage. Based on sequence similarities, they are classified into three groups. Experimental studies revealed that all trHbs display a 2-on-2 α-helical sandwich fold rather than the 3-on-3 α-helical sandwich fold of the classical hemoglobin fold. Using locally enhanced sampling (LESMD) molecular dynamics, the ligand-binding escape pathways from the distal heme binding cavity of Mt-trHbO were determined to better understand how this protein functions. The importance of specific residues, such as the group II and III invariant W(G8) residue, can be seen in terms of ligand diffusion pathways and ligand dynamics. LESMD simulations show that the wild-type Mt-trHbO has three diffusion pathways while the W(G8)F Mt-trHbO mutant has only two. The W(G8) residue plays a critical role in ligand binding and stabilization and helps regulate the rate of ligand escape from the distal heme pocket. Thus, this invariant residue is important in creating ligand diffusion pathways and possibly in the enzymatic functions of this protein.