Modeling conformational states of proteins with AlphaFold

Unveiling the influence of fastest nobel prize winner discovery: alphafold's algorithmic intelligence in medical sciences

CONTEXT

AlphaFold's advanced AI technology has transformed protein structure interpretation. By predicting three-dimensional protein structures from amino acid sequences, AlphaFold has solved the complex protein-folding problem, previously challenging for experimental methods due to numerous possible conformations. Since its inception, AlphaFold has introduced several versions, including AlphaFold2, AlphaFold DB, AlphaFold Multimer, Alpha Missense, and AlphaFold3, each further enhancing protein structure prediction. Remarkably, AlphaFold is recognized as the fastest Nobel Prize winner in science history. This technology has extensive applications, potentially transforming treatment and diagnosis in medical sciences by reducing drug design costs and time, while elucidating structural pathways of human body systems. Numerous studies have demonstrated how AlphaFold aids in understanding health conditions by providing critical information about protein mutations, abnormal protein-protein interactions, and changes in protein dynamics. Researchers have also developed new technologies and pipelines using different versions of AlphaFold to amplify its potential. However, addressing existing limitations is crucial to maximizing AlphaFold's capacity to redefine medical research. This article reviews AlphaFold's impact on five key aspects of medical sciences: protein mutation, protein-protein interaction, molecular dynamics, drug design, and immunotherapy.

METHODS

This review examines the contributions of various AlphaFold versions AlphaFold2, AlphaFold DB, AlphaFold Multimer, Alpha Missense, and AlphaFold3 to protein structure prediction. The methods include an extensive analysis of computational techniques and software used in interpreting and predicting protein structures, emphasizing advances in AI technology and its applications in medical research.

Statistical Thermodynamics of the Protein Ensemble: Mediating Function and Evolution

The growing appreciation of native state conformational fluctuations mediating protein function calls for critical reevaluation of protein evolution and adaptation. If proteins are ensembles, does nature select solely for ground state structure, or are conformational equilibria between functional states also conserved? If so, what is the mechanism and how can it be measured? Addressing these fundamental questions, we review our investigation into the role of local unfolding fluctuations in the native state ensembles of proteins. We describe the functional importance of these ubiquitous fluctuations, as revealed through studies of adenylate kinase. We then summarize elucidation of thermodynamic organizing principles, which culminate in a quantitative probe for evolutionary conservation of protein energetics. Finally, we show that these principles are predictive of sequence compatibility for multiple folds, providing a unique thermodynamic perspective on metamorphic proteins. These research areas demonstrate that the locally unfolded ensemble is an emerging, important mechanism of protein evolution.

Protein Modeling with DEER Spectroscopy

Double electron-electron resonance (DEER) combined with site-directed spin labeling can provide distance distributions between selected protein residues to investigate protein structure and conformational heterogeneity. The utilization of the full quantitative information contained in DEER data requires effective protein and spin label modeling methods. Here, we review the application of DEER data to protein modeling. First, we discuss the significance of spin label modeling for accurate extraction of protein structural information and review the most popular label modeling methods. Next, we review several important aspects of protein modeling with DEER, including site selection, how DEER restraints are applied, common artifacts, and the unique potential of DEER data for modeling structural ensembles and conformational landscapes. Finally, we discuss common applications of protein modeling with DEER data and provide an outlook.

The physics-AI dialogue in drug design

A long path has led from the determination of the first protein structure in 1960 to the recent breakthroughs in protein science. Protein structure prediction and design methodologies based on machine learning (ML) have been recognized with the 2024 Nobel prize in Chemistry, but they would not have been possible without previous work and the input of many domain scientists. Challenges remain in the application of ML tools for the prediction of structural ensembles and their usage within the software pipelines for structure determination by crystallography or cryogenic electron microscopy. In the drug discovery workflow, ML techniques are being used in diverse areas such as scoring of docked poses, or the generation of molecular descriptors. As the ML techniques become more widespread, novel applications emerge which can profit from the large amounts of data available. Nevertheless, it is essential to balance the potential advantages against the environmental costs of ML deployment to decide if and when it is best to apply it. For hit to lead optimization ML tools can efficiently interpolate between compounds in large chemical series but free energy calculations by molecular dynamics simulations seem to be superior for designing novel derivatives. Importantly, the potential complementarity and/or synergism of physics-based methods (e.g., force field-based simulation models) and data-hungry ML techniques is growing strongly. Current ML methods have evolved from decades of research. It is now necessary for biologists, physicists, and computer scientists to fully understand advantages and limitations of ML techniques to ensure that the complementarity of physics-based methods and ML tools can be fully exploited for drug design.

Characterization of conformationally heterogeneous proteins by electron paramagnetic resonance spectroscopy

The Anfinsen paradigm of representing a protein by a single conformer is challenged by the uncertainty predictions that come with AlphaFold models, which suggest a greater extent of disorder. Characterization of such conformation heterogeneity requires experimental approaches that do not depend on long-range order. Site-directed spin labeling (SDSL) coupled with electron paramagnetic resonance (EPR) spectroscopy is such an approach. The double electron-electron resonance (DEER) technique can access site-pair distance distributions in the 15-100 Å range, directly informing on ensemble width. SDSL-EPR can be applied in cellular environments, and recent work indicates that protein disorder is even more pervasive than predicted by AlphaFold. This suggests that the Anfinsen paradigm should be replaced by an ensemble paradigm.

Modeling Boltzmann-weighted structural ensembles of proteins using artificial intelligence-based methods

This review highlights recent advances in AI-driven methods for generating Boltzmann-weighted structural ensembles, which are crucial for understanding biomolecular dynamics and drug discovery. With the rise of deep learning models such as AlphaFold2, there has been a shift toward more accurate and efficient sampling of structural ensembles. The review discusses the integration of AI with traditional molecular dynamics techniques as well as experiments, the challenges of conformational sampling, and future directions for AI-driven research in structural biology, particularly in drug discovery and protein dynamics.

Molecular basis of the functional conflict between chloroquine and peptide transport in the Malaria parasite chloroquine resistance transporter PfCRT

The Plasmodium falciparum chloroquine resistance transporter (PfCRT) is a key protein contributing to resistance against the antimalarial chloroquine (CQ). Mutations such as K76T enable PfCRT to transport CQ away from its target in the parasite's digestive vacuole, but this comes at a cost to its natural peptide transport function. This creates fitness costs which can drive changes to drug susceptibility in parasite populations, but the molecular basis of this is not well understood. To investigate, here we run 130 μs of molecular dynamics simulations of CQ-sensitive and CQ-resistant PfCRT isoforms with CQ and peptide substrates. We identify the CQ binding site and characterized diverse peptide binding modes. The K76T mutation allows CQ to access the binding site but disrupts peptide binding, highlighting the importance of cavity charge in determining substrate specificity. This study provides insight into PfCRT polyspecific peptide transport and will aid in rational, structure-based inhibitor design.

Evolving concepts of the protein universe

The protein universe is the collection of all proteins on earth from all organisms both extant and extinct. Classical studies on protein folding suggested that proteins exist as a unique three-dimensional conformation that is dictated by the genetic code and is critical for function. In this perspective, we discuss ideas and developments that emerged over the past three decades regarding the protein structure-function paradigm. It is now clear that ordered (active/functional) and disordered/denatured (and hence inactive/non-functional) represent a continuum of states rather than binary states. Some proteins can switch folds without sequence change. Others exist as conformational ensembles lacking defined structure yet play critical roles in many biological processes, including forming membrane-less organelles driven by liquid-liquid phase separation. Numerous diverse proteins harbor segments with the potential to form amyloid fibrils, many of which are functional, and some possess prion-like properties enabling conformation-based transfer of heritable information. Taken together, these developments reveal the remarkable complexity of the protein universe.

TopEC: prediction of Enzyme Commission classes by 3D graph neural networks and localized 3D protein descriptor

Tools available for inferring enzyme function from general sequence, fold, or evolutionary information are generally successful. However, they can lead to misclassification if a deviation in local structural features influences the function. Here, we present TopEC, a 3D graph neural network based on a localized 3D descriptor to learn chemical reactions of enzymes from enzyme structures and predict Enzyme Commission (EC) classes. Using message-passing frameworks, we include distance and angle information to significantly improve the predictive performance for EC classification (F-score: 0.72) compared to regular 2D graph neural networks. We trained networks without fold bias that can classify enzyme structures for a vast functional space (>800 ECs). Our model is robust to uncertainties in binding site locations and similar functions in distinct binding sites. We observe that TopEC networks learn from an interplay between biochemical features and local shape-dependent features. TopEC is available as a repository on GitHub: https://github.com/IBG4-CBCLab/TopEC and https://doi.org/10.25838/d5p-66 .

Conformational ensembles for protein structure prediction

Acquisition of conformational ensembles for a protein is a challenging task, which is actually involving to the solution for protein folding problem and the study of intrinsically disordered protein. Despite AlphaFold with artificial intelligence acquired unprecedented accuracy to predict structures, its result is limited to a single state of conformation and it cannot provide multiple conformations to display protein intrinsic disorder. To overcome the barrier, a FiveFold approach was developed with a single sequence method. It applied the protein folding shape code (PFSC) uniformly to expose local folds of five amino acid residues, formed the protein folding variation matrix (PFVM) to reveal local folding variations along sequence, obtained a massive number of folding conformations in PFSC strings, and then an ensemble of multiple conformational protein structures is constructed. The P53_HUMAN as a well-known protein and LEF1_HUMAN and Q8GT36_SPIOL as typical disordered proteins are token as the benchmark to evaluate the predicted outcomes. The results demonstrated an effective algorithm and biological meaningful process well to predict protein multiple conformation structures.

Emerging approaches to investigating functional protein dynamics in modular redox enzymes: Nitric oxide synthase as a model system

Approximately 80% of eukaryotic and 65% of prokaryotic proteins are composed of multiple folding units (i.e., domains) connected by flexible linkers. These dynamic protein architectures enable diverse, essential functions such as electron transfer, respiration, and biosynthesis. This review critically assesses recent advancements in methods for studying protein dynamics, with a particular focus on modular, multidomain nitric oxide synthase (NOS) enzymes. Moving beyond traditional static "snapshots" of protein structures, current research emphasizes the dynamic nature of proteins, viewing them as flexible architectures modulated by conformational changes and interactions. In this context, the review discusses key developments in the integration of quantitative crosslinking mass spectrometry (qXL MS) with AlphaFold 2 predictions, which provides a powerful approach to disentangling NOS structural dynamics and understanding their modulation by external regulatory cues. Additionally, advances in site-specific infrared (IR) spectroscopy offer exciting potential in providing rich details about the conformational dynamics of NOSs in docked states. Moreover, optimization of genetic code expansion machinery enables the generation of genuine phosphorylated NOS enzymes, paving the way for detailed biophysical and functional analyses of phosphorylation's role in shaping NOS activity and structural flexibility; notably, this approach also empowers site-specific IR probe labeling with cyano groups. By embracing and leveraging AI-driven tools like AlphaFold 2 for structural and conformational modeling, alongside solution-based biophysical methods such as qXL MS and site-specific IR spectroscopy, researchers will gain integrative insights into functional protein dynamics. Collectively, these breakthroughs highlight the transformative potential of modern approaches in driving fundamental biological chemistry research.

Hidden Structural States of Proteins Revealed by Conformer Selection with AlphaFold-NMR

We introduce AlphaFold-NMR, a novel approach to NMR structure determination that reveals previously undetected protein conformational states. Unlike conventional NMR methods that rely on NOE-derived spatial restraints, AlphaFold-NMR combines AI-driven conformational sampling with Bayesian scoring of realistic protein models against NOESY and chemical shift data. This method uncovers alternative conformational states of the enzyme Gaussia luciferase, involving large-scale changes in the lid, binding pockets, and other surface cavities. It also identifies similar yet distinct conformational states of the human tumor suppressor Cyclin-Dependent Kinase 2-Associated Protein 1. These studies demonstrate the potential of AI-based modeling with enhanced sampling to generate diverse structural models followed by conformer selection and validation with experimental data as an alternative to traditional restraint-satisfaction protocols for protein NMR structure determination. The AlphaFold-NMR framework enables discovery of conformational heterogeneity and cryptic pockets that conventional NMR analysis methods do not distinguish, providing new insights into protein structure-function relationships.

AlphaFold prediction of structural ensembles of disordered proteins

Deep learning methods of predicting protein structures have reached an accuracy comparable to that of high-resolution experimental methods. It is thus possible to generate accurate models of the native states of hundreds of millions of proteins. An open question, however, concerns whether these advances can be translated to disordered proteins, which should be represented as structural ensembles because of their heterogeneous and dynamical nature. To address this problem, we introduce the AlphaFold-Metainference method to use AlphaFold-derived distances as structural restraints in molecular dynamics simulations to construct structural ensembles of ordered and disordered proteins. The results obtained using AlphaFold-Metainference illustrate the possibility of making predictions of the conformational properties of disordered proteins using deep learning methods trained on the large structural databases available for folded proteins.

Molecular basis of proton sensing by G protein-coupled receptors

Three proton-sensing G protein-coupled receptors (GPCRs)-GPR4, GPR65, and GPR68-respond to extracellular pH to regulate diverse physiology. How protons activate these receptors is poorly understood. We determined cryogenic-electron microscopy (cryo-EM) structures of each receptor to understand the spatial arrangement of proton-sensing residues. Using deep mutational scanning (DMS), we determined the functional importance of every residue in GPR68 activation by generating ∼9,500 mutants and measuring their effects on signaling and surface expression. Constant-pH molecular dynamics simulations provided insights into the conformational landscape and protonation patterns of key residues. This unbiased approach revealed that, unlike other proton-sensitive channels and receptors, no single site is critical for proton recognition. Instead, a network of titratable residues extends from the extracellular surface to the transmembrane region, converging on canonical motifs to activate proton-sensing GPCRs. Our approach integrating structure, simulations, and unbiased functional interrogation provides a framework for understanding GPCR signaling complexity.

An outlook on structural biology after AlphaFold: tools, limits and perspectives

AlphaFold and similar groundbreaking, AI-based tools, have revolutionized the field of structural bioinformatics, with their remarkable accuracy in ab-initio protein structure prediction. This success has catalyzed the development of new software and pipelines aimed at incorporating AlphaFold's predictions, often focusing on addressing the algorithm's remaining challenges. Here, we present the current landscape of structural bioinformatics shaped by AlphaFold, and discuss how the field is dynamically responding to this revolution, with new software, methods, and pipelines. While the excitement around AI-based tools led to their widespread application, it is essential to acknowledge that their practical success hinges on their integration into established protocols within structural bioinformatics, often neglected in the context of AI-driven advancements. Indeed, user-driven intervention is still as pivotal in the structure prediction process as in complementing state-of-the-art algorithms with functional and biological knowledge.

Proteins with alternative folds reveal blind spots in AlphaFold-based protein structure prediction

In recent years, advances in artificial intelligence (AI) have transformed structural biology, particularly protein structure prediction. Though AI-based methods, such as AlphaFold (AF), often predict single conformations of proteins with high accuracy and confidence, predictions of alternative folds are often inaccurate, low-confidence, or simply not predicted at all. Here, we review three blind spots that alternative conformations reveal about AF-based protein structure prediction. First, proteins that assume conformations distinct from their training-set homologs can be mispredicted. Second, AF overrelies on its training set to predict alternative conformations. Third, degeneracies in pairwise representations can lead to high-confidence predictions inconsistent with experiment. These weaknesses suggest approaches to predict alternative folds more reliably.

Structural basis of signaling complex inhibition by IL-6 domain-swapped dimers

Interleukin-6 (IL-6) is a multifaceted cytokine essential in many immune system processes and their regulation. It also plays a key role in hematopoiesis, and in triggering the acute phase reaction. IL-6 overproduction is critical in chronic inflammation associated with autoimmune diseases like rheumatoid arthritis and contributes to cytokine storms in COVID-19 patients. Over 20 years ago, researchers proposed that IL-6, which is typically monomeric, can also form dimers via a domain-swap mechanism, with indirect evidence supporting their existence. The physiological significance of IL-6 dimers was shown in B-cell chronic lymphocytic leukemia. However, no structures have been reported so far. Here, we present the crystal structure of an IL-6 domain-swapped dimer that computational approaches could not predict. The structure explains why the IL-6 dimer is antagonistic to the IL-6 monomer in signaling complex formation and provides insights for IL-6 targeted therapies.

Structural Biology for Target Identification and Validation

Structural biology is catalyzing a paradigm shift in drug discovery towards rational approaches in target identification and validation. Leveraging structural insights obtained through cryo-EM or X-ray crystallography not only enhances the efficiency of drug discovery projects in terms of time and cost, but also significantly improves the likelihood of achieving market approval.Initiating a successful project necessitates more than just a robust package for target credentialing; it demands a comprehensive strategy for the identification and optimization of potential drugs. The critical evaluation of target druggability is markedly enhanced when supported by experimentally derived structural information. This nuanced approach ensures a more thorough understanding of the technical feasibility of drug development from the project's inception.

An in silico framework to visualize how cancer-associated mutations influence structural plasticity of the chemokine receptor CCR3

G protein Coupled Receptors (GPCRs) are the largest family of cell surface receptors in humans. Somatic mutations in GPCRs are implicated in cancer progression and metastasis, but mechanisms are poorly understood. Emerging evidence implicates perturbation of intra-receptor activation pathway motifs whereby extracellular signals are transmitted intracellularly. Recently, sufficiently sensitive methodology was described to calculate structural strain as a function of missense mutations in AlphaFold-predicted model structures, which was extensively validated on experimental and predicted structural datasets. When paired with Molecular Dynamics (MD) simulations, these tools provide a facile approach to screen mutations in silico. We applied this framework to calculate the structural and dynamic effects of cancer-associated mutations in the chemokine receptor CCR3, a Class A GPCR involved in cancer and autoimmune disorders. Residue-residue contact scoring refined effective strain results, highlighting significant remodeling of inter- and intra-motif contacts along the highly conserved GPCR activation pathway network. We then integrated AlphaFold-derived predicted Local Distance Difference Test scores with per-residue Root Mean Square Fluctuations and activation pathway Contact Analysis (CONAN) from coarse grain MD simulations to identify statistically significant changes in receptor dynamics upon mutation. Finally, analysis of negative control mutants suggests false positive results in AlphaFold pipelines should be considered but can be mitigated with stricter control of statistical analysis. Our results indicate selected mutants influence structural plasticity of CCR3 related to ligand interaction, activation, and G protein coupling, using a framework that could be applicable to a wide range of biochemically relevant protein targets following further validation.

Quantifying Protein-Nucleic Acid Interactions for Engineering Useful CRISPR-Cas9 Genome-Editing Variants

Numerous high-specificity Cas9 variants have been engineered for precision genome editing. These variants typically harbor multiple mutations designed to alter the Cas9-single guide RNA (sgRNA)-DNA complex interactions for reduced off-target cleavage. By dissecting the contributions of individual mutations, we attempt to derive principles for designing high-specificity Cas9 variants. Here, we computationally modeled the specificity harnessing mutations of the widely used Cas9 isolated from Streptococcus pyogenes (SpCas9) and investigated their individual mutational effects. We quantified the mutational effects in terms of energy and contact changes by comparing the wild-type and mutant structures. We found that these mutations disrupt the protein-protein or protein-DNA contacts within the Cas9-sgRNA-DNA complex. We also identified additional impacted amino acid sites via energy changes that constitute the structural microenvironment encompassing the focal mutation, giving insights into how the mutations contribute to the high-specificity phenotype of SpCas9. Our method outlines a strategy to evaluate mutational effects that can facilitate rational design for Cas9 optimization.

Generating Multistate Conformations of P-type ATPases with a Conditional Diffusion Model

Understanding and predicting the diverse conformational states of membrane proteins is essential for elucidating their biological functions. Despite advancements in computational methods, accurately capturing these complex structural changes remains a significant challenge. Here, we introduce a computational approach to generate diverse and biologically relevant conformations of membrane proteins using a conditional diffusion model. Our approach integrates forward and backward diffusion processes, incorporating state classifiers and additional conditioners to control the generation gradient of conformational states. We specifically targeted the P-type ATPases, a critical family of membrane transporters, and constructed a comprehensive data set through a combination of experimental structures and molecular dynamics simulations. Our model, incorporating a graph neural network with specialized membrane constraints, demonstrates exceptional accuracy in generating a wide range of P-type ATPase conformations associated with different functional states. This approach represents a meaningful step forward in the computational generation of membrane protein conformations using AI and holds promise for studying the dynamics of other membrane proteins.

AlphaFold2-Based Characterization of Apo and Holo Protein Structures and Conformational Ensembles Using Randomized Alanine Sequence Scanning Adaptation: Capturing Shared Signature Dynamics and Ligand-Induced Conformational Changes

Proteins often exist in multiple conformational states, influenced by the binding of ligands or substrates. The study of these states, particularly the apo (unbound) and holo (ligand-bound) forms, is crucial for understanding protein function, dynamics, and interactions. In the current study, we use AlphaFold2, which combines randomized alanine sequence masking with shallow multiple sequence alignment subsampling to expand the conformational diversity of the predicted structural ensembles and capture conformational changes between apo and holo protein forms. Using several well-established datasets of structurally diverse apo-holo protein pairs, the proposed approach enables robust predictions of apo and holo structures and conformational ensembles, while also displaying notably similar dynamics distributions. These observations are consistent with the view that the intrinsic dynamics of allosteric proteins are defined by the structural topology of the fold and favor conserved conformational motions driven by soft modes. Our findings provide evidence that AlphaFold2 combined with randomized alanine sequence masking can yield accurate and consistent results in predicting moderate conformational adjustments between apo and holo states, especially for proteins with localized changes upon ligand binding. For large hinge-like domain movements, the proposed approach can predict functional conformations characteristic of both apo and ligand-bound holo ensembles in the absence of ligand information. These results are relevant for using this AlphaFold adaptation for probing conformational selection mechanisms according to which proteins can adopt multiple conformations, including those that are competent for ligand binding. The results of this study indicate that robust modeling of functional protein states may require more accurate characterization of flexible regions in functional conformations and the detection of high-energy conformations. By incorporating a wider variety of protein structures in training datasets, including both apo and holo forms, the model can learn to recognize and predict the structural changes that occur upon ligand binding.

Prediction of conformational states in a coronavirus channel using Alphafold-2 and DeepMSA2: Strengths and limitations

The envelope (E) protein is present in all coronavirus genera. This protein can form pentameric oligomers with ion channel activity which have been proposed as a possible therapeutic target. However, high resolution structures of E channels are limited to those of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), responsible for the recent COVID-19 pandemic. In the present work, we used Alphafold-2 (AF2), in ColabFold without templates, to predict the transmembrane domain (TMD) structure of six E-channels representative of genera alpha-, beta- and gamma-coronaviruses in the Coronaviridae family. High-confidence models were produced in all cases when combining multiple sequence alignments (MSAs) obtained from DeepMSA2. Overall, AF2 predicted at least two possible orientations of the α-helices in E-TMD channels: one where a conserved polar residue (Asn-15 in the SARS sequence) is oriented towards the center of the channel, 'polar-in', and one where this residue is in an interhelical orientation 'polar-inter'. For the SARS models, the comparison with the two experimental models 'closed' (PDB: 7K3G) and 'open' (PDB: 8SUZ) is described, and suggests a ∼60˚ α-helix rotation mechanism involving either the full TMD or only its N-terminal half, to allow the passage of ions. While the results obtained are not identical to the two high resolution models available, they suggest various conformational states with striking similarities to those models. We believe these results can be further optimized by means of MSA subsampling, and guide future high resolution structural studies in these and other viral channels.

Protein representations: Encoding biological information for machine learning in biocatalysis

Enzymes offer a more environmentally friendly and low-impact solution to conventional chemistry, but they often require additional engineering for their application in industrial settings, an endeavour that is challenging and laborious. To address this issue, the power of machine learning can be harnessed to produce predictive models that enable the in silico study and engineering of improved enzymatic properties. Such machine learning models, however, require the conversion of the complex biological information to a numerical input, also called protein representations. These inputs demand special attention to ensure the training of accurate and precise models, and, in this review, we therefore examine the critical step of encoding protein information to numeric representations for use in machine learning. We selected the most important approaches for encoding the three distinct biological protein representations - primary sequence, 3D structure, and dynamics - to explore their requirements for employment and inductive biases. Combined representations of proteins and substrates are also introduced as emergent tools in biocatalysis. We propose the division of fixed representations, a collection of rule-based encoding strategies, and learned representations extracted from the latent spaces of large neural networks. To select the most suitable protein representation, we propose two main factors to consider. The first one is the model setup, which is influenced by the size of the training dataset and the choice of architecture. The second factor is the model objectives such as consideration about the assayed property, the difference between wild-type models and mutant predictors, and requirements for explainability. This review is aimed at serving as a source of information and guidance for properly representing enzymes in future machine learning models for biocatalysis.

Predicting multiple conformations of ligand binding sites in proteins suggests that AlphaFold2 may remember too much

The goal of this paper is predicting the conformational distributions of ligand binding sites using the AlphaFold2 (AF2) protein structure prediction program with stochastic subsampling of the multiple sequence alignment (MSA). We explored the opening of cryptic ligand binding sites in 16 proteins, where the closed and open conformations define the expected extreme points of the conformational variation. Due to the many structures of these proteins in the Protein Data Bank (PDB), we were able to study whether the distribution of X-ray structures affects the distribution of AF2 models. We have found that AF2 generates both a cluster of open and a cluster of closed models for proteins that have comparable numbers of open and closed structures in the PDB and not too many other conformations. This was observed even with default MSA parameters, thus without further subsampling. In contrast, with the exception of a single protein, AF2 did not yield multiple clusters of conformations for proteins that had imbalanced numbers of open and closed structures in the PDB, or had substantial numbers of other structures. Subsampling improved the results only for a single protein, but very shallow MSA led to incorrect structures. The ability of generating both open and closed conformations for six out of the 16 proteins agrees with the success rates of similar studies reported in the literature. However, we showed that this partial success is due to AF2 "remembering" the conformational distributions in the PDB and that the approach fails to predict rarely seen conformations.

Challenge for Deep Learning: Protein Structure Prediction of Ligand-Induced Conformational Changes at Allosteric and Orthosteric Sites

In the realm of biomedical research, understanding the intricate structure of proteins is crucial, as these structures determine how proteins function within our bodies and interact with potential drugs. Traditionally, methods like X-ray crystallography and cryo-electron microscopy have been used to unravel these structures, but they are often challenging, time-consuming and costly. Recently, a breakthrough in computational biology has emerged with the development of deep learning algorithms capable of predicting protein structures based on their amino acid sequences (Jumper, J., et al. Nature 2021, 596, 583. Lane, T. J. Nature Methods 2023, 20, 170. Kryshtafovych, A., et al. Proteins: Structure, Function and Bioinformatics 2021, 89, 1607). This study focuses on predicting the dynamic changes that proteins undergo upon ligand binding, specifically when they bind to allosteric sites, i.e. a pocket different from the active site. Allosteric modulators are particularly important for drug discovery, as they open new avenues for designing drugs that can target proteins more effectively and with fewer side effects (Nussinov, R.; Tsai, C. J. Cell 2013, 153, 293). To study this, we curated a data set of 578 X-ray structures comprised of proteins displaying orthosteric and allosteric binding as well as a general framework to evaluate deep learning-based structure prediction methods. Our findings demonstrate the potential and current limitations of deep learning methods, such as AlphaFold2 (Jumper, J., et al. Nature 2021, 596, 583), NeuralPLexer (Qiao, Z., et al. Nat Mach Intell 2024, 6, 195), and RoseTTAFold All-Atom (Krishna, R., et al. Science 2024, 384, eadl2528) to predict not just static protein structures but also the dynamic conformational changes. Herein we show that predicting the allosteric induce-fit conformation still poses a challenge to deep learning methods as they more accurately predict the orthosteric bound conformation compared to the allosteric induce fit conformation. For AlphaFold2, we observed that conformational diversity, and sampling between the apo and holo state could be increased by modifying the MSA depth, but this did not enhance the ability to generate conformations close to the allosteric induced-fit conformation. To further support advancements in protein structure prediction field, the curated data set and evaluation framework are made publicly available.

Probing Functional Allosteric States and Conformational Ensembles of the Allosteric Protein Kinase States and Mutants: Atomistic Modeling and Comparative Analysis of AlphaFold2, OmegaFold, and AlphaFlow Approaches and Adaptations

This study reports a comprehensive analysis and comparison of several AlphaFold2 adaptations and OmegaFold and AlphaFlow approaches in predicting distinct allosteric states, conformational ensembles, and mutation-induced structural effects for a panel of state-switching allosteric ABL mutants. The results revealed that the proposed AlphaFold2 adaptation with randomized alanine sequence scanning can generate functionally relevant allosteric states and conformational ensembles of the ABL kinase that qualitatively capture a unique pattern of population shifts between the active and inactive states in the allosteric ABL mutants. Consistent with the NMR experiments, the proposed AlphaFold2 adaptation predicted that G269E/M309L/T408Y mutant could induce population changes and sample a significant fraction of the fully inactive I2 form which is a low-populated, high-energy state for the wild-type ABL protein. We also demonstrated that other ABL mutants G269E/M309L/T334I and M309L/L320I/T334I that introduce a single activating T334I mutation can reverse equilibrium and populate exclusively the active ABL form. While the precise quantitative predictions of the relative populations of the active and various hidden inactive states in the ABL mutants remain challenging, our results provide evidence that AlphaFold2 adaptation with randomized alanine sequence scanning can adequately detect a spectrum of the allosteric ABL states and capture the equilibrium redistributions between structurally distinct functional ABL conformations. We further validated the robustness of the proposed AlphaFold2 adaptation for predicting the unique inactive architecture of the BSK8 kinase and structural differences between ligand-unbound apo and ATP-bound forms of BSK8. The results of this comparative study suggested that AlpahFold2, OmegaFold, and AlphaFlow approaches may be driven by structural memorization of existing protein folds and are strongly biased toward predictions of the thermodynamically stable ground states of the protein kinases, highlighting limitations and challenges of AI-based methodologies in detecting alternative functional conformations, accurate characterization of physically significant conformational ensembles, and prediction of mutation-induced allosteric structural changes.

Assessing AF2's ability to predict structural ensembles of proteins

Recent breakthroughs in protein structure prediction have enhanced the precision and speed at which protein configurations can be determined. Additionally, molecular dynamics (MD) simulations serve as a crucial tool for capturing the conformational space of proteins, providing valuable insights into their structural fluctuations. However, the scope of MD simulations is often limited by the accessible timescales and the computational resources available, posing challenges to comprehensively exploring protein behaviors. Recently emerging approaches have focused on expanding the capability of AlphaFold2 (AF2) to predict conformational substates of protein. Here, we benchmark the performance of various workflows that have adapted AF2 for ensemble prediction and compare the obtained structures with ensembles obtained from MD simulations and NMR. We provide an overview of the levels of performance and accessible timescales that can currently be achieved with machine learning (ML) based ensemble generation. Significant minima of the free energy surfaces remain undetected.

MassiveFold: unveiling AlphaFold's hidden potential with optimized and parallelized massive sampling

Massive sampling in AlphaFold enables access to increased structural diversity. In combination with its efficient confidence ranking, this unlocks elevated modeling capabilities for monomeric structures and foremost for protein assemblies. However, the approach struggles with GPU cost and data storage. Here we introduce MassiveFold, an optimized and customizable version of AlphaFold that runs predictions in parallel, reducing the computing time from several months to hours. MassiveFold is scalable and able to run on anything from a single computer to a large GPU infrastructure, where it can fully benefit from all the computing nodes.

Develop Targeted Protein Drug Carriers through a High-Throughput Screening Platform and Rational Design

Protein-based drugs offer advantages, such as high specificity, low toxicity, and minimal side effects compared to small molecule drugs. However, delivery of proteins to target tissues or cells remains challenging due to the instability, diverse structures, charges, and molecular weights of proteins. Polymers have emerged as a leading choice for designing effective protein delivery systems, but identifying a suitable polymer for a given protein is complicated by the complexity of both proteins and polymers. To address this challenge, a fluorescence-based high-throughput screening platform called ProMatch to efficiently collect data on protein-polymer interactions, followed by in vivo and in vitro experiments with rational design is developed. Using this approach to streamline polymer selection for targeted protein delivery, candidate polymers from commercially available options are identified and a polyhexamethylene biguanide (PHMB)-based system for delivering proteins to white adipose tissue as a treatment for obesity is developed. A branched polyethylenimine (bPEI)-based system for neuron-specific protein delivery to stimulate optic nerve regeneration is also developed. The high-throughput screening methodology expedites identification of promising polymer candidates for tissue-specific protein delivery systems, thereby providing a platform to develop innovative protein-based therapeutics.

Identification and functional analysis of terpene synthases revealing the secrets of aroma formation in Chrysanthemum aromaticum

C. aromaticum is widely cultivated for its aromatic, medicinal, and tea-applicable properties, earning the nickname 'lavender in composite'. Terpenoids are the major compounds of C. aromaticum fragrance. To reveal the molecular mechanisms of terpenoid biosynthesis in C. aromaticum, NGS and SMRT sequencing were employed to identify the key terpene synthase genes. A total of 59,903 non-redundant transcripts were obtained by the transcriptome analysis. Twenty-nine terpene synthase genes (TPSs) were identified, and phylogenetic analysis showed that they belong to four subfamilies of terpene synthases. Five CaTPSs were successfully cloned. Subcellular localization showed they were present in the nucleus and cytosol. Structure models of five terpene synthases were predicted, and molecular docking results showed good binding affinities with FPP/GPP. In vitro enzymatic tests showed that CaTPS7, CaTPS8, CaTPS10 and CaTPS20 could catalyze substrates to produce terpenoids. CaTPS7 and CaTPS20 were both able to effectively convert the precursor FPP into caryophyllene. CaTPS8 could convert FPP to trans-nerolidol and nerolidyl acetate, while CaTPS10 could convert FPP to elemene and aristolochene. This study lays the groundwork for further research to depict the metabolism network of terpenoid in C. aromaticum. These identical terpene synthase genes could be introduced into the cultivated chrysanthemums to enhance their fragrance.

Click, Compute, Create: A Review of Web-based Tools for Enzyme Engineering

Enzyme engineering, though pivotal across various biotechnological domains, is often plagued by its time-consuming and labor-intensive nature. This review aims to offer an overview of supportive in silico methodologies for this demanding endeavor. Starting from methods to predict protein structures, to classification of their activity and even the discovery of new enzymes we continue with describing tools used to increase thermostability and production yields of selected targets. Subsequently, we discuss computational methods to modulate both, the activity as well as selectivity of enzymes. Last, we present recent approaches based on cutting-edge machine learning methods to redesign enzymes. With exception of the last chapter, there is a strong focus on methods easily accessible via web-interfaces or simple Python-scripts, therefore readily useable for a diverse and broad community.

AlphaFold2's training set powers its predictions of fold-switched conformations

AlphaFold2 (AF2), a deep-learning based model that predicts protein structures from their amino acid sequences, has recently been used to predict multiple protein conformations. In some cases, AF2 has successfully predicted both dominant and alternative conformations of fold-switching proteins, which remodel their secondary and tertiary structures in response to cellular stimuli. Whether AF2 has learned enough protein folding principles to reliably predict alternative conformations outside of its training set is unclear. Here, we address this question by assessing whether CFold-an implementation of the AF2 network trained on a more limited subset of experimentally determined protein structures- predicts alternative conformations of eight fold switchers from six protein families. Previous work suggests that AF2 predicted these alternative conformations by memorizing them during training. Unlike AF2, CFold's training set contains only one of these alternative conformations. Despite sampling 1300-4400 structures/protein with various sequence sampling techniques, CFold predicted only one alternative structure outside of its training set accurately and with high confidence while also generating experimentally inconsistent structures with higher confidence. Though these results indicate that AF2's current success in predicting alternative conformations of fold switchers stems largely from its training data, results from a sequence pruning technique suggest developments that could lead to a more reliable generative model in the future.

Integrative Modeling of Protein-Polypeptide Complexes by Bayesian Model Selection using AlphaFold and NMR Chemical Shift Perturbation Data

Protein-polypeptide interactions, including those involving intrinsically-disordered peptides and intrinsically-disordered regions of protein binding partners, are crucial for many biological functions. However, experimental structure determination of protein-peptide complexes can be challenging. Computational methods, while promising, generally require experimental data for validation and refinement. Here we present CSP_Rank, an integrated modeling approach to determine the structures of protein-peptide complexes. This method combines AlphaFold2 (AF2) enhanced sampling methods with a Bayesian conformational selection process based on experimental Nuclear Magnetic Resonance (NMR) Chemical Shift Perturbation (CSP) data and AF2 confidence metrics. Using a curated dataset of 108 protein-peptide complexes from the Biological Magnetic Resonance Data Bank (BMRB), we observe that while AF2 typically yields models with excellent consistency with experimental CSP data, applying enhanced sampling followed by data-guided conformational selection routinely results in ensembles of structures with improved agreement with NMR observables. For two systems, we cross-validate the CSP-selected models using independently acquired nuclear Overhauser effect (NOE) NMR data and demonstrate how CSP and NMR can be combined using our Bayesian framework for model selection. CSP_Rank is a novel method for integrative modeling of protein-peptide complexes and has broad implications for studies of protein-peptide interactions and aiding in understanding their biological functions.

Predicting Mutation-Induced Allosteric Changes in Structures and Conformational Ensembles of the ABL Kinase Using AlphaFold2 Adaptations with Alanine Sequence Scanning

Despite the success of AlphaFold2 approaches in predicting single protein structures, these methods showed intrinsic limitations in predicting multiple functional conformations of allosteric proteins and have been challenged to accurately capture the effects of single point mutations that induced significant structural changes. We examined several implementations of AlphaFold2 methods to predict conformational ensembles for state-switching mutants of the ABL kinase. The results revealed that a combination of randomized alanine sequence masking with shallow multiple sequence alignment subsampling can significantly expand the conformational diversity of the predicted structural ensembles and capture shifts in populations of the active and inactive ABL states. Consistent with the NMR experiments, the predicted conformational ensembles for M309L/L320I and M309L/H415P ABL mutants that perturb the regulatory spine networks featured the increased population of the fully closed inactive state. The proposed adaptation of AlphaFold can reproduce the experimentally observed mutation-induced redistributions in the relative populations of the active and inactive ABL states and capture the effects of regulatory mutations on allosteric structural rearrangements of the kinase domain. The ensemble-based network analysis complemented AlphaFold predictions by revealing allosteric hotspots that correspond to state-switching mutational sites which may explain the global effect of regulatory mutations on structural changes between the ABL states. This study suggested that attention-based learning of long-range dependencies between sequence positions in homologous folds and deciphering patterns of allosteric interactions may further augment the predictive abilities of AlphaFold methods for modeling of alternative protein sates, conformational ensembles and mutation-induced structural transformations.

Teaching old docks new tricks with machine learning enhanced ensemble docking

We here introduce Ensemble Optimizer (EnOpt), a machine-learning tool to improve the accuracy and interpretability of ensemble virtual screening (VS). Ensemble VS is an established method for predicting protein/small-molecule (ligand) binding. Unlike traditional VS, which focuses on a single protein conformation, ensemble VS better accounts for protein flexibility by predicting binding to multiple protein conformations. Each compound is thus associated with a spectrum of scores (one score per protein conformation) rather than a single score. To effectively rank and prioritize the molecules for further evaluation (including experimental testing), researchers must select which protein conformations to consider and how best to map each compound's spectrum of scores to a single value, decisions that are system-specific. EnOpt uses machine learning to address these challenges. We perform benchmark VS to show that for many systems, EnOpt ranking distinguishes active compounds from inactive or decoy molecules more effectively than traditional ensemble VS methods. To encourage broad adoption, we release EnOpt free of charge under the terms of the MIT license.

Enhancing the catalytic efficiency of M32 carboxypeptidase by semi-rational design and its applications in food taste improvement

BACKGROUND

Carboxypeptidase is an exopeptidase that hydrolyzes amino acids at the C-terminal end of the peptide chain and has a wide range of applications in food. However, in industrial applications, the relatively low catalytic efficiency of carboxypeptidases is one of the main limiting factors for industrialization.

RESULTS

The study has enhanced the catalytic efficiency of Bacillus megaterium M32 carboxypeptidase (BmeCPM32) through semi-rational design. Firstly, the specific activity of the optimal mutant, BmeCPM32-M2, obtained through single-site mutagenesis and combinatorial mutagenesis, was 2.2-fold higher than that of the wild type (187.9 versus 417.8 U mg-1), and the catalytic efficiency was 2.9-fold higher (110.14 versus 325.75 s-1 mmol-1). Secondly, compared to the wild type, BmeCPM32-M2 exhibited a 1.8-fold increase in half-life at 60 °C, with no significant changes in its enzymatic properties (optimal pH, optimal temperature). Finally, BmeCPM32-M2 significantly increased the umami intensity of soy protein isolate hydrolysate by 55% and reduced bitterness by 83%, indicating its potential in developing tasty protein components.

CONCLUSION

Our research has revealed that the strategy based on protein sequence evolution and computational residue mutation energy led to an improved catalytic efficiency of BmeCPM32. Molecular dynamics simulations have revealed that a smaller substrate binding pocket and increased enzyme-substrate affinity are the reasons for the enhanced catalytic efficiency. Furthermore the number of hydrogen bonds and solvent and surface area may contribute to the improvement of thermostability. Finally, the de-bittering effect of BmeCPM32-M2 in soy protein isolate hydrolysate suggests its potential in developing palatable protein components. © 2024 Society of Chemical Industry.

@TOME 3.0: Interfacing Protein Structure Modeling and Ligand Docking

Knowledge of protein-ligand complexes is essential for efficient drug design. Virtual docking can bring important information on putative complexes but it is still far from being simultaneously fast and accurate. Receptors are flexible and adapt to the incoming small molecules while docking is highly sensitive to small conformational deviations. Conformation ensemble is providing a mean to simulate protein flexibility. However, modeling multiple protein structures for many targets is seldom connected to ligand screening in an efficient and straightforward manner. @TOME-3 is an updated version of our former pipeline @TOME-2, in which protein structure modeling is now directly interfaced with flexible ligand docking. Sequence-sequence profile comparisons identify suitable PDB templates for structure modeling and ligands from these templates are used to deduce binding sites to be screened. In addition, bound ligand can be used as pharmacophoric restraint during the virtual docking. The latter is performed by PLANTS while the docking poses are analysed through multiple chemoinformatics functions. This unique combination of tools allows rapid and efficient ligand docking on multiple receptor conformations in parallel. @TOME-3 is freely available on the web at https://atome.cbs.cnrs.fr.

Exploring protein functions from structural flexibility using CABS-flex modeling

Understanding protein function often necessitates characterizing the flexibility of protein structures. However, simulating protein flexibility poses significant challenges due to the complex dynamics of protein systems, requiring extensive computational resources and accurate modeling techniques. In response to these challenges, the CABS-flex method has been developed as an efficient modeling tool that combines coarse-grained simulations with all-atom detail. Available both as a web server and a standalone package, CABS-flex is dedicated to a wide range of users. The web server version offers an accessible interface for straightforward tasks, while the standalone command-line program is designed for advanced users, providing additional features, analytical tools, and support for handling large systems. This paper examines the application of CABS-flex across various structure-function studies, facilitating investigations into the interplay among protein structure, dynamics, and function in diverse research fields. We present an overview of the current status of the CABS-flex methodology, highlighting its recent advancements, practical applications, and forthcoming challenges.

Insect-transmitted plant virus balances its vertical transmission through regulating Rab1-mediated receptor localization

Rice stripe virus (RSV) establishes infection in the ovaries of its vector insect, Laodelphax striatellus. We demonstrate that RSV infection delays ovarian maturation by inhibiting membrane localization of the vitellogenin receptor (VgR), thereby reducing the vitellogenin (Vg) accumulation essential for egg development. We identify the host protein L. striatellus Rab1 protein (LsRab1), which directly interacts with RSV nucleocapsid protein (NP) within nurse cells. LsRab1 is required for VgR surface localization and ovarian Vg accumulation. RSV inhibits LsRab1 function through two mechanisms: NP binding LsRab1 prevents GTP binding, and NP binding LsRab1-GTP complexes stimulates GTP hydrolysis, forming an inactive LsRab1 form. Through this dual inhibition, RSV infection prevents LsRab1 from facilitating VgR trafficking to the cell membrane, leading to inefficient Vg uptake. The Vg-VgR pathway is present in most oviparous animals, and the mechanisms detailed here provide insights into the vertical transmission of other insect-transmitted viruses of medical and agricultural importance.

Predicting protein conformational motions using energetic frustration analysis and AlphaFold2

Proteins perform their biological functions through motion. Although high throughput prediction of the three-dimensional static structures of proteins has proved feasible using deep-learning-based methods, predicting the conformational motions remains a challenge. Purely data-driven machine learning methods encounter difficulty for addressing such motions because available laboratory data on conformational motions are still limited. In this work, we develop a method for generating protein allosteric motions by integrating physical energy landscape information into deep-learning-based methods. We show that local energetic frustration, which represents a quantification of the local features of the energy landscape governing protein allosteric dynamics, can be utilized to empower AlphaFold2 (AF2) to predict protein conformational motions. Starting from ground state static structures, this integrative method generates alternative structures as well as pathways of protein conformational motions, using a progressive enhancement of the energetic frustration features in the input multiple sequence alignment sequences. For a model protein adenylate kinase, we show that the generated conformational motions are consistent with available experimental and molecular dynamics simulation data. Applying the method to another two proteins KaiB and ribose-binding protein, which involve large-amplitude conformational changes, can also successfully generate the alternative conformations. We also show how to extract overall features of the AF2 energy landscape topography, which has been considered by many to be black box. Incorporating physical knowledge into deep-learning-based structure prediction algorithms provides a useful strategy to address the challenges of dynamic structure prediction of allosteric proteins.

AlphaFold predictions of fold-switched conformations are driven by structure memorization

Recent work suggests that AlphaFold (AF)-a deep learning-based model that can accurately infer protein structure from sequence-may discern important features of folded protein energy landscapes, defined by the diversity and frequency of different conformations in the folded state. Here, we test the limits of its predictive power on fold-switching proteins, which assume two structures with regions of distinct secondary and/or tertiary structure. We find that (1) AF is a weak predictor of fold switching and (2) some of its successes result from memorization of training-set structures rather than learned protein energetics. Combining >280,000 models from several implementations of AF2 and AF3, a 35% success rate was achieved for fold switchers likely in AF's training sets. AF2's confidence metrics selected against models consistent with experimentally determined fold-switching structures and failed to discriminate between low and high energy conformations. Further, AF captured only one out of seven experimentally confirmed fold switchers outside of its training sets despite extensive sampling of an additional ~280,000 models. Several observations indicate that AF2 has memorized structural information during training, and AF3 misassigns coevolutionary restraints. These limitations constrain the scope of successful predictions, highlighting the need for physically based methods that readily predict multiple protein conformations.

AlphaFold two years on: Validation and impact

Two years on from the initial release of AlphaFold, we have seen its widespread adoption as a structure prediction tool. Here, we discuss some of the latest work based on AlphaFold, with a particular focus on its use within the structural biology community. This encompasses use cases like speeding up structure determination itself, enabling new computational studies, and building new tools and workflows. We also look at the ongoing validation of AlphaFold, as its predictions continue to be compared against large numbers of experimental structures to further delineate the model's capabilities and limitations.

Integrating Large-Scale Protein Structure Prediction into Human Genetics Research

The last five years have seen impressive progress in deep learning models applied to protein research. Most notably, sequence-based structure predictions have seen transformative gains in the form of AlphaFold2 and related approaches. Millions of missense protein variants in the human population lack annotations, and these computational methods are a valuable means to prioritize variants for further analysis. Here, we review the recent progress in deep learning models applied to the prediction of protein structure and protein variants, with particular emphasis on their implications for human genetics and health. Improved prediction of protein structures facilitates annotations of the impact of variants on protein stability, protein-protein interaction interfaces, and small-molecule binding pockets. Moreover, it contributes to the study of host-pathogen interactions and the characterization of protein function. As genome sequencing in large cohorts becomes increasingly prevalent, we believe that better integration of state-of-the-art protein informatics technologies into human genetics research is of paramount importance.

The power and pitfalls of AlphaFold2 for structure prediction beyond rigid globular proteins

Artificial intelligence-driven advances in protein structure prediction in recent years have raised the question: has the protein structure-prediction problem been solved? Here, with a focus on nonglobular proteins, we highlight the many strengths and potential weaknesses of DeepMind's AlphaFold2 in the context of its biological and therapeutic applications. We summarize the subtleties associated with evaluation of AlphaFold2 model quality and reliability using the predicted local distance difference test (pLDDT) and predicted aligned error (PAE) values. We highlight various classes of proteins that AlphaFold2 can be applied to and the caveats involved. Concrete examples of how AlphaFold2 models can be integrated with experimental data in the form of small-angle X-ray scattering (SAXS), solution NMR, cryo-electron microscopy (cryo-EM) and X-ray diffraction are discussed. Finally, we highlight the need to move beyond structure prediction of rigid, static structural snapshots toward conformational ensembles and alternate biologically relevant states. The overarching theme is that careful consideration is due when using AlphaFold2-generated models to generate testable hypotheses and structural models, rather than treating predicted models as de facto ground truth structures.

Allosteric modulation of fluorescence revealed by hydrogen bond dynamics in a genetically encoded maltose biosensor

Genetically encoded fluorescent biosensors (GEFBs) proved to be reliable tracers for many metabolites and cellular processes. In the simplest case, a fluorescent protein (FP) is genetically fused to a sensing protein which undergoes a conformational change upon ligand binding. This drives a rearrangement in the chromophore environment and changes the spectral properties of the FP. Structural determinants of successful biosensors are revealed only in hindsight when the crystal structures of both ligand-bound and ligand-free forms are available. This makes the development of new biosensors for desired analytes a long trial-and-error process. In the current study, we conducted μs-long all atom molecular dynamics (MD) simulations of a maltose biosensor in both the apo (dark) and holo (bright) forms. We performed detailed hydrogen bond occupancy analyses to shed light on the mechanism of ligand induced conformational change in the sensor protein and its allosteric effect on the chromophore environment. We find that two strong indicators for distinguishing bright and dark states of biosensors are due to substantial changes in hydrogen bond dynamics in the system and solvent accessibility of the chromophore.

Peptidomimetics design and characterization: Bridging experimental and computer-based approaches

Peptidomimetics, designed to mimic peptide biological activity with more drug-like properties, are increasingly pivotal in medicinal chemistry. They offer enhanced systemic delivery, cell penetration, target specificity, and protection against peptidases when compared to their native peptide counterparts. Already utilized in treating diverse diseases like neurodegenerative disorders, cancer and infectious diseases, their future in medicine seems bright, with many peptidomimetics in clinical trials or development stages. Peptidomimetics are well-suited for addressing disturbed protein-protein interactions (PPIs), which often underlie various pathologies. Structural biology and computational methods like molecular dynamics simulations facilitate rational design, whereas machine learning algorithms accelerate protein structure prediction, enabling efficient drug development. Experimental validation via various spectroscopic, biophysical, and biochemical assays confirms computational predictions and guides further optimization. Peptidomimetics, with their tailored constrained structures, represent a frontier in drug design focused on targeting PPIs. In this overview, we present a comprehensive landscape of peptidomimetics, encompassing perspectives on involvement in pathologies, chemical strategies, and methodologies for their characterization, spanning in silico, in vitro and in cell approaches. With increasing interest from pharmaceutical sectors, peptidomimetics hold promise for revolutionizing therapeutic approaches, marking a new era of precision drug discovery.

In silico studies of the open form of human tissue transglutaminase

Human tissue transglutaminase (tTG) is an intriguing multifunctional enzyme involved in various diseases, including celiac disease and neurological disorders. Although a number of tTG inhibitors have been developed, the molecular determinants governing ligand binding remain incomplete due to the lack of high-resolution structural data in the vicinity of its active site. In this study, we obtained the complete high-resolution model of tTG by in silico methods based on available PDB structures. We discovered significant differences in the active site architecture between our and known tTG models, revealing an additional loop which affects the ligand binding affinity. We assembled a library of new potential tTG inhibitors based on the obtained complete model of the enzyme. Our library substantially expands the spectrum of possible drug candidates targeting tTG and encompasses twelve molecular scaffolds, eleven of which are novel and exhibit higher binding affinity then already known ones, according to our in silico studies. The results of this study open new directions for structure-based drug design of tTG inhibitors, offering the complete protein model and suggesting a wide range of new compounds for further experimental validation.

MinActionPath2: path generation between different conformations of large macromolecular assemblies by action minimization

Recent progress in solving macromolecular structures and assemblies by cryogenic electron microscopy techniques enables sampling of their conformations in different states that are relevant to their biological function. Knowing the transition path between these conformations would provide new avenues for drug discovery. While the experimental study of transition paths is intrinsically difficult, in-silico methods can be used to generate an initial guess for those paths. The Elastic Network Model (ENM), along with a coarse-grained representation (CG) of the structures are among the most popular models to explore such possible paths. Here we propose an update to our software platform MinActionPath that generates non-linear transition paths based on ENM and CG models, using action minimization to solve the equations of motion. The new website enables the study of large structures such as ribosomes or entire virus envelopes. It provides direct visualization of the trajectories along with quantitative analyses of their behaviors at http://dynstr.pasteur.fr/servers/minactionpath/minactionpath2_submission.

APACE: AlphaFold2 and advanced computing as a service for accelerated discovery in biophysics

The prediction of protein 3D structure from amino acid sequence is a computational grand challenge in biophysics and plays a key role in robust protein structure prediction algorithms, from drug discovery to genome interpretation. The advent of AI models, such as AlphaFold, is revolutionizing applications that depend on robust protein structure prediction algorithms. To maximize the impact, and ease the usability, of these AI tools we introduce APACE, AlphaFold2 and advanced computing as a service, a computational framework that effectively handles this AI model and its TB-size database to conduct accelerated protein structure prediction analyses in modern supercomputing environments. We deployed APACE in the Delta and Polaris supercomputers and quantified its performance for accurate protein structure predictions using four exemplar proteins: 6AWO, 6OAN, 7MEZ, and 6D6U. Using up to 300 ensembles, distributed across 200 NVIDIA A100 GPUs, we found that APACE is up to two orders of magnitude faster than off-the-self AlphaFold2 implementations, reducing time-to-solution from weeks to minutes. This computational approach may be readily linked with robotics laboratories to automate and accelerate scientific discovery.