A review about nothing: are apolar cavities in proteins really empty?

Accounting for nonuniformity of bulk-solvent: A mosaic model

A flat mask-based model is almost universally used in macromolecular crystallography to account for disordered (bulk) solvent. This model assumes any voxel of the crystal unit cell that is not occupied by the atomic model is occupied by the solvent. The properties of this solvent are assumed to be exactly the same across the whole volume of the unit cell. While this is a reasonable approximation in practice, there are a number of scenarios where this model becomes suboptimal. In this work, we enumerate several of these scenarios and describe a new generalized approach to modeling the bulk-solvent which we refer to as mosaic bulk-solvent model. The mosaic bulk-solvent model allows nonuniform features of the solvent in the crystal to be accounted for in a computationally efficient way. It is implemented in the computational crystallography toolbox and the Phenix software.

Microporous water with high gas solubilities

Liquids with permanent microporosity can absorb larger quantities of gas molecules than conventional solvents¹, providing new opportunities for liquid-phase gas storage, transport and reactivity. Current approaches to designing porous liquids rely on sterically bulky solvent molecules or surface ligands and, thus, are not amenable to many important solvents, including water^2-4. Here we report a generalizable thermodynamic strategy to preserve permanent microporosity and impart high gas solubilities to liquid water. Specifically, we show how the external and internal surface chemistry of microporous zeolite and metal-organic framework (MOF) nanocrystals can be tailored to promote the formation of stable dispersions in water while maintaining dry networks of micropores that are accessible to gas molecules. As a result of their permanent microporosity, these aqueous fluids can concentrate gases, including oxygen (O₂) and carbon dioxide (CO₂), to much higher densities than are found in typical aqueous environments. When these fluids are oxygenated, record-high capacities of O₂ can be delivered to hypoxic red blood cells, highlighting one potential application of this new class of microporous liquids for physiological gas transport.

Expanding the Antiviral Potential of the Mosquito Lipid-transfer Protein AEG12 Against SARS-CoV-2 Using Hydrophobic Antiviral Ligands

The mosquito protein AEG12 encompasses a large (~ 3800 ų) hydrophobic cavity which binds and delivers unsaturated fatty acids into biological membranes, allowing it to lyse cells and neutralize a wide range of enveloped viruses. Herein, the lytic and antiviral activities are modified with non‐naturally occurring lipid ligands. We generated novel AEG12 complexes in which the endogenous fatty acid ligands were replaced with hydrophobic viral inhibitors. The resulting compounds modulated cytotoxicity and infectivity against SARS‐CoV‐2, potentially reflecting additional mechanisms of action beyond membrane destabilization. These studies provide valuable insight into the design of novel broad‐spectrum antiviral therapeutics centred on the AEG12 protein scaffold as a delivery vehicle for hydrophobic therapeutic compounds.

Hydrophobic ligands influence the structure, stability, and processing of the major cockroach allergen Bla g 1

The cockroach allergen Bla g 1 forms a novel fold consisting of 12 amphipathic alpha-helices enclosing an exceptionally large hydrophobic cavity which was previously demonstrated to bind a variety of lipids. Since lipid-dependent immunoactivity is observed in numerous allergens, understanding the structural basis of this interaction could yield insights into the molecular determinants of allergenicity. Here, we report atomic modelling of Bla g 1 bound to both fatty-acid and phospholipids ligands, with 8 acyl chains suggested to represent full stoichiometric binding. This unusually high occupancy was verified experimentally, though both modelling and circular dichroism indicate that the general alpha-helical structure is maintained regardless of cargo loading. Fatty-acid cargoes significantly enhanced thermostability while inhibiting cleavage by cathepsin S, an endosomal protease essential for antigen processing and presentation; the latter of which was found to correlate to a decreased production of known T-cell epitopes. Both effects were strongly dependent on acyl chain length, with 18-20 carbons providing the maximal increase in melting temperature (~20 °C) while completely abolishing proteolysis. Diacyl chain cargoes provided similar enhancements to thermostability, but yielded reduced levels of proteolytic resistance. This study describes how the biophysical properties of Bla g 1 ligand binding and digestion may relate to antigen processing, with potential downstream implications for immunogenicity.

How internal cavities destabilize a protein

Although many proteins possess a distinct folded structure lying at a minimum in a funneled free energy landscape, thermal energy causes any protein to continuously access lowly populated excited states. The existence of excited states is an integral part of biological function. Although transitions into the excited states may lead to protein misfolding and aggregation, little structural information is currently available for them. Here, we show how NMR spectroscopy, coupled with pressure perturbation, brings these elusive species to light. As pressure acts to favor states with lower partial molar volume, NMR follows the ensuing change in the equilibrium spectroscopically, with residue-specific resolution. For T4 lysozyme L99A, relaxation dispersion NMR was used to follow the increase in population of a previously identified "invisible" folded state with pressure, as this is driven by the reduction in cavity volume by the flipping-in of a surface aromatic group. Furthermore, multiple partly disordered excited states were detected at equilibrium using pressure-dependent H/D exchange NMR spectroscopy. Here, unfolding reduced partial molar volume by the removal of empty internal cavities and packing imperfections through subglobal and global unfolding. A close correspondence was found for the distinct pressure sensitivities of various parts of the protein and the amount of internal cavity volume that was lost in each unfolding event. The free energies and populations of excited states allowed us to determine the energetic penalty of empty internal protein cavities to be 36 cal⋅Å^-3.

Water in protein hydration and ligand recognition

This review describes selected basics of water in biomolecular recognition. We focus on a qualitative understanding of the most important physical aspects, how these change in magnitude between bulk water and protein environment, and how the roles that water plays for proteins arise from them. These roles include mechanical support, thermal coupling, dielectric screening, mass and charge transport, and the competition with a ligand for the occupation of a binding site. The presence or absence of water has ramifications that range from the thermodynamic binding signature of a single ligand up to cellular survival. The large inhomogeneity in water density, polarity and mobility around a solute is hard to assess in experiment. This is a source of many difficulties in the solvation of protein models and computational studies that attempt to elucidate or predict ligand recognition. The influence of water in a protein binding site on the experimental enthalpic and entropic signature of ligand binding is still a point of much debate. The strong water‐water interaction in enthalpic terms is counteracted by a water molecule's high mobility in entropic terms. The complete arrest of a water molecule's mobility sets a limit on the entropic contribution of a water displacement process, while the solvent environment sets limits on ligand reactivity.

Hysteresis and Allostery in Human UDP-Glucose Dehydrogenase Require a Flexible Protein Core

Human UDP-glucose dehydrogenase (hUGDH) oxidizes UDP-glucose to UDP-glucuronic acid, an essential substrate in the phase II metabolism of drugs. The activity of hUGDH is regulated by the conformation of a buried allosteric switch (T131 loop/α6 helix). Substrate binding induces the allosteric switch to slowly isomerize from an inactive E* conformation to the active E state, which can be observed as enzyme hysteresis. When the feedback inhibitor UDP-xylose binds, the allosteric switch and surrounding residues in the protein core repack, converting the hexamer into an inactive, horseshoe-shaped complex (EΩ). This allosteric transition is facilitated by large cavities and declivities in the protein core that provide the space required to accommodate the alternate packing arrangements. Here, we have used the A104L substitution to fill a cavity in the E state and sterically prevent repacking of the core into the EΩ state. Steady state analysis shows that hUGDHA104L binds UDP-xylose with lower affinity and that the inhibition is no longer cooperative. This means that the allosteric transition to the high-UDP-xylose affinity EΩ state is blocked by the substitution. The crystal structures of hUGDHA104L show that the allosteric switch still adopts the E and E* states, albeit with a more rigid protein core. However, the progress curves of hUGDHA104L do not show hysteresis, which suggests that the E* and E states are now in rapid equilibrium. Our data suggest that hysteresis in native hUGDH originates from the conformational entropy of the E* state protein core.

Crystal Structure of Chloroplastic Thioredoxin f2 from Chlamydomonas reinhardtii Reveals Distinct Surface Properties

Protein disulfide reduction by thioredoxins (TRXs) controls the conformation of enzyme active sites and their multimeric complex formation. TRXs are small oxidoreductases that are broadly conserved in all living organisms. In photosynthetic eukaryotes, TRXs form a large multigenic family, and they have been classified in different types: f, m, x, y, and z types are chloroplastic, while o and h types are located in mitochondria and cytosol. In the model unicellular alga Chlamydomonas reinhardtii, the TRX family contains seven types, with f- and h-types represented by two isozymes. Type-f TRXs interact specifically with targets in the chloroplast, controlling photosynthetic carbon fixation by the Calvin⁻Benson cycle. We solved the crystal structures of TRX f2 and TRX h1 from C. reinhardtii. The systematic comparison of their atomic features revealed a specific conserved electropositive crown around the active site of TRX f, complementary to the electronegative surface of their targets. We postulate that this surface provides specificity to each type of TRX.

Lessons from pressure denaturation of proteins

Although it is now relatively well understood how sequence defines and impacts global protein stability in specific structural contexts, the question of how sequence modulates the configurational landscape of proteins remains to be defined. Protein configurational equilibria are generally characterized by using various chemical denaturants or by changing temperature or pH. Another thermodynamic parameter which is less often used in such studies is high hydrostatic pressure. This review discusses the basis for pressure effects on protein structure and stability, and describes how the unique mechanisms of pressure-induced unfolding can provide unique insights into protein conformational landscapes.

The consequences of cavity creation on the folding landscape of a repeat protein depend upon context

The effect of introducing internal cavities on protein native structure and global stability has been well documented, but the consequences of these packing defects on folding free-energy landscapes have received less attention. We investigated the effects of cavity creation on the folding landscape of the leucine-rich repeat protein pp32 by high-pressure (HP) and urea-dependent NMR and high-pressure small-angle X-ray scattering (HPSAXS). Despite a modest global energetic perturbation, cavity creation in the N-terminal capping motif (N-cap) resulted in very strong deviation from two-state unfolding behavior. In contrast, introduction of a cavity in the most stable, C-terminal half of pp32 led to highly concerted unfolding, presumably because the decrease in stability by the mutations attenuated the N- to C-terminal stability gradient present in WT pp32. Interestingly, enlarging the central cavity of the protein led to the population under pressure of a distinct intermediate in which the N-cap and repeats 1-4 were nearly completely unfolded, while the fifth repeat and the C-terminal capping motif remained fully folded. Thus, despite modest effects on global stability, introducing internal cavities can have starkly distinct repercussions on the conformational landscape of a protein, depending on their structural and energetic context.

Tailoring the Variational Implicit Solvent Method for New Challenges: Biomolecular Recognition and Assembly

Predicting solvation free energies and describing the complex water behavior that plays an important role in essentially all biological processes is a major challenge from the computational standpoint. While an atomistic, explicit description of the solvent can turn out to be too expensive in large biomolecular systems, most implicit solvent methods fail to capture “dewetting” effects and heterogeneous hydration by relying on a pre-established (i.e., guessed) solvation interface. Here we focus on the Variational Implicit Solvent Method, an implicit solvent method that adds water “plasticity” back to the picture by formulating the solvation free energy as a functional of all possible solvation interfaces. We survey VISM's applications to the problem of molecular recognition and report some of the most recent efforts to tailor VISM for more challenging scenarios, with the ultimate goal of including thermal fluctuations into the framework. The advances reported herein pave the way to make VISM a uniquely successful approach to characterize complex solvation properties in the recognition and binding of large-scale biomolecular complexes.

Large scale analysis of protein conformational transitions from aqueous to non-aqueous media

Background

Biocatalysis in organic solvents is nowadays a common practice with a large potential in Biotechnology. Several studies report that proteins which are co-crystallized or soaked in organic solvents preserve their fold integrity showing almost identical arrangements when compared to their aqueous forms. However, it is well established that the catalytic activity of proteins in organic solvents is much lower than in water. In order to explain this diminished activity and to further characterize the behaviour of proteins in non-aqueous environments, we performed a large-scale analysis (1737 proteins) of the conformational diversity of proteins crystallized in aqueous and co-crystallized or soaked in non-aqueous media.

Results

Using proteins’ experimentally determined conformational diversity taken from CoDNaS database, we found that proteins in non-aqueous media display much lower conformational diversity when compared to the corresponding conformers obtained in water. When conformational diversity is compared between conformers obtained in different non-aqueous media, their structural differences are larger and mostly independent of the presence of cognate ligands. We also found that conformers corresponding to non-aqueous media have larger but less flexible cavities, lower number of disordered regions and lower active-site residue mobility.

Conclusions

Our results show that non-aqueous media conformers have specific structural features and that they do not adopt extreme conformations found in aqueous media. This makes them clearly different from their corresponding aqueous conformers.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2044-2) contains supplementary material, which is available to authorized users.

Contribution of buried distal amino acid residues in horse liver alcohol dehydrogenase to structure and catalysis

The dynamics of enzyme catalysis range from the slow time scale (∼ms) for substrate binding and conformational changes to the fast time (∼ps) scale for reorganization of substrates in the chemical step. The contribution of global dynamics to catalysis by alcohol dehydrogenase was tested by substituting five different, conserved amino acid residues that are distal from the active site and located in the hinge region for the conformational change or in hydrophobic clusters. X-ray crystallography shows that the structures for the G173A, V197I, I220 (V, L, or F), V222I, and F322L enzymes complexed with NAD⁺ and an analogue of benzyl alcohol are almost identical, except for small perturbations at the sites of substitution. The enzymes have very similar kinetic constants for the oxidation of benzyl alcohol and reduction of benzaldehyde as compared to the wild-type enzyme, and the rates of conformational changes are not altered. Less conservative substitutions of these amino acid residues, such as G173(V, E, K, or R), V197(G, S, or T), I220(G, S, T, or N), and V222(G, S, or T) produced unstable or poorly expressed proteins, indicating that the residues are critical for global stability. The enzyme scaffold accommodates conservative substitutions of distal residues, and there is no evidence that fast, global dynamics significantly affect the rate constants for hydride transfers. In contrast, other studies show that proximal residues significantly participate in catalysis.

Hydrogen-Deuterium Exchange of Lipoxygenase Uncovers a Relationship between Distal, Solvent Exposed Protein Motions and the Thermal Activation Barrier for Catalytic Proton-Coupled Electron Tunneling

Defining specific pathways for efficient heat transfer from protein-solvent interfaces to their active sites represents one of the compelling and timely challenges in our quest for a physical description of the origins of enzyme catalysis. Enzymatic hydrogen tunneling reactions constitute excellent systems in which to validate experimental approaches to this important question, given the inherent temperature independence of quantum mechanical wave function overlap. Herein, we present the application of hydrogen-deuterium exchange coupled to mass spectrometry toward the spatial resolution of protein motions that can be related to an enzyme's catalytic parameters. Employing the proton-coupled electron transfer reaction of soybean lipoxygenase as proof of principle, we first corroborate the impact of active site mutations on increased local flexibility and, second, uncover a solvent-exposed loop, 15-34 Å from the reactive ferric center whose temperature-dependent motions are demonstrated to mirror the enthalpic barrier for catalytic C-H bond cleavage. A network that connects this surface loop to the active site is structurally identified and supported by changes in kinetic parameters that result from site-specific mutations.

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Experimental Binding Energies in Supramolecular Complexes

On the basis of many literature measurements, a critical overview is given on essential noncovalent interactions in synthetic supramolecular complexes, accompanied by analyses with selected proteins. The methods, which can be applied to derive binding increments for single noncovalent interactions, start with the evaluation of consistency and additivity with a sufficiently large number of different host-guest complexes by applying linear free energy relations. Other strategies involve the use of double mutant cycles, of molecular balances, of dynamic combinatorial libraries, and of crystal structures. Promises and limitations of these strategies are discussed. Most of the analyses stem from solution studies, but a few also from gas phase. The empirically derived interactions are then presented on the basis of selected complexes with respect to ion pairing, hydrogen bonding, electrostatic contributions, halogen bonding, π-π-stacking, dispersive forces, cation-π and anion-π interactions, and contributions from the hydrophobic effect. Cooperativity in host-guest complexes as well as in self-assembly, and entropy factors are briefly highlighted. Tables with typical values for single noncovalent free energies and polarity parameters are in the Supporting Information.

Characterization of a novel water pocket inside the human Cx26 hemichannel structure

Connexins (Cxs) are a family of vertebrate proteins constituents of gap junction channels (GJCs) that connect the cytoplasm of adjacent cells by the end-to-end docking of two Cx hemichannels. The intercellular transfer through GJCs occurs by passive diffusion allowing the exchange of water, ions, and small molecules. Despite the broad interest to understand, at the molecular level, the functional state of Cx-based channels, there are still many unanswered questions regarding structure-function relationships, perm-selectivity, and gating mechanisms. In particular, the ordering, structure, and dynamics of water inside Cx GJCs and hemichannels remains largely unexplored. In this work, we describe the identification and characterization of a believed novel water pocket-termed the IC pocket-located in-between the four transmembrane helices of each human Cx26 (hCx26) monomer at the intracellular (IC) side. Using molecular dynamics (MD) simulations to characterize hCx26 internal water structure and dynamics, six IC pockets were identified per hemichannel. A detailed characterization of the dynamics and ordering of water including conformational variability of residues forming the IC pockets, together with multiple sequence alignments, allowed us to propose a functional role for this cavity. An in vitro assessment of tracer uptake suggests that the IC pocket residue Arg-143 plays an essential role on the modulation of the hCx26 hemichannel permeability.

Thermal green protein, an extremely stable, nonaggregating fluorescent protein created by structure-guided surface engineering

In this article, we describe the engineering and X-ray crystal structure of Thermal Green Protein (TGP), an extremely stable, highly soluble, non-aggregating green fluorescent protein. TGP is a soluble variant of the fluorescent protein eCGP123, which despite being highly stable, has proven to be aggregation-prone. The X-ray crystal structure of eCGP123, also determined within the context of this paper, was used to carry out rational surface engineering to improve its solubility, leading to TGP. The approach involved simultaneously eliminating crystal lattice contacts while increasing the overall negative charge of the protein. Despite intentional disruption of lattice contacts and introduction of high entropy glutamate side chains, TGP crystallized readily in a number of different conditions and the X-ray crystal structure of TGP was determined to 1.9 Å resolution. The structural reasons for the enhanced stability of TGP and eCGP123 are discussed. We demonstrate the utility of using TGP as a fusion partner in various assays and significantly, in amyloid assays in which the standard fluorescent protein, EGFP, is undesirable because of aberrant oligomerization.

Structure-function relationships in human testis-determining factor SRY: an aromatic buttress underlies the specific DNA-bending surface of a high mobility group (HMG) box

Human testis determination is initiated by SRY, a Y-encoded architectural transcription factor. Mutations in SRY cause 46 XY gonadal dysgenesis with female somatic phenotype (Swyer syndrome) and confer a high risk of malignancy (gonadoblastoma). Such mutations cluster in the SRY high mobility group (HMG) box, a conserved motif of specific DNA binding and bending. To explore structure-function relationships, we constructed all possible substitutions at a site of clinical mutation (W70L). Our studies thus focused on a core aromatic residue (position 15 of the consensus HMG box) that is invariant among SRY-related HMG box transcription factors (the SOX family) and conserved as aromatic (Phe or Tyr) among other sequence-specific boxes. In a yeast one-hybrid system sensitive to specific SRY-DNA binding, the variant domains exhibited reduced (Phe and Tyr) or absent activity (the remaining 17 substitutions). Representative nonpolar variants with partial or absent activity (Tyr, Phe, Leu, and Ala in order of decreasing side-chain volume) were chosen for study in vitro and in mammalian cell culture. The clinical mutation (Leu) was found to markedly impair multiple biochemical and cellular activities as respectively probed through the following: (i) in vitro assays of specific DNA binding and protein stability, and (ii) cell culture-based assays of proteosomal degradation, nuclear import, enhancer DNA occupancy, and SRY-dependent transcriptional activation. Surprisingly, however, DNA bending is robust to this or the related Ala substitution that profoundly impairs box stability. Together, our findings demonstrate that the folding, trafficking, and gene-regulatory function of SRY requires an invariant aromatic "buttress" beneath its specific DNA-bending surface.

Role of cavities and hydration in the pressure unfolding of T4 lysozyme

It is well known that high hydrostatic pressures can induce the unfolding of proteins. The physical underpinnings of this phenomenon have been investigated extensively but remain controversial. Changes in solvation energetics have been commonly proposed as a driving force for pressure-induced unfolding. Recently, the elimination of void volumes in the native folded state has been argued to be the principal determinant. Here we use the cavity-containing L99A mutant of T4 lysozyme to examine the pressure-induced destabilization of this multidomain protein by using solution NMR spectroscopy. The cavity-containing C-terminal domain completely unfolds at moderate pressures, whereas the N-terminal domain remains largely structured to pressures as high as 2.5 kbar. The sensitivity to pressure is suppressed by the binding of benzene to the hydrophobic cavity. These results contrast to the pseudo-WT protein, which has a residual cavity volume very similar to that of the L99A-benzene complex but shows extensive subglobal reorganizations with pressure. Encapsulation of the L99A mutant in the aqueous nanoscale core of a reverse micelle is used to examine the hydration of the hydrophobic cavity. The confined space effect of encapsulation suppresses the pressure-induced unfolding transition and allows observation of the filling of the cavity with water at elevated pressures. This indicates that hydration of the hydrophobic cavity is more energetically unfavorable than global unfolding. Overall, these observations point to a range of cooperativity and energetics within the T4 lysozyme molecule and illuminate the fact that small changes in physical parameters can significantly alter the pressure sensitivity of proteins.

Analysis of carotenoid isomerase activity in a prototypical carotenoid cleavage enzyme, apocarotenoid oxygenase (ACO)

Carotenoid cleavage enzymes (CCEs) constitute a group of evolutionarily related proteins that metabolize a variety of carotenoid and non-carotenoid substrates. Typically, these enzymes utilize a non-heme iron center to oxidatively cleave a carbon-carbon double bond of a carotenoid substrate. Some members also isomerize specific double bonds in their substrates to yield cis-apocarotenoid products. The apocarotenoid oxygenase from Synechocystis has been hypothesized to represent one such member of this latter category of CCEs. Here, we developed a novel expression and purification protocol that enabled production of soluble, native ACO in quantities sufficient for high resolution structural and spectroscopic investigation of its catalytic mechanism. High performance liquid chromatography and Raman spectroscopy revealed that ACO exclusively formed all-trans products. We also found that linear polyoxyethylene detergents previously used for ACO crystallization strongly inhibited the apocarotenoid oxygenase activity of the enzyme. We crystallized the native enzyme in the absence of apocarotenoid substrate and found electron density in the active site that was similar in appearance to the density previously attributed to a di-cis-apocarotenoid intermediate. Our results clearly demonstrated that ACO is in fact a non-isomerizing member of the CCE family. These results indicate that careful selection of detergent is critical for the success of structural studies aimed at elucidating structures of CCE-carotenoid/retinoid complexes.

RB69 DNA polymerase structure, kinetics, and fidelity

This review will summarize our structural and kinetic studies of RB69 DNA polymerase (RB69pol) as well as selected variants of the wild-type enzyme that were undertaken to obtain a deeper understanding of the exquisitely high fidelity of B family replicative DNA polymerases. We discuss how the structures of the various RB69pol ternary complexes can be used to rationalize the results obtained from pre-steady-state kinetic assays. Our main findings can be summarized as follows. (i) Interbase hydrogen bond interactions can increase catalytic efficiency by 5000-fold; meanwhile, base selectivity is not solely determined by the number of hydrogen bonds between the incoming dNTP and the templating base. (ii) Minor-groove hydrogen bond interactions at positions n – 1 and n – 2 of the primer strand and position n – 1 of the template strand in RB69pol ternary complexes are essential for efficient primer extension and base selectivity. (iii) Partial charge interactions among the incoming dNTP, the penultimate base pair, and the hydration shell surrounding the incoming dNTP modulate nucleotide insertion efficiency and base selectivity. (iv) Steric clashes between mismatched incoming dNTPs and templating bases with amino acid side chains in the nascent base pair binding pocket (NBP) as well as weak interactions and large gaps between the incoming dNTPs and the templating base are some of the reasons that incorrect dNTPs are incorporated so inefficiently by wild-type RB69pol. In addition, we developed a tC°–tC_nitro Förster resonance energy transfer assay to monitor partitioning of the primer terminus between the polymerase and exonuclease subdomains.

Effects of cavities at the nicotinamide binding site of liver alcohol dehydrogenase on structure, dynamics and catalysis

A role for protein dynamics in enzymatic catalysis of hydrogen transfer has received substantial scientific support, but the connections between protein structure and catalysis remain to be established. Valine residues 203 and 207 are at the binding site for the nicotinamide ring of the coenzyme in liver alcohol dehydrogenase and have been suggested to facilitate catalysis with “protein-promoting vibrations” (PPV). We find that the V207A substitution has small effects on steady-state kinetic constants and the rate of hydrogen transfer; the introduced cavity is empty and is tolerated with minimal effects on structure (determined at 1.2 Å for the complex with NAD⁺ and 2,3,4,5,6-pentafluorobenzyl alcohol). Thus, no evidence is found to support a role for Val-207 in the dynamics of catalysis. The protein structures and ligand geometries (including donor–acceptor distances) in the V203A enzyme complexed with NAD⁺ and 2,3,4,5,6-pentafluorobenzyl alcohol or 2,2,2-trifluoroethanol (determined at 1.1 Å) are very similar to those for the wild-type enzyme, except that the introduced cavity accommodates a new water molecule that contacts the nicotinamide ring. The structures of the V203A enzyme complexes suggest, in contrast to previous studies, that the diminished tunneling and decreased rate of hydride transfer (16-fold, relative to that of the wild-type enzyme) are not due to differences in ground-state ligand geometries. The V203A substitution may alter the PPV and the reorganization energy for hydrogen transfer, but the protein scaffold and equilibrium thermal motions within the Michaelis complex may be more significant for enzyme catalysis.

The designability of protein switches by chemical rescue of structure: mechanisms of inactivation and reactivation

The ability to selectively activate function of particular proteins via pharmacological agents is a longstanding goal in chemical biology. Recently, we reported an approach for designing a de novo allosteric effector site directly into the catalytic domain of an enzyme. This approach is distinct from traditional chemical rescue of enzymes in that it relies on disruption and restoration of structure, rather than active site chemistry, as a means to achieve modulate function. However, rationally identifying analogous de novo binding sites in other enzymes represents a key challenge for extending this approach to introduce allosteric control into other enzymes. Here we show that mutation sites leading to protein inactivation via tryptophan-to-glycine substitution and allowing (partial) reactivation by the subsequent addition of indole are remarkably frequent. Through a suite of methods including a cell-based reporter assay, computational structure prediction and energetic analysis, fluorescence studies, enzymology, pulse proteolysis, X-ray crystallography, and hydrogen-deuterium mass spectrometry, we find that these switchable proteins are most commonly modulated indirectly, through control of protein stability. Addition of indole in these cases rescues activity not by reverting a discrete conformational change, as we had observed in the sole previously reported example, but rather rescues activity by restoring protein stability. This important finding will dramatically impact the design of future switches and sensors built by this approach, since evaluating stability differences associated with cavity-forming mutations is a far more tractable task than predicting allosteric conformational changes. By analogy to natural signaling systems, the insights from this study further raise the exciting prospect of modulating stability to design optimal recognition properties into future de novo switches and sensors built through chemical rescue of structure.

Changing hydration level in an internal cavity modulates the proton affinity of a key glutamate in cytochrome c oxidase

Cytochrome c oxidase contributes to the transmembrane proton gradient by removing two protons from the high-pH side of the membrane each time the binuclear center active site is reduced. One proton goes to the binuclear center, whereas the other is pumped to the low-pH periplasmic space. Glutamate 286 (Glu286) has been proposed to serve as a transiently deprotonated proton donor. Using unrestrained atomistic molecular dynamics simulations, we show that the size of and water distribution in the hydrophobic cavity that holds Glu286 is controlled by the protonation state of the propionic acid of heme a3, a group on the proton outlet pathway. Protonation of the propionate disrupts hydrogen bonding to two side chains, allowing a loop to swing open. Continuum electrostatics and atomistic free-energy perturbation calculations show that the resultant changes in hydration and electrostatic interactions lower the Glu proton affinity by at least 5 kcal/mol. These changes in the internal hydration level occur in the absence of major conformational transitions and serve to stabilize needed transient intermediates in proton transport. The trigger is not the protonation of the Glu of interest, but rather the protonation of a residue ∼10 Å away. Thus, unlike local water penetration to stabilize a new charge, this finding represents a specific role for water molecules in the protein interior, mediating proton transfers and facilitating ion transport.

MP:PD--a data base of internal packing densities, internal packing defects and internal waters of helical membrane proteins

The membrane protein packing database (MP:PD) (http://proteinformatics.charite.de/mppd) is a database of helical membrane proteins featuring internal atomic packing densities, cavities and waters. Membrane proteins are not tightly packed but contain a considerable number of internal cavities that differ in volume, polarity and solvent accessibility as well as in their filling with internal water. Internal cavities are supposed to be regions of high physical compressibility. By serving as mobile hydrogen bonding donors or acceptors, internal waters likely facilitate transition between different functional states. Despite these distinct functional roles, internal cavities of helical membrane proteins are not well characterized, mainly because most internal waters are not resolved by crystal structure analysis. Here we combined various computational biophysical techniques to characterize internal cavities, reassign positions of internal waters and calculate internal packing densities of all available helical membrane protein structures and stored them in MP:PD. The database can be searched using keywords and entries can be downloaded. Each entry can be visualized in Provi, a Jmol-based protein viewer that provides an integrated display of low energy waters alongside membrane planes, internal packing density, hydrophobic cavities and hydrogen bonds.

Malachite green mediates homodimerization of antibody VL domains to form a fluorescent ternary complex with singular symmetric interfaces

We report that a symmetric small-molecule ligand mediates the assembly of antibody light chain variable domains (VLs) into a correspondent symmetric ternary complex with novel interfaces. The L5* fluorogen activating protein is a VL domain that binds malachite green (MG) dye to activate intense fluorescence. Crystallography of liganded L5* reveals a 2:1 protein:ligand complex with inclusive C2 symmetry, where MG is almost entirely encapsulated between an antiparallel arrangement of the two VL domains. Unliganded L5* VL domains crystallize as a similar antiparallel VL/VL homodimer. The complementarity-determining regions are spatially oriented to form novel VL/VL and VL/ligand interfaces that tightly constrain a propeller conformer of MG. Binding equilibrium analysis suggests highly cooperative assembly to form a very stable VL/MG/VL complex, such that MG behaves as a strong chemical inducer of dimerization. Fusion of two VL domains into a single protein tightens MG binding over 1000-fold to low picomolar affinity without altering the large binding enthalpy, suggesting that bonding interactions with ligand and restriction of domain movements make independent contributions to binding. Fluorescence activation of a symmetrical fluorogen provides a selection mechanism for the isolation and directed evolution of ternary complexes where unnatural symmetric binding interfaces are favored over canonical antibody interfaces. As exemplified by L5*, these self-reporting complexes may be useful as modulators of protein association or as high-affinity protein tags and capture reagents.

MOLE 2.0: advanced approach for analysis of biomacromolecular channels

BACKGROUND

Channels and pores in biomacromolecules (proteins, nucleic acids and their complexes) play significant biological roles, e.g., in molecular recognition and enzyme substrate specificity.

RESULTS

We present an advanced software tool entitled MOLE 2.0, which has been designed to analyze molecular channels and pores. Benchmark tests against other available software tools showed that MOLE 2.0 is by comparison quicker, more robust and more versatile. As a new feature, MOLE 2.0 estimates physicochemical properties of the identified channels, i.e., hydropathy, hydrophobicity, polarity, charge, and mutability. We also assessed the variability in physicochemical properties of eighty X-ray structures of two members of the cytochrome P450 superfamily.

CONCLUSION

Estimated physicochemical properties of the identified channels in the selected biomacromolecules corresponded well with the known functions of the respective channels. Thus, the predicted physicochemical properties may provide useful information about the potential functions of identified channels. The MOLE 2.0 software is available at http://mole.chemi.muni.cz.

Alteration in the cavity size adjacent to the active site of RB69 DNA polymerase changes its conformational dynamics

Internal cavities are a common feature of many proteins, often having profound effects on the dynamics of their interactions with substrate and binding partners. RB69 DNA polymerase (pol) has a hydrophobic cavity right below the nucleotide binding pocket at the tip of highly conserved L415 side chain. Replacement of this residue with Gly or Met in other B family pols resulted in higher mutation rates. When similar substitutions for L415 were introduced into RB69pol, only L415A and L415G had dramatic effects on pre-steady-state kinetic parameters, reducing base selectivity by several hundred fold. On the other hand, the L415M variant behaved like the wild-type. Using a novel tC^o-tC_nitro Förster Resonance Energy Transfer (FRET) assay, we were able to show that the partition of the primer terminus between pol and exonuclease (exo) domains was compromised with the L415A and L415G mutants, but not with the L415M variant. These results could be rationalized by changes in their structures as determined by high resolution X-ray crystallography.

Water mediated ligand functional group cooperativity: the contribution of a methyl group to binding affinity is enhanced by a COO(-) group through changes in the structure and thermodynamics of the hydration waters of ligand-thermolysin complexes

Ligand functional groups can modulate the contributions of one another to the ligand-protein binding thermodynamics, producing either positive or negative cooperativity. Data presented for four thermolysin phosphonamidate inhibitors demonstrate that the differential binding free energy and enthalpy caused by replacement of a H with a Me group, which binds in the well-hydrated S2' pocket, are more favorable in presence of a ligand carboxylate. The differential entropy is however less favorable. Dissection of these differential thermodynamic parameters, X-ray crystallography, and density-functional theory calculations suggest that these cooperativities are caused by variations in the thermodynamics of the complex hydration shell changes accompanying the H→Me replacement. Specifically, the COO(-) reduces both the enthalpic penalty and the entropic advantage of displacing water molecules from the S2' pocket and causes a subsequent acquisition of a more enthalpically, less entropically, favorable water network. This study contributes to understanding the important role water plays in ligand-protein binding.

Comparative void-volume analysis of psychrophilic and mesophilic enzymes: Structural bioinformatics of psychrophilic enzymes reveals sources of core flexibility

Background

Psychrophiles, cold-adapted organisms, have adapted to live at low temperatures by using a variety of mechanisms. Their enzymes are active at cold temperatures by being structurally more flexible than mesophilic enzymes. Even though, there are some indications of the possible structural mechanisms by which psychrophilic enzymes are catalytic active at cold temperatures, there is not a generalized structural property common to all psychrophilic enzymes.

Results

We examine twenty homologous enzyme pairs from psychrophiles and mesophiles to investigate flexibility as a key characteristic for cold adaptation. B-factors in protein X-ray structures are one way to measure flexibility. Comparing psychrophilic to mesophilic protein B-factors reveals that psychrophilic enzymes are more flexible in 5-turn and strand secondary structures. Enzyme cavities, identified using CASTp at various probe sizes, indicate that psychrophilic enzymes have larger average cavity sizes at probe radii of 1.4-1.5 Å, sufficient for water molecules. Furthermore, amino acid side chains lining these cavities show an increased frequency of acidic groups in psychrophilic enzymes.

Conclusions

These findings suggest that embedded water molecules may play a significant role in cavity flexibility, and therefore, overall protein flexibility. Thus, our results point to the important role enzyme flexibility plays in adaptation to cold environments.

Water in the polar and nonpolar cavities of the protein interleukin-1β

Water in the protein interior serves important structural and functional roles and is also increasingly recognized as a relevant factor in drug binding. The nonpolar cavity in the protein interleukin-1β has been reported to be filled by water on the basis of some experiments and simulations and to be empty on the basis of others. Here we study the thermodynamics of filling the central nonpolar cavity and the four polar cavities of interleukin-1β by molecular dynamics simulation. We use different water models (TIP3P and SPC/E) and protein force fields (amber94 and amber03) to calculate the semigrand partition functions term by term that quantify the hydration equilibria. We consistently find that water in the central nonpolar cavity is thermodynamically unstable, independent of force field and water model. The apparent reason is the relatively small size of the cavity, with a volume less than ∼80 Å(3). Our results are consistent with the most recent X-ray crystallographic and simulation studies but disagree with an earlier interpretation of nuclear magnetic resonance (NMR) experiments probing protein-water interactions. We show that, at least semiquantitatively, the measured nuclear Overhauser effects indicating the proximity of water to the methyl groups lining the nonpolar cavity can, in all likelihood, be attributed to interactions with buried and surface water molecules near the cavity. The same methods applied to determine the occupancy of the polar cavities show that they are filled by the same number of water molecules observed in crystallography, thereby validating the theoretical and simulation methods used to study the water occupancy in the nonpolar protein cavity.