This study synthesized and evaluated a series of benzotriazole derivatives denoted 3(a-j) and 6(a-j) for their anti-HIV-1 RT activities compared to the standard drug efavirenz. Notably, compound 3 h, followed closely by...This study synthesized and evaluated a series of benzotriazole derivatives denoted 3(a-j) and 6(a-j) for their anti-HIV-1 RT activities compared to the standard drug efavirenz. Notably, compound 3 h, followed closely by 6 h, exhibited significant anti-HIV-1 RT efficacy relative to the standard drug. In vivo oral toxicity studies were conducted for the most active compound 3 h, confirming its nontoxic nature to ascertain the safety profile. By employing molecular docking techniques, we explored the potential interactions between the synthesized compounds (ligands) and a target biomolecule (protein)(PDB ID 1RT2) at the molecular level. We undertook the molecular dynamics study of 3 h, the most active compound, within the active binding pocket of the cocrystallized structure of HIV-1 RT (PDB ID 1RT2). We aimed to learn more about how biomolecular systems behave, interact, and change at the atomic or molecular level over time. Finally, the DFT-derived HOMO and LUMO orbitals, as well as analysis of the molecular electrostatic potential map, aid in discerning the reactivity characteristics of our molecule.
Cholesterol is a negative regulator of a variety of ion channels. We have previously shown that cholesterol suppresses Kir2.2 channels via residue-residue uncoupling on the inter-subunit interfaces within the close state...Cholesterol is a negative regulator of a variety of ion channels. We have previously shown that cholesterol suppresses Kir2.2 channels via residue-residue uncoupling on the inter-subunit interfaces within the close state of the channels (3JYC). In this study, we extend this analysis to the other known structure of Kir2.2 that is closer to the open state of Kir2.2 channels (3SPI) and provide additional analysis of the residue distances between the uncoupled residues and cholesterol binding domains in the two conformation states of the channels. We found that the general phenomenon of cholesterol binding leading to uncoupling between specific residues is conserved in both channel states but the specific pattern of the uncoupling residues is distinct between the two states and implies different mechanisms. Specifically, we found that cholesterol binding in the 3SPI state results in an uncoupling of residues in three distinct regions; the transmembrane domain, membrane-cytosolic interface, and the cytosolic domain, with the first two regions forming an envelope around PI(4,5)P2 and cholesterol binding sites and the distal region overlapping with the subunit-subunit interface characterized in our previous study of the disengaged state. We also found that this uncoupling is dependent upon the number of cholesterol molecules bound to the channel. We further generated a mutant channel Kir2.2 with a single point mutation in a residue proximal to the PI(4,5)P2 binding site, which is predicted to be uncoupled from other residues in its vicinity upon cholesterol binding and found that this mutation abrogates the sensitivity of Kir2.2 to cholesterol changes in the membrane. These findings suggest that cholesterol binding to this conformation state of Kir2.2 channels may destabilize the PI(4,5)P2 interactions with the channels while in the disengaged state the destabilization occurs where the subunits interact. These findings give insight into the structural mechanistic basis for the functional effects of cholesterol binding to the Kir2.2 channel.
Aguirre-Sampieri S, Casañal A, Emsley P
… +1 more, Garza-Ramos G
J Struct Biol
· 2024 Jun · PMID 38615726
·
Full text
Many enzymes can self-assemble into higher-order structures with helical symmetry. A particularly noteworthy example is that of nitrilases, enzymes in which oligomerization of dimers into spiral homo-oligomers is a requi...Many enzymes can self-assemble into higher-order structures with helical symmetry. A particularly noteworthy example is that of nitrilases, enzymes in which oligomerization of dimers into spiral homo-oligomers is a requirement for their enzymatic function. Nitrilases are widespread in nature where they catalyze the hydrolysis of nitriles into the corresponding carboxylic acid and ammonia. Here, we present the Cryo-EM structure, at 3 Å resolution, of a C-terminal truncate nitrilase from Rhodococcus sp. V51B that assembles in helical filaments. The model comprises a complete turn of the helical arrangement with a substrate-intermediate bound to the catalytic cysteine. The structure was solved having added the substrate to the protein. The length and stability of filaments was made more substantial in the presence of the aromatic substrate, benzonitrile, but not for aliphatic nitriles or dinitriles. The overall structure maintains the topology of the nitrilase family, and the filament is formed by the association of dimers in a chain-like mechanism that stabilizes the spiral. The active site is completely buried inside each monomer, while the substrate binding pocket was observed within the oligomerization interfaces. The present structure is in a closed configuration, judging by the position of the lid, suggesting that the intermediate is one of the covalent adducts. The proximity of the active site to the dimerization and oligomerization interfaces, allows the dimer to sense structural changes once the benzonitrile was bound, and translated to the rest of the filament, stabilizing the helical structure.
Fu Z, Crooks EJ, Irizarry BA
… +4 more, Zhu X, Chowdhury S, Van Nostrand WE, Smith SO
J Struct Biol
· 2024 Jun · PMID 38615725
·
Full text
Cerebral amyloid angiopathy (CAA) is associated with the accumulation of fibrillar Aβ peptides upon and within the cerebral vasculature, which leads to loss of vascular integrity and contributes to disease progression in...Cerebral amyloid angiopathy (CAA) is associated with the accumulation of fibrillar Aβ peptides upon and within the cerebral vasculature, which leads to loss of vascular integrity and contributes to disease progression in Alzheimer's disease (AD). We investigate the structure of human-derived Aβ40 fibrils obtained from patients diagnosed with sporadic or familial Dutch-type (E22Q) CAA. Using cryo-EM, two primary structures are identified containing elements that have not been observed in in vitro Aβ40 fibril structures. One population has an ordered N-terminal fold comprised of two β-strands stabilized by electrostatic interactions involving D1, E22, D23 and K28. This charged cluster is disrupted in the second population, which exhibits a disordered N-terminus and is favored in fibrils derived from the familial Dutch-type CAA patient. These results illustrate differences between human-derived CAA and AD fibrils, and how familial CAA mutations can guide fibril formation.
Ethionamide (ETO) is a prodrug that is primarily used as a second-line agent in the treatment of tuberculosis. Among the bacterial ETO activators, the monooxygenase MymA has been recently identified, and its expression i...Ethionamide (ETO) is a prodrug that is primarily used as a second-line agent in the treatment of tuberculosis. Among the bacterial ETO activators, the monooxygenase MymA has been recently identified, and its expression is regulated by the mycobacterial regulator VirS. The discovery of VirS ligands that can enhance mymA expression and thereby increase the antimycobacterial efficacy of ETO, has led to the development of a novel therapeutic strategy against tuberculosis. This strategy involves the selection of preclinical candidates, including SMARt751. We report the first crystal structure of the AraC-like regulator VirS, in complex with SMARt751, refined at 1.69 Å resolution. Crystals were obtained via an in situ proteolysis method in the requisite presence of SMARt751. The elucidated structure corresponds to the ligand-binding domain of VirS, adopting an α/β fold with structural similarities to H-NOX domains. Within the VirS structure, SMARt751 is situated in a completely enclosed hydrophobic cavity, where it forms hydrogen bonds with Asn11 and Asn149 as well as van der Waals contacts with various hydrophobic amino acids. Comprehensive structural comparisons within the AraC family of transcriptional regulators are conducted and analyzed to figure out the effects of the SMARt751 binding on the regulatory activity of VirS.
Fusion proteins (FPs) are frequently utilized as a biotechnological tool in the determination of macromolecular structures using X-ray methods. Here, we explore the use of different protein tags in various FP, to obtain...Fusion proteins (FPs) are frequently utilized as a biotechnological tool in the determination of macromolecular structures using X-ray methods. Here, we explore the use of different protein tags in various FP, to obtain initial phases by using them in a partial molecular replacement (MR) and constructing the remaining FP structure with ARP/wARP. Usually, the tag is removed prior to crystallization, however leaving the tag on may facilitate crystal formation, and structural determination by expanding phases from known to unknown segments of the complex. In this study, the Protein Data Bank was mined for an up-to-date list of FPs with the most used protein tags, Maltose Binding Protein (MBP), Green Fluorescent Protein (GFP), Thioredoxin (TRX), Glutathione transferase (GST) and the Small Ubiquitin-like Modifier Protein (SUMO). Partial MR using the protein tag, followed by automatic model building, was tested on a subset of 116 FP. The efficiency of this method was analyzed and factors that influence the coordinate construction of a substantial portions of the fused protein were identified. Using MBP, GFP, and SUMO as phase generators it was possible to build at least 75 % of the protein of interest in 36 of the 116 cases tested. Our results reveal that tag selection has a significant impact; tags with greater structural stability, such as GFP, increase the success rate. Further statistical analysis identifies that resolution, Wilson B factor, solvent percentage, completeness, multiplicity, protein tag percentage in the FP (considering amino acids), and the linker length play pivotal roles using our approach. In cases where a structural homologous is absent, this method merits inclusion in the toolkit of protein crystallographers.
Melanin granules (melanosomes) in Asian and Caucasian black hairs were investigated by focused ion beam scanning electron microscopy (FIB-SEM). This technique facilitates a direct evaluation of the three-dimensional dist...Melanin granules (melanosomes) in Asian and Caucasian black hairs were investigated by focused ion beam scanning electron microscopy (FIB-SEM). This technique facilitates a direct evaluation of the three-dimensional distribution and morphology of melanin granules without requiring their isolation from hair. Three-dimensional reconstructed images of melanin granule distribution in hair samples were obtained using serial SEM images observed by FIB-SEM. Melanin granules in black hair tended to be three-dimensionally dense in the outer periphery of the cortex. The morphometric parameters of melanin granules were calculated using the reconstructed three-dimensional images. The results confirmed that melanin granules in Caucasian black hair were much smaller those in Asian black hair. Moreover, it was indicated that the relative frequency distribution of the volume of melanin granules was significantly different between Asians and Caucasians.
Staphylococcus aureus, an ESKAPE pathogen, is a major clinical concern due to its pathogenicity and manifold antimicrobial resistance mechanisms. The commonly used β-lactam antibiotics target bacterial penicillin-binding...Staphylococcus aureus, an ESKAPE pathogen, is a major clinical concern due to its pathogenicity and manifold antimicrobial resistance mechanisms. The commonly used β-lactam antibiotics target bacterial penicillin-binding proteins (PBPs) and inhibit crosslinking of peptidoglycan strands that comprise the bacterial cell wall mesh, initiating a cascade of effects leading to bacterial cell death. S. aureus PBP1 is involved in synthesis of the bacterial cell wall during division and its presence is essential for survival of both antibiotic susceptible and resistant S. aureus strains. Here, we present X-ray crystallographic data for S. aureus PBP1 in its apo form as well as acyl-enzyme structures with distinct classes of β-lactam antibiotics representing the penicillins, carbapenems, and cephalosporins, respectively: oxacillin, ertapenem and cephalexin. Our structural data suggest that the PBP1 active site is readily accessible for substrate, with little conformational change in key structural elements required for its covalent acylation of β-lactam inhibitors. Stopped-flow kinetic analysis and gel-based competition assays support the structural observations, with even the weakest performing β-lactams still having comparatively high acylation rates and affinities for PBP1. Our structural and kinetic analysis sheds insight into the ligand-PBP interactions that drive antibiotic efficacy against these historically useful antimicrobial targets and expands on current knowledge for future drug design and treatment of S. aureus infections.
The global spread of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) since 2019 has led to a continuous evolution of viral variants, with the latest concern being the Omicron (B.1.1.529) variant. In this stu...The global spread of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) since 2019 has led to a continuous evolution of viral variants, with the latest concern being the Omicron (B.1.1.529) variant. In this study, classical molecular dynamics simulations were conducted to elucidate the biophysical aspects of the Omicron spike protein's receptor-binding domain (RBD) in its interaction with human angiotensin-converting enzyme 2 (hACE2) and a neutralizing antibody, comparing it to the wildtype (WT). To model the Omicron variant, 15 in silico mutations were introduced in the RBD region of WT (retrieved from PDB). The simulations of WT spike-hACE2 and Omicron spike-hACE2 complexes revealed comparable binding stability and dynamics. Notably, the Q493R mutation in the Omicron spike increased interactions with hACE2, particularly with ASP38 and ASP355. Additionally, mutations such as N417K, T478K, and Y505H contributed to enhanced structural stability in the Omicron variant. Conversely, when comparing WT with Omicron in complex with a neutralizing antibody, simulation results demonstrated poorer binding dynamics and stability for the Omicron variant. The E484K mutation significantly decreased binding interactions, resulting in an overall decrease in binding energy (∼-57 kcal/mol) compared to WT (∼-84 kcal/mol). This study provides valuable molecular insights into the heightened infectivity of the Omicron variant, shedding light on the specific mutations influencing its interactions with hACE2 and neutralizing antibodies.
The goal of cryo-EM experiments in the biological sciences is to determine the atomic structure of a molecule and deduce insights into its functions and mechanisms. Despite improvements in instrumentation for data collec...The goal of cryo-EM experiments in the biological sciences is to determine the atomic structure of a molecule and deduce insights into its functions and mechanisms. Despite improvements in instrumentation for data collection and new software algorithms, in most cases, individual atoms are not resolved. Model building of proteins, nucleic acids, or molecules in general, is feasible from the experimentally determined density maps at resolutions up to the range of 3-4 Angstroms. For lower-resolution maps or parts of maps, fitting smaller structures obtained by modelling or experimental techniques with higher resolution is a way to resolve the issue. In practice, we have an atomic structure, generate its density map at a given resolution, and translate/rotate the map within a region of interest in the experimental map, computing a measure-of-fit score with the corresponding areas of the experimental map. This procedure is computationally intensive since we work in 6D space. An optimal ordered list of rotations will reduce the angular error and help to find the best-fitting positions faster for a coarse global search or a local refinement. It can be used for adaptive approaches to stop fitting algorithms earlier once the desired accuracy has been achieved. We demonstrate how the performance of some fitting algorithms can be improved by grouping sets of rotations. We present an approach to generate more efficient 3D angular sampling, and provide the computer code to generate lists of optimal orientations for single and grouped rotations and the lists themselves.
In humans, the growth pattern of the acellular extrinsic fibre cementum (AEFC) has been useful to estimate the age-at-death. However, the structural organization behind such a pattern remains poorly understood. In this s...In humans, the growth pattern of the acellular extrinsic fibre cementum (AEFC) has been useful to estimate the age-at-death. However, the structural organization behind such a pattern remains poorly understood. In this study tooth cementum from seven individuals from a Mexican modern skeletal series were analyzed with the aim of unveiling the AEFC collagenous and mineral structure using multimodal imaging approaches. The organization of collagen fibres was first determined using: light microscopy, transmission electron microscopy (TEM), electron tomography, and plasma FIB scanning electron microscopy (PFIB-SEM) tomography. The mineral properties were then investigated using: synchrotron small-angle X-ray scattering (SAXS) for T-parameter (correlation length between mineral particles); synchrotron X-ray diffraction (XRD) for L-parameter (mineral crystalline domain size estimation), alignment parameter (crystals preferred orientation) and lattice parameters a and c; as well as synchrotron X-ray fluorescence for spatial distribution of calcium, phosphorus and zinc. Results show that Sharpey's fibres branched out fibres that cover and uncover other collagen bundles forming aligned arched structures that are joined by these same fibres but in a parallel fashion. The parallel fibres are not set as a continuum on the same plane and when they are superimposed project the AEFC incremental lines due to the collagen birefringence. The orientation of the apatite crystallites is subject to the arrangement of the collagen fibres, and the obtained parameter values along with the elemental distribution maps, revealed this mineral tissue as relatively homogeneous. Therefore, no intrinsic characteristics of the mineral phase could be associated with the alternating AEFC incremental pattern.
While protein activity is traditionally studied with a major focus on the active site, the activity of enzymes has been hypothesized to be linked to the flexibility of adjacent regions, warranting more exploration into h...While protein activity is traditionally studied with a major focus on the active site, the activity of enzymes has been hypothesized to be linked to the flexibility of adjacent regions, warranting more exploration into how the dynamics in these regions affects catalytic turnover. One such enzyme is Xylanase A (XylA), which cleaves hemicellulose xylan polymers by hydrolysis at internal β-1,4-xylosidic linkages. It contains a "thumb" region whose flexibility has been suggested to affect the activity. The double mutation D11F/R122D was previously found to affect activity and potentially bias the thumb region to a more open conformation. We find that the D11F/R122D double mutation shows substrate-dependent effects, increasing activity on the non-native substrate ONPX2 but decreasing activity on its native xylan substrate. To characterize how the double mutant causes these kinetics changes, nuclear magnetic resonance (NMR) and molecular dynamics (MD) simulations were used to probe structural and flexibility changes. NMR chemical shift perturbations revealed structural changes in the double mutant relative to the wild-type, specifically in the thumb and fingers regions. Increased slow-timescale dynamics in the fingers region was observed as intermediate-exchange line broadening. Lipari-Szabo order parameters show negligible changes in flexibility in the thumb region in the presence of the double mutation. To help understand if there is increased energetic accessibility to the open state upon mutation, alchemical free energy simulations were employed that indicated thumb opening is more favorable in the double mutant. These studies aid in further characterizing how flexibility in adjacent regions affects the function of XylA.
Cryo-electron microscopy has become a powerful tool to determine three-dimensional (3D) structures of rigid biological macromolecules from noisy micrographs with single-particle reconstruction. Recently, deep neural netw...Cryo-electron microscopy has become a powerful tool to determine three-dimensional (3D) structures of rigid biological macromolecules from noisy micrographs with single-particle reconstruction. Recently, deep neural networks, e.g., CryoDRGN, have demonstrated conformational and compositional heterogeneity of complexes. However, the lack of ground-truth conformations poses a challenge to assess the performance of heterogeneity analysis methods. In this work, variational autoencoders (VAE) with three types of deep generative priors were learned for latent variable inference and heterogeneous 3D reconstruction via Bayesian inference. More specifically, VAEs with "Variational Mixture of Posteriors" priors (VampPrior-SPR), non-parametric exemplar-based priors (ExemplarPrior-SPR) and priors from latent score-based generative models (LSGM-SPR) were quantitatively compared with CryoDRGN. We built four simulated datasets composed of hypothetical continuous conformation or discrete states of the hERG K + channel. Empirical and quantitative comparisons of inferred latent representations were performed with affine-transformation-based metrics. These models with more informative priors gave better regularized, interpretable factorized latent representations with better conserved pairwise distances, less deformed latent distributions and lower within-cluster variances. They were also tested on experimental datasets to resolve compositional and conformational heterogeneity (50S ribosome assembly, cowpea chlorotic mottle virus, and pre-catalytic spliceosome) with comparable high resolution. Codes and data are available: https://github.com/benjamin3344/DGP-SPR.
Calcium carbonate is present in many biominerals, including in the exoskeletons of crustaceans and shells of mollusks. High Mg-containing calcium carbonate was synthesized by high temperatures, high pressures or high mol...Calcium carbonate is present in many biominerals, including in the exoskeletons of crustaceans and shells of mollusks. High Mg-containing calcium carbonate was synthesized by high temperatures, high pressures or high molecular organic matter. For example, biogenic high Mg-containing calcite is synthesized under strictly controlled Mg concentration at ambient temperature and pressure. The spines of sea urchins consist of calcite, which contain a high percentage of magnesium. In this study, we investigated the factors that increase the magnesium content in calcite from the spines of the sea urchin, Heliocidaris crassispina. X-ray diffraction and inductively coupled plasma mass spectrometry analyses showed that sea urchin spines contain about 4.8% Mg. The organic matrix extracted from the H. crassispina spines induced the crystallization of amorphous phase and synthesis of magnesium-containing calcite, while amorphous was synthesized without SUE (sea urchin extract). In addition, aragonite was synthesized by SUE treated with protease-K. HC tropomyosin was specifically incorporated into Mg precipitates. Recombinant HC-tropomyosin induced calcite contained 0.1-2.5% Mg synthesis. Western blotting of sea urchin spine extracts confirmed that HC tropomyosin was present in the purple sea urchin spines at a protein weight ratio of 1.5%. These results show that HC tropomyosin is one factor that increases the magnesium concentration in the calcite of H. crassispina spines.
J Struct Biol
· 2024 Jun · PMID 38431179
·
Full text
Efficient and high-accuracy filtering of cryo-electron microscopy (cryo-EM) micrographs is an emerging challenge with the growing speed of data collection and sizes of datasets. Convolutional neural networks (CNNs) are m...Efficient and high-accuracy filtering of cryo-electron microscopy (cryo-EM) micrographs is an emerging challenge with the growing speed of data collection and sizes of datasets. Convolutional neural networks (CNNs) are machine learning models that have been proven successful in many computer vision tasks, and have been previously applied to cryo-EM micrograph filtering. In this work, we demonstrate that two strategies, fine-tuning models from pretrained weights and including the power spectrum of micrographs as input, can greatly improve the attainable prediction accuracy of CNN models. The resulting software package, Miffi, is open-source and freely available for public use (https://github.com/ando-lab/miffi).
Filipek J, Chalaskiewicz K, Kosmider A
… +7 more, Nielipinski M, Michalak A, Bednarkiewicz M, Goslawski-Zeligowski M, Prucnal F, Sekula B, Pietrzyk-Brzezinska AJ
TetR family regulators (TFRs) represent a large group of one-component bacterial signal transduction systems which recognize environmental signals, like the presence of antibiotics or other bactericidal compounds, and tr...TetR family regulators (TFRs) represent a large group of one-component bacterial signal transduction systems which recognize environmental signals, like the presence of antibiotics or other bactericidal compounds, and trigger the cell response by regulating the expression of genes that secure bacterial survival in harsh environmental conditions. TFRs act as homodimers, each protomer is composed of a conserved DNA-binding N-terminal domain (NTD) and a variable ligand-binding C-terminal domain (CTD). Currently, there are about 500 structures of TFRs available in the Protein Data Bank and one-fourth of them represent the structures of TFR-ligand complexes. In this review, we summarized information on the ligands interacting with TFRs and based on structural data, we compared the CTDs of the TFR family members, as well as their ligand-binding cavities. Additionally, we divided the whole TFR family, including more than half of a million sequences, into subfamilies according to calculated multiple sequence alignment and phylogenetic tree. We also highlighted structural elements characteristic of some of the subfamilies. The presented comprehensive overview of the TFR CTDs provides good bases and future directions for further studies on TFRs that are not only important targets for battling multidrug resistance but also good candidates for many biotechnological approaches, like TFR-based biosensors.
CRISPR-Cas system is an RNA-guided adaptive immune system widespread in bacteria and archaea. Among them, type III CRISPR-Cas systems are the most ancient throughout the CRISPR-Cas family, proving anti-phage defense thro...CRISPR-Cas system is an RNA-guided adaptive immune system widespread in bacteria and archaea. Among them, type III CRISPR-Cas systems are the most ancient throughout the CRISPR-Cas family, proving anti-phage defense through a crRNA-guided RNA targeting manner and possessing multiple enzymatic activities. Type III CRISPR-Cas systems comprise four typical members (type III-A to III-D) and two atypical members (type III-E and type III-F), providing immune defense through distinct mechanisms. Here, we delve into structural studies conducted on three well-characterized members: the type III-A, III-B, and III-E systems, provide an overview of the structural insights into the crRNA-guided target RNA cleavage, self/non-self discrimination, and the target RNA-dependent regulation of enzymatic subunits in the effector complex.
Cellular cryo-electron tomography (cryo-ET) has emerged as a key method to unravel the spatial and structural complexity of cells in their near-native state at unprecedented molecular resolution. To enable quantitative a...Cellular cryo-electron tomography (cryo-ET) has emerged as a key method to unravel the spatial and structural complexity of cells in their near-native state at unprecedented molecular resolution. To enable quantitative analysis of the complex shapes and morphologies of lipid membranes, the noisy three-dimensional (3D) volumes must be segmented. Despite recent advances, this task often requires considerable user intervention to curate the resulting segmentations. Here, we present ColabSeg, a Python-based tool for processing, visualizing, editing, and fitting membrane segmentations from cryo-ET data for downstream analysis. ColabSeg makes many well-established algorithms for point-cloud processing easily available to the broad community of structural biologists for applications in cryo-ET through its graphical user interface (GUI). We demonstrate the usefulness of the tool with a range of use cases and biological examples. Finally, for a large Mycoplasma pneumoniae dataset of 50 tomograms, we show how ColabSeg enables high-throughput membrane segmentation, which can be used as valuable training data for fully automated convolutional neural network (CNN)-based segmentation.
Haynes RM, Myers J, López CS
… +3 more, Evans J, Davulcu O, Yoshioka C
J Struct Biol
· 2025 Mar · PMID 38364988
·
Full text
In recent years, cryo-electron microscopy (cryo-EM) has become a practical and effective method of determining structures at previously unattainable resolutions due to advances in detection, automation, and data processi...In recent years, cryo-electron microscopy (cryo-EM) has become a practical and effective method of determining structures at previously unattainable resolutions due to advances in detection, automation, and data processing. However, sample preparation remains a major bottleneck in the cryo-EM workflow. Even after the arduous process of biochemical sample optimization, it often takes several iterations of grid vitrification and screening to determine the optimal grid freezing parameters that yield suitable ice thickness and particle distribution for data collection. Since a high-quality sample is imperative for high-resolution structure determination, grid optimization is a vital step. For researchers who rely on cryo-EM facilities for grid screening, each iteration of this optimization process may delay research progress by a matter of months. Therefore, a more strategic and efficient approach should be taken to ensure that the grid optimization process can be completed in as few iterations as possible. Here, we present an implementation of Design of Experiments (DOE) to expedite and strategize the grid optimization process. A Fractional Factorial Design (FFD) guides the determination of a limited set of experimental conditions which can model the full parameter space of interest. Grids are frozen with these conditions and screened for particle distribution and ice thickness. Quantitative scores are assigned to each of these grid characteristics based on a qualitative rubric. Input conditions and response scores are used to generate a least-squares regression model of the parameter space in JMP, which is used to determine the conditions which should, in theory, yield optimal grids. Upon testing this approach on apoferritin and L-glutamate dehydrogenase on both the Vitrobot Mark IV and the Leica GP2 plunge freezers, the resulting grid conditions reliably yielded grids with high-quality ice and particle distribution that were suitable for collecting large overnight datasets on a Krios. We conclude that a DOE-based approach is a cost-effective and time-saving tool for cryo-EM grid preparation.