Coarse-grained and all-atom modeling of structural states and transitions in hemoglobin

Probing the Structural Dynamics of the NMDA Receptor Activation by Coarse-Grained Modeling

N-Methyl-D-aspartate (NMDA) receptors are glutamate-gated excitatory channels that play essential roles in brain functions. High-resolution structures have been solved for an allosterically inhibited and agonist-bound form of a functional NMDA receptor; however, other key functional states (particularly the active open-channel state) were only resolved at moderate resolutions by cryo-electron microscopy (cryo-EM). To decrypt the mechanism of the NMDA receptor activation, structural modeling is essential to provide presently missing information about structural dynamics. We performed systematic coarse-grained modeling using an elastic network model and related modeling/analysis tools (e.g., normal mode analysis, flexibility and hotspot analysis, cryo-EM flexible fitting, and transition pathway modeling) based on an active-state cryo-EM map. We observed extensive conformational changes that allosterically couple the extracellular regulatory and agonist-binding domains to the pore-forming trans-membrane domain (TMD), and validated these, to our knowledge, new observations against known mutational and functional studies. Our results predict two key modes of collective motions featuring shearing/twisting of the extracellular domains relative to the TMD, reveal subunit-specific flexibility profiles, and identify functional hotspot residues at key domain-domain interfaces. Finally, by examining the conformational transition pathway between the allosterically inhibited form and the active form, we predict a discrete sequence of domain motions, which propagate from the extracellular domains to the TMD. In summary, our results offer rich structural and dynamic information, which is consistent with the literature on structure-function relationships in NMDA receptors, and will guide in-depth studies on the activation dynamics of this important neurotransmitter receptor.

Probing the structural dynamics of the CRISPR-Cas9 RNA-guided DNA-cleavage system by coarse-grained modeling

In the adaptive immune systems of many bacteria and archaea, the Cas9 endonuclease forms a complex with specific guide/scaffold RNA to identify and cleave complementary target sequences in foreign DNA. This DNA targeting machinery has been exploited in numerous applications of genome editing and transcription control. However, the molecular mechanism of the Cas9 system is still obscure. Recently, high-resolution structures have been solved for Cas9 in different structural forms (e.g., unbound forms, RNA-bound binary complexes, and RNA-DNA-bound tertiary complexes, corresponding to an inactive state, a pre-target-bound state, and a cleavage-competent or product state), which offered key structural insights to the Cas9 mechanism. To further probe the structural dynamics of Cas9 interacting with RNA and DNA at the amino-acid level of details, we have performed systematic coarse-grained modeling using an elastic network model and related analyses. Our normal mode analysis predicted a few key modes of collective motions that capture the observed conformational changes featuring large domain motions triggered by binding of RNA and DNA. Our flexibility analysis identified specific regions with high or low flexibility that coincide with key functional sites (such as DNA/RNA-binding sites, nuclease cleavage sites, and key hinges). We also identified a small set of hotspot residues that control the energetics of functional motions, which overlap with known functional sites and offer promising targets for future mutagenesis efforts to improve the specificity of Cas9. Finally, we modeled the conformational transitions of Cas9 from the unbound form to the binary complex and then the tertiary complex, and predicted a distinct sequence of domain motions. In sum, our findings have offered rich structural and dynamic details relevant to the Cas9 machinery, and will guide future investigation and engineering of the Cas9 systems. Proteins 2017; 85:342-353. © 2016 Wiley Periodicals, Inc.

Comparative Study of Elastic Network Model and Protein Contact Network for Protein Complexes: The Hemoglobin Case

The overall topology and interfacial interactions play key roles in understanding structural and functional principles of protein complexes. Elastic Network Model (ENM) and Protein Contact Network (PCN) are two widely used methods for high throughput investigation of structures and interactions within protein complexes. In this work, the comparative analysis of ENM and PCN relative to hemoglobin (Hb) was taken as case study. We examine four types of structural and dynamical paradigms, namely, conformational change between different states of Hbs, modular analysis, allosteric mechanisms studies, and interface characterization of an Hb. The comparative study shows that ENM has an advantage in studying dynamical properties and protein-protein interfaces, while PCN is better for describing protein structures quantitatively both from local and from global levels. We suggest that the integration of ENM and PCN would give a potential but powerful tool in structural systems biology.

A combined coarse-grained and all-atom simulation of TRPV1 channel gating and heat activation

Coarse-grained modeling and all-atom molecular dynamics simulation provide insight into the mechanism for heat activation of TRPV1 gating.

The transient receptor potential (TRP) channels act as key sensors of various chemical and physical stimuli in eukaryotic cells. Despite years of study, the molecular mechanisms of TRP channel activation remain unclear. To elucidate the structural, dynamic, and energetic basis of gating in TRPV1 (a founding member of the TRPV subfamily), we performed coarse-grained modeling and all-atom molecular dynamics (MD) simulation based on the recently solved high resolution structures of the open and closed form of TRPV1. Our coarse-grained normal mode analysis captures two key modes of collective motions involved in the TRPV1 gating transition, featuring a quaternary twist motion of the transmembrane domains (TMDs) relative to the intracellular domains (ICDs). Our transition pathway modeling predicts a sequence of structural movements that propagate from the ICDs to the TMDs via key interface domains (including the membrane proximal domain and the C-terminal domain), leading to sequential opening of the selectivity filter followed by the lower gate in the channel pore (confirmed by modeling conformational changes induced by the activation of ICDs). The above findings of coarse-grained modeling are robust to perturbation by lipids. Finally, our MD simulation of the ICD identifies key residues that contribute differently to the nonpolar energy of the open and closed state, and these residues are predicted to control the temperature sensitivity of TRPV1 gating. These computational predictions offer new insights to the mechanism for heat activation of TRPV1 gating, and will guide our future electrophysiology and mutagenesis studies.

Concepts in bio-molecular spectroscopy: vibrational case studies on metalloenzymes

Challenges and chances in bio-molecular spectroscopy are exemplified by vibrational case studies on metalloenzymes.

Probing the flexibility of tropomyosin and its binding to filamentous actin using molecular dynamics simulations

Tropomyosin (Tm) is a coiled-coil protein that binds to filamentous actin (F-actin) and regulates its interactions with actin-binding proteins like myosin by moving between three positions on F-actin (the blocked, closed, and open positions). To elucidate the molecular details of Tm flexibility in relation to its binding to F-actin, we conducted extensive molecular dynamics simulations for both Tm alone and Tm-F-actin complex in the presence of explicit solvent (total simulation time >400 ns). Based on the simulations, we systematically analyzed the local flexibility of the Tm coiled coil using multiple parameters. We found a good correlation between the regions with high local flexibility and a number of destabilizing regions in Tm, including six clusters of core alanines. Despite the stabilization by F-actin binding, the distribution of local flexibility in Tm is largely unchanged in the absence and presence of F-actin. Our simulations showed variable fluctuations of individual Tm periods from the closed position toward the open position. In addition, we performed Tm-F-actin binding calculations based on the simulation trajectories, which support the importance of Tm flexibility to Tm-F-actin binding. We identified key residues of Tm involved in its dynamic interactions with F-actin, many of which have been found in recent mutational studies to be functionally important, and the rest of which will make promising targets for future mutational experiments.

Non-site-specific allosteric effect of oxygen on human hemoglobin under high oxygen partial pressure

Protein allostery is essential for vital activities. Allosteric regulation of human hemoglobin (HbA) with two quaternary states T and R has been a paradigm of allosteric structural regulation of proteins. It is widely accepted that oxygen molecules (O₂) act as a “site-specific” homotropic effector, or the successive O₂ binding to the heme brings about the quaternary regulation. However, here we show that the site-specific allosteric effect is not necessarily only a unique mechanism of O₂ allostery. Our simulation results revealed that the solution environment of high O₂ partial pressure enhances the quaternary change from T to R without binding to the heme, suggesting an additional “non-site-specific” allosteric effect of O₂. The latter effect should play a complementary role in the quaternary change by affecting the intersubunit contacts. This analysis must become a milestone in comprehensive understanding of the allosteric regulation of HbA from the molecular point of view.

Crystal structures of a pentameric ligand-gated ion channel provide a mechanism for activation

Pentameric ligand-gated ion channels mediate fast chemical transmission of nerve signals. The structure of a bacterial proton-gated homolog has been established in its open and locally closed conformations at acidic pH. Here we report its crystal structure at neutral pH, thereby providing the X-ray structures of the two end-points of the gating mechanism in the same pentameric ligand-gated ion channel. The large structural variability in the neutral pH structure observed in the four copies of the pentamer present in the asymmetric unit has been used to analyze the intrinsic fluctuations in this state, which are found to prefigure the transition to the open state. In the extracellular domain (ECD), a marked quaternary change is observed, involving both a twist and a blooming motion, and the pore in the transmembrane domain (TMD) is closed by an upper bend of helix M2 (as in locally closed form) and a kink of helix M1, both helices no longer interacting across adjacent subunits. On the tertiary level, detachment of inner and outer β sheets in the ECD reshapes two essential cavities at the ECD-ECD and ECD-TMD interfaces. The first one is the ligand-binding cavity; the other is close to a known divalent cation binding site in other pentameric ligand-gated ion channels. In addition, a different crystal form reveals that the locally closed and open conformations coexist as discrete ones at acidic pH. These structural results, together with site-directed mutagenesis, physiological recordings, and coarse-grained modeling, have been integrated to propose a model of the gating transition pathway.

Reaction trajectory revealed by a joint analysis of protein data bank

Structural motions along a reaction pathway hold the secret about how a biological macromolecule functions. If each static structure were considered as a snapshot of the protein molecule in action, a large collection of structures would constitute a multidimensional conformational space of an enormous size. Here I present a joint analysis of hundreds of known structures of human hemoglobin in the Protein Data Bank. By applying singular value decomposition to distance matrices of these structures, I demonstrate that this large collection of structural snapshots, derived under a wide range of experimental conditions, arrange orderly along a reaction pathway. The structural motions along this extensive trajectory, including several helical transformations, arrive at a reverse engineered mechanism of the cooperative machinery (Ren, companion article), and shed light on pathological properties of the abnormal homotetrameric hemoglobins from α-thalassemia. This method of meta-analysis provides a general approach to structural dynamics based on static protein structures in this post genomics era.

The concept of allosteric interaction and its consequences for the chemistry of the brain

Throughout this Reflections article, I have tried to follow up on the genesis in the 1960s and subsequent evolution of the concept of allosteric interaction and to examine its consequences within the past decades, essentially in the field of the neuroscience. The main conclusion is that allosteric mechanisms built on similar structural principles operate in bacterial regulatory enzymes, gene repressors (and the related nuclear receptors), rhodopsin, G-protein-coupled receptors, neurotransmitter receptors, ion channels, and so on from prokaryotes up to the human brain yet with important features of their own. Thus, future research on these basic cybernetic sensors is expected to develop in two major directions: at the elementary level, toward the atomic structure and molecular dynamics of the conformational changes involved in signal recognition and transduction, but also at a higher level of organization, the contribution of allosteric mechanisms to the modulation of brain functions.

How does hemoglobin generate such diverse functionality of physiological relevance?

The absolute values of the O2-affinities (P50, Klow, and Khigh) of hemoglobin (Hb) are regulated neither by changes in the static T-/R-quaternary and associated tertiary structures nor the ligation states. They are pre-determined and regulated by the extrinsic environmental factors such as pH, buffers, and heterotropic effectors. The effect and role of O2 on Hb are reversibly to drive the structural allosteric equilibrium between the T(deoxy)- and R(oxy)-Hb toward R(oxy)-Hb (the structural allostery). R(oxy)-Hb has a higher O2-affinity (Khigh) relative to that (Klow) of the T(deoxy)-Hb (Khigh>Klow) under any fixed environmental conditions. The apparent O2-affinity of Hb is high, as the globin matrix interferes with the dissociation process of O2, forcing the dissociated O2 geminately to re-bind to the heme Fe. This artificially increases [oxy-Hb] and concomitantly decreases [deoxy-Hb], leading to the apparent increases of the O2-affinity of Hb. The effector-linked high-frequency thermal fluctuations of the globin matrix act as a gating mechanism to modulate such physical, energetic, and kinetic barriers to enhance the dissociation process of O2, resulted in increases in [deoxy-Hb] and concomitant decrease in [oxy-Hb], leading to apparent reductions of the O2-affinity of Hb (the entropic allostery). The heme in Hb is simply a low-affinity O2-trap, the coordination structure of which is not altered by static T-/R-quaternary and associated tertiary structural changes of Hb. Thus, heterotrophic effectors are the signal molecule, which acts as a functional link between these two allosteries and generates the diverse functionality of Hb of physiological relevance. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.