Differences In Conformational Dynamics Research Articles

Conformational Differences between ERK2 and JNK3α1

Mitogen-activated protein kinases (MAPKs) are a family of protein kinases composed of extracellular-signal-regulated kinase (ERK), c-Jun N-terminal kinase (JNK), and p38. MAPKs regulate various cellular functions including cell proliferation, differentiation, responses to stresses, and cell survival. These MAPK family proteins are closely related in terms of amino acid sequences, but there are differences in their regulatory functions; for example, JNKs regulate responses to stress and inflammation while ERKs regulate events required for cell proliferation, motility, and differentiation. The high resolution X-ray crystal structures of ERK or JNK have been defined. However, differences in conformational dynamics between the MAPKs family have not clearly investigated yet. Hydrogen/deuterium exchange mass spectrometry (HDX-MS) has been focused as a powerful tool for understanding protein dynamics and determining conformational change. Here, we analyze ERK2 and JNK3α1 of MAPK family by HDX-MS to elucidate the differential conformational dynamics, which may explain different modules of MAPK signaling cascade.

DIRECT-ID: An automated method to identify and quantify conformational variations--application to β2 -adrenergic GPCR.

The conformational dynamics of a macromolecule can be modulated by a number of factors, including changes in environment, ligand binding, and interactions with other macromolecules, among others. We present a method that quantifies the differences in macromolecular conformational dynamics and automatically extracts the structural features responsible for these changes. Given a set of molecular dynamics (MD) simulations of a macromolecule, the norms of the differences in covariance matrices are calculated for each pair of trajectories. A matrix of these norms thus quantifies the differences in conformational dynamics across the set of simulations. For each pair of trajectories, covariance difference matrices are parsed to extract structural elements that undergo changes in conformational properties. As a demonstration of its applicability to biomacromolecular systems, the method, referred to as DIRECT-ID, was used to identify relevant ligand-modulated structural variations in the β2 -adrenergic (β2 AR) G-protein coupled receptor. Micro-second MD simulations of the β2 AR in an explicit lipid bilayer were run in the apo state and complexed with the ligands: BI-167107 (agonist), epinephrine (agonist), salbutamol (long-acting partial agonist), or carazolol (inverse agonist). Each ligand modulated the conformational dynamics of β2 AR differently and DIRECT-ID analysis of the inverse-agonist vs. agonist-modulated β2 AR identified residues known through previous studies to selectively propagate deactivation/activation information, along with some previously unidentified ligand-specific microswitches across the GPCR. This study demonstrates the utility of DIRECT-ID to rapidly extract functionally relevant conformational dynamics information from extended MD simulations of large and complex macromolecular systems.

Conformational Dynamics of Single HIV-1 Envelope Proteins on the Surface of Native Virions

The HIV-1 envelope (Env) mediates viral entry into host cells. While static images of Env in unliganded and ligand-bound forms define distinct conformations, direct observations of Env dynamics have yet to be realized. Here we apply single-molecule fluorescence resonance energy transfer (smFRET) imaging to elucidate the dynamics of native Env trimers on the surface of HIV-1 virions. Our observations indicated that unliganded HIV-1 Env transitions between three distinct pre-fusion conformations, which are affected by the viral receptor and co-receptors. Differences in conformational dynamics and ligand responsiveness of neutralization-sensitive and neutralization-resistant HIV-1 isolates delineated a dynamics-based mechanism of immune evasion. Broadly neutralizing antibodies stabilized one distinct pre-fusion conformation of Env, indicating the importance of the observed dynamics to HIV-1 Env function.

Conformational Preferences Underlying Reduced Activity of a Thermophilic Ribonuclease H

The conformational basis for reduced activity of the thermophilic ribonuclease HI enzyme from Thermus thermophilus, compared to its mesophilic homolog from Escherichia coli, is elucidated using a combination of NMR spectroscopy and molecular dynamics (MD) simulations. Explicit-solvent all-atom MD simulations of the two wild-type proteins and an E. coli mutant in which a glycine residue is inserted after position 80 to mimic the T. thermophilus protein reproduce the differences in conformational dynamics determined from 15N spin-relaxation NMR spectroscopy of three loop regions that surround the active site and contain functionally important residues: the glycine-rich region, the handle region, and the β5/αE loop. Examination of the MD trajectories indicates that the thermophilic protein samples conformations productive for substrate binding and activity less frequently than the mesophilic enzyme, although these differences may manifest as either increased or decreased relative flexibility of the different regions. Additional MD simulations indicate that mutations increasing activity of the T. thermophilus enzyme at mesophilic temperatures do so by reconfiguring the local environments of the mutated sites to more closely resemble active conformations. Taken together, the results show that both locally increased and decreased flexibility contribute to an overall reduction in activity of T. thermophilus ribonuclease H compared to its mesophilic E. coli homolog.

D-103 Conformational dynamics of single HIV-1 Env trimers on the surface of native virions

The HIV-1 envelope (Env) mediates viral entry into host cells. While static images of Env define distinct conformations, direct observations of Env dynamics have yet to be realized. To enable real-time imaging of conformational events in Env, we introduced fluorophores into the variable regions of the Env gp120 subunit and measured single-molecule fluorescence resonance energy transfer (smFRET) within the context of native trimers on the surface of HIV-1 virions. We observed unliganded HIV-1 Env to be intrinsically dynamic, transitioning between 3 distinct pre-fusion conformations, and receptor CD4 and co-receptor surrogate antibody 17b to remodel the conformational landscape of Env. Differences in conformational dynamics and ligand responsiveness of neutralization-sensitive and neutralization-resistant HIV-1 isolates delineated the dynamics-based mechanism of immune evasion. smFRET also revealed neutralizing antibodies VRC01, PG16, PGT128, PGT145, and 2G12, and the entry inhibitor, BMS-626529, to stabilize the ground state of Env, thereby providing dynamics-based strategies for therapeutic intervention.

Molecular dynamics simulations of the Nip7 proteins from the marine deep- and shallow-water Pyrococcus species

BackgroundThe identification of the mechanisms of adaptation of protein structures to extreme environmental conditions is a challenging task of structural biology. We performed molecular dynamics (MD) simulations of the Nip7 protein involved in RNA processing from the shallow-water (P. furiosus) and the deep-water (P. abyssi) marine hyperthermophylic archaea at different temperatures (300 and 373 K) and pressures (0.1, 50 and 100 MPa). The aim was to disclose similarities and differences between the deep- and shallow-sea protein models at different temperatures and pressures.ResultsThe current results demonstrate that the 3D models of the two proteins at all the examined values of pressures and temperatures are compact, stable and similar to the known crystal structure of the P. abyssi Nip7. The structural deviations and fluctuations in the polypeptide chain during the MD simulations were the most pronounced in the loop regions, their magnitude being larger for the C-terminal domain in both proteins. A number of highly mobile segments the protein globule presumably involved in protein-protein interactions were identified. Regions of the polypeptide chain with significant difference in conformational dynamics between the deep- and shallow-water proteins were identified.ConclusionsThe results of our analysis demonstrated that in the examined ranges of temperatures and pressures, increase in temperature has a stronger effect on change in the dynamic properties of the protein globule than the increase in pressure. The conformational changes of both the deep- and shallow-sea protein models under increasing temperature and pressure are non-uniform. Our current results indicate that amino acid substitutions between shallow- and deep-water proteins only slightly affect overall stability of two proteins. Rather, they may affect the interactions of the Nip7 protein with its protein or RNA partners.

Evolution of conformational dynamics determines the conversion of a promiscuous generalist into a specialist enzyme.

β-Lactamases are produced by many modern bacteria as a mechanism of resistance toward β-lactam antibiotics, the most common antibiotics in use. β-Lactamases, however, are ancient enzymes that originated billions of years ago. Recently, proteins corresponding to 2- to 3-Gy-old Precambrian nodes in the evolution of Class A β-lactamases have been prepared and shown to be moderately efficient promiscuous catalysts, able to degrade a variety of antibiotics with catalytic efficiency levels similar to those of an average modern enzyme. Remarkably, there are few structural differences (in particular at the active-site regions) between the resurrected enzymes and a penicillin-specialist modern β-lactamase. Here, we propose that the ancestral promiscuity originates from conformational dynamics. We investigate the differences in conformational dynamics of the ancient and extant β-lactamases through MD simulations and quantify the contribution of each position to functionally related dynamics through Dynamic Flexibility Index. The modern TEM-1 lactamase shows a comparatively rigid active-site region, likely reflecting adaptation for efficient degradation of a specific substrate (penicillin), whereas enhanced deformability at the active-site neighborhood in the ancestral resurrected proteins likely accounts for the binding and subsequent degradation of antibiotic molecules of different size and shape. Clustering of the conformational dynamics on the basis of Principal Component Analysis is in agreement with the functional divergence, as the ancient β-lactamases cluster together, separated from their modern descendant. Finally, our analysis leads to testable predictions, as sites of potential relevance for the evolution of dynamics are identified and mutations at those sites are expected to alter substrate-specificity.

Differences in conformational dynamics within the Hsp90 chaperone family reveal mechanistic insights.

The molecular chaperones of the Hsp90 family are essential in all eukaryotic cells. They assist late folding steps and maturation of many different proteins, called clients, that are not related in sequence or structure. Hsp90 interaction with its clients appears to be coupled to a series of conformational changes. Using hydrogen exchange mass spectrometry (HX-MS) we investigated the structural dynamics of human Hsp90β (hHsp90) and yeast Hsp82 (yHsp82). We found that eukaryotic Hsp90s are much more flexible than the previously studied Escherichia coli homolog (EcHtpG) and that nucleotides induce much smaller changes. More stable conformations in yHsp82 are obtained in presence of co-chaperones. The tetratricopeptide repeat (TPR) domain protein Cpr6 causes a different amide proton protection pattern in yHsp82 than the previously studied TPR-domain protein Sti1. In the simultaneous presence of Sti1 and Cpr6, protection levels are observed that are intermediate between the Sti1 and the Cpr6 induced changes. Surprisingly, no bimodal distributions of the isotope peaks are detected, suggesting that both co-chaperones affect both protomers of the Hsp90 dimer in a similar way. The cochaperones Sba1 was found previously in the crystal structure bound to the ATP hydrolysis-competent conformation of Hsp90, which did not allow to distinguish the mode of Sba1-mediated inhibition of Hsp90's ATPase activity by stabilizing the pre- or post-hydrolysis step. Our HX-MS experiments now show that Sba1 binding leads to a protection of the ATP binding lid, suggesting that it inhibits Hsp90's ATPase activity by slowing down product release. This hypothesis was verified by a single-turnover ATPase assay. Together, our data suggest that there are much smaller energy barriers between conformational states in eukaryotic Hsp90s than in EcHtpG and that co-chaperones are necessary in addition to nucleotides to stabilize defined conformational states.

Differences in Conformational Dynamics between Plasmodium falciparum and Human Hsp90 Orthologues Enable the Structure-Based Discovery of Pathogen-Selective Inhibitors

The high sequence conservation of druggable pockets of closely related proteins can make it challenging to develop selective inhibitors. We designed a new drug discovery approach that exploits both the static and dynamic differences of two orthologues. We applied it, as a proof of concept, to identify compounds that discriminate between the molecular chaperone Hsp90 of the protozoan pathogen Plasmodium falciparum (Pf) and that of its human host. We found that the ATP-binding pocket has a Pf-specific extension, whose sequence lining is identical in human Hsp90 but which differs by tertiary structure and dynamics. Using these insights for a structure-based drug screen, we discovered novel 7-azaindole compounds that exclusively bind the recombinant N-terminal domain of PfHsp90 but not of human Hsp90 nor of a PfHsp90 mutant with "human-like" dynamics. Moreover, these compounds preferentially inhibit the growth of yeast complemented by PfHsp90 and block the growth of Pf in culture.

Mechanistic Insights of β-Lactmases Evolution

β-lactmases, the enzymes produced by some bacteria, are responsible for the primary mechanism of resistance to β-lactam antibiotics, the most common antibiotics in commercial use. Earlier study has resurrected a series of ancient Class A β-lactmases in the laboratory. It showed that those 2-3 Gyr-old β-lactmases could degrade a variety of antibiotics in vitro with similar levels of catalytic efficiency, suggesting that β-lactmases evolve from substrate-promiscuous generalists to specialists. However, all those ancient β-lactmases adopt the same fold as the modern enzyme, especially at the catalytic sites. Here we propose that the enhanced substrate-promiscuity originates from protein dynamics. We investigate the differences in conformational dynamics of diverged β-lactmases through MD simulation and quantify the contribution of each residue site to functionally related dynamics through Dynamic Flexibility Index (DFI). The statistical clustering of the DFI pattern is later obtained through Principle Component Analysis (PCA). Strikingly, our result shows that those 2-3 Gyr-old β-lactmases are clustered together and far away from the modern enzyme in the principle subspace. This dynamic analysis also identifies those sites affecting dynamics the most that may play a critical in functional evolution of β-lactmases.

Light-Induced Differences in Conformational Dynamics of the Circadian Clock Regulator VIVID

The LOV (light–oxygen–voltage) domain protein VIVID (VVD) is a negative regulator of the circadian transcription factor White Collar Complex and controls light response and photoadaptation in Neurospora. Blue light converts VIVID from the dark state into the light state (VVDL) with concomitant homodimerization. Upon return to low-light conditions, VVD very slowly reverts back into the monomeric dark state (VVDD). To better understand the nature of the conformational changes that are the basis for the light–dark switch in VVD, we used hydrogen exchange mass spectrometry to probe solvent accessibility of backbone amide protons. Our data demonstrate that all structural elements of VVDD except for the N-cap region exchange according to the rare EX1 mechanism indicating a reversible unfolding with rather slow refolding rate. Interestingly, the unfolding halftimes of different elements were not identical but varied from 400 to 900s. VVDL also exchanges according to the EX1 mechanism, albeit with a halftime of 6h. Surprisingly, the dimerization interface showed very little protection suggesting a rapid dimer–monomer interconversion.

Collective Dynamics Differentiates Functional Divergence in Protein Evolution

Protein evolution is most commonly studied by analyzing related protein sequences and generating ancestral sequences through Bayesian and Maximum Likelihood methods, and/or by resurrecting ancestral proteins in the lab and performing ligand binding studies to determine function. Structural and dynamic evolution have largely been left out of molecular evolution studies. Here we incorporate both structure and dynamics to elucidate the molecular principles behind the divergence in the evolutionary path of the steroid receptor proteins. We determine the likely structure of three evolutionarily diverged ancestral steroid receptor proteins using the Zipping and Assembly Method with FRODA (ZAMF). Our predictions are within ∼2.7 Å all-atom RMSD of the respective crystal structures of the ancestral steroid receptors. Beyond static structure prediction, a particular feature of ZAMF is that it generates protein dynamics information. We investigate the differences in conformational dynamics of diverged proteins by obtaining the most collective motion through essential dynamics. Strikingly, our analysis shows that evolutionarily diverged proteins of the same family do not share the same dynamic subspace, while those sharing the same function are simultaneously clustered together and distant from those, that have functionally diverged. Dynamic analysis also enables those mutations that most affect dynamics to be identified. It correctly predicts all mutations (functional and permissive) necessary to evolve new function and ∼60% of permissive mutations necessary to recover ancestral function.

Structural Analysis of Prolyl Oligopeptidases Using Molecular Docking and Dynamics: Insights into Conformational Changes and Ligand Binding

Prolyl oligopeptidase (POP) is considered as an important pharmaceutical target for the treatment of numerous diseases. Despite enormous studies on various aspects of POPs structure and function still some of the questions are intriguing like conformational dynamics of the protein and interplay between ligand entry/egress. Here, we have used molecular modeling and docking based approaches to unravel questions like differences in ligand binding affinities in three POP species (porcine, human and A. thaliana). Despite high sequence and structural similarity, they possess different affinities for the ligands. Interestingly, human POP was found to be more specific, selective and incapable of binding to a few planar ligands which showed extrapolation of porcine POP in human context is more complicated. Possible routes for substrate entry and product egress were also investigated by detailed analyses of molecular dynamics (MD) simulations for the three proteins. Trajectory analysis of bound and unbound forms of three species showed differences in conformational dynamics, especially variations in β-propeller pore size, which was found to be hidden by five lysine residues present on blades one and seven. During simulation, β-propeller pore size was increased by ∼2 Å in porcine ligand-bound form which might act as a passage for smaller product movement as free energy barrier was reduced, while there were no significant changes in human and A. thaliana POPs. We also suggest that these differences in pore size could lead to fundamental differences in mode of product egress among three species. This analysis also showed some functionally important residues which can be used further for in vitro mutagenesis and inhibitor design. This study can help us in better understanding of the etiology of POPs in several neurodegenerative diseases.

The X‐ray crystallographic structure and specificity profile of HAD superfamily phosphohydrolase BT1666: Comparison of paralogous functions in B. thetaiotaomicron

Analysis of the haloalkanoate dehalogenase superfamily (HADSF) has uncovered homologues occurring within the same organism that are found to possess broad, overlapping substrate specificities, and low catalytic efficiencies. Here we compare the HADSF phosphatase BT1666 from Bacteroides thetaiotaomicron VPI-5482 to a homologue with high sequence identity (40%) from the same organism BT4131, a known hexose-phosphate phosphatase. The goal is to find whether these enzymes represent duplicated versus paralogous activities. The X-ray crystal structure of BT1666 was determined to 1.82 Å resolution. Superposition of the BT1666 and BT4131 structures revealed a conserved fold and identical active sites suggestive of a common physiological substrate. The steady-state kinetic constants for BT1666 were determined for a diverse panel of phosphorylated metabolites to define its substrate specificity profile and overall level of catalytic efficiency. Whereas BT1666 and BT4131 are both promiscuous, their substrate specificity profiles are distinct. The catalytic efficiency of BT1666 (k(cat) /K(m) = 4.4 × 10(2) M(-1) s(-1) for the best substrate fructose 1,6-(bis)phosphate) is an order of magnitude less than that of BT4131 (k(cat) /K(m) = 6.7 × 10(3) M(-1) s(-1) for 2-deoxyglucose 6-phosphate). The seemingly identical active-site structures point to sequence variation outside the active site causing differences in conformational dynamics or subtle catalytic positioning effects that drive the divergence in catalytic efficiency and selectivity. The overlapping substrate profiles may be understood in terms of differential regulation of expression of the two enzymes or a conferred advantage in metabolic housekeeping functions by having a larger range of possible metabolites as substrates.

Probing the Conformational Dynamics of Protein-DNA Complexes Through Hydration

The DNA-binding domain of the ETS family of transcription factors interacts with DNA in a sequence-selective manner, yet it tolerates remarkable variation in the recognition sequence. The biophysical basis for this sequence selectivity is currently unknown. Previous investigations have found significant differences in the thermodynamics of strong and weak ETS-DNA complexes. However, the ETS domain is structurally similar in both DNA-bound and unbound states. In the absence of significant coupled folding by either protein or DNA, selectivity may be linked to differences in conformational dynamics in the protein-DNA complex. Specifically, sequence selectivity may arise from differential reductions in conformational dynamics upon complex formation. We hypothesize that time-averaged differences in conformation dynamics can be resolved through their equilibrium hydration properties. The rationale for our hypothesis is that folded but conformationally mobile segments should be more solvent-accessible than segments that are relatively immobilized. To test this hypothesis, we are probing the hydration of various ETS-DNA complexes by high-resolution DNA and protein footprinting. We expect that this approach will provide a quantitative measure of a complex's overall conformational dynamics as well as identify structural elements that become conformationally stabilized upon complex formation.

Differences in conformational dynamics of [Pt3(HPTAB)]6+-DNA adducts with various cross-linking modes

We present here molecular dynamics simulations and DNA conformational dynamics for a series of trinuclear platinum [Pt3(HPTAB)]6+-DNA adducts [HPTAB = N,N,N′,N′,N′′,N′′-hexakis (2-pyridyl-methyl)-1,3,5-tris(aminomethyl) benzene], including three types of bifunctional crosslinks and four types of trifunctional crosslinks. Our simulation results reveal that binding of the trinuclear platinum compound to a DNA duplex induces the duplex unwinding in the vicinity of the platination sites, and causes the DNA to bend toward the major groove. As a consequence, this produces a DNA molecule whose minor groove is more widened and shallow compared to that of an undamaged bare-DNA molecule. Notably, for trifunctional crosslinks, we have observed extensive DNA conformational distortions, which is rarely seen for normal platinum–DNA adducts. Our findings, in this study, thus provide further support for the idea that platinum compounds with trifunctional intra-strand or long-range-inter-strand cross-linking modes can generate larger DNA conformational distortions than other types of cross-linking modes.

Molecular Dynamic Simulations of Cisplatin- and Oxaliplatin-d(GG) Intrastand Cross-links Reveal Differences in their Conformational Dynamics

Mismatch repair proteins, DNA damage-recognition proteins and translesion DNA polymerases discriminate between Pt-GG adducts containing cis-diammine ligands (formed by cisplatin (CP) and carboplatin) and trans-RR-diaminocyclohexane ligands (formed by oxaliplatin (OX)) and this discrimination is thought to be important in determining differences in the efficacy, toxicity and mutagenicity of these platinum anticancer agents. We have postulated that these proteins recognize differences in conformation and/or conformational dynamics of the DNA containing the adducts. We have previously determined the NMR solution structure of OX-DNA, CP-DNA and undamaged duplex DNA in the 5′-d(CCTCAGGCCTCC)-3′ sequence context and have shown the existence of several conformational differences in the vicinity of the Pt-GG adduct. Here we have used molecular dynamics simulations to explore differences in the conformational dynamics between OX-DNA, CP-DNA and undamaged DNA in the same sequence context. Twenty-five 10 ns unrestrained fully solvated molecular dynamics simulations were performed starting from two different DNA conformations using AMBER v8.0. All 25 simulations reached equilibrium within 4 ns, were independent of the starting structure and were in close agreement with previous crystal and NMR structures. Our data show that the cis-diammine (CP) ligand preferentially forms hydrogen bonds on the 5′ side of the Pt-GG adduct, while the trans-RR-diaminocyclohexane (OX) ligand preferentially forms hydrogen bonds on the 3′ side of the adduct. In addition, our data show that these differences in hydrogen bond formation are strongly correlated with differences in conformational dynamics, specifically the fraction of time spent in different DNA conformations in the vicinity of the adduct, for CP- and OX-DNA adducts. We postulate that differential recognition of CP- and OX-GG adducts by mismatch repair proteins, DNA damage-recognition proteins and DNA polymerases may be due, in part, to differences in the fraction of time that the adducts spend in a conformation favorable for protein binding.

A combined structural dynamics approach identifies a putative switch in factor VIIa employed by tissue factor to initiate blood coagulation

Coagulation factor VIIa (FVIIa) requires tissue factor (TF) to attain full catalytic competency and to initiate blood coagulation. In this study, the mechanism by which TF allosterically activates FVIIa is investigated by a structural dynamics approach that combines molecular dynamics (MD) simulations and hydrogen/deuterium exchange (HX) mass spectrometry on free and TF-bound FVIIa. The differences in conformational dynamics from MD simulations are shown to be confined to regions of FVIIa observed to undergo structural stabilization as judged by HX experiments, especially implicating activation loop 3 (residues 365-374{216-225}) of the so-called activation domain and the 170-loop (residues 313-322{170A-175}) succeeding the TF-binding helix. The latter finding is corroborated by experiments demonstrating rapid deglycosylation of Asn322 in free FVIIa by PNGase F but almost complete protection in the presence of TF or an active-site inhibitor. Based on MD simulations, a key switch of the TF-induced structural changes is identified as the interacting pair Leu305{163} and Phe374{225} in FVIIa, whose mutual conformations are guided by the presence of TF and observed to be closely linked to the structural stability of activation loop 3. Altogether, our findings strongly support an allosteric activation mechanism initiated by the stabilization of the Leu305{163}/Phe374{225} pair, which, in turn, stabilizes activation loop 3 and the S(1) and S(3) substrate pockets, the activation pocket, and N-terminal insertion.

A Molecular Dynamics Investigation of Mono and Dimeric States of the Outer Membrane Enzyme OMPLA

OMPLA is a phospholipase found in the outer membranes of many Gram-negative bacteria. Enzyme activation requires calcium-induced dimerisation plus bilayer perturbation. As the conformation of OMPLA in the different crystal forms (monomer versus dimer; with/without bound Ca 2+) is remarkably similar we have used multi-nanosecond molecular dynamics (MD) simulations to probe possible differences in conformational dynamics that may be related to enzyme activation. Simulations of calcium-free monomeric OMPLA, of the Ca 2+-bound dimer, and of the Ca 2+-bound dimer with a substrate analogue covalently linked to the active site serine have been performed, all with the protein embedded in a phospholipid (POPC) bilayer. All simulations were stable, but differences in the dynamic behaviour of the protein between the various states were observed. In particular, the stability of the active site and the hydrophobic substrate-binding cleft varied. Dimeric OMPLA is less flexible than monomeric OMPLA, especially around the active site. In the absence of bound substrate analogue, the hydrophobic substrate-binding cleft of dimeric OMPLA collapses. A model is proposed whereby the increased stability of the active site in dimeric OMPLA is a consequence of the local ordering of water around the nearby calcium ion. The observed collapse of the substrate-binding cleft may explain the experimentally observed occurrence of multiple dimer conformations of OMPLA, one of which is fully active while the other shows significantly reduced activity.

Conformational dynamics of a protein in the folded and the unfolded state

In a quasielastic neutron scattering experiment, the picosecond dynamics of α-amylase was investigated for the folded and the unfolded state of the protein. In order to ensure a reasonable interpretation of the internal protein dynamics, the protein was measured in D 2 O-buffer solution. The much higher structural flexibility of the pH induced unfolded state as compared to the native folded state was quantified using a simple analytical model, describing a local diffusion inside a sphere. In terms of this model the conformational volume, which is explored mainly by confined protein side-chain movements, is parameterized by the radius of a sphere (folded state, r =1.2 Å; unfolded state, 1.8 Å). Differences in conformational dynamics between the folded and the unfolded state of a protein are of fundamental interest in the field of protein science, because they are assumed to play an important role for the thermodynamics of folding/unfolding transition and for protein stability.