Search

Observing a single biomolecule performing its function is fundamental in biophysics as it provides important information for elucidating the mechanism. High-speed atomic force microscopy (HS-AFM) is a unique and powerful technique that allows the observation of biomolecular motion in a near-native environment. However, the spatial resolution of HS-AFM is limited by the physical size of the cantilever tip, which restricts the ability to obtain atomic details of molecules. In this study, we propose a novel computational algorithm designed to derive atomistic models of conformational dynamics from AFM images. Our method uses normal-mode analysis to describe the expected motions of the molecule, allowing these motions to be represented with a limited number of coordinates. This approach mitigates the problem of overinterpretation inherent in the analysis of AFM images with limited resolution. We demonstrate the effectiveness of our algorithm, NMFF-AFM, using synthetic data sets for three proteins that undergo significant conformational changes. NMFF-AFM is a fast and user-friendly program that requires minimal setup and has the potential to be a valuable tool for biophysical studies using HS-AFM.

The session "Data Science for Integrated Dynamic Structural Biology" was a notable success at the joint congress of the 21st IUPAB and the 62nd BSJ (Biophysics Society of Japan). This session included four invited talks and one contributed talk, which together delved into recent advancements in computational methods integrating physics, experimental data, and bioinformatics to study the structure and dynamic properties of molecular assemblies and their interactions., (© International Union for Pure and Applied Biophysics (IUPAB) and Springer-Verlag GmbH Germany, part of Springer Nature 2024.)

The development of near-infrared (NIR) fluorophores that have both excellent chemical stability and photostability, as well as efficient cell permeability, is highly demanded. In this study, we present phospha-rhodamine (POR) dyes which display significantly improved performance for protein labeling. This is achieved by incorporating a 2-carboxy-3-benzothiophenyl group at the 9-position of the xanthene scaffold. The resulting cis and trans isomers were successfully isolated and structurally characterized using X-ray diffraction. The HaloTag ligand conjugates of the two isomers exhibited different staining abilities in live cells. While the cis isomer showed non-specific accumulation on the organelle membranes, the trans isomer selectively labeled the HaloTag-fused proteins, enabling the long-term imaging of cell division and the 5-color imaging of cell organelles. Molecular dynamics simulations of the HaloTag ligand conjugates within the lipid membrane suggested that the cis isomer is more prone to forming oligomers in the membrane. In contrast, the oligomerization of the trans isomer is effectively suppressed by its interaction with the lipid molecules. By taking advantage of the superior labeling performance of the trans isomer and its NIR-emissive properties, multi-color time-lapse super-resolution 3D imaging, namely super-resolution 5D-imaging, of the interconnected network between the endoplasmic reticulum and microtubules was achieved in living cells., (© 2024 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

Single-particle imaging using X-ray free-electron lasers (XFELs) is a promising technique for observing nanoscale biological samples under near-physiological conditions. However, as the sample's orientation in each diffraction pattern is unknown, advanced algorithms are required to reconstruct the 3D diffraction intensity volume and subsequently the sample's density model. While most approaches perform 3D reconstruction via determining the orientation of each diffraction pattern, a correlation-based approach utilizes the averaged spatial correlations of diffraction intensities over all patterns, making it well suited for processing experimental data with a poor signal-to-noise ratio of individual patterns. Here, a method is proposed to determine the 3D structure of a sample by analyzing the double, triple and quadruple spatial correlations in diffraction patterns. This ab initio method can reconstruct the basic shape of an irregular unsymmetric 3D sample without requiring any prior knowledge of the sample. The impact of background and noise on correlations is investigated and corrected to ensure the success of reconstruction under simulated experimental conditions. Additionally, the feasibility of using the correlation-based approach to process incomplete partial diffraction patterns is demonstrated. The proposed method is a variable addition to existing algorithms for 3D reconstruction and will further promote the development and adoption of XFEL single-particle imaging techniques., (open access.)

Flexible fitting based on molecular dynamics simulation is a technique for structure modeling from cryo-EM data. It has been utilized for nearly two decades, and while cryo-EM resolution has improved significantly, it remains a powerful approach that can provide structural and dynamical insights that are not directly accessible from experimental data alone. Molecular dynamics simulations provide a means to extract atomistic details of conformational changes that are encoded in cryo-EM data and can also assist in improving the quality of structural models. Additionally, molecular dynamics simulations enable the characterization of conformational heterogeneity in cryo-EM data. We will summarize the advancements made in these techniques and highlight recent developments in this field., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

This article presents an original approach for extracting atomic-resolution landscapes of continuous conformational variability of biomolecular complexes from cryo electron microscopy (cryo-EM) single particle images. This approach is based on a new 3D-to-2D flexible fitting method, which uses molecular dynamics (MD) simulation and is embedded in an iterative conformational-landscape refinement scheme. This new approach is referred to as MDSPACE, which stands for Molecular Dynamics simulation for Single Particle Analysis of Continuous Conformational hEterogeneity. The article describes the MDSPACE approach and shows its performance using synthetic and experimental datasets., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Single-particle analysis using x-ray free-electron lasers (XFELs) is a novel method for obtaining structural information of samples in a state close to nature. In particular, it is suitable for observing the inner structure of large biomolecules by taking advantage of the high transmittance of x-rays. However, systematic studies on the resolution achievable for large molecules are lacking. In this study, the molecular size dependence of the resolution of a three-dimensional (3D) structure resulting from XFEL single-particle reconstruction is evaluated using synthetic data. Evidently, 3D structures of larger molecules can be restored with higher detail (defined relative to the molecular sizes) than smaller ones; however, reconstruction with high absolute resolution (defined in nm -1 ) is challenging. Our results provide useful information for the experimental design of 3D structure reconstruction using coherent x-ray diffraction patterns of single-particles., Competing Interests: The authors have no conflicts to disclose., (© 2023 Author(s).)

SWEET sucrose transporters play important roles in the allocation of sucrose in plants. Some SWEETs were shown to also mediate transport of the plant growth regulator gibberellin (GA). The close physiological relationship between sucrose and GA raised the questions of whether there is a functional connection and whether one or both of the substrates are physiologically relevant. To dissect these two activities, molecular dynamics were used to map the binding sites of sucrose and GA in the pore of SWEET13 and predicted binding interactions that might be selective for sucrose or GA. Transport assays confirmed these predictions. In transport assays, the N76Q mutant had 7x higher relative GA 3 activity, and the S142N mutant only transported sucrose. The impaired pollen viability and germination in sweet13;14 double mutants were complemented by the sucrose-selective SWEET13 S142N , but not by the SWEET13 N76Q mutant, indicating that sucrose is the physiologically relevant substrate and that GA transport capacity is dispensable in the context of male fertility. Therefore, GA supplementation to counter male sterility may act indirectly via stimulating sucrose supply in male sterile mutants. These findings are also relevant in the context of the role of SWEETs in pathogen susceptibility.

The mammalian cryptochrome isoforms, CRY1 and CRY2, are core circadian clock regulators that work redundantly. Recent studies revealed distinct roles of these closely related homologs in clock output pathways. Isoform-selective control of CRY1 and CRY2 is critical for further understanding their redundant and distinct roles. KL001 was the first identified small-molecule CRY modulator that activates both CRY1 and CRY2. SHP656 is an orally available KL001 derivative and has shown efficacy in blood glucose control and inhibition of glioblastoma stem cell (GSC) growth in animal models. However, CRY isoform selectivity of SHP656 was uncharacterized, limiting understanding of the roles of CRY1 and CRY2. Here, we report the elucidation of CRY2 selectivity of SHP656. SHP656 lengthened cellular circadian period in a CRY2-dependent manner and selectively interacted with CRY2. By determining the X-ray crystal structure of CRY2 in complex with SHP656 and performing molecular dynamics simulations, we elucidated compound interaction mechanisms. SHP656 binding was compatible with the intrinsic CRY2 gatekeeper W417 "in" orientation and also a close "further in" conformation. Perturbation of W417 interaction with the lid loop resulted in a reduced effect of SHP656 on CRY2, supporting an important role of gatekeeper orientation in isoform selectivity. We also identified the R form of SHP656 (called SHP1703) as the active isomer. Treatment with SHP1703 effectively reduced GSC viability. Our results suggest a direct role of CRY2 in glioblastoma antitumorigenesis and provide a rationale for the selective modulation of CRY isoforms in the therapeutic treatment of glioblastoma and other circadian clock-related diseases.

X-ray free-electron laser (XFEL) is the latest generation of the X-ray source that could become an invaluable technique in structural biology. XFEL has ultrashort pulse duration, extreme peak brilliance, and high spatial coherence, which could enable the observation of the biological molecules in near nature state at room temperature without crystallization. However, for biological systems, due to their low diffraction power and complexity of sample delivery, experiments and data analysis are not straightforward, making it extremely challenging to reconstruct three-dimensional (3D) structures from single particle XFEL data. Given the current limitations to the amount and resolution of the data from such XFEL experiments, we propose a new hybrid approach for characterizing biomolecular conformational transitions by using a single 2D low-resolution XFEL diffraction pattern in combination with another known conformation. In our method, we represent the molecular structure with a coarse-grained model, the Gaussian mixture model, to describe large conformational transitions from low-resolution XFEL data. We obtain plausible 3D structural models that are consistent with the XFEL diffraction pattern by deforming an initial structural model to maximize the similarity between the target pattern and the simulated diffraction patterns from the candidate models. We tested the proposed algorithm on two biomolecules of different sizes with different complexities of conformational transitions, adenylate kinase, and elongation factor 2, using synthetic XFEL data. The results show that, with the proposed algorithm, we can successfully describe the conformational transitions by flexibly fitting the coarse-grained model of one conformation to become consistent with an XFEL diffraction pattern simulated from another conformation. In addition, we showed that the incident beam orientation has some effect on the accuracy of the 3D structure modeling and discussed the reasons for the inaccuracies for certain orientations. The proposed method could serve as an alternative approach for retrieving information on 3D conformational transitions from the XFEL diffraction patterns to interpret experimental data. Since the molecules are represented by Gaussian kernels and no atomic structure is needed in principle, such a method could also be used as a tool to seek initial models for 3D reconstruction algorithms., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Asi, Dasgupta, Nagai, Miyashita and Tama.)

In bifidobacteria, phosphoketolase (PKT) plays a key role in the central hexose fermentation pathway called "bifid shunt." The three-dimensional structure of PKT from Bifidobacterium longum with co-enzyme thiamine diphosphate (ThDpp) was determined at 2.1 Å resolution by cryo-EM single-particle analysis using 196,147 particles to build up the structural model of a PKT octamer related by D 4 symmetry. Although the cryo-EM structure of PKT was almost identical to the X-ray crystal structure previously determined at 2.2 Å resolution, several interesting structural features were observed in the cryo-EM structure. Because this structure was solved at relatively high resolution, it was observed that several amino acid residues adopt multiple conformations. Among them, Q546-D547-H548-N549 (the QN-loop) demonstrate the largest structural change, which seems to be related to the enzymatic function of PKT. The QN-loop is at the entrance to the substrate binding pocket. The minor conformer of the QN-loop is similar to the conformation of the QN-loop in the crystal structure. The major conformer is located further from ThDpp than the minor conformer. Interestingly, the major conformer in the cryo-EM structure of PKT resembles the corresponding loop structure of substrate-bound Escherichia coli transketolase. That is, the minor and major conformers may correspond to "closed" and "open" states for substrate access, respectively. Moreover, because of the high-resolution analysis, many water molecules were observed in the cryo-EM structure of PKT. Structural features of the water molecules in the cryo-EM structure are discussed and compared with water molecules observed in the crystal structure., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Protein kinases are responsible for healthy cellular processes and signalling pathways, and their dysfunction is the basis of many pathologies. There are numerous small molecule inhibitors of protein kinases that systemically regulate dysfunctional signalling processes. However, attaining selectivity in kinase inhibition within the complex human kinome is still a challenge that inspires unconventional approaches. One of those approaches is photopharmacology, which uses light-controlled bioactive molecules to selectively activate drugs only at the intended space and time, thereby avoiding side effects outside of the irradiated area. Still, in the context of kinase inhibition, photopharmacology has thus far been rather unsuccessful in providing light-controlled drugs. Here, we present the discovery and optimisation of a photoswitchable inhibitor of casein kinase 1δ (CK1δ), important for the control of cell differentiation, circadian rhythm, DNA repair, apoptosis, and numerous other signalling processes. Varying the position at which the light-responsive azobenzene moiety has been introduced into a known CK1δ inhibitor, LH846, revealed the preferred regioisomer for efficient photo-modulation of inhibitory activity, but the photoswitchable inhibitor suffered from sub-optimal (photo)chemical properties. Replacement of the bis-phenyl azobenzene group with the arylazopyrazole moiety yielded a superior photoswitch with very high photostationary state distributions, increased solubility and a 10-fold difference in activity between irradiated and thermally adapted samples. The reasons behind those findings are explored with molecular docking and molecular dynamics simulations. Results described here show how the evaluation of privileged molecular architecture, followed by the optimisation of the photoswitchable unit, is a valuable strategy for the challenging design of the photoswitchable kinase inhibitors.

Atomic models of cryo electron microscopy (cryo-EM) maps of biomolecular conformations are often obtained by flexible fitting of the maps with available atomic structures of other conformations (e.g., obtained by X-ray crystallography). This article presents a new flexible fitting method, NMMD, which combines normal mode analysis (NMA) and molecular dynamics simulation (MD). Given an atomic structure and a cryo-EM map to fit, NMMD simultaneously estimates global atomic displacements based on NMA and local displacements based on MD. NMMD was implemented by modifying EMfit, a flexible fitting method using MD only, in GENESIS 1.4. As EMfit, NMMD can be run with replica exchange umbrella sampling procedure. The new method was tested using a variety of EM maps (synthetic and experimental, with different noise levels and resolutions). The results of the tests show that adding normal modes to MD-based fitting makes the fitting faster (40% in average) and, in the majority of cases, more accurate., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

Pannexin (PANX) family proteins form large-pore channels that mediate purinergic signaling. We analyzed the cryo-EM structures of human PANX1 in lipid nanodiscs to elucidate the gating mechanism and its regulation by the amino terminus in phospholipids. The wild-type channel has an amino-terminal funnel in the pore, but in the presence of the inhibitor probenecid, a cytoplasmically oriented amino terminus and phospholipids obstruct the pore. Functional analysis using whole-cell patch-clamp and oocyte voltage clamp showed that PANX1 lacking the amino terminus did not open and had a dominant negative effect on channel activity, thus confirming that the amino-terminal domain played an essential role in channel opening. These observations suggest that dynamic conformational changes in the amino terminus of human PANX1 are associated with lipid movement in and out of the pore. Moreover, the data provide insight into the gating mechanism of PANX1 and, more broadly, other large-pore channels.

X-ray free-electron laser (XFEL) scattering promises to probe single biomolecular complexes without crystallization, enabling the study of biomolecular structures under near-physiological conditions at room temperature. However, such structural determination of biomolecules is extremely challenging thus far. In addition to the large numbers of diffraction patterns required, the orientation of each diffraction pattern needs to be accurately estimated and the missing phase information needs to be recovered for three-dimensional (3D) structure reconstruction. Given the current limitations to the amount and resolution of the data available from single-particle XFEL scattering experiments, we propose an alternative approach to find plausible 3D biological shapes from a limited number of diffraction patterns to serve as a starting point for further analyses. In our proposed strategy, small sets of input (e.g., five) XFEL diffraction patterns were matched against a library of diffraction patterns simulated from 1628 electron microscopy (EM) models to find potential matching 3D models that are consistent with the input diffraction patterns. This approach was tested for three example cases: EMD-3457 ( Thermoplasma acidophilum 20S proteasome), EMD-5141 ( Escherichia coli 70S ribosome complex), and EMD-5152 (budding yeast Nup84 complex). We observed that choosing the best strategy to define matching regions on the diffraction patterns is critical for identifying correctly matching diffraction patterns. While increasing the number of input diffraction patterns improved the matches in some cases, we found that the resulting matches are more dependent on the uniqueness or complexity of the shape as captured in the individual input diffraction patterns and the availability of a similar 3D biological shape in the search library. The protocol could be useful for finding candidate models for a limited amount of low-resolution data, even when insufficient for reconstruction, performing a quick exploration of new data upon collection, and the analysis of the conformational heterogeneity of the particle of interest as captured within the diffraction patterns.

ClpB belongs to the cellular disaggretase machinery involved in rescuing misfolded or aggregated proteins during heat or other cellular shocks. The function of this protein relies on the interconversion between different conformations in its native condition. A recent high-speed-atomic-force-microscopy (HS-AFM) experiment on ClpB from Thermus thermophilus shows four predominant conformational classes, namely, open, closed, spiral, and half-spiral. Analyses of AFM images provide only partial structural information regarding the molecular surface, and thus computational modeling of three-dimensional (3D) structures of these conformations should help interpret dynamical events related to ClpB functions. In this study, we reconstruct 3D models of ClpB from HS-AFM images in different conformational classes. We have applied our recently developed computational method based on a low-resolution representation of 3D structure using a Gaussian mixture model, combined with a Monte-Carlo sampling algorithm to optimize the agreement with target AFM images. After conformational sampling, we obtained models that reflect conformational variety embedded within the AFM images. From these reconstructed 3D models, we described, in terms of relative domain arrangement, the different types of ClpB oligomeric conformations observed by HS-AFM experiments. In particular, we highlighted the slippage of the monomeric components around the seam. This study demonstrates that such details of information, necessary for annotating the different conformational states involved in the ClpB function, can be obtained by combining HS-AFM images, even with limited resolution, and computational modeling., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Dasgupta, Miyashita, Uchihashi and Tama.)

The circadian clock is a biological timekeeper that operates through transcription-translation feedback loops in mammals. Cryptochrome 1 (CRY1) and Cryptochrome 2 (CRY2) are highly conserved core clock components having redundant and distinct functions. We recently identified the CRY1- and CRY2-selective compounds KL101 and TH301, respectively, which provide useful tools for the exploration of isoform-selective CRY regulation. However, intrinsic differences in the compound-binding FAD (flavin adenine dinucleotide) pockets between CRY1 and CRY2 are not well understood, partly because of nonoptimal properties of previously reported apo form structures in this particular region constituted by almost identical sequences. Here, we show unliganded CRY1 and CRY2 crystal structures with well-defined electron densities that are largely free of crystal contacts at the FAD pocket and nearby lid loop. We revealed conformational isomerism in key residues. In particular, CRY1 W399 and corresponding CRY2 W417 in the FAD pocket had distinct conformations ("out" for CRY1 and "in" for CRY2) by interacting with the lid loop residues CRY1 Q407 and CRY2 F424, respectively, resulting in different overall lid loop structures. Molecular dynamics simulations supported that these conformations were energetically favorable to each isoform. Isoform-selective compounds KL101 and TH301 preferred intrinsic "out" and "in" conformations of the tryptophan residue in CRY1 and CRY2, respectively, while the nonselective compound KL001 fit to both conformations. Mutations of lid loop residues designed to perturb their isoform-specific interaction with the tryptophan resulted in reversed responses of CRY1 and CRY2 to KL101 and TH301. We propose that these intrinsic structural differences of CRY1 and CRY2 can be targeted for isoform-selective regulation., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

The circadian clock controls daily rhythms of physiological processes. The presence of the clock mechanism throughout the body is hampering its local regulation by small molecules. A photoresponsive clock modulator would enable precise and reversible regulation of circadian rhythms using light as a bio-orthogonal external stimulus. Here we show, through judicious molecular design and state-of-the-art photopharmacological tools, the development of a visible light-responsive inhibitor of casein kinase I (CKI) that controls the period and phase of cellular and tissue circadian rhythms in a reversible manner. The dark isomer of photoswitchable inhibitor 9 exhibits almost identical affinity towards the CKIα and CKIδ isoforms, while upon irradiation it becomes more selective towards CKIδ, revealing the higher importance of CKIδ in the period regulation. Our studies enable long-term regulation of CKI activity in cells for multiple days and show the reversible modulation of circadian rhythms with a several hour period and phase change through chronophotopharmacology.

CRY1 and CRY2 proteins are highly conserved components of the circadian clock that controls daily physiological rhythms. Disruption of CRY functions are related to many diseases, including circadian sleep phase disorder. Development of isoform-selective and spatiotemporally controllable tools will facilitate the understanding of shared and distinct functions of CRY1 and CRY2. Here, we developed CRY1-selective compounds that enable light-dependent manipulation of the circadian clock. From phenotypic chemical screening in human cells, we identified benzophenone derivatives that lengthened the circadian period. These compounds selectively interacted with the CRY1 photolyase homology region, resulting in activation of CRY1 but not CRY2. The benzophenone moiety rearranged a CRY1 region called the "lid loop" located outside of the compound-binding pocket and formed a unique interaction with Phe409 in the lid loop. Manipulation of this key interaction was achieved by rationally designed replacement of the benzophenone with a switchable azobenzene moiety whose cis-trans isomerization can be controlled by light. The metastable cis form exhibited sufficiently high half-life in aqueous solutions and structurally mimicked the benzophenone unit, enabling reversible period regulation over days by cellular irradiation with visible light. This study revealed an unprecedented role of the lid loop in CRY-compound interaction and paves the way for spatiotemporal regulation of CRY1 activity by photopharmacology for molecular understanding of CRY1-dependent functions in health and disease.

CRY1 and CRY2 are essential components of the circadian clock controlling daily physiological rhythms. Accumulating evidences indicate distinct roles of these highly homologous proteins, in addition to redundant functions. Therefore, the development of isoform-selective compounds represents an effective approach towards understanding the similarities and differences of CRY1 and CRY2 by controlling each isoform individually. We conducted phenotypic screenings of circadian clock modulators, and identified KL101 and TH301 that selectively stabilize CRY1 and CRY2, respectively. Crystal structures of CRY-compound complexes revealed conservation of compound-binding sites between CRY1 and CRY2. We further discovered a unique mechanism underlying compound selectivity in which the disordered C-terminal region outside the pocket was required for the differential effects of KL101 and TH301 against CRY isoforms. By using these compounds, we found a new role of CRY1 and CRY2 as enhancers of brown adipocyte differentiation, providing the basis of CRY-mediated regulation of energy expenditure.

Cryo-electron microscopy (cryo-EM) single-particle analysis has come a long way in achieving atomic-level resolution when imaging biomolecules. To obtain the best possible three-dimensional (3D) structure in cryo-EM, many parameters have to be carefully considered. Here we address the often-overlooked parameter of the pixel size, which describes the magnification of the image produced by the experiment. While efforts are made to refine and validate this parameter in the analysis of cryo-EM experimental data, there is no systematic protocol in place. Since the pixel size parameter can have an impact on the resolution and accuracy of a cryo-EM map, and the atomic resolution 3D structure models derived from it, we propose a computational protocol to estimate the appropriate pixel size parameter. In our protocol, we fit and refine atomic structures against cryo-EM maps at multiple pixel sizes. The resulting fitted and refined structures are evaluated using the GOAP (generalized orientation-dependent, all-atom statistical potential) score, which we found to perform better than other commonly used functions, such as Molprobity and the correlation coefficient from refinement. Finally, we describe the efficacy of this protocol in retrieving appropriate pixel sizes for several examples; simulated data based on yeast elongation factor 2 and experimental data from Gro-EL chaperone, beta-galactosidase, and the TRPV1 ion channel.

The structural and dynamical characterization of biomolecules holds central importance in the endeavor to understand the molecular mechanisms regulating living systems. However, owing to the inherent heterogeneity of biomolecular interactions within cells, it is often difficult to understand the overall structure and dynamics of biomolecules using any experimental method in isolation. In this regard, hybrid methods that combine data from multiple experiments to generate a comprehensive model of biomolecular complexes have gained prominence in the last few years. In this article, we discuss the advancements in hybrid methods, with a particular focus on the role of computation in their development and application. We further outline the future directions that hybrid methods are likely to take, regarding the advancements in techniques such as X-ray free-electron laser single- particle imaging, and electron cryo-tomography. Finally, we conclude the review by highlighting the future goals of broader consensus and collaboration within the integrative/hybrid structural biology community and for disseminating the data generated by hybrid modeling efforts., (Copyright © 2020 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Mammalian circadian rhythms are generated by a transcription-based feedback loop in which CLOCK:BMAL1 drives transcription of its repressors (PER1/2, CRY1/2), which ultimately interact with CLOCK:BMAL1 to close the feedback loop with ~24 hr periodicity. Here we pinpoint a key difference between CRY1 and CRY2 that underlies their differential strengths as transcriptional repressors. Both cryptochromes bind the BMAL1 transactivation domain similarly to sequester it from coactivators and repress CLOCK:BMAL1 activity. However, we find that CRY1 is recruited with much higher affinity to the PAS domain core of CLOCK:BMAL1, allowing it to serve as a stronger repressor that lengthens circadian period. We discovered a dynamic serine-rich loop adjacent to the secondary pocket in the photolyase homology region (PHR) domain that regulates differential binding of cryptochromes to the PAS domain core of CLOCK:BMAL1. Notably, binding of the co-repressor PER2 remodels the serine loop of CRY2, making it more CRY1-like and enhancing its affinity for CLOCK:BMAL1., Competing Interests: JF, AS, CS, AM, PH, CR, LN, MT, GP, ST, NZ, GL, TH, FT, CP No competing interests declared, (© 2020, Fribourgh et al.)

Background: Tim21, a subunit of a highly dynamic translocase of the inner mitochondrial membrane (TIM23) complex, translocates proteins by interacting with subunits in the translocase of the outer membrane (TOM) complex and Tim23 channel in the TIM23 complex. A loop segment in Tim21, which is in close proximity of the binding site of Tim23, has different conformations in X-ray, NMR and new crystal contact-free space (CCFS) structures. MD simulations can provide information on the structure and dynamics of the loop in solution., Methods: The conformational ensemble of the loop was characterized using loop modeling and molecular dynamics (MD) simulations., Results: MD simulations confirmed mobility of the loop. Multidimensional scaling and clustering were used to characterize the dynamic conformational ensemble of the loop. Free energy landscape showed that the CCFS crystal structure occupied a low energy region as compared to the conventional X-ray crystal structure. Analysis of crystal packing indicates that the CCFS provides larger conformational space for the motions of the loop., Conclusions: Our work reported the conformational ensemble of the loop in solution, which is in agreement with the structure obtained from CCFS approach. The combination of the experimental techniques and computational methods is beneficial for studying highly flexible regions of proteins., General Significance: Computational methods, such as loop modeling and MD simulations, have proved to be useful for studying conformational flexibility of proteins. These methods in integration with experimental techniques such as CCFS has the potential to transform the studies on flexible regions of proteins., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Background: In protein crystals, flexible loops are frequently deformed by crystal contacts, whereas in solution, the large motions result in the poor convergence of such flexible loops in NMR structure determinations. We need an experimental technique to characterize the structural and dynamic properties of intrinsically flexible loops of protein molecules., Methods: We designed an intended crystal contact-free space (CCFS) in protein crystals, and arranged the flexible loop of interest in the CCFS. The yeast Tim 21 protein was chosen as the model protein, because one of the loops (loop 2) is distorted by crystal contacts in the conventional crystal., Results: Yeast Tim21 was fused to the MBP protein by a rigid α-helical linker. The space created between the two proteins was used as the CCFS. The linker length provides adjustable freedom to arrange loop 2 in the CCFS. We re-determined the NMR structure of yeast Tim21, and conducted MD simulations for comparison. Multidimensional scaling was used to visualize the conformational similarity of loop 2. We found that the crystal contact-free conformation of loop 2 is located close to the center of the ensembles of the loop 2 conformations in the NMR and MD structures., Conclusions: Loop 2 of yeast Tim21 in the CCFS adopts a representative, dominant conformation in solution., General Significance: No single powerful technique is available for the characterization of flexible structures in protein molecules. NMR analyses and MD simulations provide useful, but incomplete information. CCFS crystallography offers a third route to this goal., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Background: Atomic Force Microscopy (AFM) is an experimental technique to study structure-function relationship of biomolecules. AFM provides images of biomolecules at nanometer resolution. High-speed AFM experiments produce a series of images following dynamics of biomolecules. To further understand biomolecular functions, information on three-dimensional (3D) structures is beneficial., Method: We aim to recover 3D information from an AFM image by computational modeling. The AFM image includes only low-resolution representation of a molecule; therefore we represent the structures by a coarse grained model (Gaussian mixture model). Using Monte-Carlo sampling, candidate models are generated to increase similarity between AFM images simulated from the models and target AFM image., Results: The algorithm was tested on two proteins to model their conformational transitions. Using a simulated AFM image as reference, the algorithm can produce a low-resolution 3D model of the target molecule. Effect of molecular orientations captured in AFM images on the 3D modeling performance was also examined and it is shown that similar accuracy can be obtained for many orientations., Conclusions: The proposed algorithm can generate 3D low-resolution protein models, from which conformational transitions observed in AFM images can be interpreted in more detail., General Significance: High-speed AFM experiments allow us to directly observe biomolecules in action, which provides insights on biomolecular function through dynamics. However, as only partial structural information can be obtained from AFM data, this new AFM based hybrid modeling method would be useful to retrieve 3D information of the entire biomolecule., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Single-particle analysis (SPA) by X-ray free electron laser (XFEL) is a novel method that can observe biomolecules and living tissue that are difficult to crystallize in a state close to nature. To reconstruct three-dimensional (3D) molecular structure from two-dimensional (2D) XFEL diffraction patterns, we have to estimate the incident beam angle to the molecule for each pattern to assemble the 3D-diffraction intensity distribution using interpolation, and retrieve the phase information. In this study, we investigated the optimal parameter sets to assemble the 3D-diffraction intensity distribution from simulated 2D-diffraction patterns of ribosome. In particular, we examined how the parameters need to be adjusted for diffraction patterns with different binning sizes and beam intensities to obtain the highest resolution of molecular structure phase retrieved from the 3D-diffraction intensity. We found that resolution of restored molecular structure is sensitive to the interpolation parameters. Using the optimal parameter set, a linear oversampling ratio of around four is found to be sufficient for correct angle estimation and phase retrieval from the diffraction patterns of SPA by XFEL., Competing Interests: Conflicts of Interest All authors declare that they have no conflict of interest., (2019 © The Biophysical Society of Japan.)

Circadian clocks, biological timekeepers that are present in almost every cell of our body, are complex systems whose disruption is connected to various diseases. Controlling cellular clock function with high temporal resolution in an inducible manner would yield an innovative approach for the circadian rhythm regulation. In the present study, we present structure-guided incorporation of photoremovable protecting groups into a circadian clock modifier, longdaysin, which inhibits casein kinase I (CKI). Using photodeprotection by UV or visible light (400 nm) as the external stimulus, we have achieved quantitative and light-inducible control over the CKI activity accompanied by an accurate regulation of circadian period in cultured human cells and mouse tissues, as well as in living zebrafish. This research paves the way for the application of photodosing in achieving precise temporal control over the biological timing and opens the door for chronophotopharmacology to deeper understand the circadian clock system.

Cell fate in eukaryotes is controlled by mitogen-activated protein kinases (MAPKs) that translate external cues into cellular responses. In plants, two MAPKs-MPK3 and MPK6-regulate diverse processes of development, environmental response and immunity. However, the mechanism that bridges these shared signalling components with a specific target remains unresolved. Focusing on the development of stomata-epidermal valves that are essential for gas exchange and transpiration-here, we report that the basic helix-loop-helix protein SCREAM functions as a scaffold that recruits MPK3/6 to downregulate SPEECHLESS, a transcription factor that initiates stomatal cell lineages. SCREAM directly binds to MPK3/6 through an evolutionarily conserved, yet unconventional, bipartite motif. Mutations in this motif abrogate association, phosphorylation and degradation of SCREAM, unmask hidden non-redundancies between MPK3 and MPK6, and result in uncontrolled stomatal differentiation. Structural analyses of MPK6 with a resolution of 2.75 Å showed bipartite binding of SCREAM to MPK6 that is distinct from an upstream MAPKK. Our findings elucidate, at the atomic resolution, the mechanism that directly links extrinsic signals to transcriptional reprogramming during the establishment of stomatal cell fate, and highlight a unique substrate-binding mode adopted by plant MAPKs.

Cryo-cooling is routinely performed before x-ray diffraction image collection to reduce the damage to crystals due to ionizing radiation. It has been suggested that although backbone structures are usually very similar between room temperature and cryo-temperature, cryo-cooling may hamper biologically relevant dynamics. In this study, the crystal of Escherichia coli dihydrofolate reductase is studied with replica-exchange molecular dynamics simulation, and the results are compared with the crystal structure determined at cryo-temperature and room temperature with the time-averaged ensemble method. Although temperature dependence of unit cell compaction and root mean-square fluctuation of Cα is found in accord with experiment, it is found that the protein structure at low temperature can be more heterogeneous than the ensemble of structures reported by using the time-averaged ensemble method, encouraging further development of the time-averaged ensemble method and indicating that data should be examined carefully to avoid overinterpretation of one average structure., (Copyright © 2018 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Compounds targeting the circadian clock have been identified as potential treatments for clock-related diseases, including cancer. Our cell-based phenotypic screen revealed uncharacterized clock-modulating compounds. Through affinity-based target deconvolution, we identified GO289, which strongly lengthened circadian period, as a potent and selective inhibitor of CK2. Phosphoproteomics identified multiple phosphorylation sites inhibited by GO289 on clock proteins, including PER2 S693. Furthermore, GO289 exhibited cell type-dependent inhibition of cancer cell growth that correlated with cellular clock function. The x-ray crystal structure of the CK2α-GO289 complex revealed critical interactions between GO289 and CK2-specific residues and no direct interaction of GO289 with the hinge region that is highly conserved among kinases. The discovery of GO289 provides a direct link between the circadian clock and cancer regulation and reveals unique design principles underlying kinase selectivity.

Flexible fitting is a powerful technique to build the 3D structures of biomolecules from cryoelectron microscopy (cryo-EM) density maps. One popular method is a cross-correlation coefficient-based approach, where the molecular dynamics (MD) simulation is carried out with the biasing potential that includes the cross-correlation coefficient between the experimental and simulated density maps. Here, we propose efficient parallelization schemes for the calculation of the cross-correlation coefficient to accelerate flexible fitting. Our schemes are tested for small, medium, and large biomolecules using CPU and hybrid CPU + GPU architectures. The scheme for the atomic decomposition MD is suitable for small proteins such as Ca 2+ -ATPase with the all-atom Go model, while that for the domain decomposition MD is better for larger systems such as ribosome with the all-atom Go or the all-atom explicit solvent models. Our methods allow flexible fitting for various biomolecules with reasonable computational cost. This approach also connects high-resolution structure refinements with investigation of protein structure-function relationship., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

The translocase of the outer membrane (TOM) mediates the membrane permeation of mitochondrial matrix proteins. Tom20 is a subunit of the TOM complex and binds to the N-terminal region (ie, presequence) in mitochondrial matrix precursor proteins. Previous experimental studies indicated that the presequence recognition by Tom20 was achieved in a dynamic-equilibrium among multiple bound states of the α-helical presequence. Accordingly, the co-crystallization of Tom20 and a presequence peptide required a disulfide-bond cross-linking. A 3-residue spacer sequence (XAG) was inserted between the presequence and the anchoring Cys residue at the C-terminus to not disturb the movement of the presequence peptide in the binding site of Tom20. Two crystalline forms were obtained according to Ala or Tyr at the X position of the spacer sequence, which may reflect the dynamic-equilibrium of the presequence. Here, we have performed replica-exchange molecular dynamics (REMD) simulations to study the effect of disulfide-bond linker and single amino acid difference in the spacer region of the linker on the conformational dynamics of Tom20-presequence complex. Free energy and network analyses of the REMD simulations were compared against previous simulations of non-tethered system. We concluded that the disulfide-bond tethering did not strongly affect the conformational ensemble of the presequence peptide in the complex. Further investigation showed that the choice of Ala or Tyr at the X position did not affect the most distributions of the conformational ensemble of the presequence. The present study provides a rational basis for the disulfide-bond tethering to study the dynamics of weakly binding complexes., (© 2018 Wiley Periodicals, Inc.)

Protein structural biology came a long way since the determination of the first three-dimensional structure of myoglobin about six decades ago. Across this period, X-ray crystallography was the most important experimental method for gaining atomic-resolution insight into protein structures. However, as the role of dynamics gained importance in the function of proteins, the limitations of X-ray crystallography in not being able to capture dynamics came to the forefront. Computational methods proved to be immensely successful in understanding protein dynamics in solution, and they continue to improve in terms of both the scale and the types of systems that can be studied. In this review, we briefly discuss the limitations of X-ray crystallography in studying protein dynamics, and then provide an overview of different computational methods that are instrumental in understanding the dynamics of proteins and biomacromolecular complexes.

We explore the advantage of Gaussian mixture model (GMM) for interpretation of single particle diffraction patterns from X-ray free electron laser (XFEL) experiments. GMM approximates a biomolecular shape by the superposition of Gaussian distributions. As the Fourier transformation of GMM can be quickly performed, we can efficiently simulate XFEL diffraction patterns from approximated structure models. We report that the resolution that GMM can accurately reproduce is proportional to the cubic root of the number of Gaussians used in the modeling. This behavior can be attributed to the correspondence between the number of adjustable parameters in GMM and the amount of sampling points in diffraction space. Furthermore, GMMs can successfully be used to perform angular assignment and to detect conformational variation. These results demonstrate that GMMs serve as useful coarse-grained models for hybrid approach in XFEL single particle experiments.

Background: Advancements in biophysical experimental techniques have pushed the limits in terms of the types of phenomena that can be characterized, the amount of data that can be produced and the resolution at which we can visualize them. Single particle techniques such as Electron Microscopy (EM) and X-ray free electron laser (XFEL) scattering require a large number of 2D images collected to resolve three-dimensional (3D) structures. In this study, we propose a quick strategy to retrieve potential 3D shapes, as low-resolution models, from a few 2D experimental images by searching a library of 2D projection images generated from existing 3D structures., Results: We developed the protocol to assemble a non-redundant set of 3D shapes for generating the 2D image library, and to retrieve potential match 3D shapes for query images, using EM data as a test. In our strategy, we disregard differences in volume size, giving previously unknown structures and conformations a greater number of 3D biological shapes as possible matches. We tested the strategy using images from three EM models as query images for searches against a library of 22750 2D projection images generated from 250 random EM models. We found that our ability to identify 3D shapes that match the query images depends on how complex the outline of the 2D shapes are and whether they are represented in the search image library., Conclusions: Through our computational method, we are able to quickly retrieve a 3D shape from a few 2D projection images. Our approach has the potential for exploring other types of 2D single particle structural data such as from XFEL scattering experiments, for providing a tool to interpret low-resolution data that may be insufficient for 3D reconstruction, and for estimating the mixing of states or conformations that could exist in such experimental data.

Three-dimensional (3D) structures of biomolecules provide insight into their functions. Using X-ray free-electron laser (XFEL) scattering experiments, it was possible to observe biomolecules that are difficult to crystallize, under conditions that are similar to their natural environment. However, resolving 3D structure from XFEL data is not without its challenges. For example, strong beam intensity is required to obtain sufficient diffraction signal and the beam incidence angles to the molecule need to be estimated for diffraction patterns with significant noise. Therefore, it is important to quantitatively assess how the experimental conditions such as the amount of data and their quality affect the expected resolution of the resulting 3D models. In this study, as an example, the restoration of 3D structure of ribosome from two-dimensional diffraction patterns created by simulation is shown. Tests are performed using the diffraction patterns simulated for different beam intensities and using different numbers of these patterns. Guidelines for selecting parameters for slice-matching 3D reconstruction procedures are established. Also, the minimum requirements for XFEL experimental conditions to obtain diffraction patterns for reconstructing molecular structures to a high-resolution of a few nanometers are discussed.

Protein kinase, casein kinase II (CK2), is ubiquitously expressed and highly conserved protein kinase that shows constitutive activity. It phosphorylates a diverse set of proteins and plays crucial role in several cellular processes. The catalytic subunit of this enzyme (CK2α) shows remarkable flexibility as evidenced in numerous crystal structures determined till now. Here, using analysis of multiple crystal structures and long timescale molecular dynamics simulations, we explore the conformational flexibility of CK2α. The enzyme shows considerably higher flexibility in the solution as compared to that observed in crystal structure ensemble. Multiple conformations of hinge region, located near the active site, were observed during the dynamics. We further observed that among these multiple conformations, the most populated conformational state was inadequately represented in the crystal structure ensemble. The catalytic spine, was found to be less dismantled in this state as compared to the "open" hinge/αD state crystal structures. The comparison of dynamics in unbound (Apo) state and inhibitor (CX4945) bound state exhibits inhibitor induced suppression in the overall dynamics of the enzyme. This is especially true for functionally important glycine-rich loop above the active site. Together, this work gives novel insights into the dynamics of CK2α in solution and relates it to the function. This work also explains the effect of inhibitor on the dynamics of CK2α and paves way for development of better inhibitors., (© 2017 Wiley Periodicals, Inc.)

Hybrid approaches for the modeling of macromolecular complexes that combine computational molecular mechanics simulations with experimental data are discussed. Experimental data for biological molecular structures are often low-resolution, and thus, do not contain enough information to determine the atomic positions of molecules. This is especially true when the dynamics of large macromolecules are the focus of the study. However, computational modeling can complement missing information. Significant increase in computational power, as well as the development of new modeling algorithms allow us to model structures of biological macromolecules reliably, using experimental data as references. We review the basics of molecular mechanics approaches, such as atomic model force field, and coarse-grained models, molecular dynamics simulation and normal mode analysis and describe how they could be used for flexible fitting hybrid modeling with experimental data, especially from cryo-EM and SAXS.

The three-dimensional (3D) structural analysis of single particles using an X-ray free-electron laser (XFEL) is a new structural biology technique that enables observations of molecules that are difficult to crystallize, such as flexible biomolecular complexes and living tissue in the state close to physiological conditions. In order to restore the 3D structure from the diffraction patterns obtained by the XFEL, computational algorithms are necessary as the orientation of the incident beam with respect to the sample needs to be estimated. A program package for XFEL single-particle analysis based on the Xmipp software package, that is commonly used for image processing in 3D cryo-electron microscopy, has been developed. The reconstruction program has been tested using diffraction patterns of an aerosol nanoparticle obtained by tomographic coherent X-ray diffraction microscopy.

Flexible fitting is a computational algorithm to derive a new conformational model that conforms to low-resolution experimental data by transforming a known structure. A common application is against data from cryo-electron microscopy to obtain conformational models in new functional states. The conventional flexible fitting algorithms cannot derive correct structures in some cases due to the complexity of conformational transitions. In this study, we show the importance of conformational ensemble in the refinement process by performing multiple fittings trials using a variety of different force constants. Application to simulated maps of Ca 2+ ATPase and diphtheria toxin as well as experimental data of release factor 2 revealed that for these systems, multiple conformations with similar agreement with the density map exist and a large number of fitting trials are necessary to generate good models. Clustering analysis can be an effective approach to avoid over-fitting models. In addition, we show that an automatic adjustment of the biasing force constants during the fitting process, implemented as replica-exchange scheme, can improve the success rate. © 2017 Wiley Periodicals, Inc., (© 2017 Wiley Periodicals, Inc.)