Dynamics of the acetylcholinesterase tetramer

Acetylcholinesterase rapidly hydrolyzes the neurotransmitter acetylcholine in cholinergic synapses, including the neuromuscular junction. The tetramer is the most important functional form of the enzyme. Two low-resolution crystal structures have been solved. One is compact with two of its four peripheral anionic sites (PAS) sterically blocked by complementary subunits. The other is a loose tetramer with all four subunits accessible to solvent. These structures lacked the C-terminal amphipathic t-peptide (WAT domain) that interacts with the proline-rich attachment domain (PRAD). A complete tetramer model (AChEt) was built based on the structure of the PRAD/WAT complex and the compact tetramer. Normal mode analysis suggested that AChEt could exist in several conformations with subunits fluctuating relative to one another. Here, a multiscale simulation involving all-atom molecular dynamics and Cα-based coarse-grained Brownian dynamics simulations was carried out to investigate the large-scale intersubunit dynamics in AChEt. We sampled the ns-μs timescale motions and found that the tetramer indeed constitutes a dynamic assembly of monomers. The intersubunit fluctuation is correlated with the occlusion of the PAS. Such motions of the subunits “gate” ligand-protein association. The gates are open more than 80% of the time on average, which suggests a small reduction in ligand-protein binding. Despite the limitations in the starting model and approximations inherent in coarse graining, these results are consistent with experiments which suggest that binding of a substrate to the PAS is only somewhat hindered by the association of the subunits.     

Open interface and large quaternary structure movements in 3D domain swapped proteins: insights from molecular dynamics simulations of the C-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two three-dimensional (3D) domain swapped dimers. Crystallographic investigations have revealed that these dimers display completely different quaternary structures: one dimer (N-dimer), which presents the swapping of the N-terminal helix, is characterized by a compact structure, whereas the other (C-dimer), which is stabilized by the exchange of the C-terminal end, shows a rather loose assembly of the two subunits. The dynamic properties of monomeric RNase A and of the N-dimer have been extensively characterized. Here, we report a molecular dynamics investigation carried out on the C-dimer. This computational experiment indicates that the quaternary structure of the C-dimer undergoes large fluctuations. These motions do not perturb the proper folding of the two subunits, which retain the dynamic properties of RNase A and the N-dimer. Indeed, the individual subunits of the C-dimer display the breathing motion of the beta-sheet structure, which is important for the enzymatic activity of pancreatic-like ribonucleases. In contrast to what has been observed for the N-dimer, the breathing motion of the two subunits of the C-dimer is not coupled. This finding suggests that the intersubunit communications in a 3D domain swapped dimer strongly rely on the extent of the interchain interface. Furthermore, the observation that the C-dimer is endowed with a high intrinsic flexibility holds interesting implications for the specific properties of 3D domain swapped dimers. Indeed, a survey of the quaternary structures of the other 3D domain swapped dimers shows that large variations are often observed when the structural determinations are conducted in different experimental conditions. The 3D domain swapping phenomenon coupled with the high flexibility of the quaternary structure may be relevant for protein-protein recognition, and in particular for the pathological aggregations.     

Weak Intra-Ring Allosteric Communications of the Archaeal Chaperonin Thermosome Revealed by Normal Mode Analysis

Chaperonins are molecular machines that use ATP-driven cycles to assist misfolded substrate proteins to reach the native state. During the functional cycle, these machines adopt distinct nucleotide-dependent conformational states, which reflect large-scale allosteric changes in individual subunits. Distinct allosteric kinetics has been described for the two chaperonin classes. Bacterial (group I) chaperonins, such as GroEL, undergo concerted subunit motions within each ring, whereas archaeal and eukaryotic chaperonins (group II) undergo sequential subunit motions. We study these distinct mechanisms through a comparative normal mode analysis of monomer and double-ring structures of the archaeal chaperonin thermosome and GroEL. We find that thermosome monomers of each type exhibit common low-frequency behavior of normal modes. The observed distinct higher-frequency modes are attributed to functional specialization of these subunit types. The thermosome double-ring structure has larger contribution from higher-frequency modes, as it is found in the GroEL case. We find that long-range intersubunit correlation of amino-acid pairs is weaker in the thermosome ring than in GroEL. Overall, our results indicate that distinct allosteric behavior of the two chaperonin classes originates from different wiring of individual subunits as well as of the intersubunit communications.     

Analysis of the dynamics of rhizomucor miehei lipase at different temperatures.

The dynamics of Rhizomucor miehei lipase has been studied by molecular dynamics simulations at temperatures ranging from 200-500K. Simulations carried out in periodic boundary conditions and using explicit water molecules were performed for 400 ps at each temperature. Our results indicate that conformational changes and internal motions in the protein are significantly influenced by the temperature increase. With increasing temperature, the number of internal hydrogen bonds decreases, while surface accessibility, radius of gyration and the number of residues in random coil conformation increase. In the temperature range studied, the motions can be described in a low dimensional subspace, whose dimensionality decreases with increasing temperature. Approximately 80% of the total motion is described by the first (i) 80 eigenvectors at T=200K, (ii) 30 eigenvectors at T=300K and (iii) 10 eigenvectors at T=400K. At high temperature, the alpha-helix covering the active site in the native Rhizomucor miehei lipase, the helix at which end the active site is located, and in particular, the loop (Gly35-Lys50) show extensive flexibility.     

Dynamic properties of the N-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two 3-dimensional domain-swapped dimers with different quaternary structures. One dimer is characterized by the swapping of the C-terminal region (C-Dimer) and presents a rather loose structure. The other dimer (N-Dimer) exhibits a very compact structure with exchange of the N-terminal helix. Here we report the results of a molecular dynamics/essential dynamics (MD/ED) study carried out on the N-Dimer. This investigation, which represents the first MD/ED analysis on a three-dimensional domain-swapped enzyme, provides information on the dynamic properties of the active site residues as well as on the global motions of the dimer subunits. In particular, the analysis of the flexibility of the active site residues agrees well with recent crystallographic and site-directed mutagenesis studies on monomeric RNase A, thus indicating that domain swapping does not affect the dynamics of the active sites. A slight but significant rearrangement of N-Dimer quaternary structure, favored by the formation of additional hydrogen bonds at subunit interface, has been observed during the MD simulation. The analysis of collective movements reveals that each subunit of the dimer retains the functional breathing motion observed for RNase A. Interestingly, the breathing motion of the two subunits is dynamically coupled, as they open and close in phase. These correlated motions indicate the presence of active site intercommunications in this dimer. On these bases, we propose a speculative mechanism that may explain negative cooperativity in systems preserving structural symmetry during the allosteric transitions.     

Dimerization affects collective dynamics of triosephosphate isomerase

Molecular dynamics simulations (30-60 ns runs) are performed on free/apo triosephosphate isomerase (TIM) to determine any correlation between collective motions and loop 6 dynamics. Native TIM is reported to be active only as a homodimer even though cooperativity has not been observed between the two identical subunits. Both dimeric and monomeric (isolated from dimer) forms of TIM are simulated in explicit water at 300 K and 1 bar to inspect any differences between the structures in terms of fluctuation dynamics and functionally important loop 6 dynamics/closure. Significant cross-correlations between residue fluctuations are observed in the dimer, which result from the global counter-rotations of the two identical subunits in the essential modes of the dimer. Specifically, the first essential mode contributing to 34% of overall motion of the dimer is strongly coupled to the loop 6's closure over the active site. In contrast, such significant correlations cannot be observed in the monomeric structure, which maintains relatively localized motions of the loops in the essential modes. Thus, the onset of collective motions at ns time scale due to dimerization has functional implications as to the coordination of loop 6 closure.     

Dynamic features of homodimer interfaces calculated by normal‐mode analysis

Knowledge of the dynamic features of protein interfaces is necessary for a deeper understanding of protein–protein interactions. We performed normal‐mode analysis (NMA) of 517 nonredundant homodimers and their protomers to characterize dimer interfaces from a dynamic perspective. The motion vector calculated by NMA for each atom of a dimer was decomposed into internal and external motion vectors in individual component subunits, followed by the averaging of time‐averaged correlations between these vectors over atom pairs in the interface. This averaged correlation coefficient (ACC) was defined for various combinations of vectors and investigated in detail. ACCs decrease exponentially with an increasing interface area and r‐value, that is, interface area divided by the entire subunit surface area. As the r‐value reflects the nature of dimer formation, the result suggests that both the interface area and the nature of dimer formation are responsible for the dynamic properties of dimer interfaces. For interfaces with small or medium r‐values and without intersubunit entanglements, ACCs are found to increase on dimer formation when compared with those in the protomer state. In contrast, ACCs do not increase on dimer formation for interfaces with large r‐values and intersubunit entanglements such as in interwinding dimers. Furthermore, relationships between ACCs for intrasubunit atom pairs and for intersubunit atom pairs are found to significantly differ between interwinding and noninterwinding dimers for external motions. External motions are considered as an important factor for characterizing dimer interfaces.     

Dynamic contact network between ribosomal subunits enables rapid large-scale rotation during spontaneous translocation

During ribosomal translation, the two ribosomal subunits remain associated through intersubunit bridges, despite rapid large-scale intersubunit rotation. The absence of large barriers hindering rotation is a prerequisite for rapid rotation. Here, we investigate how such a flat free-energy landscape is achieved, in particular considering the large shifts the bridges undergo at the periphery. The dynamics and energetics of the intersubunit contact network are studied using molecular dynamics simulations of the prokaryotic ribosome in intermediate states of spontaneous translocation. Based on observed occupancies of intersubunit contacts, residues were grouped into clusters. In addition to the central contact clusters, peripheral clusters were found to maintain strong steady interactions by changing contacts in the course of rotation. The peripheral B1 bridges are stabilized by a changing contact pattern of charged residues that adapts to the rotational state. In contrast, steady strong interactions of the B4 bridge are ensured by the flexible helix H34 following the movement of protein S15. The tRNAs which span the subunits contribute to the intersubunit binding enthalpy to an almost constant degree, despite their different positions in the ribosome. These mechanisms keep the intersubunit interaction strong and steady during rotation, thereby preventing dissociation and enabling rapid rotation.     

Intrinsic molecular properties of the protein-protein bridge facilitate ratchet-like motion of the ribosome

The ribosomal intersubunit bridges maintain the overall architecture of the ribosome and thereby play a pivotal role in the dynamics of translation. The only protein-protein bridge, b1b, is formed by the two proteins, S13 and L5 of the small and large ribosomal subunits, respectively. B1b absorbs the largest movement during ratchet-like motion, and its two proteins reorganize in different constellations during this motion of the ribosome.Our results in this study of b1b in the Escherichia coli 70S ribosome suggest that the intrinsic molecular features of the bridging proteins allow the bridge to modulate the ratchet-like motion in a controlled manner. Additionally, another large subunit protein, L31, seems to participate with S13 and L5 in the formation, dynamics, and stabilization of this bridge.     

Temperature dependence of the structure of the substrate and active site of the Thermus thermophilus chorismate mutase E x S complex

Molecular dynamics (MD) simulations of Thermus thermophilus chorismate mutase substrate complex (TtCM x S) have been carried out at 298 K, 333 K, and the temperature of optimum activity: 343 K. The enzyme exists as trimeric subunits with active sites shared between two neighboring subunits. Two features distinguish intersubunit linkages of the thermophilic and mesophilic enzyme Bacillus subtilis chorismate mutase substrate complex (BsCM x S): (i) electrostatic interactions by intersubunit ion pairs (Arg3-Glu40*/41, Arg76-Glu51* and Arg69*-Asp101, residues labeled with an asterisk are from the neighboring subunit) in the TtCM x S are not present in the structure of the BsCM x S; and (ii) replacement of polar residues with short and nonpolar residues in the interstices of the TtCM x S tighten the intersubunit hydrophobic interactions compared to BsCM x S. Concerning the active site, electrostatic interactions of the critically placed Arg6 and Arg63* with the two carboxylates of chorismate place the latter in a reactive conformation to spontaneously undergo a Claisen rearrangement. The optimum geometry at the active site has the CZ atoms of the two arginines 11 A apart. With a decrease in temperature, Arg63* moves toward Arg6 and the average conformation structure of chorismate moves further away from the reactive ground state conformation. This movement is due to the decrease in distance separating the electrostatic (in the main) and hydrophobic interacting pairs holding the two subunits together.     

Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.     

Essential spaces defined by NMR structure ensembles and molecular dynamics simulation show significant overlap

Large concerted motions of proteins which span its “essential space,” are an important component of protein dynamics. We investigate to what extent structure ensembles generated with standard structure calculation techniques such as simulated annealing can capture these motions by comparing them to long-time molecular dynamics (MD) trajectories. The motions are analyzed by principal component analysis and compared using inner products of eigenvectors of the respective covariance matrices. Two very different systems are studied, the β-spectrin PH domain and the single-stranded DNA binding protein (ssDBP) from the filamentous phage Pf3. A comparison of the ensembles from NMR and MD shows significant overlap of the essential spaces, which in the case of ssDBP is extraordinarily high. The influence of variations in the specifications of distance restraints is investigated. We also study the influence of the selection criterion for the final structure ensemble on the definition of mobility. The results suggest a modified criterion that improves conformational sampling in terms of amplitudes of correlated motion. Proteins 31:370–382, 1998. © 1998 Wiley-Liss, Inc.     

Normal‐mode‐based modeling of allosteric couplings that underlie cyclic conformational transition in F1 ATPase

F₁ ATPase, a rotary motor comprised of a central stalk ( γ subunit) enclosed by three α and β subunits alternately arranged in a hexamer, features highly cooperative binding and hydrolysis of ATP. Despite steady progress in biophysical, biochemical, and computational studies of this fascinating motor, the structural basis for cooperative ATPase involving its three catalytic sites remains not fully understood. To illuminate this key mechanistic puzzle, we have employed a coarse‐grained elastic network model to probe the allosteric couplings underlying the cyclic conformational transition in F₁ ATPase at a residue level of detail. We will elucidate how ATP binding and product (ADP and phosphate) release at two catalytic sites are coupled with the rotation of γ subunit via various domain motions in α ₃ β ₃ hexamer (including intrasubunit hinge‐bending motions in β subunits and intersubunit rigid‐body rotations between adjacent α and β subunits). To this end, we have used a normal‐mode‐based correlation analysis to quantify the allosteric couplings of these domain motions to local motions at catalytic sites and the rotation of γ subunit. We have then identified key amino acid residues involved in the above couplings, some of which have been validated against past studies of mutated and γ ‐truncated F₁ ATPase. Our finding strongly supports a binding change mechanism where ATP binding to the empty catalytic site triggers a series of intra‐ and intersubunit domain motions leading to ATP hydrolysis and product release at the other two closed catalytic sites. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Isoflurane alters the structure and dynamics of GLIC

Pentameric ligand-gated ion channels are targets of general anesthetics. Although the search for discrete anesthetic binding sites has achieved some degree of success, little is known regarding how anesthetics work after the events of binding. Using the crystal structures of the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), which is sensitive to a variety of general anesthetics, we performed multiple molecular dynamics simulations in the presence and absence of the general anesthetic isoflurane. Isoflurane bound to several locations within GLIC, including the transmembrane pocket identified crystallographically, the extracellular (EC) domain, and the interface of the EC and transmembrane domains. Isoflurane also entered the channel after the pore was dehydrated in one of the simulations. Isoflurane disrupted the quaternary structure of GLIC, as evidenced in a striking association between the binding and breakage of intersubunit salt bridges in the EC domain. The pore-lining helix experienced lateral and inward radial tilting motion that contributed to the channel closure. Isoflurane binding introduced strong anticorrelated motions between different subunits of GLIC. The demonstrated structural and dynamical modulations by isoflurane aid in the understanding of the underlying mechanism of anesthetic inhibition of GLIC and possibly other homologous pentameric ligand-gated ion channels.     

Global ribosome motions revealed with elastic network model

The motions of large systems such as the ribosome are not fully accessible with conventional molecular simulations. A coarse-grained, less-than-atomic-detail model such as the anisotropic network model (ANM) is a convenient informative tool to study the cooperative motions of the ribosome. The motions of the small 30S subunit, the larger 50S subunit, and the entire 70S assembly of the two subunits have been analyzed using ANM. The lowest frequency collective modes predicted by ANM show that the 50S subunit and 30S subunit are strongly anti-correlated in the motion of the 70S assembly. A ratchet-like motion is observed that corresponds well to the experimentally reported ratchet motion. Other slow modes are also examined because of their potential links to the translocation steps in the ribosome. We identify several modes that may facilitate the E-tRNA exiting from the assembly. The A-site t-RNA and P-site t-RNA are found to be strongly coupled and positively correlated in these slow modes, suggesting that the translocations of these two t-RNAs occur simultaneously, while the motions of the E-site t-RNA are less correlated, and thus less likely to occur simultaneously. Overall the t-RNAs exhibit relatively large deformations. Animations of these slow modes of motion can be viewed at.     

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.     

Diffusion coefficients of segmentally flexible macromolecules with two subunits: a study of broken rods

A general treatment for the solution dynamics of segmentally flexible macromolecules having two subunits is presented. Bead modeling allows for a complete inclusion of hydrodynamic interactions in this treatment. The finite size of the beads is also considered, so that it is therefore possible to account properly for torsional motions of the subunits. Expressions for the components of the resistance matrix are derived. From them, the translational and rotational diffusion coefficients can be calculated. Distinction is made between hinged macromolecules, whose only internal motion is bending, and swivel-jointed macromolecules, for which torsions of the subunits are also allowed. Numerical results are presented for broken rods with the two types of flexibility. The effects of hydrodynamic interaction between arms of broken rods are about 25% for translation and under 10% for rotation. These findings give support to the treatments of Harvey, Wegener, and co-workers in which interactions were neglected. The rotational dynamics of hinged and swivel-jointed rods are compared. Although there are differences in the short-time behavior, the longest relaxation time is the same for the two cases. Finally, the validity of Wegener's rotational diffusion constants is discussed.     

Lactate dehydrogenase from the tetraploid weatherfish Misgurnus fossilis during temperature adaptation: determination of structural differences in two forms of the enzyme by molecular modeling methods

A structural analysis of two lactate dehydrogenase M4 protein forms has been performed. These structures are the protein products of two lactate dehydrogenase gene (LDH-A) copies in the weatherfish Misgurnus fossilis genome after thermal adaptation (acclimation) to 5 degrees C and 18 degrees C. The localization of three earlier identified amino acid substitutions (Gly214Val, Leu304Ile, Asp312Glu) has been determined, and the molecular dynamics simulation and computer modeling of two forms of the enzyme from skeletal muscles LDH-M4 have been carried out. After molecular dynamics trajectory calculations carried out at 5, 18, and 25 degrees C, the intersubunit distances for all structures used in calculations have been determined. It has been found that the Gly214Val substitution localized in the intersubunit region leads to a new intersubunit interaction, which plays a role in the stabilization of tetrameric enzyme structure after the adaptation to 18 degrees C.     

Conformational flexibility of neurophysin as investigated by local motions of fluorophores. Relationships with neurohypophyseal hormone binding

Flexibility of various structural domains of neurophysin and neurophysin-neurohypophyseal hormone complexes has been investigated through the fast rotational motion of fluorophores in highly viscous medium. Despite seven intrachain disulfide links, it is shown that some domains of neurophysin remain highly flexible. Dimerization of neurophysin does not affect the structural integrity of the individual subunits, each subdomain being conformationally equivalent within each protomer of the unliganded dimer. The absence of heterogeneous fluorescence anisotropy precludes the existence of a dimer tautomerization equilibrium. Binding of the hormonal ligands to neurophysin dimer promotes a large conformational change over the whole protein structure as assessed by differential alterations of the flexibility-rigidity and intrasegmental interaction properties of domains that do not participate directly to the dimerization/binding areas. The order of free-energy coupling between ligand binding and protein subunit association has been evaluated. Data are consistent with a model in which the first mole of bound ligand stabilizes the dimer by increasing the intersubunit contacts while the second mole of ligand induces most of the described conformational change. Accordingly, the positive cooperativity between the two dimeric binding sites is linked mainly to the binding of the second ligand. The induced structural change is perceived differently by each subunit as assessed by opposite local motions of Tyr49 in each liganded protomer and leads to the formation of a dimeric complex with a global pseudospherical symmetry although containing domains of local asymmetry.     

Nanosecond-timescale conformational dynamics of the human alpha7 nicotinic acetylcholine receptor

We explore the conformational dynamics of a homology model of the human alpha7 nicotinic acetylcholine receptor using molecular dynamics simulation and analyses of root mean-square fluctuations, block partitioning of segmental motion, and principal component analysis. The results reveal flexible regions and concerted global motions of the subunits encompassing extracellular and transmembrane domains of the subunits. The most relevant motions comprise a bending, hinged at the beta10-M1 region, accompanied by concerted tilting of the M2 helices that widens the intracellular end of the channel. Despite the nanosecond timescale, the observations suggest that tilting of the M2 helices may initiate opening of the pore. The results also reveal direct coupling between a twisting motion of the extracellular domain and dynamic changes of M2. Covariance analysis of interresidue motions shows that this coupling arises through a network of residues within the Cys and M2-M3 loops where Phe135 is stabilized within a hydrophobic pocket formed by Leu270 and Ile271. The resulting concerted motion causes a downward shift of the M2 helices that disrupts a hydrophobic girdle formed by 9' and 13' residues.