Conformational Changes Research Articles

Mapping the Intersubunit Interdomain FMN-Heme Interactions in Neuronal Nitric Oxide Synthase by Targeted Quantitative Cross-Linking Mass Spectrometry.

Nitric oxide synthase (NOS) in mammals is a family of multidomain proteins in which interdomain electron transfer (IET) is controlled by domain-domain interactions. Calmodulin (CaM) binds to the canonical CaM-binding site in the linker region between the FMN and heme domains of NOS and allows tethered FMN domain motions, enabling an intersubunit FMN-heme IET in the output state for NO production. Our previous cross-linking mass spectrometric (XL MS) results demonstrated site-specific protein dynamics in the CaM-responsive regions of rat neuronal NOS (nNOS) reductase construct, a monomeric protein [Jiang et al., Biochemistry, 2023, 62, 2232-2237]. In this work, we have extended our combined approach of XL MS structural mapping and AlphaFold structural prediction to examine the homodimeric nNOS oxygenase/FMN (oxyFMN) construct, an established model of the NOS output state. We employed parallel reaction monitoring (PRM) based quantitative XL MS (qXL MS) to assess the CaM-induced changes in interdomain dynamics and interactions. Intersubunit cross-links were identified by mapping the cross-links onto top AlphaFold structural models, which was complemented by comparing their relative abundances in the cross-linked dimeric and monomeric bands. Furthermore, contrasting the CaM-free and CaM-bound nNOS samples shows that CaM enables the formation of the intersubunit FMN-heme docking complex and that CaM binding induces extensive, allosteric conformational changes across the NOS regions. Moreover, the observed cross-links sites specifically respond to changes in ionic strength. This indicates that interdomain salt bridges are responsible for stabilizing and orienting the output state for efficient FMN-heme IET. Taken together, our targeted qXL MS results have revealed that CaM and ionic strength modulate specific dynamic changes in the CaM/FMN/heme complexes, particularly in the context of intersubunit interdomain FMN-heme interactions.

Molecular Insights into the Effects of F16L and F19L Substitutions on the Conformation and Aggregation Dynamics of Human Calcitonin.

Human calcitonin (hCT) regulates calcium-phosphorus metabolism, but its amyloid aggregation disrupts physiological activity, increases thyroid carcinoma risk, and hampers its clinical use for bone-related diseases like osteoporosis and Paget's disease. Improving hCT with targeted modifications to mitigate amyloid formation while maintaining its function holds promise as a strategy. Understanding how each residue in hCT's amyloidogenic core affects its structure and aggregation dynamics is crucial for designing effective analogues. Mutants F16L-hCT and F19L-hCT, where Phe residues in the core are replaced with Leu as in nonamyloidogenic salmon calcitonin, showed different aggregation kinetics. However, the molecular effects of these substitutions in hCT are still unclear. Here, we systematically investigated the folding and self-assembly conformational dynamics of hCT, F16L-hCT, and F19L-hCT through multiple long-time scale independent atomistic discrete molecular dynamics (DMD) simulations. Our results indicated that the hCT monomer primarily assumed unstructured conformations with dynamic helices around residues 4-12 and 14-21. During self-assembly, the amyloidogenic core of hCT14-21 converted from dynamic helices to β-sheets. However, substituting F16L did not induce significant conformational changes, as F16L-hCT exhibited characteristics similar to those of wild-type hCT in both monomeric and oligomeric states. In contrast, F19L-hCT exhibited substantially more helices and fewer β-sheets than did hCT, irrespective of their monomers or oligomers. The substitution of F19L significantly enhanced the stability of the helical conformation for hCT14-21, thereby suppressing the helix-to-β-sheet conformational conversion. Overall, our findings elucidate the molecular mechanisms underlying hCT aggregation and the effects of F16L and F19L substitutions on the conformational dynamics of hCT, highlighting the critical role of F19 as an important target in the design of amyloid-resistant hCT analogs for future clinical applications.

Adaptive Workflows of Machine Learning Illuminate the Sequential Operation Mechanism of the TAK1's Allosteric Network.

Allostery is a fundamental mechanism driving biomolecular processes that holds significant therapeutic concern. Our study rigorously investigates how two distinct machine-learning algorithms uniquely classify two already close-to-active DFG-in states of TAK1, differing just by the presence or absence of its allosteric activator TAB1, from an ensemble mixture of conformations (obtained from 2.4 μs molecular dynamics (MD) simulations). The novelty, however, lies in understanding the deeper algorithmic potentials to systematically derive a diverse set of differential residue connectivity features that reconstruct the essential mechanistic architecture for TAK1-TAB1 allostery in such a close-to-active biochemical scenario. While the recursive, random forest-based workflow displays the potential of conducting discretized, hierarchical derivation of allosteric features, a multilayer perceptron-based approach gains considerable efficacy in revealing fluid connected patterns of features when hybridized with mutual information scoring. Interestingly, both pipelines benchmark similar directions of functional conformational changes for TAK1's activation. The findings significantly advance the depth of mechanistic understanding by highlighting crucial activation signatures along a directed C-lobe → activation loop → ATP pocket channel of information flow, including (1) the αF-αE biterminal alignments and (2) the "catalytic" drift of the activation loop toward kinase active site. Besides, some novel allosteric hotspots (K253, Y206, N189, etc.) are further recognized as TAB1 sensors, transducers, and responders, including a benchmark E70 mutation site, precisely mapping the important structural segments for sequential allosteric execution. Hence, our work demonstrates how to navigate through greater structural depths and dimensions of dynamic allosteric machineries just by leveraging standard ML methods in suitable streamlined workflows adaptive to the specific system and objectives.

Chemical manipulation of an activation/inhibition switch in the nuclear receptor PXR

Nuclear receptors are ligand-activated transcription factors that can often be useful drug targets. Unfortunately, ligand promiscuity leads to two-thirds of receptors remaining clinically untargeted. PXR is a nuclear receptor that can be activated by diverse compounds to elevate metabolism, negatively impacting drug efficacy and safety. This presents a barrier to drug development because compounds designed to target other proteins must avoid PXR activation while retaining potency for the desired target. This problem could be avoided by using PXR antagonists, but these compounds are rare, and their molecular mechanisms remain unknown. Here, we report structurally related PXR-selective agonists and antagonists and their corresponding co-crystal structures to describe mechanisms of antagonism and selectivity. Structural and computational approaches show that antagonists induce PXR conformational changes incompatible with transcriptional coactivator recruitment. These results guide the design of compounds with predictable agonist/antagonist activities and bolster efforts to generate antagonists to prevent PXR activation interfering with other drugs.

Structural basis for double-stranded RNA recognition by SID1.

The nucleic acid transport properties of the systemic RNAi-defective (SID) 1 family make them attractive targets for developing RNA-based therapeutics and drugs. However, the molecular basis for double-stranded (ds) RNA recognition by SID1 family remains elusive. Here, we report the cryo-EM structures of Caenorhabditis elegans (c) SID1 alone and in complex with dsRNA, both at a resolution of 2.2 Å. The dimeric cSID1 interacts with two dsRNA molecules simultaneously. The dsRNA is located at the interface between β-strand rich domain (BRD)1 and BRD2 and nearly parallel to the membrane plane. In addition to extensive ionic interactions between basic residues and phosphate backbone, several hydrogen bonds are formed between 2'-hydroxyl group of dsRNA and the contact residues. Additionally, the electrostatic potential surface shows three basic regions are fitted perfectly into three major grooves of dsRNA. These structural characteristics enable cSID1 to bind dsRNA in a sequence-independent manner and to distinguish between DNA and RNA. The cSID1 exhibits no conformational changes upon binding dsRNA, with the exception of a few binding surfaces. Structural mapping of dozens of loss-of-function mutations allows potential interpretation of their diverse functional mechanisms. Our study marks an important step toward mechanistic understanding of the SID1 family-mediated dsRNA uptake.

(Thermo)mechanical and chemical characteristics of photochemically crosslinked acrylates from bio‐ and fossil‐based origin

AbstractThe acrylates with oligomers and monomers from (partially) bio‐based feedstock become available at (semi‐)industrial scale, which can be processed through photochemical crosslinking for applications in coatings, additive manufacturing, electronics, or inks. Although fossil‐ and bio‐based acrylates may have a similar chemical composition, it requires good understanding of processing and structure–property relationships as minor changes in microstructure may strongly alter the performance. A comparative study on mechanical properties and chemical structure of bio‐ and fossil‐based acrylates with different functionalities and backbone structures reveals higher ductility of bio‐based acrylates, in relation with a more complex organization of the intrinsic molecular structure. The latter is confirmed by mechanical testing and visco‐elastic characteristics (dynamic mechanical analysis) yielding lower stiffness and higher dampening of bio‐based acrylates, in parallel with a lower glass transition temperature (differential scanning calorimetry). The complex molecular arrangements include a nanoscale morphology with ordered structure (X‐ray diffraction), conformational changes (infrared spectroscopy), and a residual high‐molecular weight fraction (size exclusion chromatography). The visco‐elastic calculations indicate only 4% to 5% lower crosslinking density and around 10% higher mean molar mass of the polymer chains segments between chemical crosslinks and trapped chain entanglements, which explain the unique structure and performance of bio‐based acrylates.

Mechanical Stability and Unfolding Pathways of Parallel Tetrameric G-Quadruplexes Probed by Pulling Simulations.

Guanine quadruplex (GQ) is a noncanonical nucleic acid structure formed by guanine-rich DNA and RNA sequences. Folding of GQs is a complex process, where several aspects remain elusive, despite being important for understanding structure formation and biological functions of GQs. Pulling experiments are a common tool for acquiring insights into the folding landscape of GQs. Herein, we applied a computational pulling strategy─steered molecular dynamics (SMD) simulations─in combination with standard molecular dynamics (MD) simulations to explore the unfolding landscapes of tetrameric parallel GQs. We identified anisotropic properties of elastic conformational changes, unfolding transitions, and GQ mechanical stabilities. Using a special set of structural parameters, we found that the vertical component of pulling force (perpendicular to the average G-quartet plane) plays a significant role in disrupting GQ structures and weakening their mechanical stabilities. We demonstrated that the magnitude of the vertical force component depends on the pulling anchor positions and the number of G-quartets. Typical unfolding transitions for tetrameric parallel GQs involve base unzipping, opening of the G-stem, strand slippage, and rotation to cross-like structures. The unzipping was detected as the first and dominant unfolding event, and it usually started at the 3'-end. Furthermore, results from both SMD and standard MD simulations indicate that partial spiral conformations serve as a transient ensemble during the (un)folding of GQs.

Folding Kinetics and Volume Variation of the β-Hairpin Peptide Chignolin upon Ultrafast pH-Jumps.

In-depth characterization of fundamental folding steps of small model peptides is crucial for a better understanding of the folding mechanisms of more complex biomacromolecules. We have previously reported on the folding/unfolding kinetics of a model α-helix. Here, we study folding transitions in chignolin (GYDPETGTWG), a short β-hairpin peptide previously used as a model to study conformational changes in β-sheet proteins. Although previously suggested, until now, the role of the Tyr2-Trp9 interaction in the folding mechanism of chignolin was not clear. In the present work, pH-dependent conformational changes of chignolin were characterized by circular dichroism (CD), nuclear magnetic resonance (NMR), ultrafast pH-jump coupled with time-resolved photoacoustic calorimetry (TR-PAC), and molecular dynamics (MD) simulations. Taken together, our results present a comprehensive view of chignolin's folding kinetics upon local pH changes and the role of the Tyr2-Trp9 interaction in the folding process. CD data show that chignolin's β-hairpin formation displays a pH-dependent skew bell-shaped curve, with a maximum close to pH 6, and a large decrease in β-sheet content at alkaline pH. The β-hairpin structure is mainly stabilized by aromatic interactions between Tyr2 and Trp9 and CH-π interactions between Tyr2 and Pro4. Unfolding of chignolin at high pH demonstrates that protonation of Tyr2 is essential for the stability of the β-hairpin. Refolding studies were triggered by laser-induced pH-jumps and detected by TR-PAC. The refolding of chignolin from high pH, mainly due to the protonation of Tyr2, is characterized by a volume expansion (10.4 mL mol-1), independent of peptide concentration, in the microsecond time range (lifetime of 1.15 μs). At high pH, the presence of the deprotonated hydroxyl (tyrosinate) hinders the formation of the aromatic interaction between Tyr2 and Trp9 resulting in a more disorganized and dynamic tridimensional structure of the peptide. This was also confirmed by comparing MD simulations of chignolin under conditions mimicking neutral and high pH.

Normal-to-Supercooled Liquid Transition in Molecular Glass-Formers: A Hidden Structural Transformation Fuelled by Conformational Interconversion.

Molecular dynamics and transport coefficients change significantly around the so-called Arrhenius crossover in glass-forming systems. In this article, we revisit the dynamic processes occurring in a glass-forming macrocyclic crown thiaether MeBzS2O above its glass transition, revealing two crossover temperatures: TB at 309 and TA at 333 K. We identify the second one as the Arrhenius crossover that is closely related to the normal-to-supercooled liquid transition in this compound. We show that the transformation occurring at this point goes far beyond molecular dynamics (where the temperature dependence of structural relaxation times changes its character from activation-like to super-Arrhenius), being reflected also in the internal structure and diffraction pattern. In this respect, we found a twofold local organization of the nearest-neighbor molecules via weak van der Waals forces, without the formation of any medium-range order or mesophases. The nearest surrounding of each molecule evolves structurally in time due to the ongoing fast conformational changes. We identify several conformers of MeBzS2O, demonstrating that its lowest-energy conformation is preferred mainly at lower temperatures, i.e., in the supercooled liquid state. Its increased prevalence modifies locally the short-range intermolecular order and promotes vitrification. Consequently, we indicate that the Arrhenius transition is fuelled rather by conformational changes in this glass-forming macrocyclic crown thiaether, which is a different scenario from the so-far existing concepts. Our studies combine broadband dielectric spectroscopy (BDS), X-ray diffraction, Fourier transform infrared (FTIR) spectroscopy, molecular dynamics (MD) simulations, and density functional theory (DFT) calculations.

In vitro susceptibility of grass carp fish (Ctenopharyngodon idella) to oxidized myosin on structural change of protein secondary structure associated with reactive oxygen species (ROS) treatment

This study aimed to investigate the in vitro susceptibility of grass carp (Ctenopharyngodon idella) oxidized myosin, with the treatment of H2O2 levels being 0, 25, 50, and 100 μmol/L, respectively. The results showed that in vitro oxidation of myosin, reactive oxidative species promoted protein modifications that propagated the conformation changes on myosin secondary structure, leading to a significant increase (P < 0.05) in the protein carbonyl content in the high oxidation treatment, when compared to the low and medium treatments. The results of the Western blot analysis and sodium-polyacrylamide amide-dodecyl sulfate gel electrophoresis demonstrated that in vitro oxidation actively promoted the degradation of protein molecular weight and Desmin (53 kDa) deterioration in the sample with the highest concentration treatment of H2O2. Circular dichroism spectra showed a decrease in α-helix content from 70.1% to 29.8% and an increase in β-sheet percentage from 5.6% to 17.4%. Meanwhile, cysteine and methionine as the free amino acids were most susceptible to oxidative change due to their high sensitivity to the reaction of the sulfur group of amino acids. Protein degradation induced by oxidation resulted in significant changes to the protein secondary structure, as evidenced by the observation of oxidized myosin under different H2O2 levels.

Unveiling the transition of tumbling and tank-treading dynamics of star-shaped polyelectrolytes in shear flow by modulating the solution's dielectric properties.

The effects of the solution's dielectric properties on the conformation and dynamics of star-shaped polyelectrolytes in shear flow are investigated using a hybrid simulation method coupling multi-particle collision dynamics and molecular dynamics. The simulation results showed that by modulating the dielectric properties of the solution, star-shaped polyelectrolytes showed a three-step dynamic behavior transition from tumbling to tank-treading to tumbling dynamics under shear flow. The analysis indicated that this distinct transition in dynamics could be attributed to the uneven distribution of counterions induced by shear on the chain, resulting in a change in the polyelectrolyte conformation and degree of segmental alignment in arms. These findings contribute to a comprehensive understanding of the non-equilibrium dynamics of star-shaped polyelectrolytes in shear flow and offer a viable approach for controlling the dynamic behavior of star-shaped polyelectrolytes by adjusting the dielectric properties of the solution.

Stimulus-Induced Dynamic Behavior Regulation of Flexible Crystals through the Tuning of Module Rigidity.

Introducing dynamic behavior into periodic frameworks has borne fruit in the form of flexible porous crystals. The detailed molecular design of frameworks in order to control their collective dynamics is of particular interest, for example, to achieve stimulus-induced behavior. Herein, by varying the degree of rigidity of ditopic pillar linkers, two isostructural flexible metal-organic frameworks (MOFs) with common rigid supermolecular building bilayers were constructed. The subtle substitution of single (in bibenzyl-4,4'-dicarboxylic acid; H2BBDC) with double (in 4,4'-stilbenedicarboxylic acid; H2SDC) C-C bonds in pillared linkers led to markedly different flexible behavior of these two MOFs. Upon the removal of guest molecules, both frameworks clearly show reversible single-crystal-to-single-crystal transformations involving the cis-trans conformation change and a resulting swing of the corresponding pillar linkers, which gives rise to Flex-Cd-MOF-1a and Flex-Cd-MOF-2a, respectively. Strikingly, a more favorable gas-induced dynamic behavior in Flex-Cd-MOF-2a was verified in detail by stepwise C3H6/C3H8 sorption isotherms and the corresponding in situ powder X-ray diffraction experiments. These insights are strongly supported by molecular modeling studies on the sorption mechanism that explores the sorption landscape. Furthermore, a consistency between the macroscopic elasticity and microscopic flexibility of Flex-Cd-MOF-2 was observed. This work fuels a growing interest in developing MOFs with desired chemomechanical functions and presents detailed insights into the origins of flexible MOFs.

Impact of protein conformations on binding free energy calculations in the beta-secretase 1 system.

In binding free energy calculations, simulations must sample all relevant conformations of the system in order to obtain unbiased results. For instance, different ligands can bind to different metastable states of a protein, and if these protein conformational changes are not sampled in relative binding free energy calculations, the contribution of these states to binding is not accounted for and thus calculated binding free energies are inaccurate. In this work, we investigate the impact of different beta-sectretase 1 (BACE1) protein conformations obtained from x-ray crystallography on the binding of BACE1 inhibitors. We highlight how these conformational changes are not adequately sampled in typical molecular dynamics simulations. Furthermore, we show that insufficient sampling of relevant conformations induces substantial error in relative binding free energy calculations, as judged by a variation in calculated relative binding free energies up to 2kcal/mol depending on the starting protein conformation. These results emphasize the importance of protein conformational sampling and pose this BACE1 system as a challenge case for further method development in the area of enhanced protein conformational sampling, either in combination with binding calculations or as an endpoint correction.

Molecular Dynamics Simulations of the Mutated Proton-Transferring a-Subunit of E. coli FoF1-ATP Synthase

The membrane Fo factor of ATP synthase is highly sensitive to mutations in the proton half-channel leading to the functional blocking of the entire protein. To identify functionally important amino acids for the proton transport, we performed molecular dynamic simulations on the selected mutants of the membrane part of the bacterial FoF1-ATP synthase embedded in a native lipid bilayer: there were nine different mutations of a-subunit residues (aE219, aH245, aN214, aQ252) in the inlet half-channel. The structure proved to be stable to these mutations, although some of them (aH245Y and aQ252L) resulted in minor conformational changes. aH245 and aN214 were crucial for proton transport as they directly facilitated H+ transfer. The substitutions with nonpolar amino acids disrupted the transfer chain and water molecules or neighboring polar side chains could not replace them effectively. aE219 and aQ252 appeared not to be determinative for proton translocation, since an alternative pathway involving a chain of water molecules could compensate the ability of H+ transmembrane movement when they were substituted. Thus, mutations of conserved polar residues significantly affected hydration levels, leading to drastic changes in the occupancy and capacity of the structural water molecule clusters (W1–W3), up to their complete disappearance and consequently to the proton transfer chain disruption.

Skipping events impose repeated binding attempts: profound kinetic implications of protein-DNA conformational changes.

The kinetics of protein-DNA recognition, along with its thermodynamic properties, including affinity and specificity, play a central role in shaping biological function. Protein-DNA recognition kinetics are characterized by two key elements: the time taken to locate the target site amid various nonspecific alternatives; and the kinetics involved in the recognition process, which may necessitate overcoming an energetic barrier. In this study, we developed a coarse-grained (CG) model to investigate interactions between a transcription factor called the sex-determining region Y (SRY) protein and DNA, in order to probe how DNA conformational changes affect SRY-DNA recognition and binding kinetics. We find that, not only does a requirement for such a conformational DNA transition correspond to a higher energetic barrier for binding and therefore slower kinetics, it may further impede the recognition kinetics by increasing unsuccessful binding events (skipping events) where the protein partially binds its DNA target site but fails to form the specific protein-DNA complex. Such skipping events impose the need for additional cycles protein search of nonspecific DNA sites, thus significantly extending the overall recognition time. Our results highlight a trade-off between the speed with which the protein scans nonspecific DNA and the rate at which the protein recognizes its specific target site. Finally, we examine molecular approaches potentially adopted by natural systems to enhance protein-DNA recognition despite its intrinsically slow kinetics.

Mechanistic basis of the dynamic response of TWIK1 ionic selectivity to pH

Highly selective for K+ at neutral pH, the TWIK1 channel becomes permeable to Na+ upon acidification. Using molecular dynamics simulations, we identify a network of residues involved in this unique property. Between the open and closed states previously observed by electron microscopy, molecular dynamics simulations show that the channel undergoes conformational changes between pH 7.5–6 involving residues His122, Glu235, Lys246 and Phe109. A complex network of interactions surrounding the selectivity filter at high pH transforms into a simple set of stronger interactions at low pH. In particular, His122 protonated by acidification moves away from Lys246 and engages in a salt bridge with Glu235. In addition, stacking interactions between Phe109 and His122, which stabilize the selectivity filter in its K+-selective state at high pH, disappear upon acidification. This leads to dissociation of the Phe109 aromatic side chain from this network, resulting in the Na+-permeable conformation of the channel.

Hybrid Small-Molecule/Protein Fluorescent Probes.

Hybrid small-molecule/protein fluorescent probes are powerful tools for visualizing protein localization and function in living cells. These hybrid probes are constructed by diverse site-specific chemical protein labeling approaches through chemical reactions to exogenous peptide/small protein tags, enzymatic post-translational modifications, bioorthogonal reactions for genetically incorporated unnatural amino acids, and ligand-directed chemical reactions. The hybrid small-molecule/protein fluorescent probes are employed for imaging protein trafficking, conformational changes, and bioanalytes surrounding proteins. In addition, fluorescent hybrid probes facilitate visualization of protein dynamics at the single-molecule level and the defined structure with super-resolution imaging. In this review, we discuss development and the bioimaging applications of fluorescent probes based on small-molecule/protein hybrids.

Nonalternant Nanographenes Containing N-Centered Cyclopenta[ef]heptalene and Aza[7]Helicene Units.

Introducing helical subunits into negatively curved π-systems has a significant effect on both the molecular geometry and photophysical properties; however, the synthesis of these helical π-systems embedded with nonbenzenoid subunits remains challenging due to the high strain deriving from both the curvature and helix. Here, we report a family of nonalternant nanographenes containing a nitrogen (N)-doped cyclopenta[ef]heptalene unit. Among them, CPH-2 and CPH-3 can be viewed as hybrids of benzoannulated cyclopenta[ef]heptalene and aza[7]helicene. The crystal structures revealed a saddle geometry for CPH-1, a saddle-helix hybrid for CPH-2, and a twist-helix hybrid for CPH-3. Experimental measurements and theoretical calculations indicate that the saddle moieties in CPHs undergo flexible conformational changes at room temperature, while the aza[7]helicene subunit exhibits a dramatically increased racemization energy barrier (78.2 kcal mol-1 for CPH-2, 143.2 kcal mol-1 for CPH-3). The combination of the nitrogen lone electron pairs of the N-doped cyclopenta[ef]heptalene unit with the twisted helix fragments results in rich photophysics with distinctive fluorescence and phosphorescence in CPH-1 and CPH-2 and the similar energy fluorescence and phosphorescence in CPH-3. Both enantiopure CPH-2 and CPH-3 display distinct circular dichroism (CD) signals in the UV-vis range. Notably, compared to the reported fully π-extended helical nanographenes, CPH-3 exhibits excellent chiroptical properties with a |gabs| value of 1.0 × 10-2 and a |glum| value of 7.0 × 10-3; these values are among the highest for helical nanographenes.

Significance of Zn2+ in RyR1 for Structural Integrity and Ligand Binding: Insight from Molecular Dynamics.

Ryanodine receptor type 1 (RyR1) is a Ca2+-release channel central to skeletal muscle excitation-contraction (EC) coupling. RyR1's cryo-EM structures reveal a zinc-finger motif positioned within the cytoplasmic C-terminal domain (CTD). Yet, owing to limitations in cryo-EM resolution, RyR1 structures lack precision in detailing the metal coordination structure, prompting the need for an accurate model. In this study, we employed molecular dynamics (MD) simulations and the density functional theory (DFT) method to refine the binding characteristics of Zn2+ in the zinc-finger site of the RyR1 channel. Our findings also highlight substantial conformational changes in simulations conducted in the absence of Zn2+. Notably, we observed a loss of contact at the interface between protein domains proximal to the zinc-finger site, indicating a crucial role of Zn2+ in maintaining structural integrity and interdomain interactions within RyR1. Furthermore, this study provides valuable insights into the modulation of ATP, Ca2+, and caffeine binding, shedding light on the intricate relationship between Zn2+ coordination and the dynamic behavior of RyR1. Our integrative approach combining MD simulations and DFT calculations enhances our understanding of the molecular mechanisms governing ligand binding in RyR1.

Cadmium binding by the F-box domain induces p97-mediated SCF complex disassembly to activate stress response programs

The F-box domain is a highly conserved structural motif that defines the largest class of ubiquitin ligases, Skp1/Cullin1/F-box protein (SCF) complexes. The only known function of the F-box motif is to form the protein interaction surface with Skp1. Here we show that the F-box domain can function as an environmental sensor. We demonstrate that the F-box domain of Met30 is a cadmium sensor that blocks the activity of the SCFMet30 ubiquitin ligase during cadmium stress. Several highly conserved cysteine residues within the Met30 F-box contribute to binding of cadmium with a KD of 8 µM. Binding induces a conformational change that allows for Met30 autoubiquitylation, which in turn leads to recruitment of the segregase Cdc48/p97/VCP followed by active SCFMet30 disassembly. The resulting inactivation of SCFMet30 protects cells from cadmium stress. Our results show that F-box domains participate in regulation of SCF ligases beyond formation of the Skp1 binding interface.