Protein dynamics

Network of Conformational Transitions Revealed by Molecular Dynamics Simulations of the Carbonic Anhydrase II Apo-Enzyme

Human carbonic anhydrase II (HCA II) is an enzyme that catalyzes the reversible hydration of CO2 into bicarbonate (HCO3-) and a proton (H+) as well as other reactions at an extremely high rate. This enzyme plays fundamental roles in human physiology/pathology, such as controlling the pH level in cells and so on. However, the binding mechanism between apo-HCA II and CO2 or other ligands as well as related conformational changes remains poorly understood, and atomic investigation into it could promote our understanding of related internal physiological/pathological mechanisms. In this study, long-time atomic molecular dynamics simulations as well as the clustering and free-energy analysis were performed to reveal the dynamics of apo-HCA II as well as the mechanism upon ligand binding. Our simulations indicate that the crystallographic B-factors considerably underestimate the loop dynamics: multiple conformations can be adopted by loops 1 and 2, especially for loop 1 because loop 1 is one side of the binding pocket, and its left-to-right movement can compress or extend the binding pocket, leading to one inactive (closed) state, three intermediate (semiopen) states, and one active (open) state; CO2 cannot get into the binding pocket of the inactive state but can get into those of intermediate and active states. The coexistence of multiple conformational states proposes a possible conformational selection model for the binding mechanism between apo-HCA II and CO2 or other ligands, revising our previous view of its functional mechanism of conformational change upon ligand binding and offering valuable structural insights into the workings of HCA II.

A Network of Conformational Transitions in the Apo Form of NDM-1 Enzyme Revealed by MD Simulation and a Markov State Model

New Delhi metallo-β-lactamase-1 (NDM-1) is a novel β-lactamase enzyme that confers enteric bacteria with nearly complete resistance to all β-lactam antibiotics, so it raises a formidable and global threat to human health. However, the binding mechanism between apo-NDM-1 and antibiotics as well as related conformational changes remains poorly understood, which largely hinders the overcoming of its antibiotic resistance. In our study, long-time conventional molecular dynamics simulation and Markov state models were applied to reveal both the dynamical and conformational landscape of apo-NDM-1: the MD simulation demonstrates that loop L3, which is responsible for antibiotic binding, is the most flexible and undergoes dramatic conformational changes; moreover, the Markov state model built from the simulation maps four metastable states including open, semiopen, and closed conformations of loop L3 as well as frequent transitions between the states. Our findings propose a possible conformational selection model for the binding mechanism between apo-NDM-1 and antibiotics, which facilitates the design of novel inhibitors and antibiotics.

A Network of Conformational Transitions Revealed by Molecular Dynamics Simulations of the Binary Complex of Escherichia coli 6-Hydroxymethyl-7,8-dihydropterin Pyrophosphokinase with MgATP

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the first reaction in the folate biosynthetic pathway. Comparison of its X-ray and nuclear magnetic resonance structures suggests that the enzyme undergoes significant conformational change upon binding to its substrates, especially in three catalytic loops. Experimental research has shown that, in its binary form, even bound by analogues of MgATP, loops 2 and 3 remain rather flexible; this raises questions about the putative large-scale induced-fit conformational change of the HPPK-MgATP binary complex. In this work, long-time all-atomic molecular dynamics simulations were conducted to investigate the loop dynamics in this complex. Our simulations show that, with loop 3 closed, multiple conformations of loop 2, including the open, semiopen, and closed forms, are all accessible to the binary complex. These results provide valuable structural insights into the details of conformational changes upon 6-hydroxymethyl-7,8-dihydropterin (HP) binding and biological activities of HPPK. Conformational network analysis and principal component analysis related to the loops are also discussed.

Molecular dynamics simulations of the Escherichia coli HPPK apo-enzyme reveal a network of conformational transitions

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the first reaction in the folate biosynthetic pathway. Comparison of its X-ray and nuclear magnetic resonance structures suggests that the enzyme undergoes significant conformational change upon binding to its substrates, especially in three catalytic loops. Experimental research has shown that even when confined by crystal contacts, loops 2 and 3 remain rather flexible when the enzyme is in its apo form, raising questions about the putative large-scale induced-fit conformational change of HPPK. To investigate the loop dynamics in a crystal-free environment, we performed conventional molecular dynamics simulations of the apo-enzyme at two different temperatures (300 and 350 K). Our simulations show that the crystallographic B-factors considerably underestimate the loop dynamics; multiple conformations of loops 2 and 3, including the open, semi-open, and closed conformations that an enzyme must adopt throughout its catalytic cycle, are all accessible to the apo-enzyme. These results revise our previous view of the functional mechanism of conformational change upon MgATP binding and offer valuable structural insights into the workings of HPPK. In this paper, conformational network analysis and principal component analysis related to the loops are discussed to support the presented conclusions.

Flexibility and molecular recognition in the immune system

Photon echo spectroscopy has been used to measure the response of three antibody-binding sites to perturbation from electronic excitation of a bound antigen, fluorescein. The three antibodies show motions that range in time scale from tens of femtoseconds to nanoseconds. Relative to the others, one antibody, 4-4-20, possesses a rigid binding site that likely results from a short and inflexible heavy chain complementarity-determining region 3 (HCDR3) loop and a critical Tyr that acts as a "molecular splint," rigidifying the antigen across its most flexible internal degree of freedom. The remaining two antibodies, 34F10 and 40G4, despite being generated against the same antigen, possess binding sites that are considerably more flexible. The more flexible combining sites likely result from longer HCDR3 loops and a deletion in the light chain complementarity-determining region 1 (LCDR1) that removes the critical Tyr residue. The binding site flexibilities may result in varying mechanisms of antigen recognition including lock-and-key, induced-fit, and conformational selection.

Investigation of the structural stability of the human acidic fibroblast growth factor by hydrogen-deuterium exchange

The conformational stability of the human acidic fibroblast growth factor (hFGF-1) is investigated using amide proton exchange and temperature-dependent chemical shifts, monitored by two-dimensional NMR spectroscopy. The change in free energy of unfolding (DeltaG(u)) of hFGF-1 is estimated to be 5.00 +/- 0.09 kcal.mol(-)(1). Amide proton-exchange rates of 74 residues (in hFGF-1) have been unambiguously measured, and the exchange process occurs predominately according to the conditions of the EX2 limit. The exchange rates of the fast-exchanging amide protons exposed to the solvent have been measured using the clean SEA-HSQC technique. The amide proton protection factor and temperature coefficient estimates show reasonably good correlation. Residues in beta-strands II and VI appear to constitute the stability core of the protein. Among the 12 beta-strands constituting the beta-barrel architecture of hFGF-1, beta-strand XI, located in the heparin binding domain, exhibits the lowest average protection factor value. Amide protons involved in the putative folding nucleation site in hFGF-1, identified by quench-flow NMR studies, do not represent the slow-exchanging core. Residues in portions of hFGF-1 experiencing high conformational flexibility mostly correspond to those involved in receptor recognition and binding.

Hydrogen/deuterium exchange and mass spectrometric analysis of a protein containing multiple disulfide bonds: Solution structure of recombinant macrophage colony stimulating factor-beta (rhM-CSFbeta)

Studies with the homodimeric recombinant human macrophage colony-stimulating factor beta (rhM-CSFbeta), show for the first time that a large number (9) of disulfide linkages can be reduced after amide hydrogen/deuterium (H/D) exchange, and the protein digested and analyzed successfully for the isotopic composition by electrospray mass spectrometry. Analysis of amide H/D after exchange-in shows that in solution the conserved four-helix bundle of (rhM-CSFbeta) has fast and moderately fast exchangeable sections of amide hydrogens in the alphaA helix, and mostly slow exchanging sections of amide hydrogens in the alphaB, alphaC, and alphaD helices. Most of the amide hydrogens in the loop between the beta1 and beta4 sheets exhibited fast or moderately fast exchange, whereas in the amino acid 63-67 loop, located at the interface of the two subunits, the exchange was slow. Solvent accessibility as measured by H/D exchange showed a better correlation with the average depth of amide residues calculated from reported X-ray crystallographic data for rhM-CSFalpha than with the average B-factor. The rates of H/D exchange in rhM-CSFbeta appear to correlate well with the exposed surface calculated for each amino acid residue in the crystal structure except for the alphaD helix. Fast hydrogen isotope exchange throughout the segment amino acids 150-221 present in rhM-CSFbeta, but not rhM-CSFalpha, provides evidence that the carboxy-terminal region is unstructured. It is, therefore, proposed that the anomalous behavior of the alphaD helix is due to interaction of the carboxy-terminal tail with this helical segment.

A view at the millennium: the efficiency of enzymatic catalysis

Binding TS in preference to S and increasing TDeltaS++by freezing out motions in E X S and E X TS have been accepted as the driving forces in enzymatic catalysis; however, the smaller value of DeltaG++ for a one-substrate enzymatic reaction, as compared to its nonenzymatic counterpart, is generally the result of a smaller value of DeltaH++. Ground-state conformers (E X NACs) are formed in enzymatic reactions that structurally resemble E X TS. E X NACs are in thermal equilibrium with all other E X S conformers and are turnstiles through which substrate molecules must pass to arrive at the lowest-energy TS. TS in E X TS may or may not be bound tighter than NAC in E X NAC.

A theoretical examination of the acid-catalyzed and noncatalyzed ring-opening reaction of an oxirane by nucleophilic addition of acetate. Implications to epoxide hydrolases

Ab initio and density functional calculations have been performed to gain a better understanding of the epoxide ring-opening reaction catalyzed by epoxide hydrolase. The S(N)2 reaction of acetate with 1S,2S-trans-2-methylstyrene oxide to provide the corresponding diol acetate ester was studied with and without general-acid catalysis. MP2 and DFT (B3LYP) calculations predict, for the noncatalyzed reaction, a central barrier of approximately 20-21 kcal/mol separating the reactants from products depending on which carbon center in the epoxide is undergoing attack. From these gas-phase reactions the immediate alkoxide products are not energetically far below their associated transition states such that the reaction is predicted to be endothermic. Inclusion of aqueous solvation effects via a polarizable continuum model predicts the activation barrier to increase by almost 10 kcal/mol due to the solvation of the acetate ion nucleophile. The activation barrier for the epoxide ring-opening reaction is reduced to approximately 10 kcal/mol when phenol, as the general-acid catalyst, is included in the gas-phase calculations. This is due to the immediate product being the neutral ester rather than the corresponding alkoxide. The transition state in the general-acid-catalyzed reaction is earlier than that for the noncatalyzed reaction and the reaction is highly exothermic. Molecular mechanics calculations of 1S,2S-trans-2-methylstyrene oxide in the active site of murine epoxide hydrolase show two possible binding conformations. Both conformers have the epoxide oxygen forming hydrogen bonds with the acidic hydrogens of the catalytic tyrosines (Tyr381 and Tyr465). These two conformations likely lead to different products since the nucleophile (Asp333-CO(2)(-)) is positioned to react with either carbon center in the epoxide.

The effect of matrix metalloproteinase complex formation on the conformational mobility of tissue inhibitor of metalloproteinases-2 (TIMP-2)

The backbone mobility of the N-terminal domain of tissue inhibitor of metalloproteinases-2 (N-TIMP-2) was determined both for the free protein and when bound to the catalytic domain of matrix metalloproteinase-3 (N-MMP-3). Regions of the protein with internal motion were identified by comparison of the T(1) and T(2) relaxation times and (1)H-(15)N nuclear Overhauser effect values for the backbone amide (15)N signals for each residue in the sequence. This analysis revealed rapid internal motion on the picosecond to nanosecond time scale for several regions of free N-TIMP-2, including the extended beta-hairpin between beta-strands A and B, which forms part of the MMP binding site. Evidence of relatively slow motion indicative of exchange between two or more local conformations on a microsecond to millisecond time scale was also found in the free protein, including two other regions of the MMP binding site (the CD and EF loops). On formation of a tight N-TIMP-2. N-MMP-3 complex, the rapid internal motion of the AB beta-hairpin was largely abolished, a change consistent with tight binding of this region to the MMP-3 catalytic domain. The extended AB beta-hairpin is not a feature of all members of the TIMP family; therefore, the binding of this highly mobile region to a site distant from the catalytic cleft of the MMPs suggests a key role in TIMP-2 binding specificity.

Molecular aspects of magnetic resonance imaging and spectroscopy

Magnetic resonance imaging (MRI) is a well known diagnostic tool in radiology that produces unsurpassed images of the human body, in particular of soft tissue. However, the medical community is often not aware that MRI is an important yet limited segment of magnetic resonance (MR) or nuclear magnetic resonance (NMR) as this method is called in basic science. The tremendous morphological information of MR images sometimes conceal the fact that MR signals in general contain much more information, especially on processes on the molecular level. NMR is successfully used in physics, chemistry, and biology to explore and characterize chemical reactions, molecular conformations, biochemical pathways, solid state material, and many other applications that elucidate invisible characteristics of matter and tissue. In medical applications, knowledge of the molecular background of MRI and in particular MR spectroscopy (MRS) is an inevitable basis to understand molecular phenomenon leading to macroscopic effects visible in diagnostic images or spectra. This review shall provide the necessary background to comprehend molecular aspects of magnetic resonance applications in medicine. An introduction into the physical basics aims at an understanding of some of the molecular mechanisms without extended mathematical treatment. The MR typical terminology is explained such that reading of original MR publications could be facilitated for non-MR experts. Applications in MRI and MRS are intended to illustrate the consequences of molecular effects on images and spectra.

Determinants of backbone dynamics in native BPTI: cooperative influence of the 14-38 disulfide and the Tyr35 side-chain

15Nitrogen relaxation experiments were used to characterize the backbone dynamics of two modified forms of bovine pancreatic trypsin inhibitor (BPTI). In one form, the disulfide between Cys14 and Cys38 in the wild-type protein was selectively reduced and methylated to generate an analog of the final intermediate in the disulfide-coupled folding pathway. The second form was generated by similarly modifying a mutant protein in which Tyr35 was replaced with Gly (Y35G). For both selectively reduced proteins, the overall conformation of native BPTI was retained, and the relaxation data for these proteins were compared with those obtained previously with the native wild-type and Y35G proteins. Removing the disulfide from either protein had only small effects on the observed longitudinal relaxation rates (R1) or heteronuclear cross relaxation rates (nuclear Overhauser effect), suggesting that the 14-38 disulfide has little influence on the fast (ps to ns) backbone dynamics of either protein. In the wild-type protein, the pattern of residues undergoing slower (micros to ms) internal motions, reflected in unusually large transverse relaxation rates (R2), was also largely unaffected by the removal of this disulfide. It thus appears that the large R2 rates previously observed in native wild-type protein are not a direct consequence of isomerization of the 14-38 disulfide. In contrast with the wild-type protein, reducing the disulfide in Y35G BPTI significantly decreased the number of backbone amides displaying large R2 rates. In addition, the frequencies of the backbone motions in the modified protein, estimated from R2 values measured at multiple refocusing delays, appear to span a wider range than those seen in native Y35G BPTI. Together, these observations suggest that the slow internal motions in Y35G BPTI are more independent in the absence of the 14-38 disulfide and that formation of this bond may lead to a substantial loss of conformational entropy. These effects may account for the previous observation that the Y35G substitution greatly destabilizes the disulfide. The results also demonstrate that the disulfide and the buried side-chain influence the dynamics of the folded protein in a highly cooperative fashion, with the effects of removing either being much greater in the absence of the other.

Molecular dynamics simulation of the rare amino acid LL-dityrosine and a dityrosine-containing peptide: comparison with time-resolved fluorescence

The fluorescence of the rare amino acid LL-dityrosine, which is found in insoluble biological materials with structural features, was recently shown to decay non-exponentially (Kungl et al. (1992) J. Fluorescence 2, 63-74). Here we investigated the time-resolved fluorescence of a dityrosine-containing peptide (DCP) to study the influence of side chains on the fluorescence decay of the chromophore. The fluorescence decay of DCP was best fitted by three exponential terms including a sub-nanosecond rise term, the values of which are quite similar to the parameters obtained for the decay of free dityrosine. They were found to depend on the pH of the aqueous solution but not on the temperature. Analysis by an exponential series method revealed broad fluorescence lifetime distributions for DCP. Compared to the corresponding analysis of dityrosine transients, similar lifetime centers were found whereas the widths of the distributions were found broader for DCP. Molecular dyamics (MD) simulations of dityrosine at 300 K show that chi 1 and chi 2 side chain conformers (rotamers) of both tyrosine subunits interconvert on a picosecond timescale. The rates of interconversion were shown to depend critically upon the MD technique applied: in vacuo simulations yielded lower interconversion rates compared to stochastic dynamics (SD) and full MD (water explicitly included). However, MD simulations of the dityrosine-containing peptide revealed no interconversions of the chi 1 and chi 2 side chain rotamers of both tyrosine subunits within a 400 ps trajectory. Interconversions could be induced by raising the temperature of the system (DCP plus solvent) to 340 K. Side chain rotamers of dityrosine are not stable on a fluorescence time scale but are stable when a dityrosine-containing peptide is regarded. Nevertheless both molecules yield similar fluorescence decay patterns. We therefore conclude that the rotamer model proposed for the fluorescence decay of tyrosine and tryptophan cannot be applied to the fluorescence decay of dityrosine and peptides containing this chromophore. This should be of future interest when dityrosine is used as an intrinsic sensor to study complex dityrosine-containing macromolecules by fluorescence spectroscopy.