Importance of hydrophobic cavities in allosteric regulation of formylglycinamide synthetase: insight from xenon trapping and statistical coupling analysis

Studying the Effects of Dissolved Noble Gases and High Hydrostatic Pressure on the Spherical DOPC Bilayer Using Molecular Dynamic Simulations

Fine-grained molecular dynamics simulations have been conducted to depict lipid objects enclosed in water and interacting with a series of noble gases dissolved in the medium. The simple point-charge (SPC) water system, featuring a boundary composed of 1,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPC) molecules, maintained stability throughout the simulation under standard conditions. This allowed for the accurate modeling of the effects of hydrostatic pressure at an ambient pressure of 25 bar. The chosen pressure references the 240 m depth of seawater: the horizon frequently used by commercial divers, who comprise the primary patient population of the neurological complication of inert gas narcosis and the consequences of high-pressure neurological syndrome. To quantify and validate the neurological effects of noble gases and discriminate them from high hydrostatic pressure, we reduced the dissolved gas molar concentration to 1.5%, three times smaller than what we previously tested for the planar bilayer (3.5%). The nucleation and growth of xenon, argon and neon nanobubbles proved consistent with the data from the planar bilayer simulations. On the other hand, hyperbaric helium induces only a residual distorting effect on the liposome, with no significant condensed gas fraction observed within the hydrophobic core. The bubbles were distributed over a large volume-both in the bulk solvent and in the lipid phase-thereby causing substantial membrane distortion. This finding serves as evidence of the validity of the multisite distortion hypothesis for the neurological effect of inert gases at high pressure.

Deciphering protein microenvironment by using a cysteine specific switch-ON fluorescent probe

Fluorescent probes provide an unparalleled opportunity to visualize and quantify dynamic events. Here, we employ a medium-size, cysteine specific coumarin based switch-ON fluorescent probe 'L' to track protein unfolding profiles and accessibility of cysteine residues in proteins. It was established that 'L' is highly selective and exhibits no artifact due to interaction with other bystander species. 'L' is able to gauge subtle changes in protein microenvironment and proved to be effective in delineating early unfolding events that are difficult to otherwise discern by classic techniques such as circular dichroism. By solving the X-ray structure of TadA and probing the temperature dependent fluorescence-ON response with native TadA and its cysteine mutants, it was revealed that unfolding occurs in a stage-wise manner and the regions that are functionally important form compact sub-domains and unfold at later stages. Our results assert that probe 'L' serves as an efficient tool to monitor subtle changes in protein structure and can be employed as a generic dye to study processes such as protein unfolding.

An Ultra-Microporous Metal-Organic Framework with Exceptional Xe Capacity

Molecular confinement plays a significant effect on trapped gas and solvent molecules. A fundamental understanding of gas adsorption within the porous confinement provides information necessary to design a material with improved selectivity. In this regard, metal-organic framework (MOF) adsorbents are ideal candidate materials to study confinement effects for weakly interacting gas molecules, such as noble gases. Among the noble gases, xenon (Xe) has practical applications in the medical, automotive and aerospace industries. In this Communication, we report an ultra-microporous nickel-isonicotinate MOF with exceptional Xe uptake and selectivity compared to all benchmark MOF and porous organic cage materials. The selectivity arises because of the near perfect fit of the atomic Xe inside the porous confinement. Notably, at low partial pressure, the Ni-MOF interacts very strongly with Xe compared to the closely related Krypton gas (Kr) and more polarizable CO₂ . Further ¹²⁹ Xe NMR suggests a broad isotropic chemical shift due to the reduced motion as a result of confinement.

Small design from big alignment: engineering proteins with multiple sequence alignment as the starting point

Multiple sequence alignment (MSA) is a fundamental way to gain information that cannot be obtained from the analysis of any individual sequence included in the alignment. It provides ways to investigate the relationship between sequence and function from a perspective of evolution. Thus, the MSA of proteins can be employed as a reference for protein engineering. In this paper, we reviewed the recent advances to highlight how protein engineering was benefited from the MSA of proteins. These methods include (1) engineering the thermostability or solubility of proteins by making it closer to the consensus sequence of the alignment through introducing site mutations; (2) structure-based engineering proteins with comparative modeling; (3) creating paleoenzymes featured with high thermostability and promiscuity by constructing the ancestral sequences derived from multiple sequence alignment; and (4) incorporating site-mutations targeting the evolutionarily coupled sites identified from multiple sequence alignment.

Role of allosteric switches and adaptor domains in long-distance cross-talk and transient tunnel formation

Transient tunnels that assemble and disassemble to facilitate passage of unstable intermediates in enzymes containing multiple reaction centers are controlled by allosteric cues. Using the 140-kDa purine biosynthetic enzyme PurL as a model system and a combination of biochemical and x-ray crystallographic studies, we show that long-distance communication between ~25-Å distal active sites is initiated by an allosteric switch, residing in a conserved catalytic loop, adjacent to the synthetase active site. Further, combinatory experiments seeded from molecular dynamics simulations help to delineate transient states that bring out the central role of nonfunctional adaptor domains. We show that carefully orchestrated conformational changes, facilitated by interplay of dynamic interactions at the allosteric switch and adaptor-domain interface, control reactivity and concomitant formation of the ammonia tunnel. This study asserts that substrate channeling is modulated by allosteric hotspots that alter protein energy landscape, thereby allowing the protein to adopt transient conformations paramount to function.

Deep Analysis of Residue Constraints (DARC): identifying determinants of protein functional specificity

Protein functional constraints are manifest as superfamily and functional-subgroup conserved residues, and as pairwise correlations. Deep Analysis of Residue Constraints (DARC) aids the visualization of these constraints, characterizes how they correlate with each other and with structure, and estimates statistical significance. This can identify determinants of protein functional specificity, as we illustrate for bacterial DNA clamp loader ATPases. These load ring-shaped sliding clamps onto DNA to keep polymerase attached during replication and contain one δ, three γ, and one δ’ AAA+ subunits semi-circularly arranged in the order δ-γ₁-γ₂-γ₃-δ’. Only γ is active, though both γ and δ’ functionally influence an adjacent γ subunit. DARC identifies, as functionally-congruent features linking allosterically the ATP, DNA, and clamp binding sites: residues distinctive of γ and of γ/δ’ that mutually interact in trans, centered on the catalytic base; several γ/δ’-residues and six γ/δ’-covariant residue pairs within the DNA binding N-termini of helices α2 and α3; and γ/δ’-residues associated with the α2 C-terminus and the clamp-binding loop. Most notable is a trans-acting γ/δ’ hydroxyl group that 99% of other AAA+ proteins lack. Mutation of this hydroxyl to a methyl group impedes clamp binding and opening, DNA binding, and ATP hydrolysis—implying a remarkably clamp-loader-specific function.

Structural Insights of Stereospecific Inhibition of Human Acetylcholinesterase by VX and Subsequent Reactivation by HI-6

Over 50 years ago, the toxicity of irreversible organophosphate inhibitors targeting human acetylcholinesterase (hAChE) was observed to be stereospecific. The therapeutic reversal of hAChE inhibition by reactivators has also been shown to depend on the stereochemistry of the inhibitor. To gain clarity on the mechanism of stereospecific inhibition, the X-ray crystallographic structures of hAChE inhibited by a racemic mixture of VX (P _R/S) and its enantiomers were obtained. Beyond identifying hAChE structural features that lend themselves to stereospecific inhibition, structures of the reactivator HI-6 bound to hAChE inhibited by VX enantiomers of varying toxicity, or in its uninhibited state, were obtained. Comparison of hAChE in these pre-reactivation and post-reactivation states along with enzymatic data reveals the potential influence of unproductive reactivator poses on the efficacy of these types of therapeutics. The recognition of structural features related to hAChE's stereospecificity toward VX shed light on the molecular influences of toxicity and their effect on reactivators. In addition to providing a better understanding of the innate issues with current reactivators, an avenue for improvement of reactivators is envisioned.

Inferring joint sequence-structural determinants of protein functional specificity

Residues responsible for allostery, cooperativity, and other subtle but functionally important interactions remain difficult to detect. To aid such detection, we employ statistical inference based on the assumption that residues distinguishing a protein subgroup from evolutionarily divergent subgroups often constitute an interacting functional network. We identify such networks with the aid of two measures of statistical significance. One measure aids identification of divergent subgroups based on distinguishing residue patterns. For each subgroup, a second measure identifies structural interactions involving pattern residues. Such interactions are derived either from atomic coordinates or from Direct Coupling Analysis scores, used as surrogates for structural distances. Applying this approach to N-acetyltransferases, P-loop GTPases, RNA helicases, synaptojanin-superfamily phosphatases and nucleases, and thymine/uracil DNA glycosylases yielded results congruent with biochemical understanding of these proteins, and also revealed striking sequence-structural features overlooked by other methods. These and similar analyses can aid the design of drugs targeting allosteric sites.

V67L Mutation Fills an Internal Cavity To Stabilize RecA Mtu Intein

Inteins mediate protein splicing, which has found extensive applications in protein science and biotechnology. In the Mycobacterium tuberculosis RecA mini-mini intein (ΔΔIhh), a single valine to leucine substitution at position 67 (V67L) dramatically increases intein stability and activity. However, crystal structures show that the V67L mutation causes minimal structural rearrangements, with a root-mean-square deviation of 0.2 Å between ΔΔIhh-V67 and ΔΔIhh-L67. Thus, the structural mechanisms for V67L stabilization and activation remain poorly understood. In this study, we used intrinsic tryptophan fluorescence, high-pressure nuclear magnetic resonance (NMR), and molecular dynamics (MD) simulations to probe the structural basis of V67L stabilization of the intein fold. Guanidine hydrochloride denaturation monitored by fluorescence yielded free energy changes (ΔG_f°) of -4.4 and -6.9 kcal mol^-1 for ΔΔIhh-V67 and ΔΔIhh-L67, respectively. High-pressure NMR showed that ΔΔIhh-L67 is more resistant to pressure-induced unfolding than ΔΔIhh-V67 is. The change in the volume of folding (ΔV_f) was significantly larger for V67 (71 ± 2 mL mol^-1) than for L67 (58 ± 3 mL mol^-1) inteins. The measured difference in ΔV_f (13 ± 3 mL mol^-1) roughly corresponds to the volume of the additional methylene group for Leu, supporting the notion that the V67L mutation fills a nearby cavity to enhance intein stability. In addition, we performed MD simulations to show that V67L decreases side chain dynamics and conformational entropy at the active site. It is plausible that changes in cavities in V67L can also mediate allosteric effects to change active site dynamics and enhance intein activity.

Higher order structures in purine and pyrimidine metabolism

The recent discovery of several forms of higher order protein structures in cells has shifted the paradigm of how we think about protein organization and metabolic regulation. These dynamic and controllable protein assemblies, which are composed of dozens or hundreds of copies of an enzyme or related enzymes, have emerged as important players in myriad cellular processes. We are only beginning to appreciate the breadth of function of these types of macromolecular assemblies. These higher order structures, which can be assembled in response to varied cellular stimuli including changing metabolite concentrations or signaling cascades, give the cell the capacity to modulate levels of biomolecules both temporally and spatially. This provides an added level of control with distinct kinetics and unique features that can be harnessed as a subtle, yet powerful regulatory mechanism. Due, in large part, to advances in structural methods, such as crystallography and cryo-electron microscopy, and the advent of super-resolution microscopy techniques, a rapidly increasing number of these higher order structures are being identified and characterized. In this review, we detail what is known about the structure, function and control mechanisms of these mesoscale protein assemblies, with a particular focus on those involved in purine and pyrimidine metabolism. These structures have important implications both for our understanding of fundamental cellular processes and as fertile ground for new targets for drug discovery and development.

Gleaning structural and functional information from correlations in protein multiple sequence alignments

The availability of vast amounts of protein sequence data facilitates detection of subtle statistical correlations due to imposed structural and functional constraints. Recent breakthroughs using Direct Coupling Analysis (DCA) and related approaches have tapped into correlations believed to be due to compensatory mutations. This has yielded some remarkable results, including substantially improved prediction of protein intra- and inter-domain 3D contacts, of membrane and globular protein structures, of substrate binding sites, and of protein conformational heterogeneity. A complementary approach is Bayesian Partitioning with Pattern Selection (BPPS), which partitions related proteins into hierarchically-arranged subgroups based on correlated residue patterns. These correlated patterns are presumably due to structural and functional constraints associated with evolutionary divergence rather than to compensatory mutations. Hence joint application of DCA- and BPPS-based approaches should help sort out the structural and functional constraints contributing to sequence correlations.

Active site coupling in Plasmodium falciparum GMP synthetase is triggered by domain rotation

GMP synthetase (GMPS), a key enzyme in the purine biosynthetic pathway performs catalysis through a coordinated process across two catalytic pockets for which the mechanism remains unclear. Crystal structures of Plasmodium falciparum GMPS in conjunction with mutational and enzyme kinetic studies reported here provide evidence that an 85° rotation of the GATase domain is required for ammonia channelling and thus for the catalytic activity of this two-domain enzyme. We suggest that conformational changes in helix 371–375 holding catalytic residues and in loop 376–401 along the rotation trajectory trigger the different steps of catalysis, and establish the central role of Glu374 in allostery and inter-domain crosstalk. These studies reveal the mechanism of domain rotation and inter-domain communication, providing a molecular framework for the function of all single polypeptide GMPSs and form a solid basis for rational drug design targeting this therapeutically important enzyme.

GMP synthetase, a key enzyme in purine biosynthesis, is of interest for understanding purine metabolism processes and for developing therapeutic applications. Here, the authors propose a molecular mechanism and the structural basis for the catalytic activity of this enzyme.

Determination of the formylglycinamide ribonucleotide amidotransferase ammonia pathway by combining 3D-RISM theory with experiment

Molecular tunnels in enzyme systems possess variable architecture and are therefore difficult to predict. In this work, we design and apply an algorithm to resolve the pathway followed by ammonia using the bifunctional enzyme formylglycinamide ribonucleotide amidotransferase (FGAR-AT) as a model system. Though its crystal structure has been determined, an ammonia pathway connecting the glutaminase domain to the 30 Å distal FGAR/ATP binding site remains elusive. Crystallography suggested two purported paths: an N-terminal-adjacent path (path 1) and an auxiliary ADP-adjacent path (path 2). The algorithm presented here, RismPath, which enables fast and accurate determination of solvent distribution inside a protein channel, predicted path 2 as the preferred mode of ammonia transfer. Supporting experimental studies validate the identity of the path, and results lead to the conclusion that the residues in the middle of the channel do not partake in catalytic coupling and serve only as channel walls facilitating ammonia transfer.