Search NASA⌕ Search

SEARCH · Search NASA

Results for “conformational equilibrium”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Dual-action kinase inhibitors influence p38α MAP kinase dephosphorylation

Reversible protein phosphorylation directs essential cellular processes including cell division, cell growth, cell death, inflammation, and differentiation. Because protein phosphorylation drives diverse diseases, kinases and phosphatases have been targets for drug discovery, with some achieving remarkable clinical success. Most protein kinases are activated by phosphorylation of their activation loops, which shifts the conformational equilibrium of the kinase toward the active state. To turn off the kinase, protein phosphatases dephosphorylate these sites, but how the conformation of the dynamic activation loop contributes to dephosphorylation was not known. To answer this, we modulated the activation loop conformational equilibrium of human p38α ΜΑP kinase with existing kinase inhibitors that bind and stabilize specific inactive activation loop conformations. From this, we identified three inhibitors that increase the rate of dephosphorylation of the activation loop phospho-threonine by the PPM serine/threonine phosphatase WIP1. Hence, these compounds are “dual-action” inhibitors that simultaneously block the active site and promote p38α dephosphorylation. Our X-ray crystal structures of phosphorylated p38α bound to the dual-action inhibitors reveal a shared flipped conformation of the activation loop with a fully accessible phospho-threonine. In contrast, our X-ray crystal structure of phosphorylated apo human p38α reveals a different activation loop conformation with an inaccessible phospho-threonine, thereby explaining the increased rate of dephosphorylation upon inhibitor binding. These findings reveal a conformational preference of phosphatases for their targets and suggest a unique approach to achieving improved potency and specificity for therapeutic kinase inhibitors.

Science & Technology - Other Topics↗

An autoinhibited state of 53BP1 revealed by small molecule antagonists and protein engineering

The recruitment of 53BP1 to chromatin, mediated by its recognition of histone H4 dimethylated at lysine 20 (H4K20me2), is important for DNA double-strand break repair. Using a series of small molecule antagonists, we demonstrate a conformational equilibrium between an open and a pre-existing lowly populated closed state of 53BP1 in which the H4K20me2 binding surface is buried at the interface between two interacting 53BP1 molecules. In cells, these antagonists inhibit the chromatin recruitment of wild type 53BP1, but do not affect 53BP1 variants unable to access the closed conformation despite preservation of the H4K20me2 binding site. Thus, this inhibition operates by shifting the conformational equilibrium toward the closed state. Our work therefore identifies an auto-associated form of 53BP1—autoinhibited for chromatin binding—that can be stabilized by small molecule ligands encapsulated between two 53BP1 protomers. Such ligands are valuable research tools to study the function of 53BP1 and have the potential to facilitate the development of new drugs for cancer therapy.

59 BASIC BIOLOGICAL SCIENCES↗

Conformational influence of the ribose 2'-hydroxyl group: crystal structures of DNA-RNA chimeric duplexes

We have crystallized three double-helical DNA-RNA chimeric duplexes and determined their structures by X-ray crystallography at resolutions between 2 and 2.25 A. The two self-complementary duplexes [r(G)d(CGTATACGC)]2 and [d(GCGT)r(A)d(TACGC)]2, as well as the Okazaki fragment d(GGGTATACGC).r(GCG)d(TATACCC), were found to adopt A-type conformations. The crystal structures are non-isomorphous, and the crystallographic environments for the three chimeras are different. A number of intramolecular interactions of the ribose 2'-hydroxyl groups contribute to the stabilization of the A-conformation. Hydrogen bonds between 2'-hydroxyls and 5'-oxygens or phosphate oxygens, in addition to the previously observed hydrogen bonds to 1'-oxygens of adjacent riboses and deoxyriboses, are observed in the DNA-RNA chimeric duplexes. The crystalline chimeric duplexes do not show a transition between the DNA A- and B-conformations. CD spectra suggest that the Okazaki fragment assumes an A-conformation in solution as well. In this molecule the three RNA residues may therefore lock the complete decamer in the A-conformation. Crystals of an all-DNA strand with the same sequence as the self-complementary chimeras show a morphology which is different from those of the chimera crystals. Moreover, the oligonucleotide does not match any of the sequence characteristics of DNAs usually adopting the A-conformation in the crystalline state (e.g., octamers with short alternating stretches of purines and pyrimidines). In DNA-RNA chimeric duplexes, it is therefore possible that a single RNA residue can drive the conformational equilibrium toward the A-conformation.

NASA Discipline Number 52-20↗

Revealing the mechanism of action of a first-in-class covalent inhibitor of $\mathrm{KRASG12C}$ $\mathrm{(ON)}$ and other functional properties of oncogenic $\text{KRAS}$ by 31 P $\mathrm{NMR}$

Individual oncogenic KRAS mutants confer distinct differences in biochemical properties and signaling for reasons that are not well understood. KRAS activity is closely coupled to protein dynamics and is regulated through two interconverting conformations: state 1 (inactive, effector binding deficient) and state 2 (active, effector binding enabled). Here, we use 31 P NMR to delineate the differences in state 1 and state 2 populations present in WT and common KRAS oncogenic mutants (G12C, G12D, G12V, G13D, and Q61L) bound to its natural substrate GTP or a commonly used nonhydrolyzable analog GppNHp (guanosine-5'-[(β,γ)-imido] triphosphate). Our results show that GppNHp-bound proteins exhibit significant state 1 population, whereas GTP-bound KRAS is primarily (90% or more) in state 2 conformation. This observation suggests that the predominance of state 1 shown here and in other studies is related to GppNHp and is most likely nonexistent in cells. We characterize the impact of this differential conformational equilibrium of oncogenic KRAS on RAF1 kinase effector RAS-binding domain and intrinsic hydrolysis. Through a KRAS G12C drug discovery, we have identified a novel small-molecule inhibitor, BBO-8956, which is effective against both GDP- and GTP-bound KRAS G12C. We show that binding of this inhibitor significantly perturbs state 1–state 2 equilibrium and induces an inactive state 1 conformation in GTP-bound KRAS G12C. In the presence of BBO-8956, RAF1–RAS-binding domain is unable to induce a signaling competent state 2 conformation within the ternary complex, demonstrating the mechanism of action for this novel and active-conformation inhibitor.

31P NMR↗

A biophysical framework for double-drugging kinases

Selective orthosteric inhibition of kinases has been challenging due to the conserved active site architecture of kinases and emergence of resistance mutants. Simultaneous inhibition of distant orthosteric and allosteric sites, which we refer to as “double-drugging”, has recently been shown to be effective in overcoming drug resistance. However, detailed biophysical characterization of the cooperative nature between orthosteric and allosteric modulators has not been undertaken. Here, we provide a quantitative framework for double-drugging of kinases employing isothermal titration calorimetry, Förster resonance energy transfer, coupled-enzyme assays, and X-ray crystallography. We discern positive and negative cooperativity for Aurora A kinase (AurA) and Abelson kinase (Abl) with different combinations of orthosteric and allosteric modulators. We find that a conformational equilibrium shift is the main principle governing cooperativity. Notably, for both kinases, we find a synergistic decrease of the required orthosteric and allosteric drug dosages when used in combination to inhibit kinase activities to clinically relevant inhibition levels. X-ray crystal structures of the double-drugged kinase complexes reveal the molecular principles underlying the cooperative nature of double-drugging AurA and Abl with orthosteric and allosteric inhibitors. Finally, we observe a fully closed conformation of Abl when bound to a pair of positively cooperative orthosteric and allosteric modulators, shedding light on the puzzling abnormality of previously solved closed Abl structures. Collectively, our data provide mechanistic and structural insights into rational design and evaluation of double-drugging strategies.

59 BASIC BIOLOGICAL SCIENCES↗

Interactions of anesthetics with the membrane-water interface

Although the potency of conventional anesthetics correlates with lipophilicity, an affinity to water also is essential. It was recently found that compounds with very low affinities to water do not produce anesthesia regardless of their lipophilicity. This finding implies that clinical anesthesia might arise because of interactions at molecular sites near the interface of neuronal membranes with the aqueous environment and, therefore, might require increased concentrations of anesthetic molecules at membrane interfaces. As an initial test of this hypothesis, we calculated in molecular dynamics simulations the free energy profiles for the transfer of anesthetic 1,1,2-trifluoroethane and nonanesthetic perfluoroethane across water-membrane and water-hexane interfaces. Consistent with the hypothesis, it was found that trifluoroethane, but not perfluoroethane, exhibits a free energy minimum and, therefore, increased concentrations at both interfaces. The transfer of trifluoroethane from water to the nonpolar hexane or interior of the membrane is accompanied by a considerable, solvent-induced shift in the conformational equilibrium around the C-C bond.

NASA Center ARC↗

Data-driven equation-free dynamics applied to many-protein complexes: The microtubule tip relaxation

Microtubules (MTs) constitute the largest components of the eukaryotic cytoskeleton and play crucial roles in various cellular processes, including mitosis and intracellular transport. The property allowing MTs to cater to such diverse roles is attributed to dynamic instability, which is coupled to the hydrolysis of GTP (guanosine-5'-triphosphate) to GDP (guanosine-5'-diphosphate) within the β-tubulin monomers. Understanding the equilibrium dynamics and the structural features of both GDP- and GTP-complexed MT tips, especially at an all-atom level, remains challenging for both experimental and computational methods because of their dynamic nature and the prohibitive computational demands of simulating large, many-protein systems. This study employs the “equation-free” multiscale computational method to accelerate the relaxation of all-atom simulations of MT tips toward their putative equilibrium conformation. Using large MT lattice systems (14 protofilaments × 8 heterodimers) comprising ~21-38 million atoms, we applied this multiscale approach to leapfrog through time and nearly double the computational efficiency in realizing relaxed all-atom conformations of GDP- and GTP-complexed MT tips. Commencing from an initial 4 μs unbiased all-atom simulation, we interleave coarse projective “equation-free” jumps with short bursts of all-atom molecular dynamics simulation to realize an additional effective simulation time of 1.875 μs. Our 5.875 μs of effective simulation trajectories for each system expose the subtle yet essential differences in the structures of MT tips as a function of whether β-tubulin monomer is complexed with GDP or GTP, as well as the lateral interactions within the MT tip, offering a refined understanding of features underlying MT dynamic instability. Furthermore, the approach presents a robust and generalizable framework for future explorations of large biomolecular systems at atomic resolution.

Wu, Jiangbo [University of Chicago, IL (United Sta↗

Discovery of Multiple Light-Harvesting States of the Photosynthetic Protein PE545

Cryptophytes are photosynthetic microalga that flourish in a remarkable diversity of natural environments by using pigment-containing proteins with absorption maxima tuned to each ecological niche. While this diversity in the absorption has been well established, the subsequent photophysics is highly sensitive to the local protein environment and so may exhibit similar variation. Thermal fluctuations of the protein conformation are expected to introduce photophysical heterogeneity of the pigments that may have evolved important functional properties in a manner similar to that of the absorption. However, such heterogeneity is averaged out in ensemble measurements and, therefore, has not yet been probed. Here, we report single-molecule measurements of phycoerythrin 545 (PE545), the prototypical cryptophyte antenna protein, in its native dimeric form. A conformational ensemble was resolved consisting of distinct photophysical states with different light-harvesting properties. Proteins that did not quench, partially quenched, or fully quenched absorbed light were observed. Light intensity increased the quenched-state population of the dimer, potentially as a mechanism to deal with the extreme light intensities found in aqueous environments. Cross-linking, which mimics local interactions, introduces this light-dependent functionality while also suppressing other conformational dynamics. The cellular organization can, therefore, actively modulate the protein conformation and dynamics, selecting for distinct levels of light harvesting. Furthermore, the complex conformational equilibrium provides an additional mechanism for cryptophytes and likely other photosynthetic organisms to optimize solar energy capture and conversion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reduced dynamic complexity allows structure elucidation of an excited state of KRAS G13D

Localized dynamics of RAS, including regions distal to the nucleotide-binding site, is of high interest for elucidating the mechanisms by which RAS proteins interact with effectors and regulators and for designing inhibitors. Among several oncogenic mutants, methyl relaxation dispersion experiments reveal highly synchronized conformational dynamics in the active (GMPPNP-bound) KRAS G13D , which suggests an exchange between two conformational states in solution. Methyl and 31 P NMR spectra of active KRAS G13D in solution confirm a two-state ensemble interconverting on the millisecond timescale, with a major Pγ atom peak corresponding to the dominant State 1 conformation and a secondary peak indicating an intermediate state different from the known State 2 conformation recognized by RAS effectors. High-resolution crystal structures of active KRAS G13D and KRAS G13D -RAF1 RBD complex provide snapshots of the State 1 and 2 conformations, respectively. We use residual dipolar couplings to solve and cross-validate the structure of the intermediate state of active KRAS G13D , showing a conformation distinct from those of States 1 and 2 outside the known flexible switch regions. The dynamic coupling between the conformational exchange in the effector lobe and the breathing motion in the allosteric lobe is further validated by a secondary mutation in the allosteric lobe, which affects the conformational population equilibrium.

59 BASIC BIOLOGICAL SCIENCES↗

Towards modelling AR Sco: generalized particle dynamics and strong radiation-reaction regimes

Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is especially highlighted in pulsar modelling due to the extreme associated electromagnetic fields and particle Lorentz factors. Including the radiation-reaction force in the particle dynamics adds even more complexity to the problem, but is crucial for such extreme astrophysical sources. We have also realized the need for such modelling concerning magnetic mirroring and particle injection models proposed for AR Sco, the first white dwarf pulsar. This paper demonstrates the benefits of using higher-order explicit numerical integrators with adaptive time-step methods to solve the full particle dynamics with radiation-reaction forces included. We show that for standard test scenarios, namely various combinations of uniform E- and B-fields and a static dipole B-field, the schemes we use are equivalent to and in extreme field cases outperform standard symplectic integrators in accuracy. We show that the higher-order schemes have massive computational time improvements due to the adaptive time-steps we implement, especially in non-uniform field scenarios and included radiation reaction where the particle gyro-radius rapidly changes. When balancing accuracy and computational time, we identified the adaptive Dormand–Prince eighth-order scheme to be ideal for our use cases. The schemes we use maintain accuracy and stability in describing the particle dynamics and we indicate how a charged particle enters radiation-reaction equilibrium and conforms to the analytical Aristotelian Electrodynamics expectations.

79 ASTRONOMY AND ASTROPHYSICS↗

Computer modeling of polymers

A Polymer Molecular Analysis Display System (p-MADS) was developed for computer modeling of polymers. This method of modeling allows for the theoretical calculation of molecular properties such as equilibrium geometries, conformational energies, heats of formations, crystal packing arrangements, and other properties. Furthermore, p-MADS has the following capabilities: constructing molecules from internal coordinates (bonds length, angles, and dihedral angles), Cartesian coordinates (such as X-ray structures), or from stick drawings; manipulating molecules using graphics and making hard copy representation of the molecules on a graphics printer; and performing geometry optimization calculations on molecules using the methods of molecular mechanics or molecular orbital theory.

Green, Terry J.↗

Boronate affinity adsorption of RNA: possible role of conformational changes

Batch equilibrium adsorption isotherm determination is used to characterize the adsorption of mixed yeast RNA on agarose-immobilized m-aminophenylboronic acid. It is shown that the affinity-enhancing influence of divalent cations depends strongly on the precise nature of the cation used, with barium being far more effective than the conventionally-used magnesium. This adsorption-promoting influence of barium is suggested to arise primarily from ionic influences on the structure and rigidity of the RNA molecule, as the adsorption of ribose-based small molecules is not similarly affected. The substitution of barium for the standard magnesium counterion does not greatly promote the adsorption of DNA, implying that the effect is specific to RNA and may be useful in boronate-based RNA separations. RNA adsorption isotherms exhibit sharp transitions as functions of temperature, and these transitions occur at different temperatures with Mg2+ and with Ba2+. Adsorption affinity and capacity were found to increase markedly at lower temperatures, suggestive of an enthalpically favored interaction process. The stoichiometric displacement parameter, Z, in Ba2+ buffer is three times the value in Mg2+ buffer, and is close to unity.

NASA Discipline Environmental Health↗

Independent Mutation of Two Bridging Carboxylate Ligands Stabilizes Alternate Conformers of the Photosynthetic O 2 -Evolving Mn 4 CaO 5 Cluster in Photosystem II

The O 2 -evolving Mn 4 CaO 5 cluster in Photosystem II is ligated by six carboxylate residues. One of these is D170 of the D1 subunit. This carboxylate bridges between one Mn ion (Mn4) and the Ca ion. A second carboxylate ligand is D342 of the D1 subunit. This carboxylate bridges between two Mn ions (Mn1 and Mn2). D170 and D342 are located on opposite sides of the Mn 4 CaO 5 cluster. Recently, it was shown that the D170E mutation perturbs both the intricate networks of H-bonds that surround the Mn 4 CaO 5 cluster and the equilibrium between different conformers of the cluster in two of its lower oxidation states, S 1 and S 2 , while still supporting O 2 evolution at approx. 50% the rate of wild-type. In this study, we show that the D342E mutation produces much the same alterations to the cluster’s FTIR and EPR spectra as D170E, while still supporting O 2 evolution at approx. 20% the rate of wild-type. Furthermore, the double mutation, D170E + D342E, behaves similarly to the two single mutations. We conclude that D342E alters the equilibrium between different conformers of the cluster in its S 1 and S 2 states in the same manner as D170E and perturbs the H-bond networks in a similar fashion. This is the second identification of a Mn 4 CaO 5 metal ligand whose mutation influences the equilibrium between the different conformers of the S 1 and S 2 states without eliminating O 2 evolution. Finally, this finding has implications for our understanding of the mechanism of O 2 formation in terms of catalytically active/inactive conformations of the Mn 4 CaO 5 cluster in its lower oxidation states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prehydrodynamic evolution in large and small systems

We extend our previous investigation of the effects of prehydrodynamic evolution on final-state observables in heavy-ion collisions [38] to smaller systems. We use a state-of-the-art hybrid model for the numerical simulations with optimal parameters obtained from a previous Bayesian study. By studying p–Pb collisions, we find that the effects due to the assumption of a conformal evolution in the prehydrodynamical stage are even more important in small systems. We also show that this effect depends on the time duration of the pre-equilibrium stage, which is further enhanced in small systems. Finally, we show that the recent proposal of a free-streaming with subluminal velocity for the pre-equilibrium stage, thus effectively breaking conformal invariance, can alleviate the contamination of final-state observables. Furthermore, our study further reinforces the need for moving beyond conformal approaches in pre-equilibrium dynamics modeling, especially when extracting transport coefficients from hybrid models in the high-precision era of heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Visualizing the Active Site Oxyanion Loop Transition Upon Ensitrelvir Binding and Transient Dimerization of SARS-CoV-2 Main Protease

N-terminal autoprocessing from its polyprotein precursor enables creating the mature-like stable dimer interface of SARS-CoV-2 main protease (MPro), concomitant with the active site oxyanion loop equilibrium transitioning to the active conformation (E*) and onset of catalytic activity. Here, through mutagenesis of critical interface residues and evaluating noncovalent inhibitor (ensitrelvir, ESV) facilitated dimerization through its binding to MPro, we demonstrate that residues extending from Ser1 through Glu14 are critical for dimerization. Combined mutations G11A, E290A and R298A (MPro™) restrict dimerization even upon binding of ESV to monomeric MPro™ with an inhibitor dissociation constant of 7.4 ± 1.6 µM. Contrasting the covalent inhibitor NMV or GC373 binding to monomeric MPro, ESV binding enabled capturing the transition of the oxyanion loop conformations in the absence of a reactive warhead and independent of dimerization. Characterization of complexes by room-temperature X-ray crystallography reveals ESV bound to the E* state of monomeric MPro as well as an intermediate approaching the inactive state (E). It appears that the E* to E equilibrium shift occurs initially from G138-F140 residues, leading to the unwinding of the loop and formation of the 3 10 -helix. Finally, we describe a transient dimer structure of the MPro precursor held together through interactions of residues A5-G11 with distinct states of the active sites, E and E*, likely representing an intermediate in the autoprocessing pathway.

59 BASIC BIOLOGICAL SCIENCES↗

Genetic switching by the Lac repressor is based on two-state Monod–Wyman–Changeux allostery

High-resolution NMR spectroscopy enabled us to characterize allosteric transitions between various functional states of the dimeric Escherichia coli Lac repressor. In the absence of ligands, the dimer exists in a dynamic equilibrium between DNA-bound and inducer-bound conformations. Binding of either effector shifts this equilibrium toward either bound state. Analysis of the ternary complex between repressor, operator DNA, and inducer shows how adding the inducer results in allosteric changes that disrupt the interdomain contacts between the inducer binding and DNA binding domains and how this in turn leads to destabilization of the hinge helices and release of the Lac repressor from the operator. Based on our data, the allosteric mechanism of the induction process is in full agreement with the well-known Monod–Wyman–Changeux model.

59 BASIC BIOLOGICAL SCIENCES↗

Ultrafast Correlation Energy Estimator

A virtually no-cost method is proposed that can compute the correlation energies of general, covalently bonded, organic, and inorganic molecules (including conjugated π-electron systems) with a well-defined dominant Lewis structure at the accuracy of 99.5% of the near-exact values determined by the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] in the complete-basis-set (CBS) limit. This Correlation Energy Per Bond (CEPB) method assigns a partial correlation energy to each bond type (characterized by the identities of the two atoms forming the bond and its integer bond order) and to a lone pair, regardless of the bond length, bond angle, sp-hybridization, π-electron conjugation, ionicity, noncovalent interactions, etc. At its current stage, the method is mainly suitable for near-equilibrium geometries. The correlation energies per bond are determined by a fit to the CCSD(T)/CBS benchmarks. It can neither improve the equilibrium structures nor discern conformers or positional isomers, yet its accuracy for reaction energies rivals that of the second-order Møller–Plesset perturbation theory, which is far more expensive. Its promising performance underscores the possibility that surprisingly compact, chemically intuitive molecular fragments exist into which correlation energies can be partitioned, leading to various ultrafast correlation-energy estimators tailored to different purposes.

Binding energy↗