Illuminating ligand-induced dynamics in nuclear receptors through MD simulations

Nuclear receptors (NRs) regulate gene expression in critical physiological processes, with their functionality finely tuned by ligand-induced conformational changes. While NRs may sometimes undergo significant conformational motions in response to ligand-binding, these effects are more commonly subtle and challenging to study by traditional structural or biophysical methods. Molecular dynamics (MD) simulations are a powerful tool to bridge the gap between static protein-ligand structures and dynamical changes that govern NR function. Here, we summarize a handful of recent studies that apply MD simulations to study NRs. We present diverse methodologies for analyzing simulation data with a detailed examination of the information each method can yield. By delving into the strengths, limitations and unique contributions of these tools, this review provides guidance for extracting meaningful data from MD simulations to advance the goal of understanding the intricate mechanisms by which ligands orchestrate a range of functional outcomes in NRs.

Four Parallel Pathways in T4 Ligase-Catalyzed Repair of Nicked DNA with Diverse Bending Angles

The structural diversity of biological macromolecules in different environments contributes complexity to enzymological processes vital for cellular functions. Fluorescence resonance energy transfer and electron microscopy are used to investigate the enzymatic reaction of T4 DNA ligase catalyzing the ligation of nicked DNA. The data show that both the ligase-AMP complex and the ligase-AMP-DNA complex can have four conformations. This finding suggests the parallel occurrence of four ligation reaction pathways, each characterized by specific conformations of the ligase-AMP complex that persist in the ligase-AMP-DNA complex. Notably, these complexes have DNA bending angles of ≈0°, 20°, 60°, or 100°. The mechanism of parallel reactions challenges the conventional notion of simple sequential reaction steps occurring among multiple conformations. The results provide insights into the dynamic conformational changes and the versatile attributes of T4 DNA ligase and suggest that the parallel multiple reaction pathways may correspond to diverse T4 DNA ligase functions. This mechanism may potentially have evolved as an adaptive strategy across evolutionary history to navigate complex environments.

Peptide-based allosteric inhibitor targets TNFR1 conformationally active region and disables receptor-ligand signaling complex

Tumor necrosis factor (TNF) receptor 1 (TNFR1) plays a pivotal role in mediating TNF induced downstream signaling and regulating inflammatory response. Recent studies have suggested that TNFR1 activation involves conformational rearrangements of preligand assembled receptor dimers and targeting receptor conformational dynamics is a viable strategy to modulate TNFR1 signaling. Here, we used a combination of biophysical, biochemical, and cellular assays, as well as molecular dynamics simulation to show that an anti-inflammatory peptide (FKCRRWQWRMKK), which we termed FKC, inhibits TNFR1 activation allosterically by altering the conformational states of the receptor dimer without blocking receptor-ligand interaction or disrupting receptor dimerization. We also demonstrated the efficacy of FKC by showing that the peptide inhibits TNFR1 signaling in HEK293 cells and attenuates inflammation in mice with intraperitoneal TNF injection. Mechanistically, we found that FKC binds to TNFR1 cysteine-rich domains (CRD2/3) and perturbs the conformational dynamics required for receptor activation. Importantly, FKC increases the frequency in the opening of both CRD2/3 and CRD4 in the receptor dimer, as well as induces a conformational opening in the cytosolic regions of the receptor. This results in an inhibitory conformational state that impedes the recruitment of downstream signaling molecules. Together, these data provide evidence on the feasibility of targeting TNFR1 conformationally active region and open new avenues for receptor-specific inhibition of TNFR1 signaling.

Molecular dynamics simulations reveal phosphorylation-induced conformational dynamics of the fibroblast growth factor receptor 1 kinase

The Fibroblast Growth Factor Receptor1 (FGFR1) kinase wields exquisite control on cell fate, proliferation, differentiation, and homeostasis. An imbalance of FGFR1 signaling leads to several pathogeneses of diseases ranging from multiple cancers to allergic and neurodegenerative disorders. In this study, we investigated the phosphorylation-induced conformational dynamics of FGFR1 in apo and ATP-bound states via all-atom molecular dynamics simulations. All simulations were performed for 2 × 2 µs. We have also investigated the energetics of the binding of ATP to FGFR1 using the molecular mechanics Poisson-Boltzmann scheme. Our study reveals that the FGFR1 kinase can reach a fully active configuration through phosphorylation and ATP binding. A 3-10 helix formation in the activation loop signifies its rearrangement leading to stability upon ATP binding. The interaction of phosphorylated tyrosine (pTyr654) with positively charged residues forms strong salt-bridge interactions, driving the compactness of the structure. The dynamic cross-correlation map reveals phosphorylation enhances correlated motions and reduces anti-correlated motions between different domains. We believe that the mechanistic understanding of large-conformational changes upon the activation of the FGFR1 kinase will aid the development of novel targeted therapeutics.Communicated by Ramaswamy H. Sarma.

Mutations influence the conformational dynamics of the GDP/KRAS complex

Mutations near allosteric sites can have a significant impact on the function of KRAS. Three specific mutations, K104Q, G12D/K104Q, and G12D/G75A, which are located near allosteric positions, were selected to investigate the molecular mechanisms behind mutation-induced influences on the activity of KRAS. Gaussian accelerated molecular dynamics (GaMD) simulations followed by the principal component analysis (PCA) were performed to improve the sampling of conformational states. The results revealed that these mutations significantly alter the structural flexibility, correlated motions, and dynamic behavior of the switch regions that are essential for KRAS binding to effectors or regulators. Furthermore, the mutations have a significant impact on the hydrogen bonding interactions between GDP and the switch regions, as well as on the electrostatic interactions of magnesium ions (Mg2+) with these regions. Our results verified that these mutations strongly influence the binding of KRAS to its effectors or regulators and allosterically regulate the activity. We believe that this work can provide valuable theoretical insights into a deeper understanding of KRAS function.Communicated by Ramaswamy H. Sarma.

Mispacking of the F87 sidechain drives aggregation-promoting conformational fluctuations in the subunit interfaces of the transthyretin tetramer

Aberrant formation and deposition of human transthyretin (TTR) aggregates causes transthyretin amyloidosis. To initialize aggregation, transthyretin tetramers must first dissociate into monomers that partially unfold to promote entry into the aggregation pathway. The native TTR tetramer (T) is stabilized by docking of the F87 sidechain into an interfacial cavity enclosed by several hydrophobic residues including A120. We have previously shown that an alternative tetramer (T*) with mispacked F87 sidechains is more prone to dissociation and aggregation than the native T state. However, the molecular basis for the reduced stability in T* remains unclear. Here we report characterization of the A120L mutant, where steric hindrance is introduced into the F87 binding site. The X-ray structure of A120L shows that the F87 sidechain is displaced from its docking site across the subunit interface. In A120S, a naturally occurring pathogenic mutant that is less aggregation-prone than A120L, the F87 sidechain is correctly docked, as in the native TTR tetramer. Nevertheless, 19F-NMR aggregation assays show an elevated population of a monomeric aggregation intermediate in A120S relative to a control containing the native A120, due to accelerated tetramer dissociation and slowed monomer tetramerization. The mispacking of the F87 sidechain is associated with enhanced exchange dynamics for interfacial residues. At 298 K, the T* populations of various naturally occurring mutants fall between 4-7% (ΔG ~ 1.5-1.9 kcal/mol), consistent with the free energy change expected for undocking and solvent exposure of one of the four F87 sidechains in the tetramer (ΔG ~ 1.6 kcal/mol). Our data provide a molecular-level picture of the likely universal F87 sidechain mispacking in tetrameric TTR that promotes interfacial conformational dynamics and increases aggregation propensity.

The conformational dynamics of Hepatitis C Virus E2 glycoprotein with the increasing number of N-glycosylation unraveled by molecular dynamics simulations

The Hepatitis C Virus (HCV), responsible for causing hepatitis and a significant contributor to liver disorders, presents a challenge for treatment due to its high genetic variability. Despite efforts, there is still no effective medication available for this virus. One of the promising targets for drug development involves targeting glycoprotein E2. However, our understanding of the dynamic behavior of E2 and its associated glycans remains limited. In this study, we investigated the dynamic characteristics of E2 with varying degrees of glycosylation using all-atom molecular dynamics simulations. We also explored glycan's interactions with the protein and among themselves. An overall increase in correlation between the vital protein regions was observed with an increase in glycan number. The protein dynamics is followed by the analysis of glycan dynamics, where the flexibility of the individual glycans was analyzed in their free and bound state, which revealed a decrease in their fluctuation in some cases. Furthermore, we generated the free energy landscape of individual N-glycan linkages in both free and bound states and observed both increases and decreases in flexibility, which can be attributed to the formation and breakage of hydrogen bonds with amino acids. Finally, we found that for a high glycosylation system, glycans interact with glycoprotein and form hydrogen bonds among themselves. Moreover, the hydrogen bond profiles of a given glycan can vary when influenced by other glycans. In summary, our study provides valuable insights into the dynamics of the core region of HCV E2 glycoprotein and its associated glycans.Communicated by Ramaswamy H. Sarma.

Studying Conformational Properties of Transmembrane Domain of KCNE3 in a Lipid Bilayer Membrane Using Molecular Dynamics Simulations

KCNE3 is a single-pass integral membrane protein that regulates numerous voltage-gated potassium channel functions such as KCNQ1. Previous solution NMR studies suggested a moderate degree of curved α-helical structure in the transmembrane domain (TMD) of KCNE3 in lyso-myristoylphosphatidylcholine (LMPC) micelles and isotropic bicelles with the residues T71, S74 and G78 situated along the concave face of the curved helix. During the interaction of KCNE3 and KCNQ1, KCNE3 pushes its transmembrane domain against KCNQ1 to lock the voltage sensor in its depolarized conformation. A cryo-EM study of KCNE3 complexed with KCNQ1 in nanodiscs suggested a deviation of the KCNE3 structure from its independent structure in isotropic bicelles. Despite the biological significance of KCNE3 TMD, the conformational properties of KCNE3 are poorly understood. Here, all atom molecular dynamics (MD) simulations were utilized to investigate the conformational dynamics of the transmembrane domain of KCNE3 in a lipid bilayer containing a mixture of POPC and POPG lipids (3:1). Further, the effect of the interaction impairing mutations (V72A, I76A and F68A) on the conformational properties of the KCNE3 TMD in lipid bilayers was investigated. Our MD simulation results suggest that the KCNE3 TMD adopts a nearly linear α helical structural conformation in POPC-POPG lipid bilayers. Additionally, the results showed no significant change in the nearly linear α-helical conformation of KCNE3 TMD in the presence of interaction impairing mutations within the sampled time frame. The KCNE3 TMD is more stable with lower flexibility in comparison to the N-terminal and C-terminal of KCNE3 in lipid bilayers. The overall conformational flexibility of KCNE3 also varies in the presence of the interaction-impairing mutations. The MD simulation data further suggest that the membrane bilayer width is similar for wild-type KCNE3 and KCNE3 containing mutations. The Z-distance measurement data revealed that the TMD residue site A69 is close to the lipid bilayer center, and residue sites S57 and S82 are close to the surfaces of the lipid bilayer membrane for wild-type KCNE3 and KCNE3 containing interaction-impairing mutations. These results agree with earlier KCNE3 biophysical studies. The results of these MD simulations will provide complementary data to the experimental outcomes of KCNE3 to help understand its conformational dynamic properties in a more native lipid bilayer environment.

Conformational dynamics of α-1 acid glycoprotein (AGP) in cancer: A comparative study of glycosylated and unglycosylated AGP

α-1 acid glycoprotein (AGP) is one of the most abundant plasma proteins. It fulfills two important functions: immunomodulation, and binding to various drugs and receptors. These different functions are closely associated and modulated via changes in glycosylation and cancer missense mutations. From a structural point of view, glycans alter the local biophysical properties of the protein leading to a diverse ligand-binding spectrum. However, glycans can typically not be observed in the resolved X-ray crystallography structure of AGP due to their high flexibility and microheterogeneity, so limiting our understanding of AGP's conformational dynamics 70 years after its discovery. We here investigate how mutations and glycosylation interfere with AGP's conformational dynamics changing its biophysical behavior, by using molecular dynamics (MD) simulations and sequence-based dynamics predictions. The MD trajectories show that glycosylation decreases the local backbone flexibility of AGP and increases the flexibility of distant regions through allosteric effects. We observe that mutations near the glycosylation site affect glycan's conformational preferences. Thus, we conclude that mutations control glycan dynamics which modulates the protein's backbone flexibility directly affecting its accessibility. These findings may assist in the drug design targeting AGP's glycosylation and mutations in cancer.

Lateral gating mechanism and plasticity of the β-barrel assembly machinery complex in micelles and Escherichia coli

The β-barrel assembly machinery (BAM) mediates the folding and insertion of the majority of outer membrane proteins (OMPs) in gram-negative bacteria. BAM is a penta-heterooligomeric complex consisting of the central β-barrel BamA and four interacting lipoproteins BamB, C, D, and E. The conformational switching of BamA between inward-open (IO) and lateral-open (LO) conformations is required for substrate recognition and folding. However, the mechanism for the lateral gating or how the structural details observed in vitro correspond with the cellular environment remains elusive. In this study, we addressed these questions by characterizing the conformational heterogeneity of BamAB, BamACDE, and BamABCDE complexes in detergent micelles and/or Escherichia coli using pulsed dipolar electron spin resonance spectroscopy (PDS). We show that the binding of BamB does not induce any visible changes in BamA, and the BamAB complex exists in the IO conformation. The BamCDE complex induces an IO to LO transition through a coordinated movement along the BamA barrel. However, the extracellular loop 6 (L6) is unaffected by the presence of lipoproteins and exhibits large segmental dynamics extending to the exit pore. PDS experiments with the BamABCDE complex in intact E. coli confirmed the dynamic behavior of both the lateral gate and the L6 in the native environment. Our results demonstrate that the BamCDE complex plays a key role in the function by regulating lateral gating in BamA.

The Role of Key Amino Acids of the Human Fe(II)/2OG-Dependent Dioxygenase ALKBH3 in Structural Dynamics and Repair Activity toward Methylated DNA

Non-heme dioxygenases of the AlkB family hold a unique position among enzymes that repair alkyl lesions in nucleic acids. These enzymes activate the Fe(II) ion and molecular oxygen through the coupled decarboxylation of the 2-oxoglutarate co-substrate to subsequently oxidize the substrate. ALKBH3 is a human homolog of E. coli AlkB, which displays a specific activity toward N1-methyladenine and N3-methylcytosine bases in single-stranded DNA. Due to the lack of a DNA-bound structure of ALKBH3, the basis of its substrate specificity and structure-function relationships requires further exploration. Here we have combined biochemical and biophysical approaches with site-directed mutational analysis to elucidate the role of key amino acids in maintaining the secondary structure and catalytic activity of ALKBH3. Using stopped-flow fluorescence spectroscopy we have shown that conformational dynamics play a crucial role in the catalytic repair process catalyzed by ALKBH3. A transient kinetic mechanism, which comprises the steps of the specific substrate binding, eversion, and anchoring within the DNA-binding cleft, has been described quantitatively by rate and equilibrium constants. Through CD spectroscopy, we demonstrated that replacing side chains of Tyr143, Leu177, and His191 with alanine results in significant alterations in the secondary structure content of ALKBH3 and decreases the stability of mutant proteins. The bulky side chain of Tyr143 is critical for binding the methylated base and stabilizing its flipped-out conformation, while its hydroxyl group is likely involved in facilitating the product release. The removal of the Leu177 and His191 side chains substantially affects the secondary structure content and conformational flexibility, leading to the complete inactivation of the protein. The mutants lacking enzymatic activity exhibit a marked decrease in antiparallel β-strands, offset by an increase in the helical component.

Dissection of integrated readout reveals the structural thermodynamics of DNA selection by transcription factors

Nucleobases such as inosine have been extensively utilized to map direct contacts by proteins in the DNA groove. Their deployment as targeted probes of dynamics and hydration, which are dominant thermodynamic drivers of affinity and specificity, has been limited by a paucity of suitable experimental models. We report a joint crystallographic, thermodynamic, and computational study of the bidentate complex of the arginine side chain with a Watson-Crick guanine (Arg×GC), a highly specific configuration adopted by major transcription factors throughout the eukaryotic branches in the Tree of Life. Using the ETS-family factor PU.1 as a high-resolution structural framework, inosine substitution for guanine resulted in a sharp dissection of conformational dynamics and hydration and elucidated their role in the DNA specificity of PU.1. Our work suggests an under-exploited utility of modified nucleobases in untangling the structural thermodynamics of interactions, such as the Arg×GC motif, where direct and indirect readout are tightly integrated.

SAXS/MC studies of the mixed-folded protein Cdt1 reveal monomeric, folded over conformations

Cdt1 is a protein critical for DNA replication licensing and is well-established to be a binding partner of the minichromosome maintenance (MCM) complex. Cdt1 has also been demonstrated to have an emerging, "moonlighting" role at the kinetochore via direct binding to microtubules and to the Ndc80 complex. However, it is not known how the structure and conformations of Cdt1 could allow for these multiple, completely unique sets of protein complexes. And while there exist multiple robust methods to study entirely folded or entirely unfolded proteins, structure-function studies of combined, mixed folded/disordered proteins remain challenging. It this work, we employ multiple orthogonal biophysical and computational techniques to provide a detailed structural characterization of human Cdt1 92-546. DSF and DSCD show both folded winged helix (WH) domains of Cdt1 are relatively unstable. CD and NMR show the N-terminal and the linker regions are intrinsically disordered. Using DLS and SEC-MALS, we show that Cdt1 is polydisperse, monomeric at high concentrations, and without any apparent inter-molecular self-association. SEC-SAXS of the monomer in solution enabled computational modeling of the protein in silico. Using the program SASSIE, we performed rigid body Monte Carlo simulations to generate a conformational ensemble. Using experimental SAXS data, we filtered for conformations which did and did not fit our data. We observe that neither fully extended nor extremely compact Cdt1 conformations are consistent with our SAXS data. The best fit models have the N-terminal and linker regions extended into solution and the two folded domains close to each other in apparent "folded over" conformations. The best fit Cdt1 conformations are consistent with a function as a scaffold protein which may be sterically blocked without the presence of binding partners. Our studies also provide a template for combining experimental and computational biophysical techniques to study mixed-folded proteins.

Dynamic conformational changes of a tardigrade group-3 late embryogenesis abundant protein modulate membrane biophysical properties

A number of intrinsically disordered proteins (IDPs) encoded in stress-tolerant organisms, such as tardigrade, can confer fitness advantage and abiotic stress tolerance when heterologously expressed. Tardigrade-specific disordered proteins including the cytosolic-abundant heat-soluble proteins are proposed to confer stress tolerance through vitrification or gelation, whereas evolutionarily conserved IDPs in tardigrades may contribute to stress tolerance through other biophysical mechanisms. In this study, we characterized the mechanism of action of an evolutionarily conserved, tardigrade IDP, HeLEA1, which belongs to the group-3 late embryogenesis abundant (LEA) protein family. HeLEA1 homologs are found across different kingdoms of life. HeLEA1 is intrinsically disordered in solution but shows a propensity for helical structure across its entire sequence. HeLEA1 interacts with negatively charged membranes via dynamic disorder-to-helical transition, mainly driven by electrostatic interactions. Membrane interaction of HeLEA1 is shown to ameliorate excess surface tension and lipid packing defects. HeLEA1 localizes to the mitochondrial matrix when expressed in yeast and interacts with model membranes mimicking inner mitochondrial membrane. Yeast expressing HeLEA1 shows enhanced tolerance to hyperosmotic stress under nonfermentative growth and increased mitochondrial membrane potential. Evolutionary analysis suggests that although HeLEA1 homologs have diverged their sequences to localize to different subcellular organelles, all homologs maintain a weak hydrophobic moment that is characteristic of weak and reversible membrane interaction. We suggest that such dynamic and weak protein-membrane interaction buffering alterations in lipid packing could be a conserved strategy for regulating membrane properties and represent a general biophysical solution for stress tolerance across the domains of life.

Integrated spectroscopic and MD simulation approach to decipher the effect of pH on the structure function of Staphylococcus aureus thymidine kinase

Staphylococcus aureus is a major human pathogen responsible for a variety of clinical infections, becoming increasingly resistant to antibiotics. To address this challenge, there is a need to identify new cellular targets and innovative approaches to expand treatment options. One such target is thymidine kinase (TK), a crucial enzyme in the pyrimidine salvage pathway, which plays a key role in the phosphorylation of thymidine, an essential component in DNA synthesis and repair. In this study, we have successfully cloned, expressed, and purified the TK protein. A comprehensive investigation into how different pH levels affect the structure and functional activity of TK, using a combination of spectroscopy, classical molecular dynamics simulations, and enzyme activity assays was conducted. Our study revealed that variation in pH disrupts secondary and tertiary structures of TK with noticeable aggregate formation at pH 5.0. Enzyme activity studies demonstrated that TK exhibited its maximum kinase activity within the physiological pH range. These findings strongly suggest a connection between structural changes and enzymatic activity, which was further supported by the agreement between the spectroscopic features we measured and the results of our MD simulations. Our study provides a deeper insight into the structural features of TK, which could potentially be harnessed for the development of therapeutic strategies aimed at combatting infectious diseases. Conformational dynamics plays an essential role in the design and development of effective inhibitors. Considering the effects of pH on the conformational dynamics of TK, our findings may be implicated in the development of potent and selective inhibitors.Communicated by Ramaswamy H. Sarma.

Decoding molecular mechanism of species-selective targeting of fungal versus human HSP90 using multiple replica molecular dynamics simulations and binding free energy calculations

As a highly evolutionarily conserved molecular chaperone, heat shock protein (HSP90), plays an important role in virulence traits, representing a therapeutic target for the treatment of fungal infections. The close evolutionary relationship between fungi and their human hosts poses a key challenge for the development of selective antifungal agents. In this work, molecular docking, multiple replica microsecond-based molecular dynamics (MD) simulations, and binding free energy calculations were performed to decode molecular mechanism of species-selective targeting of fungal versus human HSP90 triggered by the compound A11. MD simulations reveal that binding of compound A11 to human HSP90 nucleotide-binding domain (NBD) leads to obvious conformational changes relative to fungal HSP90 NBD. Binding free energy calculations show that the binding of compound A11 to fungal HSP90 NBD is stronger than that to human HSP90 NBD. Per residue-based free energy decomposition analysis was used to evaluate the inhibitor - residue interaction profile. The results efficiently identify the hot spot residues that play vital roles in favorable binding of compound A11 to fungal HSP90 NBD. This study is expected to provide a useful guidance for the development of selective inhibitors toward fungal HSP90.Communicated by Ramaswamy H. Sarma.

Magnesium Ion-Driven Folding and Conformational Switching Kinetics of Tetracycline Binding Aptamer: Implications for in vivo Riboswitch Engineering

Engineering in vitro selected RNA aptamers into in vivo functional riboswitches represents a long-standing challenge in molecular biology. The highly specific aptamer domain of the riboswitch undergoes a conformational adjustment in response to ligand sensing, which in turn exerts the regulatory function. Besides essential factors like structural complexity and ligand binding kinetics, the active role of magnesium ions in stabilizing RNA tertiary structures and assisting in ligand binding can be a vital criterion. We present spectroscopic studies on the magnesium ion-driven folding of the Tetracycline binding aptamer. Using fluorescent labels, the aptamer pre-folding and subsequent ligand binding is monitored by magnesium titration experiments and time-resolved stopped-flow measurements. A minimum concentration of 0.5 mM magnesium is required to fold into a magnesium ion-stabilized binding-competent state with a preformed binding pocket. Tetracycline binding causes a pronounced conformational change that results in the establishment of the triple helix core motif, and that further propagates towards the closing stem. By a dynamic acquisition of magnesium ions, a kink motif is formed at the intersection of the triple helix and closing stem regions. This ultimately entails a stabilization of the closing stem which is discussed as a key element in the regulatory function of the Tetracycline aptamer.

Insights into the solution structure and transcriptional regulation of the MazE9 antitoxin in Mycobacterium tuberculosis

The present study endeavors to decode the details of the transcriptional autoregulation effected by the MazE9 antitoxin of the Mycobacterium tuberculosis MazEF9 toxin-antitoxin system. Regulation of this bicistronic operon at the level of transcription is a critical biochemical process that is key for the organism's stress adaptation and virulence. Here, we have reported the solution structure of the DNA binding domain of MazE9 and scrutinized the thermodynamic and kinetic parameters operational in its interaction with the promoter/operator region, specific to the mazEF9 operon. A HADDOCK model of MazE9 bound to its operator DNA has been calculated based on the information on interacting residues obtained from these studies. The thermodynamics and kinetics of the interaction of MazE9 with the functionally related mazEF6 operon indicate that the potential for intracellular cross-regulation is unlikely. An interesting feature of MazE9 is the cis ⇌ trans conformational isomerization of proline residues in the intrinsically disordered C-terminal domain of this antitoxin.

Permissive Conformations of a Transmembrane Helix Allow Intramembrane Proteolysis by γ-Secretase

The intramembrane protease γ-secretase activates important signaling molecules, such as Notch receptors. It is still unclear, however, how different elements within the primary structure of substrate transmembrane domains (TMDs) contribute to their cleavability. Using a newly developed yeast-based cleavage assay, we identified three crucial regions within the TMDs of the paralogs Notch1 and Notch3 by mutational and gain-of-function approaches. The AAAA or AGAV motifs within the N-terminal half of the TMDs were found to confer strong conformational flexibility to these TMD helices, as determined by mutagenesis coupled to deuterium/hydrogen exchange. Crucial amino acids within the C-terminal half may support substrate docking into the catalytic cleft of presenilin, the enzymatic subunit of γ-secretase. Further, residues close to the C-termini of the TMDs may stabilize a tripartite β-sheet in the substrate/enzyme complex. NMR structures reveal different extents of helix bending as well as an ability to adopt widely differing conformational substates, depending on the sequence of the N-terminal half. The difference in cleavability between Notch1 and Notch3 TMDs is jointly determined by the conformational repertoires of the TMD helices and the sequences of the C-terminal half, as suggested by mutagenesis and building molecular models. In sum, cleavability of a γ-secretase substrate is enabled by different functions of cooperating TMD regions, which deepens our mechanistic understanding of substrate/non-substrate discrimination in intramembrane proteolysis.

Solution structure of the HOIL-1L NZF domain reveals a conformational switch regulating linear ubiquitin affinity

Attachment of polyubiquitin (poly-Ub) chains to proteins is a major posttranslational modification in eukaryotes. Linear ubiquitin chain assembly complex, consisting of HOIP (HOIL-1-interacting protein), HOIL-1L (heme-oxidized IRP2 Ub ligase 1), and SHARPIN (Shank-associated RH domain-interacting protein), specifically synthesizes "head-to-tail" poly-Ub chains, which are linked via the N-terminal methionine α-amino and C-terminal carboxylate of adjacent Ub units and are thus commonly called "linear" poly-Ub chains. Linear ubiquitin chain assembly complex-assembled linear poly-Ub chains play key roles in immune signaling and suppression of cell death and have been associated with immune diseases and cancer; HOIL-1L is one of the proteins known to selectively bind linear poly-Ub via its Npl4 zinc finger (NZF) domain. Although the structure of the bound form of the HOIL-1L NZF domain with linear di-Ub is known, several aspects of the recognition specificity remain unexplained. Here, we show using NMR and orthogonal biophysical methods, how the NZF domain evolves from a free to the specific linear di-Ub-bound state while rejecting other potential Ub species after weak initial binding. The solution structure of the free NZF domain revealed changes in conformational stability upon linear Ub binding, and interactions between the NZF core and tail revealed conserved electrostatic contacts, which were sensitive to charge modulation at a reported phosphorylation site: threonine-207. Phosphomimetic mutations reduced linear Ub affinity by weakening the integrity of the linear di-Ub-bound conformation. The described molecular determinants of linear di-Ub binding provide insight into the dynamic aspects of the Ub code and the NZF domain's role in full-length HOIL-1L.

Single-molecule analysis reveals that a glucagon-bound extracellular domain of the glucagon receptor is dynamic

Dynamic information is vital to understanding the activation mechanism of G protein-coupled receptors (GPCRs). Despite the availability of high-resolution structures of different conformational states, the dynamics of those states at the molecular level are poorly understood. Here, we used total internal reflection fluorescence microscopy to study the extracellular domain (ECD) of the glucagon receptor (GCGR), a class B family GPCR that controls glucose homeostasis. Single-molecule fluorescence resonance energy transfer was used to observe the ECD dynamics of GCGR molecules expressed and purified from mammalian cells. We observed that for apo-GCGR, the ECD is dynamic and spent time predominantly in a closed conformation. In the presence of glucagon, the ECD is wide open and also shows more dynamic behavior than apo-GCGR, a finding that was not previously reported. These results suggest that both apo-GCGR and glucagon-bound GCGRs show reversible opening and closing of the ECD with respect to the seven-transmembrane (7TM) domain. This work demonstrates a molecular approach to visualizing the dynamics of the GCGR ECD and provides a foundation for understanding the conformational changes underlying GPCR activation, which is critical in the development of new therapeutics.

Understanding IP3R channels: From structural underpinnings to ligand-dependent conformational landscape

Inositol 1,4,5-trisphosphate receptors (IP3Rs) are ubiquitously expressed large-conductance Ca2+-permeable channels predominantly localized to the endoplasmic reticulum (ER) membranes of virtually all eukaryotic cell types. IP3Rs work as Ca2+ signaling hubs through which diverse extracellular stimuli and intracellular inputs are processed and then integrated to result in delivery of Ca2+ from the ER lumen to generate cytosolic Ca2+ signals with precise temporal and spatial properties. IP3R-mediated Ca2+ signals control a vast repertoire of cellular functions ranging from gene transcription and secretion to the more enigmatic brain activities such as learning and memory. IP3Rs open and release Ca2+ when they bind both IP3 and Ca2+, the primary channel agonists. Despite overwhelming evidence supporting functional interplay between IP3 and Ca2+ in activation and inhibition of IP3Rs, the mechanistic understanding of how IP3R channels convey their gating through the interplay of two primary agonists remains one of the major puzzles in the field. The last decade has seen much progress in the use of cryogenic electron microscopy to elucidate the molecular mechanisms of ligand binding, ion permeation, ion selectivity and gating of the IP3R channels. The results of these studies, summarized in this review, provide a prospective view of what the future holds in structural and functional research of IP3Rs.

Partial agonism in heteromeric GLUK2/GLUK5 kainate receptor

Kainate receptors are a subtype of ionotropic glutamate receptors that form transmembrane channels upon binding glutamate. Here, we have investigated the mechanism of partial agonism in heteromeric GluK2/K5 receptors, where the GluK2 and GluK5 subunits have distinct agonist binding profiles. Using single-molecule Förster resonance energy transfer, we found that at the bi-lobed agonist-binding domain, the partial agonist AMPA-bound receptor occupied intermediate cleft closure conformational states at the GluK2 cleft, compared to the more open cleft conformations in apo form and more closed cleft conformations in the full agonist glutamate-bound form. In contrast, there is no significant difference in cleft closure states at the GluK5 agonist-binding domain between the partial agonist AMPA- and full agonist glutamate-bound states. Additionally, unlike the glutamate-bound state, the dimer interface at the agonist-binding domain is not decoupled in the AMPA-bound state. Our findings suggest that partial agonism observed with AMPA binding is mediated primarily due to differences in the GluK2 subunit, highlighting the distinct contributions of the subunits towards activation.

Inner Amino Acid Contacts Are Key Factors of Multistage Structural Rearrangements of DNA and Affect Substrate Specificity of Apurinic/Apyrimidinic Endonuclease APE1

Apurinic/apyrimidinic endonuclease 1 (APE1) is one of the most important enzymes in base excision repair. Studies on this enzyme have been conducted for a long time, but some aspects of its activity remain poorly understood. One such question concerns the mechanism of damaged-nucleotide recognition by the enzyme, and the answer could shed light on substrate specificity control in all enzymes of this class. In the present study, by pulsed electron-electron double resonance (DEER, also known as PELDOR) spectroscopy and pre-steady-state kinetic analysis along with wild-type (WT) APE1 from Danio rerio (zAPE1) or three mutants (carrying substitution N253G, A254G, or E260A), we aimed to elucidate the molecular events in the process of damage recognition. The data revealed that the zAPE1 mutant E260A has much higher activity toward DNA substrates containing 5,6-dihydro-2'-deoxyuridine (DHU), 2'-deoxyuridine (dU), alpha-2'-deoxyadenosine (αA), or 1,N6-ethenoadenosine (εA). Examination of conformational changes in DNA clearly revealed multistep DNA rearrangements during the formation of the catalytic complex. These structural rearrangements of DNA are directly associated with the capacity of damaged DNA for enzyme-induced bending and unwinding, which are required for eversion of the damaged nucleotide from the DNA duplex and for its placement into the active site of the enzyme. Taken together, the results experimentally prove the factors that control substrate specificity of the AP endonuclease zAPE1.

Individual Contributions of Amido Acid Residues Tyr122, Ile168, and Asp173 to the Activity and Substrate Specificity of Human DNA Dioxygenase ABH2

Human Fe(II)/α-ketoglutarate-dependent dioxygenase ABH2 plays a crucial role in the direct reversal repair of nonbulky alkyl lesions in DNA nucleobases, e.g., N¹-methyladenine (m¹A), N³-methylcytosine (m³C), and some etheno derivatives. Moreover, ABH2 is capable of a less efficient oxidation of an epigenetic DNA mark called 5-methylcytosine (m⁵C), which typically is a specific target of DNA dioxygenases from the TET family. In this study, to elucidate the mechanism of the substrate specificity of ABH2, we investigated the role of several active-site amino acid residues. Functional mapping of the lesion-binding pocket was performed through the analysis of the functions of Tyr122, Ile168, and Asp173 in the damaged base recognition mechanism. Interactions of wild-type ABH2, or its mutants Y122A, I168A, or D173A, with damaged DNA containing the methylated base m¹A or m³C or the epigenetic marker m⁵C were analyzed by molecular dynamics simulations and kinetic assays. Comparative analysis of the enzymes revealed an effect of the substitutions on DNA binding and on catalytic activity. Obtained data clearly demonstrate the effect of the tested amino acid residues on the catalytic activity of the enzymes rather than the DNA-binding ability. Taken together, these data shed light on the molecular and kinetic consequences of the substitution of active-site residues for the mechanism of the substrate recognition.

Coupling of Conformational Dynamics and Inhibitor Binding in the Phosphodiesterase-5 Family

Phosphodiesterase-5 (PDE5) is responsible for regulating the concentration of the second messenger molecule cGMP by hydrolyzing it into 5'-GMP. PDE5 is implicated in erectile dysfunction and cardiovascular diseases. The substrate binding site in the catalytic domain of PDE5 is surrounded by several dynamic structural motifs (including the α14 helix, M-loop, and H-loop) that are known to switch between inactive and active conformational states via currently unresolved structural intermediates. We evaluated the conformational dynamics of these structural motifs in the apo state and upon binding of an allosteric inhibitor (evodamine) or avanafil, a competitive inhibitor. We employed enhanced sampling based replica exchange solute scaling (REST2) method, principal component analysis (PCA), time-lagged independent component analysis (tICA), molecular dynamics (MD) simulations, and well-tempered metadynamics simulations to probe the conformational changes in these structural motifs. Our results support a regulatory mechanism for PDE5, where the α14 helix alternates between an inward (lower activity) conformation and an outward (higher activity) conformation that is accompanied by the folding/unfolding of the α 8' and α 8″ helices of the H-loop. When the allosteric inhibitor evodamine is bound to PDE5, the inward (inactive) state of the α 14 helix is preferred, thus preventing substrate access to the catalytic site. In contrast, competitive inhibitors of PDE5 block catalysis by occupying the active site accompanied by stabilization of the outward conformation of the α 14 helix. Defining the conformational dynamics underlying regulation of PDE5 activation will be helpful in rational design of next-generation small molecules modulators of PDE5 activity. This article is protected by copyright. All rights reserved.

Conformational changes in the activation loop of a bacterial PASTA kinase

Many bacterial genomes encode a transmembrane protein kinase belonging to the PASTA kinase family, which control numerous processes in diverse bacterial pathogens, including antibiotic resistance, cell division, stress resistance, toxin production, and virulence. PASTA kinases share a conserved three-part domain architecture, consisting of an extracellular PASTA domain, proposed to sense the peptidoglycan layer status, a single transmembrane helix, and an intracellular Ser/Thr kinase domain. The crystal structures of the kinase domain from two homologous PASTA kinases reveal a characteristic two-lobed structure typical of eukaryotic protein kinases with a centrally located, but unresolved, activation loop that becomes phosphorylated and regulates downstream signaling pathways. We previously identified three sites of phosphorylation on the activation loop (T163, T166, T168) of IreK, a PASTA kinase from the pathogen Enterococcus faecalis, as well as a distal phosphorylation site (T218) that each influence IreK activity in vivo. Still, the mechanism by which loop phosphorylation regulates PASTA kinase function is yet unknown. Therefore, we utilized site-directed spin labeling (SDSL) and continuous wave (CW) electron paramagnetic resonance (EPR) spectroscopy to assess the E. faecalis IreK kinase activation loop dynamics, including the effects of phosphorylation on activation loop motion, and the IreK-IreB interaction. Our results reveal that the IreK activation loop occupies a more immobile state when dephosphorylated, and that loop autophosphorylation shifts the loop to a more mobile state that can then enable interaction with IreB, a known substrate. This article is protected by copyright. All rights reserved.

Arachidonic acid reverses cholesterol and zinc inhibition of human voltage-gated proton channels

Unlike other members of the voltage-gated ion channel superfamily, voltage-gated proton (Hv) channels are solely composed of voltage sensor domains without separate ion-conducting pores. Due to their unique dependence on both voltage and transmembrane pH gradients, Hv channels normally open to mediate proton efflux. Multiple cellular ligands were also found to regulate the function of Hv channels, including Zn²⁺, cholesterol, polyunsaturated arachidonic acid, and albumin. Our previous work showed that Zn²⁺ and cholesterol inhibit the human voltage-gated proton channel hHv1 by stabilizing its S4 segment at resting state conformations. Released from phospholipids by phospholipase A2 in cells upon infection or injury, arachidonic acid regulates the function of many ion channels, including hHv1. In the present work, we examined the effects of arachidonic acid on purified hHv1 channels using liposome flux assays and revealed underlying structural mechanisms using single-molecule Fluorescence Resonance Energy Transfer (smFRET). Our data indicated that arachidonic acid strongly activates hHv1 channels by promoting transitions of the S4 segment towards opening or 'pre-opening' conformations. Moreover, we found that arachidonic acid even activates hHv1 channels inhibited by Zn²⁺ and cholesterol, providing a biophysical mechanism to activate hHv1 channels in non-excitable cells upon infection or injury.

Changes in Active Site Loop Conformation Relate to the Transition toward a Novel Enzymatic Activity

Enzymatic promiscuity, the ability of enzymes to catalyze multiple, distinct chemical reactions, has been well documented and is hypothesized to be a major driver for the emergence of new enzymatic functions. Yet, the molecular mechanisms involved in the transition from one activity to another remain debated and elusive. Here, we evaluated the redesign of the active site binding cleft of the lactonase SsoPox using structure-based design and combinatorial libraries. We created variants with largely improved catalytic abilities against phosphotriesters, the best ones being > 1,000-fold better compared to the wild-type enzyme. The observed shifts in activity specificity are large, ~1,000,000-fold and beyond, since some variants completely lost their initial activity. The selected combinations of mutations have considerably reshaped the active site cavity via side chain changes but mostly through large rearrangements of the active site loops, as revealed by a suite of crystal structures. This suggests that specific active site loop configuration is critical to the lactonase activity. Interestingly, analysis of high-resolution structures hints at the potential role of conformational sampling and its directionality in defining an enzyme activity profile.

Exploiting conformational dynamics to modulate the function of designed proteins

With the recent success in calculating protein structures from amino acid sequences using artificial intelligence-based algorithms, an important next step is to decipher how dynamics is encoded by the primary protein sequence so as to better predict function. Such dynamics information is critical for protein design, where strategies could then focus not only on sequences that fold into particular structures that perform a given task, but would also include low-lying excited protein states that could influence the function of the designed protein. Herein, we illustrate the importance of dynamics in modulating the function of C34, a designed α/β protein that captures β-strands of target ligands and is a member of a family of proteins designed to sequester β-strands and β hairpins of aggregation-prone molecules that lead to a variety of pathologies. Using a strategy to "see" regions of apo C34 that are invisible to NMR spectroscopy as a result of pervasive conformational exchange, as well as a mutagenesis approach whereby C34 molecules are stabilized into a single conformer, we determine the structures of the predominant conformations that are sampled by C34 and show that these attenuate the affinity for cognate peptide. Subsequently, the observed motion is exploited to develop an allosterically regulated peptide binder whose binding affinity can be controlled through the addition of a second molecule. Our study emphasizes the unique role that NMR can play in directing the design process and in the construction of new molecules with more complex functionality.

Activation loop phosphorylation tunes conformational dynamics underlying Pyk2 tyrosine kinase activation

Pyk2 is a multidomain non-receptor tyrosine kinase that undergoes a multistage activation mechanism. Activation is instigated by conformational rearrangements relieving autoinhibitory FERM domain interactions. The kinase autophosphorylates a central linker residue to recruit Src kinase. Pyk2 and Src mutually phosphorylate activation loops to confer full activation. While the mechanisms of autoinhibition are established, the conformational dynamics associated with autophosphorylation and Src recruitment remain unclear. We employ hydrogen/deuterium exchange mass spectrometry and kinase activity profiling to map the conformational dynamics associated with substrate binding and Src-mediated activation loop phosphorylation. Nucleotide engagement stabilizes the autoinhibitory interface, while phosphorylation deprotects both FERM and kinase regulatory surfaces. Phosphorylation organizes active site motifs linking catalytic loop with activation segment. Dynamics of the activation segment anchor propagate to EF/G helices to prevent reversion of the autoinhibitory FERM interaction. We employ targeted mutagenesis to dissect how phosphorylation-induced conformational rearrangements elevate kinase activity above the basal autophosphorylation rate.

Slow conformational changes in the rigid and highly stable chymotrypsin inhibitor 2

Slow conformational changes are often directly linked to protein function. It is however less clear how such processes may perturb the overall folding stability of a protein. We previously found that the stabilizing double mutant L49I/I57V in the small protein chymotrypsin inhibitor 2 from barley led to distributed increased nanosecond and faster dynamics. Here we asked what effects the L49I and I57V substitutions, either individually or together, have on the slow conformational dynamics of CI2. We used ¹⁵ N CPMG spin relaxation dispersion experiments to measure the kinetics, thermodynamics and structural changes associated with slow conformational change in CI2. These changes result in an excited state that is populated to 4.3% at 1 °C. As the temperature is increased the population of the excited state decreases. Structural changes in the excited state are associated with residues that interact with water molecules that have well defined positions and are found at these positions in all crystal structures of CI2. The substitutions in CI2 have only little effect on the structure of the excited state whereas the stability of the excited state to some extent follows the stability of the main state. The minor state is thus most populated for the most stable CI2 variant and least populated for the least stable variant. We hypothesize that the interactions between the substituted residues and the well-ordered water molecules links subtle structural changes around the substituted residues to the region in the protein that experience slow conformational changes. This article is protected by copyright. All rights reserved.

Controlled modulation of the dynamics of the Deinococcus grandis Dps N-terminal tails by divalent metals

DNA-binding proteins from starved cells (Dps) are small multifunctional nanocages expressed by prokaryotes in acute oxidative stress conditions or during the starvation-induced stationary phase, as a bacterial defense mechanism. Dps proteins protect bacterial DNA from damage by either direct binding or by removing precursors of reactive oxygen species from solution. The DNA-binding properties of most Dps proteins studied so far are related to their unordered, flexible, N- and C-terminal extensions. In a previous work, we revealed that the N-terminal tails of Deinoccocus grandis Dps shift from an extended to a compact conformation depending on the ionic strength of the buffer and detected a novel high-spin ferrous iron center in the proximal ends of those tails. In this work, we further explore the conformational dynamics of the protein by probing the effect of divalent metals binding to the tail by comparing the metal-binding properties of the wild-type protein with a binding site-impaired D34A variant using size exclusion chromatography, dynamic light scattering, synchrotron radiation circular dichroism, and small-angle X-ray scattering. The N-terminal ferrous species was also characterized by Mössbauer spectroscopy. The results herein presented reveal that the conformation of the N-terminal tails is altered upon metal binding in a gradual, reversible, and specific manner. These observations may point towards the existence of a regulatory process for the DNA-binding properties of Dps proteins through metal binding to their N- and/or C-terminal extensions.

Conformational Dynamics and Stability of Bilayers Formed by Mycolic Acids from the Mycobacterium tuberculosis Outer Membrane

Bilayers of mycolic acids (MAs) form the outer membrane of Mycobacterium tuberculosis that has high strength and extremely low permeability for external molecules (including antibiotics). For the first time, we were able to study them using the all-atom long-term molecular dynamic simulations (from 300 ns up to 1.2 μs) in order to investigate the conformational changes and most favorable structures of the mycobacterial membranes. The structure and properties of the membranes are crucially dependent on the initial packing of the α-mycolic acid (AMA) molecules, as well as on the presence of the secondary membrane components, keto- and methoxy mycolic acids (KMAs and MMAs). In the case of AMA-based membranes, the most labile conformation is W while other types of conformations (sU as well as sZ, eU, and eZ) are much more stable. In the multicomponent membranes, the presence of the KMA and MMA components (in the W conformation) additionally stabilizes both the W and eU conformations of AMA. The membrane in which AMA prevails in the eU conformation is much thicker and, at the same time, much denser. Such a packing of the MA molecules promotes the formation of a significantly stronger outer mycobacterial membrane that should be much more resistant to the threatening external factors.

A protein engineering approach toward understanding FKBP51 conformational dynamics and mechanisms of ligand binding

Most proteins are flexible molecules that coexist in an ensemble of several conformations. Point mutations in the amino acid sequence of a protein can trigger structural changes that drive the protein population to a conformation distinct from the native state. Here, we report a protein engineering approach to better understand protein dynamics and ligand binding of the FK506-binding protein 51 (FKBP51), a prospective target for stress-related diseases, metabolic disorders, some types of cancers and chronic pain. By randomizing selected regions of its ligand-binding domain and sorting yeast display libraries expressing these variants, mutants with high affinity to conformation-specific FKBP51 selective ligands were identified. These improved mutants are valuable tools for the discovery of novel selective ligands that preferentially and specifically bind the FKBP51 active site in its open conformation state. Moreover, they will help us understand the conformational dynamics and ligand binding mechanics of the FKBP51 binding pocket.

Developments in rapid hydrogen-deuterium exchange methods

Biological macromolecules, such as proteins, nucleic acids, and carbohydrates, contain heteroatom-bonded hydrogens that undergo exchange with solvent hydrogens on timescales ranging from microseconds to hours. In hydrogen-deuterium exchange mass spectrometry (HDX-MS), this exchange process is used to extract information about biomolecular structure and dynamics. This minireview focuses on millisecond timescale HDX-MS measurements, which, while less common than 'conventional' timescale (seconds to hours) HDX-MS, provide a unique window into weakly structured species, weak (or fast cycling) binding interactions, and subtle shifts in conformational dynamics. This includes intrinsically disordered proteins and regions (IDPs/IDRs) that are associated with cancer and amyloidotic neurodegenerative disease. For nucleic acids and carbohydrates, structures such as isomers, stems, and loops, can be elucidated and overall structural rigidity can be assessed. We will provide a brief overview of technical developments in rapid HDX followed by highlights of various applications, emphasising the importance of broadening the HDX timescale to improve throughput and to capture a wider range of function-relevant dynamic and structural shifts.

Ligands-induced open-close conformational change during DapE catalysis: Insights from molecular dynamics simulations

The microbial enzyme DapE plays a critical role in the lysine biosynthetic pathway and is considered as a potentially safe antibiotic target. In this study, atomistic simulations are employed to identify the modes of essential dynamics that define the conformational response of the enzyme to ligand binding and unbinding. The binding modes and the binding affinities of the products to the DapE enzyme are estimated from the MM-PBSA method, and the residues contributing to the ligand binding are identified. Various structural analyses and the principal component analysis of the molecular dynamics trajectories reveal that the removal of products from the active site causes a significant change in the overall enzyme structure. Both Cartesian and dihedral principal component analyses are used to characterize the structural changes in terms of domain unfolding and domain twisting motions. In the most dominant mode, that is, the domain unfolding motion, the two catalytic domains move away from the two dimerization domains of the dimeric enzyme, representing a closed-to-open conformational change. The conformational changes are initiated by the coordinated movement of three loops (Asp75-Pro82, Gly240-Asn244, and Thr347-Glu353) that trigger a domain-level movement. From multiple short trajectories, the time constant associated with the domain opening motion is estimated as 43.6 ns. Physiologically, this close-to-open conformational change is essential for the regeneration of the initial state of the enzyme for the subsequent cycle of catalytic action and provides the apo enzyme enough flexibility for efficient substrate binding.

Single-Molecule Study of DNAzyme Reveals Its Intrinsic Conformational Dynamics

DNAzyme is a class of DNA molecules that can perform catalytic functions with high selectivity towards specific metal ions. Due to its potential applications for biosensors and medical therapeutics, DNAzyme has been extensively studied to characterize the relationships between its biochemical properties and functions. Similar to protein enzymes and ribozymes, DNAzymes have been found to undergo conformational changes in a metal-ion-dependent manner for catalysis. Despite the important role the conformation plays in the catalysis process, such structural and dynamic information might not be revealed by conventional approaches. Here, by using the single-molecule fluorescence resonance energy transfer (smFRET) technique, we were able to investigate the detailed conformational dynamics of a uranyl-specific DNAzyme 39E. We observed conformation switches of 39E to a folded state with the addition of Mg²⁺ and to an extended state with the addition of UO₂²⁺. Furthermore, 39E can switch to a more compact configuration with or without divalent metal ions. Our findings reveal that 39E can undergo conformational changes spontaneously between different configurations.

Peptide Triazole Inhibitors of HIV-1: Hijackers of Env Metastability

With 1.5 million new infections and 690,000 AIDS-related deaths globally each year, HIV- 1 remains a pathogen of significant public health concern. Although a wide array of effective antiretroviral drugs have been discovered, these largely target intracellular stages of the viral infectious cycle, and inhibitors that act at or before the point of viral entry still require further advancement. A unique class of HIV-1 entry inhibitors, called peptide triazoles (PTs), has been developed, which irreversibly inactivates Env trimers by exploiting the protein structure's innate metastable nature. PTs, and a related group of inhibitors called peptide triazole thiols (PTTs), are peptide compounds that dually engage the CD4 receptor and coreceptor binding sites of Env's gp120 subunit. This triggers dramatic conformational rearrangements of Env, including the shedding of gp120 (PTs and PTTs) and lytic transformation of the gp41 subunit to a post-fusion-like arrangement (PTTs). Due to the nature of their dual receptor site engagement, PT/PTT-induced conformational changes may elucidate mechanisms behind the native fusion program of Env trimers following receptor and coreceptor engagement, including the role of thiols in fusion. In addition to inactivating Env, PTT-induced structural transformation enhances the exposure of important and conserved neutralizable regions of gp41, such as the membrane proximal external region (MPER). PTT-transformed Env could present an intriguing potential vaccine immunogen prototype. In this review, we discuss the origins of the PT class of peptide inhibitors, our current understanding of PT/PTT-induced structural perturbations and viral inhibition, and prospects for using these antagonists for investigating Env structural mechanisms and for vaccine development.

Conformational dynamics of the Hsp70 chaperone throughout key steps of its ATPase cycle

The 70 kDa heat shock proteins (Hsp70s) are highly versatile molecular chaperones that assist in a wide variety of protein-folding processes. They exert their functions by continuously cycling between states of low and high affinity for client polypeptides, driven by ATP-binding and hydrolysis. This cycling is tuned by cochaperones and clients. Although structures for the high and low client affinity conformations of Hsp70 and Hsp70 domains in complex with various cochaperones and peptide clients are available, it is unclear how structural rearrangements in the presence of cochaperones and clients are orchestrated in space and time. Here, we report insights into the conformational dynamics of the prokaryotic model Hsp70 DnaK throughout its adenosine-5'-triphosphate hydrolysis (ATPase) cycle using proximity-induced fluorescence quenching. Our data suggest that ATP and cochaperone-induced structural rearrangements in DnaK occur in a sequential manner and resolve hitherto unpredicted cochaperone and client-induced structural rearrangements. Peptides induce large conformational changes in DnaK·ATP prior to ATP hydrolysis, whereas a protein client induces significantly smaller changes but is much more effective in stimulating ATP hydrolysis. Analysis of the enthalpies of activation for the ATP-induced opening of the DnaK lid in the presence of clients indicates that the lid does not exert an enthalpic pulling force onto bound clients, suggesting entropic pulling as a major mechanism for client unfolding. Our data reveal important insights into the mechanics, allostery, and dynamics of Hsp70 chaperones. We established a methodology for understanding the link between dynamics and function, Hsp70 diversity, and activity modulation.

The Optimized Conformation Dynamics of the KcsA Filter as a Probe for Lateral Membrane Effects: A First Principle Based Femto-Sec Resolution MD Study

We provide a high resolution, all-atom, femto-second molecular dynamics (MD) simulation of the passage of K⁺ ions and H₂O molecules through the selectivity filter of the KcsA potassium ion channel, based on first principle physical methods. Our results show that a change in the length of the selectivity filter of as little as 3%, regardless of whether the filter is made longer or shorter, will reduce the K⁺ ion current by around 50%. In addition, further squeezing or stretching by about 9% can effectively stop the current. Our results demonstrate optimized conformational dynamics that associate an increased mobility of parts in the filter linings with a standard configuration, leading to maximized conduction rates that are highly sensitive to geometrical distortions. We discuss this latter aspect in relation to lateral membrane effects on the filter region of ion channels and the 'force from lipids' hypothesis.

Membrane phosphoinositides regulate GPCR-β-arrestin complex assembly and dynamics

Binding of arrestin to phosphorylated G protein-coupled receptors (GPCRs) is crucial for modulating signaling. Once internalized, some GPCRs remain complexed with β-arrestins, while others interact only transiently; this difference affects GPCR signaling and recycling. Cell-based and in vitro biophysical assays reveal the role of membrane phosphoinositides (PIPs) in β-arrestin recruitment and GPCR-β-arrestin complex dynamics. We find that GPCRs broadly stratify into two groups, one that requires PIP binding for β-arrestin recruitment and one that does not. Plasma membrane PIPs potentiate an active conformation of β-arrestin and stabilize GPCR-β-arrestin complexes by promoting a fully engaged state of the complex. As allosteric modulators of GPCR-β-arrestin complex dynamics, membrane PIPs allow for additional conformational diversity beyond that imposed by GPCR phosphorylation alone. For GPCRs that require membrane PIP binding for β-arrestin recruitment, this provides a mechanism for β-arrestin release upon translocation of the GPCR to endosomes, allowing for its rapid recycling.

Encoded Conformational Dynamics of the HIV Splice Site A3 Regulatory Locus: Implications for Differential Binding of hnRNP Splicing Auxiliary Factors

Alternative splicing of the HIV transcriptome is controlled through cis regulatory elements functioning as enhancers or silencers depending on their context and the type of host RNA binding proteins they recruit. Splice site acceptor A3 (ssA3) is one of the least used acceptor sites in the HIV transcriptome and its activity determines the levels of tat mRNA. Splice acceptor 3 is regulated by a combination of cis regulatory sequences, auxiliary splicing factors, and presumably RNA structure. The mechanisms by which these multiple regulatory components coordinate to determine the frequency in which ssA3 is utilized is poorly understood. By NMR spectroscopy and phylogenetic analysis, we show that the ssA3 regulatory locus is conformationally heterogeneous and that the sequences that encompass the locus are conserved across most HIV isolates. Despite the conformational heterogeneity, the major stem loop (A3SL1) observed in vitro folds to base pair the Polypyrimdine Tract (PPyT) to the Exon Splicing Silencer 2p (ESS2p) element and to a conserved downstream linker. The 3D structure as determined by NMR spectroscopy further reveals that the A3 consensus cleavage site is embedded within a unique stereochemical environment within the apical loop, where it is surrounded by alternating base-base interactions. Despite being described as a receptor for hnRNP H, the ESS2p element is sequestered by base pairing to the 3' end of the PPyT and within this context it cannot form a stable complex with hnRNP H. By comparison, hnRNP A1 directly binds to the A3 consensus cleavage site located within the apical loop, suggesting that it can directly modulate U2AF assembly. Sequence mutations designed to destabilize the PPyT:ESS2p helix results in an increase usage of ssA3 within HIV-infected cells, consistent with the PPyT becoming more accessible for U2AF recognition. Additional mutations introduced into the downstream ESS2 element synergize with ESS2p to cause further increases in ssA3 usage. When taken together, our work provides a unifying picture by which cis regulatory sequences, splicing auxiliary factors and RNA structure cooperate to provide stringent control over ssA3. We describe this as the pair-and-lock mechanism to restrict access of the PPyT, and posit that it operates to regulate a subset of the heterogenous structures encompassing the ssA3 regulatory locus.

Can Second Coordination Sphere and Long-Range Interactions Modulate Hydrogen Atom Transfer in a Non-Heme Fe(II)-Dependent Histone Demethylase?

Fe(II)-dependent oxygenases employ hydrogen atom transfer (HAT) to produce a myriad of products. Understanding how such enzymes use dynamic processes beyond the immediate vicinity of the active site to control the selectivity and efficiency of HAT is important for metalloenzyme engineering; however, obtaining such knowledge by experiments is challenging. This study develops a computational framework for identifying second coordination sphere (SCS) and especially long-range (LR) residues relevant for catalysis through dynamic cross-correlation analysis (DCCA) using the human histone demethylase PHF8 (KDM7B) as a model oxygenase. Furthermore, the study explores the mechanistic pathways of influence of the SCS and LR residues on the HAT reaction. To demonstrate the plausibility of the approach, we investigated the effect of a PHF8 F279S clinical mutation associated with X-linked intellectual disability, which has been experimentally shown to ablate PHF8-catalyzed demethylation. In agreement, the molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) studies showed a change in the H3_1-14K9me2 substrate orientation and an increased HAT barrier. We systematically analyzed the pathways by which the identified SCS and LR residues may influence HAT by exploring changes in H3K9me2 substrate orientation, interdomain correlated motions, HAT transition state stabilization, reaction energetics, electron transfer mechanism, and alterations in the intrinsic electric field of PHF8. Importantly, SCS and LR variations decrease key motions of α9-α12 of the JmjC domain toward the Fe(IV)-center that are associated with tighter binding of the H3_1-14K9me2 substrate. SCS and LR residues alter the intrinsic electric field of the enzyme along the reaction coordinate and change the individual energetic contributions of residues toward TS stabilization. The overall results suggest that DCCA can indeed identify non-active-site residues relevant for catalysis. The substitutions of such dynamically correlated residues might be used as a tool to tune HAT in non-heme Fe(II)- and 2OG-dependent enzymes.

Identification of 11-Hydroxytephrosin and Torosaflavone A as Potential Inhibitors of 3-Phosphoinositide-Dependent Protein Kinase 1 (PDPK1): Toward Anticancer Drug Discovery

The 3-phosphoinositide-dependent protein kinase 1 (PDPK1) has a significant role in cancer progression and metastasis as well as other inflammatory disorders, and has been proposed as a promising therapeutic target for several malignancies. In this work, we conducted a systematic virtual screening of natural compounds from the IMPPAT database to identify possible PDPK1 inhibitors. Primarily, the Lipinski rules, ADMET, and PAINS filter were applied and then the binding affinities, docking scores, and selectivity were carried out to find effective hits against PDPK1. Finally, we identified two natural compounds, 11-Hydroxytephrosin and Torosaflavone A, bearing substantial affinity with PDPK1. Both compounds showed drug-likeness as predicted by the ADMET analysis and their physicochemical parameters. These compounds preferentially bind to the ATP-binding pocket of PDPK1 and interact with functionally significant residues. The conformational dynamics and complex stability of PDPK1 with the selected compounds were then studied using interaction analysis and molecular dynamics (MD) simulations for 100 ns. The simulation results revealed that PDPK1 forms stable docked complexes with the elucidated compounds. The findings show that the newly discovered 11-Hydroxytephrosin and Torosaflavone A bind to PDPK1 in an ATP-competitive manner, suggesting that they could one day be used as therapeutic scaffolds against PDPK1-associated diseases including cancer.

Conformation Dynamics of Single Polymer Strands in Solution

Conformational changes in macromolecules significantly affect their functions and assembly into high-level structures. Despite advances in theoretical and experimental studies, investigations into the intrinsic conformational variations and dynamic motions of single macromolecules remain challenging. Here, liquid-phase transmission electron microscopy enables the real-time tracking of single-chain polymers. Imaging linear polymers, synthetically dendronized with conjugated aromatic groups, in organic solvent confined within graphene liquid cells, directly exhibits chain-resolved conformational dynamics of individual semiflexible polymers. These experimental and theoretical analyses reveal that the dynamic conformational transitions of the single-chain polymer originate from the degree of intrachain interactions. In situ observations also show that such dynamics of the single-chain polymer are significantly affected by environmental factors, including surfaces and interfaces.

Principles of Alternating Access in LeuT-fold Transporters: Commonalities and Divergences

Found in all domains of life, transporters belonging to the LeuT-fold class mediate the import and exchange of hydrophilic and charged compounds such as amino acids, metals, and sugar molecules. Nearly two decades of investigations on the eponymous bacterial transporter LeuT have yielded a library of high-resolution snapshots of its conformational cycle linked by solution-state experimental data obtained from multiple techniques. In parallel, its topology has been observed in symporters and antiporters characterized by a spectrum of substrate specificities and coupled to gradients of distinct ions. Here we review and compare mechanistic models of transport for LeuT, its well-studied homologs, as well as functionally distant members of the fold, emphasizing the commonalities and divergences in alternating access and the corresponding energy landscapes. Our integrated summary illustrates how fold conservation, a hallmark of the LeuT fold, coincides with divergent choreographies of alternating access that nevertheless capitalize on recurrent structural motifs. In addition, it highlights the knowledge gap that hinders the leveraging of the current body of research into detailed mechanisms of transport for this important class of membrane proteins.

Integrating Conformational Dynamics and Perturbation-Based Network Modeling for Mutational Profiling of Binding and Allostery in the SARS-CoV-2 Spike Variant Complexes with Antibodies: Balancing Local and Global Determinants of Mutational Escape Mechanisms

In this study, we combined all-atom MD simulations, the ensemble-based mutational scanning of protein stability and binding, and perturbation-based network profiling of allosteric interactions in the SARS-CoV-2 spike complexes with a panel of cross-reactive and ultra-potent single antibodies (B1-182.1 and A23-58.1) as well as antibody combinations (A19-61.1/B1-182.1 and A19-46.1/B1-182.1). Using this approach, we quantify the local and global effects of mutations in the complexes, identify protein stability centers, characterize binding energy hotspots, and predict the allosteric control points of long-range interactions and communications. Conformational dynamics and distance fluctuation analysis revealed the antibody-specific signatures of protein stability and flexibility of the spike complexes that can affect the pattern of mutational escape. A network-based perturbation approach for mutational profiling of allosteric residue potentials revealed how antibody binding can modulate allosteric interactions and identified allosteric control points that can form vulnerable sites for mutational escape. The results show that the protein stability and binding energetics of the SARS-CoV-2 spike complexes with the panel of ultrapotent antibodies are tolerant to the effect of Omicron mutations, which may be related to their neutralization efficiency. By employing an integrated analysis of conformational dynamics, binding energetics, and allosteric interactions, we found that the antibodies that neutralize the Omicron spike variant mediate the dominant binding energy hotpots in the conserved stability centers and allosteric control points in which mutations may be restricted by the requirements of the protein folding stability and binding to the host receptor. This study suggested a mechanism in which the patterns of escape mutants for the ultrapotent antibodies may not be solely determined by the binding interaction changes but are associated with the balance and tradeoffs of multiple local and global factors, including protein stability, binding affinity, and long-range interactions.

Metalloid-Cluster Ligands Enabling Stable and Active FeN4 -Ten Motifs for the Oxygen Reduction Reaction

In nature, the oxygen reduction reaction (ORR) is catalyzed by cytochrome P450 (CYP) enzymes containing heme iron centers with an axial thiolate ligand (FeN₄ -S), which are among the most finely developed catalysts by natural selection. However, the exceptional ORR activity and selectivity of CYP enzymes originate from their non-rigid and self-adaptive coordination network with molecular ligands, which sacrifices the stability of the active motifs under electrochemical reaction conditions. Here, a design strategy to circumvent this dilemma by incorporating Fe-N₄ motifs into carbon matrices instead of the protein scaffold and replacing the axial molecular thiolate ligand with a stable tellurium cluster (Te_n ) is demonstrated. Theoretical calculations indicate a moderate interaction between Fe 3d and Te 5p orbitals once n > 2, allowing the FeTe bond to dynamically change its strength to adaptively facilitate the intermediate steps during the ORR process, which renders FeN₄ -Te_n active sites with superior ORR activity. This adaptive behavior mimics the conformational dynamics of an enzyme during the reaction, but retains the stability nature as a heterogeneous catalyst. The experiments validate that the as-designed catalyst with a characterized FeN₄ -Te_n structure outperforms the commercial Pt/C catalyst both on activity and stability.