Large-scale Conformational Changes Research Articles

Precise Steric Features Control Aminoacyl-tRNA Accommodation on the Ribosome.

Protein synthesis involves a complex series of large-scale conformational changes in the ribosome. While long-lived intermediate states of these processes can be characterized by experiments, computational methods can be used to identify the interactions that contribute to the rate-limiting free-energy barriers. To this end, we use a simplified energetic model to perform molecular dynamics (MD) simulations of aminoacyl-tRNA (aa-tRNA) accommodation on the ribosome. While numerous studies have probed the energetics of the early stages of accommodation, we focus on the final stage of accommodation, where the 3'-CCA tail of aa-tRNA enters the peptidyl transferase center (PTC). These simulations show how a distinct intermediate is induced by steric confinement of the tail, immediately before it completes accommodation. Multiple pathways for 3'-CCA tail accommodation can be quantitatively distinguished, where the tail enters the PTC by moving past a pocket enclosed by Helix 89, 90, and 92, or through an alternate route formed by Helix 93 and the P-site tRNA. C2573, located within Helix 90, is shown to provide the largest contribution to this late-accommodation steric barrier, such that sub-Å perturbations to this residue can alter the time scale of tail accommodation by nearly an order of magnitude. In terms of biological function, these calculations suggest how this late-stage sterically induced barrier may contribute to tRNA proofreading by the ribosome.

Single-molecule experiments reveal the elbow as an essential folding guide in SMC coiled-coil arms

Structural maintenance of chromosome (SMC) complexes form ring-like structures through exceptional elongated coiled-coils (CCs). Recent studies found that variable CC conformations, including open and collapsed forms, which might result from discontinuities in the CC, facilitate the diverse functions of SMCs in DNA organization. However, a detailed description of the SMC CC architecture is still missing. Here, we study the structural composition and mechanical properties of SMC proteins with optical tweezers unfolding experiments using the isolated Psm3 CC as a model system. We find a comparatively unstable protein with three unzipping intermediates, which we could directly assign to CC features by crosslinking experiments and state-of-the-art prediction software. Particularly, the CC elbow is shown to be a flexible, potentially non-structured feature, which divides the CC into sections, induces a pairing shift from one CC strand to the other and could facilitate large-scale conformational changes, most likely via thermal fluctuations of the flanking CC sections. A replacement of the elbow amino acids hinders folding of the consecutive CC region and frequently leads to non-native misalignments, revealing the elbow as a guide for proper folding. Additional in vivo manipulation of the elbow flexibility resulted in impaired cohesin complexes, which directly link the sensitive CC architecture to the biological function of cohesin.

Domain-Based Protein Docking with Extremely Large Conformational Changes

Proteins are key components in many processes in living cells, and physical interactions with other proteins and nucleic acids often form key parts of their functions. In many cases, large flexibility of proteins as they interact is key to their function. To understand the mechanisms of these processes, it is necessary to consider the 3D structures of such protein complexes. When such structures are not yet experimentally determined, protein docking has long been present to computationally generate useful structure models. However, protein docking has long had the limitation that the consideration of flexibility is usually limited to very small movements or very small structures. Methods have been developed which handle minor flexibility via normal mode or other structure sampling, but new methods are required to model ordered proteins which undergo large-scale conformational changes to elucidate their function at the molecular level. Here, we present Flex-LZerD, a framework for docking such complexes. Via partial assembly multidomain docking and an iterative normal mode analysis admitting curvilinear motions, we demonstrate the ability to model the assembly of a variety of protein–protein and protein-nucleic acid complexes.

Protein folding intermediates on the dimensionality reduced landscape with UMAP and native contact likelihood

To understand protein folding mechanisms from molecular dynamics (MD) simulations, it is important to explore not only folded/unfolded states but also representative intermediate structures on the conformational landscape. Here, we propose a novel approach to construct the landscape using the uniform manifold approximation and projection (UMAP) method, which reduces the dimensionality without losing data-point proximity. In the approach, native contact likelihood is used as feature variables rather than the conventional Cartesian coordinates or dihedral angles of protein structures. We tested the performance of UMAP for coarse-grained MD simulation trajectories of B1 domain in protein G and observed on-pathway transient structures and other metastable states on the UMAP conformational landscape. In contrast, these structures were not clearly distinguished on the dimensionality reduced landscape using principal component analysis or time-lagged independent component analysis. This approach is also useful to obtain dynamical information through Markov state modeling and would be applicable to large-scale conformational changes in many other biomacromolecules.

Visualizing the Conformational Dynamics of Membrane Receptors Using Single-Molecule FRET.

The ability of cells to respond to external signals is essential for cellular development, growth, and survival. To respond to a signal from the environment, a cell must be able to recognize and process it. This task mainly relies on the function of membrane receptors, whose role is to convert signals into the biochemical language of the cell. G protein-coupled receptors (GPCRs) constitute the largest family of membrane receptor proteins in humans. Among GPCRs, metabotropic glutamate receptors (mGluRs) are a unique subclass that function as obligate dimers and possess a large extracellular domain that contains the ligand-binding site. Recent advances in structural studies of mGluRs have improved the understanding of their activation process. However, the propagation of large-scale conformational changes through mGluRs during activation and modulation is poorly understood. Single-molecule fluorescence resonance energy transfer (smFRET) is a powerful technique to visualize and quantify the structural dynamics of biomolecules at the single-protein level. To visualize the dynamic process of mGluR2 activation, fluorescent conformational sensors based on unnatural amino acid (UAA) incorporation were developed that allowed site-specific protein labeling without perturbation of the native structure of receptors. The protocol described here explains how to perform these experiments, including the novel UAA labeling approach, sample preparation, and smFRET data acquisition and analysis. These strategies are generalizable and can be extended to investigate the conformational dynamics of a variety of membrane proteins.

Cross-validation of distance measurements in proteins by PELDOR/DEER and single-molecule FRET

Pulsed electron-electron double resonance spectroscopy (PELDOR/DEER) and single-molecule Förster resonance energy transfer spectroscopy (smFRET) are frequently used to determine conformational changes, structural heterogeneity, and inter probe distances in biological macromolecules. They provide qualitative information that facilitates mechanistic understanding of biochemical processes and quantitative data for structural modelling. To provide a comprehensive comparison of the accuracy of PELDOR/DEER and smFRET, we use a library of double cysteine variants of four proteins that undergo large-scale conformational changes upon ligand binding. With either method, we use established standard experimental protocols and data analysis routines to determine inter-probe distances in the presence and absence of ligands. The results are compared to distance predictions from structural models. Despite an overall satisfying and similar distance accuracy, some inconsistencies are identified, which we attribute to the use of cryoprotectants for PELDOR/DEER and label-protein interactions for smFRET. This large-scale cross-validation of PELDOR/DEER and smFRET highlights the strengths, weaknesses, and synergies of these two important and complementary tools in integrative structural biology.

Architecture of the linker-scaffold in the nuclear pore.

INTRODUCTION In eukaryotic cells, the selective bidirectional transport of macromolecules between the nucleus and cytoplasm occurs through the nuclear pore complex (NPC). Embedded in nuclear envelope pores, the ~110-MDa human NPC is an ~1200-Å-wide and ~750-Å-tall assembly of ~1000 proteins, collectively termed nucleoporins. Because of the NPC's eightfold rotational symmetry along the nucleocytoplasmic axis, each of the ~34 different nucleoporins occurs in multiples of eight. Architecturally, the NPC's symmetric core is composed of an inner ring encircling the central transport channel and two outer rings anchored on both sides of the nuclear envelope. Because of its central role in the flow of genetic information from DNA to RNA to protein, the NPC is commonly targeted in viral infections and its nucleoporin constituents are associated with a plethora of diseases. RATIONALE Although the arrangement of most scaffold nucleoporins in the NPC's symmetric core was determined by quantitative docking of crystal structures into cryo-electron tomographic (cryo-ET) maps of intact NPCs, the topology and molecular details of their cohesion by multivalent linker nucleoporins have remained elusive. Recently, in situ cryo-ET reconstructions of NPCs from various species have indicated that the NPC's inner ring is capable of reversible constriction and dilation in response to variations in nuclear envelope membrane tension, thereby modulating the diameter of the central transport channel by ~200 Å. We combined biochemical reconstitution, high-resolution crystal and single-particle cryo-electron microscopy (cryo-EM) structure determination, docking into cryo-ET maps, and physiological validation to elucidate the molecular architecture of the linker-scaffold interaction network that not only is essential for the NPC's integrity but also confers the plasticity and robustness necessary to allow and withstand such large-scale conformational changes. RESULTS By biochemically mapping scaffold-binding regions of all fungal and human linker nucleoporins and determining crystal and single-particle cryo-EM structures of linker-scaffold complexes, we completed the characterization of the biochemically tractable linker-scaffold network and established its evolutionary conservation, despite considerable sequence divergence. We determined a series of crystal and single-particle cryo-EM structures of the intact Nup188 and Nup192 scaffold hubs bound to their Nic96, Nup145N, and Nup53 linker nucleoporin binding regions, revealing that both proteins form distinct question mark-shaped keystones of two evolutionarily conserved hetero‑octameric inner ring complexes. Linkers bind to scaffold surface pockets through short defined motifs, with flanking regions commonly forming additional disperse interactions that reinforce the binding. Using a structure‑guided functional analysis in Saccharomycescerevisiae, we confirmed the robustness of linker‑scaffold interactions and established the physiological relevance of our biochemical and structural findings. The near-atomic composite structures resulting from quantitative docking of experimental structures into human and S.cerevisiae cryo-ET maps of constricted and dilated NPCs structurally disambiguated the positioning of the Nup188 and Nup192 hubs in the intact fungal and human NPC and revealed the topology of the linker-scaffold network. The linker-scaffold gives rise to eight relatively rigid inner ring spokes that are flexibly interconnected to allow for the formation of lateral channels. Unexpectedly, we uncovered that linker‑scaffold interactions play an opposing role in the outer rings by forming tight cross-link staples between the eight nuclear and cytoplasmic outer ring spokes, thereby limiting the dilatory movements to the inner ring. CONCLUSION We have substantially advanced the structural and biochemical characterization of the symmetric core of the S.cerevisiae and human NPCs and determined near-atomic composite structures. The composite structures uncover the molecular mechanism by which the evolutionarily conserved linker‑scaffold establishes the NPC's integrity while simultaneously allowing for the observed plasticity of the central transport channel. The composite structures are roadmaps for the mechanistic dissection of NPC assembly and disassembly, the etiology of NPC‑associated diseases, the role of NPC dilation in nucleocytoplasmic transport of soluble and integral membrane protein cargos, and the anchoring of asymmetric nucleoporins. [Figure: see text].

Extended conformational states dominate the Hsp90 chaperone dynamics

The heat shock protein 90 (Hsp90) is a molecular chaperone central to client protein folding and maturation in eukaryotic cells. During its chaperone cycle, Hsp90 undergoes ATPase-coupled large-scale conformational changes between open and closed states, where the N-terminal and middle domains of the protein form a compact dimerized conformation. However, the molecular principles of the switching motion between the open and closed states remain poorly understood. Here we show by integrating atomistic and coarse-grained molecular simulations with small-angle X-ray scattering experiments and NMR spectroscopy data that Hsp90 exhibits rich conformational dynamics modulated by the charged linker, which connects the N-terminal with the middle domain of the protein. We show that the dissociation of these domains is crucial for the conformational flexibility of the open state, with the separation distance controlled by a β-sheet motif next to the linker region. Taken together, our results suggest that the conformational ensemble of Hsp90 comprises highly extended states, which could be functionally crucial for client processing.

Molecular characterization of linker and loop-mediated structural modulation and hinge motion in the C4-C5 domains of cMyBPC

IntroductionThe central C4 and C5 domains (C4C5) of cardiac myosin binding protein C (cMyBPC) contain a flexible interdomain linker and a cardiac-isoform specific loop. However, their importance in the functional regulation of cMyBPC has not been extensively studied. Methods and resultsWe expressed recombinant C4C5 proteins with deleted linker and loop regions and performed biophysical experiments to determine each of their structural and dynamic roles. We show that the linker and C5 loop regions modulate the secondary structure and thermal stability of C4C5. Furthermore, we provide evidence through extended molecular dynamics simulations and principle component analyses that C4C5 can adopt a completely bent or latched conformation. The simulation trajectory and interaction network analyses reveal that the completely bent conformation of C4C5 exhibits a specific pattern of residue-level interactions. Therefore, we propose a “hinge-and-latch” mechanism where the linker allows a great degree of flexibility and bending, while the loop aids in achieving a completely bent and latched conformation. Although this may be one of many bent positions that C4C5 can adopt, we illustrate for the first time in molecular detail that this type of large scale conformational change can occur in the central domains of cMyBPC. ConclusionsOur hinge-and-latch mechanism demonstrates that the linker and loop regions participate in dynamic modulation of cMyBPC’s motion and global conformation. These structural and dynamic features may contribute to muscle isoform-specific regulation of actomyosin activity, and have potential implications regarding its ability to propagate or retract cMyBPC’s regulatory N-terminal domains.

Impact of structural dynamics on biological functions of flaviviruses.

Flaviviruses comprise a number of mosquito- or tick-transmitted human pathogens of global public health importance. Advances in structural biology techniques have contributed substantially to our current understanding of the life cycle of these small enveloped RNA viruses and led to deep insights into details of virus assembly, maturation and cell entry. In addition to large-scale conformational changes and oligomeric rearrangements of envelope proteins during these processes, there is increasing evidence that smaller-scale protein dynamics (referred to as virus "breathing") can confer extra flexibility to these viruses for the fine-tuning of their interactions with the immune system and possibly with cellular factors they encounter in their complex ecological cycles in arthropod and vertebrate hosts. In this review, we discuss how work with tick-borne encephalitis virus has extended our view on flavivirus breathing, leading to the identification of a novel mechanism of antibody-mediated infection enhancement and demonstrating breathing intermediates of the envelope protein in the process of membrane fusion. These data are discussed in the context of other flaviviruses and the perspective of a potential role of virus breathing to cope with the requirements of adaptation and replication in evolutionarily very different hosts.

Extended sister-chromosome catenation leads to massive reorganization of the E. coli genome.

In bacteria, chromosome segregation occurs progressively from the origin to terminus within minutes of replication of each locus. Between replication and segregation, sister loci are held in an apparent cohesive state by topological links. The decatenation activity of topoisomerase IV (Topo IV) is required for segregation of replicated loci, yet little is known about the structuring of the chromosome maintained in a cohesive state. In this work, we investigated chromosome folding in cells with altered decatenation activities. Within minutes after Topo IV inactivation, massive chromosome reorganization occurs, associated with increased in contacts between nearby loci, likely trans-contacts between sister chromatids, and in long-range contacts between the terminus and distant loci. We deciphered the respective roles of Topo III, MatP and MukB when TopoIV activity becomes limiting. Topo III reduces short-range inter-sister contacts suggesting its activity near replication forks. MatP, the terminus macrodomain organizing system, and MukB, the Escherichia coli SMC, promote long-range contacts with the terminus. We propose that the large-scale conformational changes observed under these conditions reveal defective decatenation attempts involving the terminus area. Our results support a model of spatial and temporal partitioning of the tasks required for sister chromosome segregation.

Microscopic Characterization of the Chloride Permeation Pathway in the Human Excitatory Amino Acid Transporter 1 (EAAT1).

Excitatory amino acid transporters (EAATs) are glutamate transporters that belong to the solute carrier 1A (SLC1A) family. They couple glutamate transport to the cotransport of three sodium (Na+) ions and one proton (H+) and the counter-transport of one potassium (K+) ion. In addition to this coupled transport, binding of cotransported species to EAATs activates a thermodynamically uncoupled chloride (Cl-) conductance. Structures of SLC1A family members have revealed that these transporters use a twisting elevator mechanism of transport, where a mobile transport domain carries substrate and coupled ions across the membrane, while a static scaffold domain anchors the transporter in the membrane. We recently demonstrated that the uncoupled Cl- conductance is activated by the formation of an aqueous pore at the domain interface during the transport cycle in archaeal GltPh. However, a pathway for the uncoupled Cl- conductance has not been reported for the EAATs, and it is unclear if such a pathway is conserved. Here, we employ all-atom molecular dynamics (MD) simulations combined with enhanced sampling, free-energy calculations, and experimental mutagenesis to approximate large-scale conformational changes during the transport process and identified a Cl--conducting conformation in human EAAT1 (hEAAT1). Sampling the large-scale structural transitions in hEAAT1 allowed us to capture an intermediate conformation formed during the transport cycle with a continuous aqueous pore at the domain interface. The free-energy calculations performed for the conduction of Cl- and Na+ ions through the captured conformation highlight the presence of two hydrophobic gates that control low-barrier movement of Cl- through the aqueous pathway. Overall, our findings provide insights into the mechanism by which a human neurotransmitter transporter supports functional duality of active transport and passive Cl- permeation and confirm the commonality of this mechanism in different members of the SLC1A family.

Time-resolved functional expansion of visual rhodopsin detected by x-ray scattering

Visual rhodopsin is an excellent model system to understand spatio-temporal dynamics of class A G-protein-coupled receptors (GPCRs). We aim to understand the mechanistic details of how the local changes in the retinal cofactor propagate to become large scale conformational changes leading to visual signaling. Our hypothesis was that a volumetric expansion occurs in the light-activation process [1]. Here we show the first time-resolved X-ray scattering-based detection of the expansion of rhodopsin leading to the signaling state.

Simulations and geometric analysis of the ribosome reveal precise single-molecule strategies for monitoring head dynamics during translocation

The ribosome is a complex ribonucleoprotein machine that translates mRNA chains into polypeptide chains, which subsequently fold into proteins in all biological cells. The ribosome consists of a large subunit and a small subunit. The small subunit is further subdivided into a body domain and a head domain. The ribosome undergoes large-scale conformational changes during the elongation step of the process of translation. Structural studies and simulations reveal that the head domain undergoes coupled rotatory motions around multiple axes during mRNA-tRNA translocation.

Conformational dynamics of nuclear pore complexes

Nuclear pore complexes (NPCs) fuse the inner and outer nuclear membranes and mediatenucleocytoplasmic exchange. They are made of 30 different nucleoporins that form an intricatecylindrical architecture around an aqueous central channel. This architecture is highly dynamic inspace and time. Variations in NPC diameter were reported, but the physiological circumstances andthe molecular details remain unknown. Here we combined cryo-electron tomography andsubtomogram averaging with integrative structural modeling to capture a molecular movie of therespective large-scale conformational changes in cellulo.

Structure of the Core Postfusion Porcine Endogenous Retrovirus Fusion Protein.

ABSTRACTRetroviral elements from endogenous retroviruses have functions in mammalian physiology. The best-known examples are the envelope proteins that function in placenta development and immune suppression. Porcine endogenous retroviruses (PERVs) are an understudied class of endogenous retroviruses that infect cultured human cells, raising concern regarding porcine xenografts. The PERV envelope glycoprotein has also been proposed as a possible swine syncytin with a role in placental development. Despite the growing interest in PERVs, their envelope glycoproteins remain poorly characterized. Here, we successfully determined the postfusion crystal structure of the PERV core fusion ectodomain. The PERV fusion protein structure reveals a conserved class I viral fusion protein six-helix bundle. Biophysical experiments demonstrated that the thermodynamic stability of the PERV fusion protein secondary structure was the same at physiological and acidic pHs. A conserved surface analysis highlights the high degree of sequence conservation among retroviral fusogens in the chain reversal region that facilitates the large-scale conformational change required for membrane fusion. Further structural alignment of class I viral fusogens revealed a phylogenetic clustering that shows evolution into various lineages that correlate with virus mechanisms of cell entry. Our work indicates that structural dendrograms can be used to qualitatively infer insights into the fusion mechanisms of newly discovered class I viral fusogen structures.

Extended-ensemble docking to probe dynamic variation of ligand binding sites during large-scale structural changes of proteins.

Proteins can sample a broad landscape as they undergo conformational transition between different functional states. At the same time, as key players in almost all cellular processes, proteins are important drug targets. Considering the different conformational states of a protein is therefore central for a successful drug-design strategy. Here we introduce a novel docking protocol, termed extended-ensemble docking, pertaining to proteins that undergo large-scale (global) conformational changes during their function. In its application to multidrug ABC-transporter P-glycoprotein (Pgp), extensive non-equilibrium molecular dynamics simulations employing system-specific collective variables are first used to describe the transition cycle of the transporter. An extended set of conformations (extended ensemble) representing the full transition cycle between the inward- and the outward-facing states is then used to seed high-throughput docking calculations of known substrates, non-substrates, and modulators of the transporter. Large differences are predicted in the binding affinities to different conformations, with compounds showing stronger binding affinities to intermediate conformations compared to the starting crystal structure. Hierarchical clustering of the binding modes shows all ligands preferably bind to the large central cavity of the protein, formed at the apex of the transmembrane domain (TMD), whereas only small binding populations are observed in the previously described R and H sites present within the individual TMD leaflets. Based on the results, the central cavity is further divided into two major subsites, first preferably binding smaller substrates and high-affinity inhibitors, whereas the second one shows preference for larger substrates and low-affinity modulators. These central subsites along with the low-affinity interaction sites present within the individual TMD leaflets may respectively correspond to the proposed high- and low-affinity binding sites in Pgp. We propose further an optimization strategy for developing more potent inhibitors of Pgp, based on increasing its specificity to the extended ensemble of the protein, instead of using a single protein structure, as well as its selectivity for the high-affinity binding site. In contrast to earlier in silico studies using single static structures of Pgp, our results show better agreement with experimental studies, pointing to the importance of incorporating the global conformational flexibility of proteins in future drug-discovery endeavors.

Ceramide-1-phosphate transfer protein promotes sphingolipid reorientation needed for binding during membrane interaction

Lipid transfer proteins acquire and release their lipid cargoes by interacting transiently with source and destination biomembranes. In the GlycoLipid Transfer Protein (GLTP) superfamily, the two-layer all-α-helical GLTP-fold defines proteins that specifically target sphingolipids (SLs) containing either sugar or phosphate headgroups via their conserved but evolutionarily-modified SL recognitions centers. Despite comprehensive structural insights provided by X-ray crystallography, the conformational dynamics associated with membrane interaction and SL uptake/release by GLTP superfamily members have remained unknown. Herein, we report insights gained from molecular dynamics (MD) simulations into the conformational dynamics that enable ceramide-1-phosphate transfer proteins (CPTPs) to acquire and deliver ceramide-1-phosphate (C1P) during interaction with 1-palmitoyl-2-oleoyl phosphatidylcholine bilayers. The focus on CPTP reflects this protein’s involvement in regulating pro-inflammatory eicosanoid production and autophagy-dependent inflammasome assembly that drives interleukin (IL-1β and IL-18) production and release by surveillance cells. We found that membrane penetration by CPTP involved α-6 helix and the α-2 helix N-terminal region, was confined to one bilayer leaflet, and was relatively shallow. Large-scale dynamic conformational changes were minimal for CPTP during membrane interaction or C1P uptake except for the α-3/α-4 helices connecting loop, which is located near the membrane interface and interacts with certain phosphoinositide headgroups. Apart from functioning as a shallow membrane-docking element, α-6 helix was found to adeptly reorient membrane lipids to help guide C1P hydrocarbon chain insertion into the interior hydrophobic pocket of the SL binding site.These findings support a proposed ‘hydrocarbon chain-first’ mechanism for C1P uptake, in contrast to the ‘lipid polar headgroup-first’ uptake used by most lipid-transfer proteins.

Cryo-EM reveals local and global structural rearrangements in RYR mutants

Single-point mutations in ryanodine receptors (RYRs), large intracellular Ca2+ channels that play a critical role in EC coupling, are linked to debilitating and lethal disorders such as central core disease, malignant hyperthermia (for the skeletal isoform, RYR1), catecholaminergic polymorphic ventricular tachycardia, and ARVD2 (for the cardiac isoform, RYR2). Mutant RYRs result in elevated [Ca2+]cyto due to steady leak from the sarcoplasmic reticulum. To explore the nature of long-range allosteric mechanisms of malfunction, we determined the structure of two N-terminal domain mutants of RYR1, situated far away from the pore. Cryo-electron microscopy of the N-terminal subdomain A (NTDA) and subdomain C (NTDC) full-length mutants, RYR1 R163C (determined to 3.5 Å resolution), and RYR1 Y522S (determined to 4.0 Å resolution), respectively, reveal large-scale conformational changes in the cytoplasmic assembly under closed-state conditions (i.e., absence of activating Ca2+). The multidomain changes suggest that the mutations induce a preactivated state of the channel in R164C by altering the NTDA+/CD interface, and in Y522S by rearrangement of the α-helical bundle in NTDC. However, the extent of preactivation is considerably higher in Y522S as compared with R163C, which agrees with the increased severity of the Y522S mutation as established by various functional studies. The Y522S mutation represents loss of a spacer residue that is crucial for maintaining optimal orientation of α helices in NTDC, alteration of which has long-range effects felt as far away as ∼100 Å. Additionally, the structure of the Y522S mutant channel under open-state conditions also differs from RYR1 WT open channels. Our developing work with RYR mutants exhibits the diverse mechanisms by which these single-point mutations exert an effect on the channel's function and highlight the complexity of the multidomain channel, as well as the need for targeted therapies.

Modeling Catalysis in Allosteric Enzymes: Capturing Conformational Consequences.

Greater understanding of enzymatic mechanisms aids the discovery of new targets for biologics, the development of biocatalytic transformations, and de novo enzyme design. Methods using quantum mechanical (QM) potentials, such as Density Functional Theory (DFT), have enabled complex multistep enzymatic mechanisms to be studied, often in quantitative detail. Nevertheless, the dynamic interconversion of enzyme conformations between active and inactive catalytic forms, involving length- and timescales inaccessible to QM treatments, presents a formidable challenge for the development of computational models for allosterically modulated enzymes. We present an overview of the key concepts underlying multistate models of enzyme catalysis, enzyme allostery, and the challenge that large-scale conformational changes pose for methods using QM, QM/MM, and MM potentials. Structural clustering is highlighted as a valuable approach to bridge molecular dynamics conformational sampling of MM potentials and quantum chemical cluster models of catalysis. Particularly relevant to this discussion is structural allostery, which serves as the exemplar of conformational consequences. Here, a well-characterized allosteric enzyme, Imidazole Glycerol Phosphate Synthase (IGPS), is used to showcase the importance of multiple conformations and guide a new direction for qualitative understanding and quantitative modeling in enzyme catalysis.