An official website of the United States government
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.
more »
« less
Conformational dynamics of adenylate kinase in crystals
Adenylate kinase is a ubiquitous enzyme in living systems and undergoes dramatic conformational changes during its catalytic cycle. For these reasons, it is widely studied by genetic, biochemical, and biophysical methods, both experimental and theoretical. We have determined the basic crystal structures of three differently liganded states of adenylate kinase from Methanotorrus igneus, a hyperthermophilic organism whose adenylate kinase is a homotrimeric oligomer. The multiple copies of each protomer in the asymmetric unit of the crystal provide a unique opportunity to study the variation in the structure and were further analyzed using advanced crystallographic refinement methods and analysis tools to reveal conformational heterogeneity and, thus, implied dynamic behaviors in the catalytic cycle.
more »
« less
- Award ID(s):
- 1231306
- PAR ID:
- 10589096
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- Structural Dynamics
- Volume:
- 11
- Issue:
- 1
- ISSN:
- 2329-7778
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
BY-kinases constitute a protein tyrosine kinase family that encodes unique catalytic domains that deviate from those of eukaryotic kinases resembling P-loop nucleotide triphosphatases (NTPases) instead. We have used computational and supporting biochemical approaches using the catalytic domain of the Escherichia coli BY-kinase, Wzc, to illustrate mechanistic divergences between BY-kinases and NTPases despite their deployment of similar catalytic motifs. In NTPases, the “arginine finger” drives the reactive conformation of ATP while also displacing its solvation shell, thereby making favorable enthalpic and entropic contributions toward βγ-bond cleavage. In BY-kinases, the reactive state of ATP is enabled by ATP·Mg 2+ -induced global conformational transitions coupled to the conformation of the Walker-A lysine. While the BY-kinase arginine finger does promote the desolvation of ATP, it does so indirectly by generating an ordered active site in combination with other structural elements. Bacteria, using these mechanistic variations, have thus repurposed an ancient fold to phosphorylate on tyrosine.more » « less
-
How enzymes achieve their enormous rate enhancements remains a central question in biology, and our understanding to date has impacted drug development, influenced enzyme design, and deepened our appreciation of evolutionary processes. While enzymes position catalytic and reactant groups in active sites, physics requires that atoms undergo constant motion. Numerous proposals have invoked positioning or motions as central for enzyme function, but a scarcity of experimental data has limited our understanding of positioning and motion, their relative importance, and their changes through the enzyme’s reaction cycle. To examine positioning and motions and test catalytic proposals, we collected “room temperature” X-ray crystallography data forPseudomonas putidaketosteroid isomerase (KSI), and we obtained conformational ensembles for this and a homologous KSI from multiple PDB crystal structures. Ensemble analyses indicated limited change through KSI’s reaction cycle. Active site positioning was on the 1- to 1.5-Å scale, and was not exceptional compared to noncatalytic groups. The KSI ensembles provided evidence against catalytic proposals invoking oxyanion hole geometric discrimination between the ground state and transition state or highly precise general base positioning. Instead, increasing or decreasing positioning of KSI’s general base reduced catalysis, suggesting optimized Ångstrom-scale conformational heterogeneity that allows KSI to efficiently catalyze multiple reaction steps. Ensemble analyses of surrounding groups for WT and mutant KSIs provided insights into the forces and interactions that allow and limit active-site motions. Most generally, this ensemble perspective extends traditional structure–function relationships, providing the basis for a new era of “ensemble–function” interrogation of enzymes.more » « less
-
Ever since the first structure of an enzyme, lysozyme, was solved, scientists have been eager to explore how these molecules perform their catalytic function. There has been an overwhelmingly large body of publications that report the X-ray structures of enzymes determined after substrate and ligand binding. None of them truly show the structures of an enzyme working freely through a sequence of events that range from the formation of the enzyme–substrate complex to the dissociation of the product. The technical difficulties were too severe. By 1969, Sluyterman and de Graaf had pointed out that there might be a way to start a reaction in an enzyme crystal by diffusion and following its catalytic cycle in its entirety with crystallographic methods. The crystal only has to be thin enough so that the diffusion is not rate limiting. Of course, the key questions are as follows: How thin should the crystal be? Will the existing X-ray sources be able to collect data from a thin enough crystal fast enough? This review shines light on these questions.more » « less
-
Accurate modeling of conformational energies is key to the crystal structure prediction of conformational polymorphs. Focusing on molecules XXXI and XXXII from the seventh blind test of crystal structure prediction, this study employs various electronic structure methods up to the level of domain-local pair natural orbital coupled cluster singles and doubles with perturbative triples [DLPNO-CCSD(T1)] to benchmark the conformational energies and to assess their impact on the crystal energy landscapes. Molecule XXXI proves to be a relatively straightforward case, with the conformational energies from generalized gradient approximation (GGA) functional B86bPBE-XDM changing only modestly when using more advanced density functionals such as PBE0-D4, ωB97M-V, and revDSD-PBEP86-D4, dispersion-corrected second-order Møller–Plesset perturbation theory (SCS-MP2D), or DLPNO-CCSD(T1). In contrast, the conformational energies of molecule XXXII prove difficult to determine reliably, and variations in the computed conformational energies appreciably impact the crystal energy landscape. Even high-level methods such as revDSD-PBEP86-D4 and SCS-MP2D exhibit significant disagreements with the DLPNO-CCSD(T1) benchmarks for molecule XXXII, highlighting the difficulty of predicting conformational energies for complex, drug-like molecules. The best-converged predicted crystal energy landscape obtained here for molecule XXXII disagrees significantly with what has been inferred about the solid-form landscape experimentally. The identified limitations of the calculations are probably insufficient to account for the discrepancies between theory and experiment on molecule XXXII, and further investigation of the experimental solid-form landscape would be valuable. Finally, assessment of several semi-empirical methods findsr2SCAN-3c to be the most promising, with conformational energy accuracy intermediate between the GGA and hybrid functionals and a low computational cost.more » « less
