Artificial Intelligence: Exploring the conformational diversity of proteins
The human body contains a vast number of different proteins that carry out distinct roles. Proteins are made up of combinations of 20 amino acids, each with different physicochemical properties. The number and sequence of amino acids in a protein determine how it will fold into the specific three-dimensional structure or shape that the protein needs to perform its role.
Several proteins, including membrane proteins, do not simply fold into a single conformation. Instead, they need to be able to ‘flip’ between different conformations to do their job. Innovations in the experimental techniques used to determine protein structures – such as cryo-electron microscopy, nuclear magnetic resonance spectroscopy or X-ray crystallography –have provided valuable insights into the different conformations of many membrane proteins. However, these methods are costly and time consuming.
Using computational methods to predict the structures of proteins could allow scientists to fill the gap between protein sequence and structural knowledge, without having to rely on expensive experimental methods (Baker and Sali, 2001). Recently, an artificial intelligence-based method to predict protein structures, called AlphaFold2 (AF2), has taken structural biology by storm (Jumper et al., 2021).
AF2 emerged as a valuable tool for predicting the structures of proteins from their sequences with an accuracy comparable to that obtained by experimental techniques at a fraction of their time and costs, as shown for various biological problems (Evans et al., 2021; Mosalaganti et al., 2021; Tunyasuvunakool et al., 2021; McCoy et al., 2022). Now, in eLife, Diego del Alamo, Davide Sala, Hassane Mchaourab and Jens Meiler report how AF2 can also predict different conformations of membrane proteins (Del Alamo et al., 2022).
To do so, the researchers – based at Vanderbilt University and Leipzig University – used a set of eight membrane proteins representing different structural classes and mechanisms of action. This included five unique transporters (LAT1, ZnT8, MCT1, STP10, and ASCT2), whose structures had been previously experimentally determined in both inward- and outward-facing conformations (Figure 1), and three representative G-protein-coupled receptors (CGRPR, PTH1R, and FZD7), whose structures had been solved experimentally in active and inactive states. None of these proteins were part of the original AF2 training set, which included structures located in the Protein Data Bank (PDB).
The researchers then tested the ability of AF2 to model distinct conformational states by varying different parameters of the predictor, such as the number of models generated, and by using known structures of protein homologues as templates. One key feature of AF2 is its ability to generate a multiple sequence alignment (MSA) of evolutionarily related sequences, which is critical for accurate modeling. These MSAs, which can include thousands of sequences, are used by AF2 to identify residues that have co-evolved, thereby highlighting the contacts critical for defining the three-dimensional fold of the protein.
Remarkably, del Alamo et al. demonstrated that by reducing the size of the input MSA or alignment depth from 5,120 to as few as 16 sequences, the conformational diversity explored by AF2 increased, thereby capturing the structures that were experimentally determined in different conformations. This procedure also generated misfolded or outlier models, which were identified by the lack of structural similarity to other models and excluded from further analysis. This provides an important step to distinguish structural models that represent biologically relevant states. Moreover, in some cases using templates as input increased the conformational diversity of the generated models when MSAs with reduced number of aligned sequences (shallow MSAs) were used. Taken together, these results suggest that minor modifications to the input parameters allow AF2 to explore a larger area of the conformational space of proteins to capture distinct, biologically relevant states.
However, the analysis was performed on a relatively small benchmark set of proteins, due to the limited number of membrane protein structures not included in the AF2 training set and resolved in multiple states. Furthermore, del Alamo et al. did not identify a one-size-fits-all protocol that could accurately model the conformational diversity of all the membrane proteins in their benchmark set. A more generalized approach would be useful to study a larger variety of proteins that adopt different conformations, including enzymes and transcription factors. Finally, given the significant role of membrane proteins as drug targets, it will be crucial to assess whether the models generated with the proposed approach can be used for rational drug design, which typically requires accurate modeling of the protein’s amino acid sidechains.
In conclusion, the work by del Alamo et al. extends the scope of AF2 beyond structure prediction of a single state to the exploration of the conformational diversity of proteins. Even though determining the populations of alternative conformations and the interconversion pathways between them still appears to be out of reach, this work represents a crucial step towards describing the dynamic nature of proteins with modern artificial intelligence-based structure predictors.
References
-
Implications of AlphaFold2 for crystallographic phasing by molecular replacementActa Crystallographica. Section D, Structural Biology 78:1–13.https://doi.org/10.1107/S2059798321012122
Article and author information
Author details
Publication history
Copyright
© 2022, Schlessinger and Bonomi
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 5,048
- views
-
- 465
- downloads
-
- 12
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Structural Biology and Molecular Biophysics
The M2 proton channel aids in the exit of mature influenza viral particles from the host plasma membrane through its ability to stabilize regions of high negative Gaussian curvature (NGC) that occur at the neck of budding virions. The channels are homo-tetramers that contain a cytoplasm-facing amphipathic helix (AH) that is necessary and sufficient for NGC generation; however, constructs containing the transmembrane spanning helix, which facilitates tetramerization, exhibit enhanced curvature generation. Here, we used all-atom molecular dynamics (MD) simulations to explore the conformational dynamics of M2 channels in lipid bilayers revealing that the AH is dynamic, quickly breaking the fourfold symmetry observed in most structures. Next, we carried out MD simulations with the protein restrained in four- and twofold symmetric conformations to determine the impact on the membrane shape. While each pattern was distinct, all configurations induced pronounced curvature in the outer leaflet, while conversely, the inner leaflets showed minimal curvature and significant lipid tilt around the AHs. The MD-generated profiles at the protein–membrane interface were then extracted and used as boundary conditions in a continuum elastic membrane model to calculate the membrane-bending energy of each conformation embedded in different membrane surfaces characteristic of a budding virus. The calculations show that all three M2 conformations are stabilized in inward-budding, concave spherical caps and destabilized in outward-budding, convex spherical caps, the latter reminiscent of a budding virus. One of the C2-broken symmetry conformations is stabilized by 4 kT in NGC surfaces with the minimum energy conformation occurring at a curvature corresponding to 33 nm radii. In total, our work provides atomistic insight into the curvature sensing capabilities of M2 channels and how enrichment in the nascent viral particle depends on protein shape and membrane geometry.
-
- Structural Biology and Molecular Biophysics
Small-molecule drug design hinges on obtaining co-crystallized ligand-protein structures. Despite AlphaFold2’s strides in protein native structure prediction, its focus on apo structures overlooks ligands and associated holo structures. Moreover, designing selective drugs often benefits from the targeting of diverse metastable conformations. Therefore, direct application of AlphaFold2 models in virtual screening and drug discovery remains tentative. Here, we demonstrate an AlphaFold2-based framework combined with all-atom enhanced sampling molecular dynamics and Induced Fit docking, named AF2RAVE-Glide, to conduct computational model-based small-molecule binding of metastable protein kinase conformations, initiated from protein sequences. We demonstrate the AF2RAVE-Glide workflow on three different mammalian protein kinases and their type I and II inhibitors, with special emphasis on binding of known type II kinase inhibitors which target the metastable classical DFG-out state. These states are not easy to sample from AlphaFold2. Here, we demonstrate how with AF2RAVE these metastable conformations can be sampled for different kinases with high enough accuracy to enable subsequent docking of known type II kinase inhibitors with more than 50% success rates across docking calculations. We believe the protocol should be deployable for other kinases and more proteins generally.