Deep learning-based mixed-dimensional Gaussian mixture model for characterizing variability in cryo-EM.

Integrating Molecular Models Into CryoEM Heterogeneity Analysis Using Scalable High-resolution Deep Gaussian Mixture Models

The motor enzymes that belong to the family of RNA helicases catalyze the strand separation of duplex RNA via ATP hydrolysis. Among these enzymes, Escherichia coli DbpA is a unique RNA helicase because it possesses ATPase-specific activity toward the peptidyl transferase center in 23 S ribosomal RNA. For this reason, it has been the subject of numerous biochemical and structure-function studies. The ATP-stimulated unwinding activity of DbpA toward specific and nonspecific RNA duplexes has been demonstrated. However, the underlying molecular and structural basis, which facilitates its helicase activities, is presently not known. We combined time-dependent limited proteolysis digestion, fluorescence spectroscopy, and three-dimensional structural homology modeling techniques to study the structural conformations of DbpA with respect to its binding to stoichiometric ratios of RNA and cofactors. We show that the conformational state of DbpA is markedly different in the ADP-bound state than in any other state (ATP- or RNA-bound). These results, together with structural homology studies, suggest that a hinge region located in the core domain of DbpA mediates such conformational changes.

In P(VDF-TrFE) random copolymers with compositions from 60/40 to 80/20 mol%, the ferroelectric phase transition exhibits a strong first order character and is accompanied by a large expansion of the crystal lattice. In the ferroelectric structure, the polymer chains have an all-trans conformation and are close-packed parallel to one another in a pseudo-hexagonal lattice. According to the m2m symmetry, the spontaneous polarization Py is perpendicular to the chain axis 2. In the paraelectric phase, the structure is disordered : the chains undergo continuous conformational changes, they have a statistical Cm symmetry and the expanded lattice is hexagonal (6/m point group)

Resolving structural heterogeneity in situ through cryogenic electron tomography

PCNA is a master coordinator of many DNA-metabolic events. During DNA replication, the maturation of Okazaki fragments involves at least four DNA enzymes, all of which contain PCNA-interacting motifs. However, the temporal relationships and functional modulations between these PCNA-binding proteins are unclear. Here, we developed a strategy to purify endogenous PCNA-containing complexes from native chromatin, and characterized their structures using cryo-EM. Two structurally resolved classes (PCNA-FEN1 and PCNA-FEN1-RNaseH2 complexes) have captured a series of 3D snapshots for the primer-removal steps of Okazaki fragment maturation. These structures show that product release from FEN1 is a rate-liming step. Furthermore, both FEN1 and RNaseH2 undergo continuous conformational changes on PCNA that result in constant fluctuations in the bending angle of substrate DNA at the nick site, implying that these enzymes could regulate each other through conformational modulation of the bound DNA. The structures of the PCNA-FEN1-RNaseH2 complex confirm the toolbelt function of PCNA and suggests a potential unrecognized role of RNaseH2, as a dsDNA binding protein, in promoting the 5′-flap cleaving activity of FEN1.

RNAs play critical roles in most biological processes. Although the three-dimensional (3D) structures of RNAs primarily determine their functions, it remains challenging to experimentally determine these 3D structures due to their conformational heterogeneity and intrinsic dynamics. Cryogenic electron microscopy (cryo-EM) has recently played an emerging role in resolving dynamic conformational changes and understanding structure-function relationships of RNAs including ribozymes, riboswitches and bacterial and viral noncoding RNAs. A variety of methods and pipelines have been developed to facilitate cryo-EM structure determination of challenging RNA targets with small molecular weights at subnanometer to near-atomic resolutions. While a wide range of conditions have been used to prepare RNAs for cryo-EM analysis, correlations between the variables in these conditions and cryo-EM visualizations and reconstructions remain underexplored, which continue to hinder optimizations of RNA samples for high-resolution cryo-EM structure determination. Here we present a protocol that describes rigorous screenings and iterative optimizations of RNA preparation conditions that facilitate cryo-EM structure determination, supplemented by cryo-EM data processing pipelines that resolve RNA dynamics and conformational changes and RNA modeling algorithms that generate atomic coordinates based on moderate- to high-resolution cryo-EM density maps. The current protocol is designed for users with basic skills and experience in RNA biochemistry, cryo-EM and RNA modeling. The expected time to carry out this protocol may range from 3 days to more than 3 weeks, depending on the many variables described in the protocol. For particularly challenging RNA targets, this protocol could also serve as a starting point for further optimizations.

Domain-Based Protein Docking with Extremely Large Conformational Changes

In the last decade, various coarse-grained elastic network models have been developed to study the large-scale motions of proteins and protein complexes where computer simulations using detailed all-atom models are not feasible. Among these models, the Gaussian Network Model (GNM) and Anisotropic Network Model (ANM) have been widely used. Both models have strengths and limitations. GNM can predict the relative magnitudes of the fluctuations well, but due to its isotropic assumption, it can not be applied to predict the directions of the fluctuations. In contrast, ANM adds the ability to do the latter, but it loses a significant amount of precision in the prediction of the magnitudes. In this article, we develop a generalized spring tensor model (STeM) that is able to predict well both the magnitudes and the directions of the fluctuations. STeM also outperforms ANM in explaining protein conformation changes. All of these are accomplished without sacrificing the essential features that have made ANM and GNM attractive.

Rotationally averaged collision cross section (CCS) values for a series of proteins and protein complexes ranging in size from 8.6 to 810 kDa are reported. The CCSs were obtained using a native electrospray ionization drift tube ion mobility-Orbitrap mass spectrometer specifically designed to enhance sensitivity while having high-resolution ion mobility and mass capabilities. Periodic focusing (PF)-drift tube (DT)-ion mobility (IM) provides first-principles determination of the CCS of large biomolecules that can then be used as CCS calibrants. The experimental, first-principles CCS values are compared to previously reported experimentally determined and computationally calculated CCS using projected superposition approximation (PSA), the Ion Mobility Projection Approximation Calculation Tool (IMPACT), and Collidoscope. Experimental CCS values are generally in agreement with previously reported CCSs, with values falling within ∼5.5%. In addition, an ion mobility resolution (CCS centroid divided by CCS fwhm) of ∼60 is obtained for pyruvate kinase (MW ∼ 233 kDa); however, ion mobility resolution for bovine serum albumin (MW ∼ 68 kDa) is less than ∼20, which arises from sample impurities and underscores the importance of sample quality. The high resolution afforded by the ion mobility-Orbitrap mass analyzer provides new opportunities to understand the intricate details of protein complexes such as the impact of post-translational modifications (PTMs), stoichiometry, and conformational changes induced by ligand binding.

Allosteric changes modulate the enzymatic activity, leading to activation or inhibition of the molecular target. Understanding the induced fit accommodation mechanism of a ligand in its lowest-free energy state and the subsequent conformational changes induced in the protein are important questions for drug design. In the present study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were applied to analyze the glycerol-3-phosphate dehydrogenase of Leishmania mexicana (LmGPDH) conformational changes induced by its cofactor and substrate binding. GPDH is a nicotinamide adenine dinucleotide (NAD)-dependent enzyme, which has been reported as an interesting target for drug discovery and development against leishmaniasis. Despite its relevance for glycolysis and pentose phosphate pathways, the structural flexibility and conformational motions of LmGPDH in complex with NADH and dihydroxyacetone phosphate (DHAP) remain unexplored. Here, we analyzed the conformational dynamics of the enzyme-NADH complex (cofactor), and the enzyme-NADH-DHAP complex (adduct), mapped the hydrogen-bond interactions for the complexes and pointed some structural determinants of the enzyme that emerge from these contacts to NADH and DHAP. Finally, we proposed a consistent mechanism for the conformational changes on the first step of the reversible redox conversion of dihydroxyacetone phosphate to glycerol 3-phosphate, indicating key residues and interactions that could be further explored in drug discovery.

Dynamic conformational and structural changes in proteins and protein complexes play a central and ubiquitous role in the regulation of protein function, yet it is very challenging to study these changes, especially for large protein complexes, under physiological conditions. Here, we introduce a novel isobaric crosslinker, Qlinker, for studying conformational and structural changes in proteins and protein complexes using quantitative crosslinking mass spectrometry. Qlinkers are small and simple, amine-reactive molecules with an optimal extended distance of ~10 Å, which use MS2 reporter ions for relative quantification of Qlinker-modified peptides derived from different samples. We synthesized the 2-plex Q2linker and showed that the Q2linker can provide quantitative crosslinking data that pinpoints key conformational and structural changes in biosensors, binary and ternary complexes composed of the general transcription factors TBP, TFIIA, and TFIIB, and RNA polymerase II complexes.

Dynamic conformational and structural changes in proteins and protein complexes play a central and ubiquitous role in the regulation of protein function, yet it is very challenging to study these changes, especially for large protein complexes, under physiological conditions. Here, we introduce a novel isobaric crosslinker, Qlinker, for studying conformational and structural changes in proteins and protein complexes using quantitative crosslinking mass spectrometry. Qlinkers are small and simple, amine-reactive molecules with an optimal extended distance of ~10 Å, which use MS2 reporter ions for relative quantification of Qlinker-modified peptides derived from different samples. We synthesized the 2-plex Q2linker and showed that the Q2linker can provide quantitative crosslinking data that pinpoints key conformational and structural changes in biosensors, binary and ternary complexes composed of the general transcription factors TBP, TFIIA, and TFIIB, and RNA polymerase II complexes.

Dynamic conformational and structural changes in proteins and protein complexes play a central and ubiquitous role in the regulation of protein function, yet it is very challenging to study these changes, especially for large protein complexes, under physiological conditions. Here we introduce a novel isobaric crosslinker, Qlinker, for studying conformational and structural changes in proteins and protein complexes using quantitative crosslinking mass spectrometry (qCLMS). Qlinkers are small and simple, amine-reactive molecules with an optimal extended distance of ∼10 Å which use MS2 reporter ions for relative quantification of Qlinker-modified peptides derived from different samples. We synthesized the 2-plex Q2linker and showed that the Q2linker can provide quantitative crosslinking data that pinpoints key conformational and structural changes in biosensors, binary and ternary complexes composed of the general transcription factors TBP, TFIIA, and TFIIB, and RNA polymerase II (pol II) complexes.

The karyopherin family of nuclear transport receptors is composed of a long array of amphiphilic α-helices and undergoes flexible conformational changes to pass through the hydrophobic crowding barrier of the nuclear pore. Here, we focused on the characteristic enrichment of prolines in the middle of the outer α-helices of importin-β. When these prolines were substituted with alanine, nuclear transport activity was reduced drastically in vivo and in vitro, and caused a severe defect in mitotic progression. These mutations did not alter the overall folding of the helical repeat or affect its interaction with cargo or the regulatory factor Ran. However, in vitro and in silico analyses revealed that the mutant lost structural flexibility and could not undergo rapid conformational changes when transferring from a hydrophilic to hydrophobic environment or vice versa. These findings reveal the essential roles of prolines in ensuring the structural flexibility and functional integrity of karyopherins.

1. Three floods (July 2000, August 2002, September 2003) and a hurricane (October 2001) that occurred in a lowland forest in the southern Maya Mountains of Belize presented an opportunity to evaluate the influence of these disturbances on the structure of a small mammal assemblage. 2. Four terrestrial and four primarily scansorial/arboreal species were trapped July 2000-March 2005 in six grids over 14 irregularly spaced trapping periods. 3. Community dynamics were characterized more by changes in species composition than changes in diversity. The dynamics were driven by species-specific variation in abundance, with changes in composition generally, but not exclusively, due to the occurrence or disappearance of species at low abundance. Despite the disturbances, species richness remained relatively constant. Evenness within the assemblage was consistently low, primarily as a result of dominance by one species, Heteromys desmarestianus. 4. Effects of flooding on community structure were direct but relatively brief (< 1 year), and varied with the duration and intensity of flooding. Effects from the hurricane were indirect but long-lasting and strongly related to severely reduced food resources. 5. This study suggests that long-term dynamics in the structure of many animal communities in the tropics often results from interactions between direct and indirect effects of disturbance. It also suggests that community resistance will depend on variation in disturbance type and regime, but resilience will be determined by the life-history characteristics of each species.