From cryo-EM, multiple protein structures in one shot

Single particle electron cryomicroscopy (cryo-EM) has recently made significant progress in high-resolution structure determination of macromolecular complexes due to improvements in electron microscopic instrumentation and computational image analysis. However, cryo-EM structures can be highly non-uniform in local resolution and all structures available to date have been limited to resolutions above 3Å. Here we present the cryo-EM structure of the 70S ribosome from Escherichia coli in complex with elongation factor Tu, aminoacyl-tRNA and the antibiotic kirromycin at 2.65-2.9Å resolution using spherical aberration (Cs)-corrected cryo-EM. Overall, the cryo-EM reconstruction at 2.9Å resolution is comparable to the best-resolved X-ray structure of the E. coli 70S ribosome (2.8Å), but provides more detailed information (2.65Å) at the functionally important ribosomal core. The cryo-EM map elucidates for the first time the structure of all 35 rRNA modifications in the bacterial ribosome, explaining their roles in fine-tuning ribosome structure and function and modulating the action of antibiotics. We also obtained atomic models for flexible parts of the ribosome such as ribosomal proteins L9 and L31. The refined cryo-EM-based model presents the currently most complete high-resolution structure of the E. coli ribosome, which demonstrates the power of cryo-EM in structure determination of large and dynamic macromolecular complexes.

Cryogenic electron microscopy (cryo-EM) has become in the past 10 years one of the major tools for the structure determination of proteins. Nowadays, the structure prediction field is experiencing the same revolution and, using AlphaFold2, it is possible to have high-confidence atomic models for virtually any polypeptide chain, smaller than 4000 amino acids, in a simple click. Even in a scenario where all polypeptide chain folding were to be known, cryo-EM retains specific characteristics that make it a unique tool for the structure determination of macromolecular complexes. Using cryo-EM, it is possible to obtain near-atomic structures of large and flexible mega-complexes, describe conformational panoramas, and potentially develop a structural proteomic approach from fully ex vivo specimens.

Volta phase plate data collection facilitates image processing and cryo-EM structure determination.

Cryo-electron microscopy (cryoEM) is a powerful technique for structure determination of macromolecular complexes, via single particle analysis (SPA). The overall process involves i) vitrifying the specimen in a thin film supported on a cryoEM grid; ii) screening the specimen to assess particle distribution and ice quality; iii) if the grid is suitable, collecting a single particle dataset for analysis; and iv) image processing to yield an EM density map. In this protocol, an overview for each of these steps is provided, with a focus on the variables which a user can modify during the workflow and the troubleshooting of common issues. With remote microscope operation becoming standard in many facilities, variations on imaging protocols to assist users in efficient operation and imaging when physical access to the microscope is limited will be described.

Cryogenic electron tomography (cryo-ET) combined with subtomogram averaging, allows in situ visualization and structure determination of macromolecular complexes at subnanometre resolution. Cryogenic focused ion beam (cryo-FIB) micromachining is used to prepare a thin lamella-shaped sample out of a frozen-hydrated cell for cryo-ET imaging, but standard cryo-FIB fabrication is blind to the precise location of the structure or proteins of interest. Fluorescence-guided focused ion beam (FIB) milling at target locations requires multiple sample transfers prone to contamination, and relocation and registration accuracy is often insufficient for 3D targeting. Here, we present in situ fluorescence microscopy-guided FIB fabrication of a frozen-hydrated lamella to address this problem: we built a coincident three-beam cryogenic correlative microscope by retrofitting a compact cryogenic microcooler, custom positioning stage, and an inverted widefield fluorescence microscope (FM) on an existing FIB scanning electron microscope. We show FM controlled targeting at every milling step in the lamella fabrication process, validated with transmission electron microscope tomogram reconstructions of the target regions. The ability to check the lamella during and after the milling process results in a higher success rate in the fabrication process and will increase the throughput of fabrication for lamellae suitable for high-resolution imaging.

Nowadays, more and more attention has been attracted to develop selective PI3Kγ inhibitors, but the unique structural features of PI3Kγ protein make it a very big challenge. In the present study, a virtual screening strategy based on machine learning with multiple PI3Kγ protein structures was developed to screen novel PI3Kγ inhibitors. First, six mainstream docking programs were chosen to evaluate their scoring power and screening power; CDOCKER and Glide show satisfactory reliability and accuracy against the PI3Kγ system. Next, virtual screening integrating multiple PI3Kγ protein structures was demonstrated to significantly improve the screening enrichment rate comparing to that with an individual protein structure. Last, a multi-conformational Naïve Bayesian Classification model with the optimal docking programs was constructed, and it performed a true capability in the screening of PI3Kγ inhibitors. Taken together, the current study could provide some guidance for the docking-based virtual screening to discover novel PI3Kγ inhibitors.

We have developed a receptor-based pharmacophore method which utilizes a collection of protein structures to account for inherent protein flexibility in structure-based drug design. Several procedures were systematically evaluated to derive the most general protocol for using multiple protein structures. Most notably, incorporating more protein flexibility improved the performance of the method. The pharmacophore models successfully discriminate known inhibitors from drug-like non-inhibitors. Furthermore, the models correctly identify the bound conformations of some ligands. We used unliganded HIV-1 protease to develop and validate this method. Drug design is always initiated with a protein-ligand structure, and such success with unbound protein structures is remarkable - particularly in the case of HIV-1 protease, which has a large conformational change upon binding. This technique holds the promise of successful computer-based drug design before bound crystal structures are even discovered, which can mean a jump-start of 1-3 years in tackling some medically relevant systems with computational methods.

Pose prediction and virtual screening performance of a molecular docking method depend on the choice of protein structures used for docking. Multiple structures for a target protein are often used to take into account the receptor flexibility and problems associated with a single receptor structure. However, the use of multiple receptor structures is computationally expensive when docking a large library of small molecules. Here, we propose a new cross-docking pipeline suitable to dock a large library of molecules while taking advantage of multiple target protein structures. Our method involves the selection of a suitable receptor for each ligand in a screening library utilizing ligand 3D shape similarity with crystallographic ligands. We have prospectively evaluated our method in D3R Grand Challenge 2 and demonstrated that our cross-docking pipeline can achieve similar or better performance than using either single or multiple-receptor structures. Moreover, our method displayed not only decent pose prediction performance but also better virtual screening performance over several other methods.

Automatic comparison and classification of protein structures

Developing methods to incorporate protein flexibility into structure-based drug design is an important challenge. Our approach uses multiple protein structures (MPS) to create a receptor-based pharmacophore model of the desired target. We have previously demonstrated the success of the method by applying it to human immunodeficiency virus-1 protease (HIV-1p). Our models, based on an apo structure, discriminated known HIV-1p inhibitors from druglike inactive compounds and also accurately identified bound conformations of known inhibitors. Here, we test the method by applying it to all three unbound crystal structures of HIV-1p. We have also improved our method with denser probe mapping of the binding site and refined our selection criteria for pharmacophore elements. Our improved protocol has led to the development of a consistent 8-site pharmacophore model for HIV-1p, which is independent of starting structure, and a robust MPS pharmacophore method that is more amenable to automation.

Anaplastic lymphoma kinase (ALK) has been regarded as a promising target for the therapy of various cancers. A large number of ALK inhibitors with diverse scaffolds have been discovered, and most of them belong to Type-I inhibitors that only occupy the ATP-binding pocket. Recently, we reported a series of novel and potent Type-I1/2 inhibitors of ALK with the 1-purine-3-piperidinecarboxamide scaffold, which can bind to both the ATP-binding site of ALK and the adjacent hydrophobic allosteric pocket. In this study, the binding mechanisms of these Type-I1/2 ALK inhibitors were elucidated by multiple molecular modeling techniques. The calculation results demonstrate that the ensemble docking based on multiple protein structures and the MM/PB(GB)SA calculations based on molecular dynamics (MD) simulations yield better predictions than conventional rigid receptor docking (Glide, Surflex-Dock, and Autodock Vina), highlighting the importance of incorporating receptor flexibility in the predictions of binding poses and binding affinities of Type-I1/2 ALK inhibitors. Furthermore, the umbrella sampling (US) simulations and MM/GBSA binding free energy decomposition analyses indicate that Leu1122, Leu1198, Gly1202 and Glu1210 in the hinge region and Glu1197, Ile1171, Phe1174, Ile1179, His1247, Ile1268, Asp1270 and Phe1271 in the allosteric pocket of ALK are the key residues for determining the relative binding strength of the studied inhibitors. Besides, we found that the most potent inhibitor (001-017) tends to form stronger transient interactions with residues along the dissociation channel due to the high electronegativity of its bulky 4-(trifluoromethoxy) phenylamine tail. As a whole, both the stronger binding affinity and the higher energetic barrier (which may prolong the drug-target residence time) of 001-017 contribute to its excellent anti-proliferation activity against ALK-positive cancer cells.

As an extension to a previous published study (McGaughey et al., J Chem Inf Model 47:1504-1519, 2007) comparing 2D and 3D similarity methods to docking, we apply a subset of those virtual screening methods (TOPOSIM, SQW, ROCS-color, and Glide) to a set of protein/ligand pairs where the protein is the target for docking and the cocrystallized ligand is the target for the similarity methods. Each protein is represented by a maximum of five crystal structures. We search a diverse subset of the MDDR as well as a diverse small subset of the MCIDB, Merck's proprietary database. It is seen that the relative effectiveness of virtual screening methods, as measured by the enrichment factor, is highly dependent on the particular crystal structure or ligand, and on the database being searched. 2D similarity methods appear very good for the MDDR, but poor for the MCIDB. However, ROCS-color (a 3D similarity method) does well for both databases.

Enzymes are highly efficient bio-catalysts interesting for industries and medicine. Therefore, a goal of utmost importance in biochemical research is to understand how an enzyme catalyzes a chemical reaction. Here, the computational identification of functionally or structurally important residue positions can be of tremendous help. The datasets that are most informative for the algorithms are the 3D structure of a protein and a multiple sequence alignment (MSA) composed of homologous sequences. For example, an MSA allows for the quantification of residue conservation. Residue conservation at a given position indicates that only one type of amino acid fulfills all constraints imposed by protein structure or function. Furthermore, a detailed analysis of less strictly conserved residue positions may identify pairs, whose orchestration is mutually dependent and induces correlated mutations. Both of these conservation signals are indicative of functionally or structurally important positions. In the first part of this thesis, methods of machine learning were used to identify and classify these residue positions. It was the aim to predict in a mutually exclusively manner a role in catalysis, ligand-binding or protein stability for each residue position of a protein. Unfortunately, for many proteins the 3D structure is unknown. For other proteins, the number of known homologs is not sufficient to compile a meaningful MSA. Therefore, three variants of a classifier were designed and implemented, named CLIPS-1D, CLIPS-3D, and CLIPS-4D. These multi-class support vector machines allow for a classification based on an MSA (CLIPS-1D), a 3D structure (CLIPS-3D), and a combination of both (CLIPS-4D). CLIPS-1D exploits seven sequence-based features, whereas CLIPS-3D utilizes seven structure-based features. CLIPS-4D combines the seven sequence-based features of CLIPS-1D with those two structure-based features that increased its classification performance. A comparison with existing methods and a detailed analysis on a well-studied enzyme confirmed state-of-the-art prediction quality for CLIPS-1D and CLIPS-4D. In the second part of this thesis an algorithm for the identification of correlated mutations was improved. A common method for the identification of correlated mutations is to deduce the mutual information (MI) of a pair of residue positions from an MSA. The classical MI is based on Shannon’s information theory that utilizes probabilities only. Consequently, these approaches do not consider the similarity of residue pairs, which is a severe limitation. In order to improve these algorithms, H2rs was developed for this thesis. Thus, the MIvalues originate from the von Neumann entropy (vNE), which takes into account amino acid similarities modeled by means of a substitution matrix. To further improve the specificity of H2rs, the significance of MIvNE-values was assessed with a bootstrapping approach. The analysis of a large in silico testbed and the detailed assessment of five well-studied enzymes demonstrated state-of-the-art performance.

Revealing protein binding sites with other molecules, such as nucleic acids, peptides, or small ligands, sheds light on disease mechanism elucidation and novel drug design. With the explosive growth of proteins in sequence databases, how to accurately and efficiently identify these binding sites from sequences becomes essential. However, current methods mostly rely on expensive multiple sequence alignments or experimental protein structures, limiting their genome-scale applications. Besides, these methods haven’t fully explored the geometry of the protein structures. Here, we propose GPSite, a multi-task network for simultaneously predicting binding residues of DNA, RNA, peptide, protein, ATP, HEM, and metal ions on proteins. GPSite was trained on informative sequence embeddings and predicted structures from protein language models, while comprehensively extracting residual and relational geometric contexts in an end-to-end manner. Experiments demonstrate that GPSite substantially surpasses state-of-the-art sequence-based and structure-based approaches on various benchmark datasets, even when the structures are not well-predicted. The low computational cost of GPSite enables rapid genome-scale binding residue annotations for over 568,000 sequences, providing opportunities to unveil unexplored associations of binding sites with molecular functions, biological processes, and genetic variants. The GPSite webserver and annotation database can be freely accessed at https://bio-web1.nscc-gz.cn/app/GPSite.

We present a new structurally derived pair-to-pair substitution matrix (P2PMAT). This matrix is constructed from a very large amount of integrated high quality multiple sequence alignments (Blocks) and protein structures. It evaluates the likelihoods of all 160,000 pair-to-pair substitutions. P2PMAT matrix implicitly accounts for evolutionary conservation, correlated mutations, and residue-residue contact potentials. The usefulness of the matrix for structural predictions is shown in this article. Predicting protein residue-residue contacts from sequence information alone, by our method (P2PConPred) is particularly accurate in the protein cores, where it performs better than other basic contact prediction methods (increasing accuracy by 25-60%). The method mean accuracy for protein cores is 24% for 59 diverse families and 34% for a subset of proteins shorter than 100 residues. This is above the level that was recently shown to be sufficient to significantly improve ab initio protein structure prediction. We also demonstrate the ability of our approach to identify native structures within large sets of (300-2000) protein decoys. On the basis of evolutionary information alone our method ranks the native structure in the top 0.3% of the decoys in 4/10 of the sets, and in 8/10 of sets the native structure is ranked in the top 10% of the decoys. The method can, thus, be used to assist filtering wrong models, complementing traditional scoring functions.