Elucidating structural variability in p53 conformers using combinatorial refinement strategies and molecular dynamics

Low molecular weight proteins and protein assemblies can now be investigated using cryo-electron microscopy (EM) as a complement to traditional structural biology techniques. It is important, however, to not lose sight of the dynamic information inherent in macromolecules that give rise to their exquisite functionality. As computational methods continue to advance the field of biomedical imaging, so must strategies to resolve the minute details of disease-related entities. Here, we employed combinatorial modeling approaches to assess flexible properties among low molecular weight proteins (~100 kDa or less). Through a blend of rigid body refinement and simulated annealing, we determined new hidden conformations for wild type p53 monomer and dimer forms. Structures for both states converged to yield new conformers, each revealing good stereochemistry and dynamic information about the protein. Based on these insights, we identified fluid parts of p53 that complement the stable central core of the protein responsible for engaging DNA. Molecular dynamics simulations corroborated the modeling results and helped pinpoint the more flexible residues in wild type p53. Overall, the new computational methods may be used to shed light on other small protein features in a vast ensemble of structural data that cannot be easily delineated by other algorithms.

Delineating Conformational Variability in Small Protein Structures Using Combinatorial Refinement Strategies

As small protein assemblies and even small proteins are becoming more amenable to cryo-Electron Microscopy (EM) structural studies, it is important to consider the complementary dynamic information present in the data. Current computational strategies are limited in their ability to resolve minute differences among low molecular weight entities. Here, we demonstrate a new combinatorial approach to delineate flexible conformations among small proteins using real-space refinement applications. We performed a meta-analysis of structural data for the SARS CoV-2 Nucleocapsid (N) protein using a combination of rigid-body refinement and simulated annealing methods. For the N protein monomer, we determined three new flexible conformers with good stereochemistry and quantitative comparisons provided new evidence of their dynamic properties. A similar analysis performed for the N protein dimer showed only minor structural differences among the flexible models. These results suggested a more stable view of the N protein dimer than the monomer structure. Taken together, the new computational strategies can delineate conformational changes in low molecular weight proteins that may go unnoticed by conventional assessments. The results also suggest that small proteins may be further stabilized for structural studies through the use of solution components that limit the movement of external flexible regions.

Analytic modeling of inhomogeneous-resolution maps in cryo-electron microscopy and crystallography

Refinement of macromolecular atomic models versus experimental maps in crystallography and cryo-electron microscopy is a critical step in structure solution. For an appropriate comparison, model maps should mimic the imperfections in the experimental maps, mainly atomic disorder and limited resolution, which are often inhomogeneous over the molecular region. In the suggested method, these model maps are calculated as the sum of atomic contributions expressed through a specifically designed function describing a solitary spherical wave. Thanks to this function, atomic contributions are analytically expressed through their atomic displacement parameter and local resolution, a value now associated with each atom. Such a full analytic dependence of inhomogeneous-resolution map values on model parameters permits the refinement of all of these parameters together.

Making ripples in the comparison of calculated and experimental maps for real-space refinement and assessment by analytic modelling of local resolution

Commentary is given on a paper [Urzhumtsev & Lunin (2022). IUCrJ, 9, 728-734] proposing a new method for the analytic modelling of inhomogeneous resolution in electrostatic potential volumes and electron density maps for improved real-space refinement.

Macromolecular refinement of X-ray and cryoelectron microscopy structures with Phenix/OPLS3e for improved structure and ligand quality

With the advent of the resolution revolution in cryoelectron microscopy (cryo-EM), low-resolution refinement is common, and likewise increases the need for a reliable force field. Here, we report on the incorporation of the OPLS3e force field with the VSGB2.1 solvation model in the structure determination package Phenix. Our results show significantly improved structure quality and reduced ligand strain at lower resolution for X-ray refinement. For refinement of cryo-EM-based structures, we find comparable quality structures, goodness-of-fit, and reduced ligand strain. We also show how structure quality and ligand strain are related to the map-model cross-correlation as a function of data weight, and how that can detect overfitting. Signs of overfitting are found in over half of our cryo-EM dataset, which can be remedied by a re-refinement at a lower data weight. Finally, a start-to-end script for refining structures with Phenix/OPLS3e is available in the Schrödinger 2020-3 distribution.

Real-space quantum-based refinement for cryo-EM: Q|R#3

Electron cryo-microscopy (cryo-EM) is rapidly becoming a major competitor to X-ray crystallography, especially for large structures that are difficult or impossible to crystallize. While recent spectacular technological improvements have led to significantly higher resolution three-dimensional reconstructions, the average quality of cryo-EM maps is still at the low-resolution end of the range compared with crystallography. A long-standing challenge for atomic model refinement has been the production of stereochemically meaningful models for this resolution regime. Here, it is demonstrated that including accurate model geometry restraints derived from ab initio quantum-chemical calculations (HF-D3/6-31G) can improve the refinement of an example structure (chain A of PDB entry 3j63). The robustness of the procedure is tested for additional structures with up to 7000 atoms (PDB entry 3a5x and chain C of PDB entry 5fn5) using the less expensive semi-empirical (GFN1-xTB) model. The necessary algorithms enabling real-space quantum refinement have been implemented in the latest version of qr.refine and are described here.

Real-space refinement in PHENIX for cryo-EM and crystallography

This article describes the implementation of real-space refinement in the phenix.real_space_refine program from the PHENIX suite. The use of a simplified refinement target function enables very fast calculation, which in turn makes it possible to identify optimal data-restraint weights as part of routine refinements with little runtime cost. Refinement of atomic models against low-resolution data benefits from the inclusion of as much additional information as is available. In addition to standard restraints on covalent geometry, phenix.real_space_refine makes use of extra information such as secondary-structure and rotamer-specific restraints, as well as restraints or constraints on internal molecular symmetry. The re-refinement of 385 cryo-EM-derived models available in the Protein Data Bank at resolutions of 6 Å or better shows significant improvement of the models and of the fit of these models to the target maps.

ISOLDE: a physically realistic environment for model building into low-resolution electron-density maps

This paper introduces ISOLDE, a new software package designed to provide an intuitive environment for high-fidelity interactive remodelling/refinement of macromolecular models into electron-density maps. ISOLDE combines interactive molecular-dynamics flexible fitting with modern molecular-graphics visualization and established structural biology libraries to provide an immersive interface wherein the model constantly acts to maintain physically realistic conformations as the user interacts with it by directly tugging atoms with a mouse or haptic interface or applying/removing restraints. In addition, common validation tasks are accelerated and visualized in real time. Using the recently described 3.8 Å resolution cryo-EM structure of the eukaryotic minichromosome maintenance (MCM) helicase complex as a case study, it is demonstrated how ISOLDE can be used alongside other modern refinement tools to avoid common pitfalls of low-resolution modelling and improve the quality of the final model. A detailed analysis of changes between the initial and final model provides a somewhat sobering insight into the dangers of relying on a small number of validation metrics to judge the quality of a low-resolution model.

Taura syndrome virus IRES initiates translation by binding its tRNA-mRNA-like structural element in the ribosomal decoding center

In cap-dependent translation initiation, the open reading frame (ORF) of mRNA is established by the placement of the AUG start codon and initiator tRNA in the ribosomal peptidyl (P) site. Internal ribosome entry sites (IRESs) promote translation of mRNAs in a cap-independent manner. We report two structures of the ribosome-bound Taura syndrome virus (TSV) IRES belonging to the family of Dicistroviridae intergenic IRESs. Intersubunit rotational states differ in these structures, suggesting that ribosome dynamics play a role in IRES translocation. Pseudoknot I of the IRES occupies the ribosomal decoding center at the aminoacyl (A) site in a manner resembling that of the tRNA anticodon-mRNA codon. The structures reveal that the TSV IRES initiates translation by a previously unseen mechanism, which is conceptually distinct from initiator tRNA-dependent mechanisms. Specifically, the ORF of the IRES-driven mRNA is established by the placement of the preceding tRNA-mRNA-like structure in the A site, whereas the 40S P site remains unoccupied during this initial step.

Structure of the ribosome with elongation factor G trapped in the pretranslocation state

During protein synthesis, tRNAs and their associated mRNA codons move sequentially on the ribosome from the A (aminoacyl) site to the P (peptidyl) site to the E (exit) site in a process catalyzed by a universally conserved ribosome-dependent GTPase [elongation factor G (EF-G) in prokaryotes and elongation factor 2 (EF-2) in eukaryotes]. Although the high-resolution structure of EF-G bound to the posttranslocation ribosome has been determined, the pretranslocation conformation of the ribosome bound with EF-G and A-site tRNA has evaded visualization owing to the transient nature of this state. Here we use electron cryomicroscopy to determine the structure of the 70S ribosome with EF-G, which is trapped in the pretranslocation state using antibiotic viomycin. Comparison with the posttranslocation ribosome shows that the small subunit of the pretranslocation ribosome is rotated by ∼12° relative to the large subunit. Domain IV of EF-G is positioned in the cleft between the body and head of the small subunit outwardly of the A site and contacts the A-site tRNA. Our findings suggest a model in which domain IV of EF-G promotes the translocation of tRNA from the A to the P site as the small ribosome subunit spontaneously rotates back from the hybrid, rotated state into the nonrotated posttranslocation state.

New model-fitting and model-completion programs for automated iterative nucleic acid refinement

In the past decade many structures of nucleic acids have been determined, which have contributed to our understanding of their biological functions. However, crystals containing nucleic acids often diffract X-rays poorly. This makes electron-density interpretation difficult and requires a great deal of expertise in crystallography and knowledge of nucleic acid structure. Here, new programs called NAFIT and NABUILD for fitting and extending nucleic acid models are presented. These programs can be used as modules in the automated refinement system LAFIRE, as well as acting as independent programs. NAFIT performs sequential grouped fitting with empirical torsion-angle restraints and antibumping restraints including H atoms. NABUILD extends the model using a skeletonized map in a coarse-grained manner. It has been shown that NAFIT greatly improves electron-density fit and geometric quality and that iterative refinement with NABUILD significantly reduces the Rfree factor.

Application of DEN refinement and automated model building to a difficult case of molecular-replacement phasing: the structure of a putative succinyl-diaminopimelate desuccinylase from Corynebacterium glutamicum

Phasing by molecular replacement remains difficult for targets that are far from the search model or in situations where the crystal diffracts only weakly or to low resolution. Here, the process of determining and refining the structure of Cgl1109, a putative succinyl-diaminopimelate desuccinylase from Corynebacterium glutamicum, at ∼3 Å resolution is described using a combination of homology modeling with MODELLER, molecular-replacement phasing with Phaser, deformable elastic network (DEN) refinement and automated model building using AutoBuild in a semi-automated fashion, followed by final refinement cycles with phenix.refine and Coot. This difficult molecular-replacement case illustrates the power of including DEN restraints derived from a starting model to guide the movements of the model during refinement. The resulting improved model phases provide better starting points for automated model building and produce more significant difference peaks in anomalous difference Fourier maps to locate anomalous scatterers than does standard refinement. This example also illustrates a current limitation of automated procedures that require manual adjustment of local sequence misalignments between the homology model and the target sequence.

Protein structure fitting and refinement guided by cryo-EM density

For many macromolecular assemblies, both a cryo-electron microscopy map and atomic structures of its component proteins are available. Here we describe a method for fitting and refining a component structure within its map at intermediate resolution (<15 A). The atomic positions are optimized with respect to a scoring function that includes the crosscorrelation coefficient between the structure and the map as well as stereochemical and nonbonded interaction terms. A heuristic optimization that relies on a Monte Carlo search, a conjugate-gradients minimization, and simulated annealing molecular dynamics is applied to a series of subdivisions of the structure into progressively smaller rigid bodies. The method was tested on 15 proteins of known structure with 13 simulated maps and 3 experimentally determined maps. At approximately 10 A resolution, Calpha rmsd between the initial and final structures was reduced on average by approximately 53%. The method is automated and can refine both experimental and predicted atomic structures.