Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, San Francisco, United States
Ashraya Ravikumar
Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, San Francisco, United States
Shivani Sharma
Structural Biology Initiative, CUNY Advanced Science Research Center, New York, United States; Ph.D. Program in Biology, The Graduate Center, City University of New York, New York, United States
Blake Riley
Structural Biology Initiative, CUNY Advanced Science Research Center, New York, United States
Akshay Raju
Structural Biology Initiative, CUNY Advanced Science Research Center, New York, United States
Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, San Francisco, United States
Jessica Flowers
Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, San Francisco, United States
Henry van den Bedem
Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, San Francisco, United States; Atomwise Inc, San Francisco, United States
Structural Biology Initiative, CUNY Advanced Science Research Center, New York, United States; Department of Chemistry and Biochemistry, City College of New York, New York, United States; Ph.D. Programs in Biochemistry, Biology and Chemistry, The Graduate Center, City University of New York, New York, United States
In their folded state, biomolecules exchange between multiple conformational states that are crucial for their function. Traditional structural biology methods, such as X-ray crystallography and cryogenic electron microscopy (cryo-EM), produce density maps that are ensemble averages, reflecting molecules in various conformations. Yet, most models derived from these maps explicitly represent only a single conformation, overlooking the complexity of biomolecular structures. To accurately reflect the diversity of biomolecular forms, there is a pressing need to shift toward modeling structural ensembles that mirror the experimental data. However, the challenge of distinguishing signal from noise complicates manual efforts to create these models. In response, we introduce the latest enhancements to qFit, an automated computational strategy designed to incorporate protein conformational heterogeneity into models built into density maps. These algorithmic improvements in qFit are substantiated by superior Rfree and geometry metrics across a wide range of proteins. Importantly, unlike more complex multicopy ensemble models, the multiconformer models produced by qFit can be manually modified in most major model building software (e.g., Coot) and fit can be further improved by refinement using standard pipelines (e.g., Phenix, Refmac, Buster). By reducing the barrier of creating multiconformer models, qFit can foster the development of new hypotheses about the relationship between macromolecular conformational dynamics and function.