Large-scale Conformational Changes Research Articles

Nanomechanical collective vibration of SARS-CoV-2 spike proteins.

The development of effective therapeutics against COVID-19 requires a thorough understanding of the receptor recognition mechanism of the SARS-CoV-2 spike (S) protein. Here the multidomain collective dynamics on the trimer of the spike protein has been analyzed using normal mode analysis (NMA). A common nanomechanical profile was identified in the spike proteins of SARS-CoV-2 and its variants. The profile involves collective vibrations of the receptor-binding domain (RBD) and the N-terminal domain (NTD), which may mediate the physical interaction process. Quantitative analysis of the collective modes suggests a nanomechanical property involving large-scale conformational changes, which explains the difference in receptor binding affinity among different variants. These results support the use of intrinsic global dynamics as a valuable perspective for studying the allosteric and functional mechanisms of the S protein. This approach also provides a low-cost theoretical toolkit for screening potential pathogenic mutations and drug targets.

Time-Resolved Ion Mobility Mass Spectrometry to Solve Conformational Changes in a Cryptochrome.

Many biological mechanisms rely on the precise control of conformational changes in proteins. Understanding such dynamic processes requires methods for determining structures and their temporal evolution. In this study, we introduce a novel approach to time-resolved ion mobility mass spectrometry. We validated the method on a simple photoreceptor model and applied it to a more complex system, the animal-like cryptochrome from Chlamydomonas reinhardtii (CraCRY), to determine the role of specific amino acids affecting the conformational dynamics as reaction to blue light activation. In our setup, using a high-power LED mounted in the source region of an ion mobility mass spectrometer, we allow a time-resolved evaluation of mass and ion mobility spectra. Cryptochromes like CraCRY are a widespread type of blue light photoreceptors and mediate various light-triggered biological functions upon excitation of their inbuilt flavin chromophore. Another hallmark of cryptochromes is their flexible carboxy-terminal extension (CTE), whose structure and function as well as the details of its interaction with the photolyase homology region are not yet fully understood and differ among different cryptochromes types. Here, we addressed the highly conserved C-terminal domain of CraCRY, to study the effects of single mutations on the structural transition of the C-terminal helix α22 and the attached CTE upon lit-state formation. We show that D321, the putative proton acceptor of the terminal proton-coupled electron transfer event from Y373, is essential for triggering the large-scale conformational changes of helix α22 and the CTE in the lit state, while D323 influences the timing.

Testing the feasibility of targeting a conserved region on the S2 domain of the SARS-CoV-2 spike protein

The efficacy of vaccines against the SARS-CoV-2 virus significantly declines with the emergence of mutant strains, prompting investigation into the feasibility of targeting highly conserved but often cryptic regions on the S2 domain of spike protein. Using tools from molecular dynamics, we find that exposure of a conserved S2 epitope located in the central helices below the receptor binding domains would require large-scale motion beyond receptor binding domain up-down motion, but, along the reaction coordinates we explored, it is unlikely to be exposed by such large-scale dynamic fluctuations of the S1 domain without any external facilitating factors, despite some previous computational evidence suggesting transient exposure of this region. Furthermore, glycans, particularly those on N165 and N234, hinder S2-exposing opening dynamics, and thus stabilize spike in addition to immunologically shielding the protein surface. Although the S2 epitope region examined here is central to large-scale conformational changes during viral entry, free energy landscape analysis obtained using the path coordinate formalism reveals no inherent “loaded spring” effect, suggesting that a vaccine immunogen would tend to present the epitope in a prefusion-like conformation and may be effective in neutralization. These findings contribute to a deeper understanding of the dynamic origins of the function of the spike protein, as well as further characterizing the feasibility of the S2 epitope as a therapeutic target.

Structural transitions modulate the chaperone activities of Grp94

Hsp90s are ATP-dependent chaperones that collaborate with co-chaperones and Hsp70s to remodel client proteins. Grp94 is the ER Hsp90 homolog essential for folding multiple secretory and membrane proteins. Grp94 interacts with the ER Hsp70, BiP, although the collaboration of the ER chaperones in protein remodeling is not well understood. Grp94 undergoes large-scale conformational changes that are coupled to chaperone activity. Within Grp94, a region called the pre-N domain suppresses ATP hydrolysis and conformational transitions to the active chaperone conformation. In this work, we combined in vivo and in vitro functional assays and structural studies to characterize the chaperone mechanism of Grp94. We show that Grp94 directly collaborates with the BiP chaperone system to fold clients. Grp94's pre-N domain is not necessary for Grp94-client interactions. The folding of some Grp94 clients does not require direct interactions between Grp94 and BiP in vivo, suggesting that the canonical collaboration may not be a general chaperone mechanism for Grp94. The BiP co-chaperone DnaJB11 promotes the interaction between Grp94 and BiP, relieving the pre-N domain suppression of Grp94's ATP hydrolysis activity. In structural studies, we find that ATP binding by Grp94 alters the ATP lid conformation, while BiP binding stabilizes a partially closed Grp94 intermediate. Together, BiP and ATP push Grp94 into the active closed conformation for client folding. We also find that nucleotide binding reduces Grp94's affinity for clients, which is important for productive client folding. Alteration of client affinity by nucleotide binding may be a conserved chaperone mechanism for a subset of ER chaperones.

Archaeal GPN-loop GTPases involve a lock-switch-rock mechanism for GTP hydrolysis.

GPN-loop GTPases have been found to be crucial for eukaryotic RNA polymerase II assembly and nuclear trafficking. Despite their ubiquitous occurrence in eukaryotes and archaea, the mechanism by which these GTPases mediate their function is unknown. Our study on an archaeal representative from Sulfolobus acidocaldarius showed that these dimeric GTPases undergo large-scale conformational changes upon GTP hydrolysis, which can be summarized as a lock-switch-rock mechanism. The observed requirement of SaGPN for motility appears to be due to its large footprint on the archaeal proteome.

Glycosylation and Crowded Membrane Effects on Influenza Neuraminidase Stability and Dynamics.

All protein simulations are conducted with varying degrees of simplification, oftentimes with unknown ramifications about how these simplifications affect the interpretability of the results. In this work, we investigated how protein glycosylation and lateral crowding effects modulate an array of properties characterizing the stability and dynamics of influenza neuraminidase. We constructed three systems: (1) glycosylated neuraminidase in a whole virion (i.e., crowded membrane) environment, (2) glycosylated neuraminidase in its own lipid bilayer, and (3) unglycosylated neuraminidase in its own lipid bilayer. We saw that glycans tend to stabilize the protein structure and reduce its conformational flexibility while restricting the solvent movement. Conversely, a crowded membrane environment encouraged exploration of the free energy landscape and a large-scale conformational change, while making the protein structure more compact. Understanding these effects informs what factors one must consider in attempting to recapture the desired level of physical accuracy.

Crystal structures of human and mouse ketohexokinase provide a structural basis for species- and isoform-selective inhibitor design.

A molecular understanding of the proteins involved in fructose metabolism is essential for controlling the current spread of fructose-related obesity, diabetes and related adverse metabolic states in Western populations. Fructose catabolism starts with the phosphorylation of D-fructose to fructose 1-phosphate by ketohexokinase (KHK). KHK exists in two alternatively spliced isoforms: the hepatic and intestinal isoform KHK-C and the peripheral isoform KHK-A. Here, the structure of apo murine KHK (mKHK), which differs from structures of human KHK in overall conformation, is reported. An isoform-selective ligand, which offers a 50-fold higher potency on mKHK and human KHK-A compared with KHK-C, is further characterized. In mKHK, large-scale conformational changes are observed upon ligand binding. The structures suggest a combined strategy for the design of species- and isoform-selective KHK inhibitors.

Redox driven B12-ligand switch drives CarH photoresponse

CarH is a coenzyme B12-dependent photoreceptor involved in regulating carotenoid biosynthesis. How light-triggered cleavage of the B12 Co-C bond culminates in CarH tetramer dissociation to initiate transcription remains unclear. Here, a series of crystal structures of the CarH B12-binding domain after illumination suggest formation of unforeseen intermediate states prior to tetramer dissociation. Unexpectedly, in the absence of oxygen, Co-C bond cleavage is followed by reorientation of the corrin ring and a switch from a lower to upper histidine-Co ligation, corresponding to a pentacoordinate state. Under aerobic conditions, rapid flash-cooling of crystals prior to deterioration upon illumination confirm a similar B12-ligand switch occurs. Removal of the upper His-ligating residue prevents monomer formation upon illumination. Combined with detailed solution spectroscopy and computational studies, these data demonstrate the CarH photoresponse integrates B12 photo- and redox-chemistry to drive large-scale conformational changes through stepwise Co-ligation changes.

Quinone Catalysis Modulates Proton Transfer Reactions in the Membrane Domain of Respiratory Complex I.

Complex I is a redox-driven proton pump that drives electron transport chains and powers oxidative phosphorylation across all domains of life. Yet, despite recently resolved structures from multiple organisms, it still remains unclear how the redox reactions in Complex I trigger proton pumping up to 200 Å away from the active site. Here, we show that the proton-coupled electron transfer reactions during quinone reduction drive long-range conformational changes of conserved loops and trans-membrane (TM) helices in the membrane domain of Complex I from Yarrowia lipolytica. We find that the conformational switching triggers a π → α transition in a TM helix (TM3ND6) and establishes a proton pathway between the quinone chamber and the antiporter-like subunits, responsible for proton pumping. Our large-scale (>20 μs) atomistic molecular dynamics (MD) simulations in combination with quantum/classical (QM/MM) free energy calculations show that the helix transition controls the barrier for proton transfer reactions by wetting transitions and electrostatic effects. The conformational switching is enabled by re-arrangements of ion pairs that propagate from the quinone binding site to the membrane domain via an extended network of conserved residues. We find that these redox-driven changes create a conserved coupling network within the Complex I superfamily, with point mutations leading to drastic activity changes and mitochondrial disorders. On a general level, our findings illustrate how catalysis controls large-scale protein conformational changes and enables ion transport across biological membranes.

Structural basis of the T4 bacteriophage primosome assembly and primer synthesis

The T4 bacteriophage gp41 helicase and gp61 primase assemble into a primosome to couple DNA unwinding with RNA primer synthesis for DNA replication. How the primosome is assembled and how the primer length is defined are unclear. Here we report a series of cryo-EM structures of T4 primosome assembly intermediates. We show that gp41 alone is an open spiral, and ssDNA binding triggers a large-scale scissor-like conformational change that drives the ring closure and activates the helicase. Helicase activation exposes a cryptic hydrophobic surface to recruit the gp61 primase. The primase binds the helicase in a bipartite mode in which the N-terminal Zn-binding domain and the C-terminal RNA polymerase domain each contain a helicase-interacting motif that bind to separate gp41 N-terminal hairpin dimers, leading to the assembly of one primase on the helicase hexamer. Our study reveals the T4 primosome assembly process and sheds light on the RNA primer synthesis mechanism.

Principles of human pre-60S biogenesis.

During the early stages of human large ribosomal subunit (60S) biogenesis, an ensemble of assembly factors establishes and fine-tunes the essential RNA functional centers of pre-60S particles by an unknown mechanism. Here, we report a series of cryo-electron microscopy structures of human nucleolar and nuclear pre-60S assembly intermediates at resolutions of 2.5 to 3.2 angstroms. These structures show how protein interaction hubs tether assembly factor complexes to nucleolar particles and how guanosine triphosphatases and adenosine triphosphatase couple irreversible nucleotide hydrolysis steps to the installation of functional centers. Nuclear stages highlight how a conserved RNA-processing complex, the rixosome, couples large-scale RNA conformational changes with pre-ribosomal RNA processing by the RNA degradation machinery. Our ensemble of human pre-60S particles provides a rich foundation with which to elucidate the molecular principles of ribosome formation.

Insight into the Nucleotide Based Modulation of the Grp94 Molecular Chaperone Using Multiscale Dynamics.

Grp94, an ER-localized molecular chaperone, is required for the folding and activation of many membrane and secretory proteins. Client activation by Grp94 is mediated by nucleotide and conformational changes. In this work, we aim to understand how microscopic changes from nucleotide hydrolysis can potentiate large-scale conformational changes of Grp94. We performed all-atom molecular dynamics simulations on the ATP-hydrolysis competent state of the Grp94 dimer in four different nucleotide bound states. We found that Grp94 was the most rigid when ATP was bound. ATP hydrolysis or nucleotide removal enhanced mobility of the N-terminal domain and ATP lid, resulting in suppression of interdomain communication. In an asymmetric conformation with one hydrolyzed nucleotide, we identified a more compact state, similar to experimental observations. We also identified a potential regulatory role of the flexible linker, as it formed electrostatic interactions with the Grp94 M-domain helix near the region where BiP is known to bind. These studies were complemented with normal-mode analysis of an elastic network model to investigate Grp94's large-scale conformational changes. SPM analysis identified residues that are important in signaling conformational change, many of which have known functional relevance in ATP coordination and catalysis, client binding, and BiP binding. Our findings suggest that ATP hydrolysis in Grp94 alters allosteric wiring and facilitates conformational changes.

Comprehensive Approach to Simulating Large Scale Conformational Changes in Biological Systems Utilizing a Path Collective Variable and New Barrier Restraint.

Conformational sampling of complex biomolecules is an emerging frontier in drug discovery. Advances in lab-based structural biology and related computational approaches like AlphaFold have made great strides in obtaining static protein structures for biologically relevant targets. However, biology is in constant motion, and many important biological processes rely on conformationally driven events. Conventional molecular dynamics (MD) simulations run on standard hardware are impractical for many drug design projects, where conformationally driven biological events can take microseconds to milliseconds or longer. An alternative approach is to focus the search on a limited region of conformational space defined by a putative reaction coordinate (i.e., path collective variable). The search space is typically limited by applying restraints, which can be guided by insights about the underlying biological process of interest. The challenge is striking a balance between the degree to which the system is constrained and still allowing for natural motions along the path. A plethora of restraints exist to limit the size of conformational search space, although each has drawbacks when simulating complex biological motions. In this work, we present a three-stage procedure to construct realistic path collective variables (PCVs) and introduce a new kind of barrier restraint that is particularly well suited for complex conformationally driven biological events, such as allosteric modulations and conformational signaling. The PCV presented here is all-atom (as opposed to C-alpha or backbone only) and is derived from all-atom MD trajectory frames. The new restraint relies on a barrier function (specifically, the scaled reciprocal function), which we show is particularly beneficial in the context of molecular dynamics, where near-hard-wall restraints are needed with zero tolerance to restraint violation. We have implemented our PCV and barrier restraint within a hybrid sampling framework that combines well-tempered metadynamics and extended-Lagrangian adaptive biasing force (meta-eABF). We use three particular examples of high pharmaceutical interest to demonstrate the value of this approach: (1) sampling the distance from ubiquitin to a protein of interest within the supramolecular cullin-RING ligase complex, (2) stabilizing the wild-type conformation of the oncogenic mutant JAK2-V617F pseudokinase domain, and (3) inducing an activated state of the stimulator of interferon genes (STING) protein observed upon ligand binding. For examples 2 and 3, we present statistical analysis of meta-eABF free energy estimates and, for each case, code for reproducing this work.

Solution Structure Ensembles of the Open and Closed Forms of the ∼130 kDa Enzyme I via AlphaFold Modeling, Coarse Grained Simulations, and NMR.

Large-scale interdomain rearrangements are essential to protein function, governing the activity of large enzymes and molecular machineries. Yet, obtaining an atomic-resolution understanding of how the relative domain positioning is affected by external stimuli is a hard task in modern structural biology. Here, we show that combining structural modeling by AlphaFold2 with coarse-grained molecular dynamics simulations and NMR residual dipolar coupling data is sufficient to characterize the spatial domain organization of bacterial enzyme I (EI), a ∼130 kDa multidomain oligomeric protein that undergoes large-scale conformational changes during its catalytic cycle. In particular, we solve conformational ensembles for EI at two different experimental temperatures and demonstrate that a lower temperature favors sampling of the catalytically competent closed state of the enzyme. These results suggest a role for conformational entropy in the activation of EI and demonstrate the ability of our protocol to detect and characterize the effect of external stimuli (such as mutations, ligand binding, and post-translational modifications) on the interdomain organization of multidomain proteins. We expect the ensemble refinement protocol described here to be easily transferrable to the investigation of the structure and dynamics of other uncharted multidomain systems and have assembled a Google Colab page (https://potoyangroup.github.io/Seq2Ensemble/) to facilitate implementation of the presented methodology elsewhere.

Unusual substructure conformations observed in crystal structures of a dicistrovirus RNA-dependent RNA polymerase suggest contribution of the N-terminal extension in proper folding

The Dicistroviridae is a virus family that includes many insect pathogens. These viruses contain a positive-sense RNA genome that is replicated by the virally encoded RNA-dependent RNA polymerase (RdRP) also named 3Dpol. Compared with the Picornaviridae RdRPs such as poliovirus (PV) 3Dpol, the Dicistroviridae representative Israeli acute paralysis virus (IAPV) 3Dpol has an additional N-terminal extension (NE) region that is about 40-residue in length. To date, both the structure and catalytic mechanism of the Dicistroviridae RdRP have remain elusive. Here we reported crystal structures of two truncated forms of IAPV 3Dpol, namely Δ85 and Δ40, both missing the NE region, and the 3Dpol protein in these structures exhibited three conformational states. The palm and thumb domains of these IAPV 3Dpol structures are largely consistent with those of the PV 3Dpol structures. However, in all structures, the RdRP fingers domain is partially disordered, while different conformations of RdRP substructures and interactions between them are also present. In particular, a large-scale conformational change occurred in the motif B-middle finger region in one protein chain of the Δ40 structure, while a previously documented alternative conformation of motif A was observed in all IAPV structures. These experimental data on one hand show intrinsic conformational variances of RdRP substructures, and on the other hand suggest possible contribution of the NE region in proper RdRP folding in IAPV.

Role of the Coiled-Coil Domain in Allosteric Activity Regulation in Soluble Guanylate Cyclase.

Soluble guanylate cyclase (sGC) is the primary nitric oxide (NO) receptor in higher eukaryotes, including humans. NO-dependent signaling via sGC is associated with important physiological effects in the vascular, pulmonary, and neurological systems, and sGC itself is an established drug target for the treatment of pulmonary hypertension due to its central role in vasodilation. Despite isolation in the late 1970s, high-resolution structural information on full-length sGC remained elusive until recent cryo-electron microscopy structures were determined of the protein in both the basal unactivated state and the NO-activated state. These structures revealed large-scale conformational changes upon activation that appear to be centered on rearrangements within the coiled-coil (CC) domains in the enzyme. Here, a structure-guided approach was used to engineer constitutively unactivated and constitutively activated sGC variants through mutagenesis of the CC domains. These results demonstrate that the activation-induced conformational change in the CC domains is necessary and sufficient for determining the level of sGC activity.

A conserved oligomerization domain in the disordered linker of coronavirus nucleocapsid proteins.

The nucleocapsid (N-)protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a key role in viral assembly and scaffolding of the viral RNA. It promotes liquid-liquid phase separation (LLPS), forming dense droplets that support the assembly of ribonucleoprotein particles with as-of-yet unknown macromolecular architecture. Combining biophysical experiments, molecular dynamics simulations, and analysis of the mutational landscape, we describe a heretofore unknown oligomerization site that contributes to LLPS, is required for the assembly of higher-order protein-nucleic acid complexes, and is coupled to large-scale conformational changes of N-protein upon nucleic acid binding. The self-association interface is located in a leucine-rich sequence of the intrinsically disordered linker between N-protein folded domains and formed by transient helices assembling into trimeric coiled-coils. Critical residues stabilizing hydrophobic and electrostatic interactions between adjacent helices are highly protected against mutations in viable SARS-CoV-2 genomes, and the oligomerization motif is conserved across related coronaviruses, thus presenting a target for antiviral therapeutics.

Highly connected 3D chromatin networks established by an oncogenic fusion protein shape tumor cell identity

Cell fate transitions observed in embryonic development involve changes in three-dimensional genomic organization that provide proper lineage specification. Whether similar events occur within tumor cells and contribute to cancer evolution remains largely unexplored. We modeled this process in the pediatric cancer Ewing sarcoma and investigated high-resolution looping and large-scale nuclear conformation changes associated with the oncogenic fusion protein EWS-FLI1. We show that chromatin interactions in tumor cells are dominated by highly connected looping hubs centered on EWS-FLI1 binding sites, which directly control the activity of linked enhancers and promoters to establish oncogenic expression programs. Conversely, EWS-FLI1 depletion led to the disassembly of these looping networks and a widespread nuclear reorganization through the establishment of new looping patterns and large-scale compartment configuration matching those observed in mesenchymal stem cells, a candidate Ewing sarcoma progenitor. Our data demonstrate that major architectural features of nuclear organization in cancer cells can depend on single oncogenes and are readily reversed to reestablish latent differentiation programs.

Allosteric communication in the gating mechanism for controlled protein degradation by the bacterial ClpP peptidase.

Proteolysis is essential for the control of metabolic pathways and the cell cycle. Bacterial caseinolytic proteases (Clp) use peptidase components, such as ClpP, to degrade defective substrate proteins and to regulate cellular levels of stress-response proteins. To ensure selective degradation, access to the proteolytic chamber of the double-ring ClpP tetradecamer is controlled by a critical gating mechanism of the two axial pores. The binding of conserved loops of the Clp ATPase component of the protease or small molecules, such as acyldepsipeptide (ADEP), at peripheral ClpP ring sites, triggers axial pore opening through dramatic conformational transitions of flexible N-terminal loops between disordered conformations in the "closed" pore state and ordered hairpins in the "open" pore state. In this study, we probe the allosteric communication underlying these conformational changes by comparing residue-residue couplings in molecular dynamics simulations of each configuration. Both principal component and normal mode analyses highlight large-scale conformational changes in the N-terminal loop regions and smaller amplitude motions of the peptidase core. Community network analysis reveals a switch between intra- and inter-protomer coupling in the open-closed pore transition. Allosteric pathways that connect the ADEP binding sites to N-terminal loops are rewired in this transition, with shorter network paths in the open pore configuration supporting stronger intra- and inter-ring coupling. Structural perturbations, either through the removal of ADEP molecules or point mutations, alter the allosteric network to weaken the coupling.

Indoline CD4-mimetic compounds mediate potent and broad HIV-1 inhibition and sensitization to antibody-dependent cellular cytotoxicity

Binding to the host cell receptors, CD4 and CCR5/CXCR4, triggers large-scale conformational changes in the HIV-1 envelope glycoprotein (Env) trimer [(gp120/gp41)3] that promote virus entry into the cell. CD4-mimetic compounds (CD4mcs) comprise small organic molecules that bind in the highly conserved CD4-binding site of gp120 and prematurely induce inactivating Env conformational changes, including shedding of gp120 from the Env trimer. By inducing more "open," antibody-susceptible Env conformations, CD4mcs also sensitize HIV-1 virions to neutralization by antibodies and infected cells to antibody-dependent cellular cytotoxicity (ADCC). Here, we report the design, synthesis, and evaluation of novel CD4mcs based on an indoline scaffold. Compared with our current lead indane scaffold CD4mc, BNM-III-170, several indoline CD4mcs exhibit increased potency and breadth against HIV-1 variants from different geographic clades. Viruses that were selected for resistance to the lead indane CD4mc, BNM-III-170, are susceptible to inhibition by the indoline CD4mcs. The indoline CD4mcs also potently sensitize HIV-1-infected cells to ADCC mediated by plasma from HIV-1-infected individuals. Crystal structures indicate that the indoline CD4mcs gain potency compared to the indane CD4mcs through more favorable π-π overlap from the indoline pose and by making favorable contacts with the vestibule of the CD4-binding pocket on gp120. The rational design of indoline CD4mcs thus holds promise for further improvements in antiviral activity, potentially contributing to efforts to treat and prevent HIV-1 infection.