Molecular dynamics simulation reveals a surface salt bridge forming a kinetic trap in unfolding of truncated Staphylococcal nuclease

Interaction strength of osmolytes with the anion of a salt-bridge determines its stability

In order to understand the role of osmolytes in regulating physicochemical behavior of proteins, we investigated the influence of protein destabilizing (urea and guanidinium chloride) and stabilizing osmolytes (TMAO, glycerol, and betaine) on a model salt-bridge (SB) formed between structural analogues of arginine and glutamate/aspartate sidechains in a solvent continuum using first-principles quantum chemical calculations based on DFT and MP2 methods. The binding strength of the osmolyte with the SB is found to be in the order of betaine > TMAO > Gdm+ > glycerol > urea. The osmolytes (TMAO and betaine) that preferentially bind to the SB cation have a marginal influence on SB stability. Also, pure π-π stacking interaction between Gdm+ and the SB cation plays an insignificant role in destabilizing the SB. In fact, the interaction strength of osmolytes with the SB anion mainly determines the stability of SB. For instance, a competition between Gdm+ and the SB cation to bind with the SB anion is responsible for instability and subsequent dissociation of the SB. The competition provided by other osmolytes is too weak to break the SB. Exploiting this information, we designed three structural derivatives of Gdm+, all having a stronger interaction with SB anion, and thereby show a stronger SB dissociation potential. Furthermore, we find an excellent linear anti-correlation between SB interaction energy and the energy of interaction between osmolyte and the SB anion, which suggests that by knowing only the strength of osmolyteacetate interaction, one can predict the influence of osmolytes on the salt-bridge instability. This information is useful in fine-tuning the SB dissociation power of Gdm+, which has a practical significance in obtaining the mechanistic insight into the influence of GdmCl on protein stability. Our results also provide a basis for understanding the chemistry of other ion-pairs formed between a cationic hydrogen donor and an anionic acceptor.

Enhanced Thermostability and Enzymatic Activity of Cel6A Variants from Thermobifida fusca by Empirical Domain Engineering (Short Title: Domain Engineering of Cel6A)

Cellulases are a set of lignocellulolytic enzymes, capable of producing eco-friendly low-cost renewable bioethanol. However, low stability and hydrolytic activity limit their wide-scale applicability at the industrial scale. In this work, we report the domain engineering of endoglucanase (Cel6A) of Thermobifida fusca to improve their catalytic activity and thermal stability. Later, enzymatic activity and thermostability of the most efficient variant named as Cel6A.CBC was analyzed by molecular dynamics simulations. This variant demonstrated profound activity against soluble and insoluble cellulosic substrates like filter paper, alkali-treated bagasse, regenerated amorphous cellulose (RAC), and bacterial microcrystalline cellulose. The variant Cel6A.CBC showed the highest catalysis of carboxymethyl cellulose (CMC) and other related insoluble substrates at a pH of 6.0 and a temperature of 60 °C. Furthermore, a sound rationale was observed between experimental findings and molecular modeling of Cel6A.CBC which revealed thermostability of Cel6A.CBC at 26.85, 60.85, and 74.85 °C as well as structural flexibility at 126.85 °C. Therefore, a thermostable derivative of Cel6A engineered in the present work has enhanced biological performance and can be a useful construct for the mass production of bioethanol from plant biomass.

Surfactant Charge Modulates Structure and Stability of Lipase-Embedded Monolayers at Marine-Relevant Aerosol Surfaces

Lipases, as well as other enzymes, are present and active within the sea surface microlayer (SSML). Upon bubble bursting, lipases partition into sea spray aerosol (SSA) along with surface-active molecules such as lipids. Lipases are likely to be embedded in the lipid monolayer at the SSA surface and thus have the potential to influence SSA interfacial structure and chemistry. Elucidating the structure of the lipid monolayer at SSA interfaces and how this structure is altered upon interaction with a protein system like lipase is of interest, given the importance of how aerosols interact with sunlight, influence cloud formation, and provide surfaces for chemical reactions. Herein, we report an integrated experimental and computational study of Burkholderia cepacia lipase (BCL) embedded in a lipid monolayer and highlight the important role of electrostatic, rather than hydrophobic, interactions as a driver for monolayer stability. Specifically, we combine Langmuir film experiments and molecular dynamics (MD) simulations to examine the detailed interactions between the zwitterionic dipalmitoylphosphatidylcholine (DPPC) monolayer and BCL. Upon insertion of BCL from the underlying subphase into the lipid monolayer, it is shown that BCL permeates and largely disorders the monolayer while strongly interacting with zwitterionic DPPC molecules, as experimentally observed by Langmuir adsorption curves and infrared reflectance absorbance spectroscopy. Explicitly solvated, all-atom MD is then used to provide insights into inter- and intramolecular interactions that drive these observations, with specific attention to the formation of salt bridges or ionic-bonding interactions. We show that after insertion into the DPPC monolayer, lipase is maintained at high surface pressures and in large BCL concentrations by forming a salt-bridge-stabilized lipase-DPPC complex. In comparison, when embedded in an anionic monolayer at low surface pressures, BCL preferentially forms intramolecular salt bridges, reducing its total favorable interactions with the surfactant and partitioning out of the monolayer shortly after injection. Overall, this study shows that the structure and dynamics of lipase-embedded SSA surfaces vary based on surface charge and pressure and that these variations have the potential to differentially modulate the properties of marine aerosols.

Influence of Glu/Arg, Asp/Arg, and Glu/Lys Salt Bridges on α-Helical Stability and Folding Kinetics

Using a combination of ultraviolet circular dichroism, temperature-jump transient-infrared spectroscopy, and molecular dynamics simulations, we investigate the effect of salt bridges between different types of charged amino-acid residue pairs on α-helix folding. We determine the stability and the folding and unfolding rates of 12 alanine-based α-helical peptides, each of which has a nearly identical composition containing three pairs of positively and negatively charged residues (either Glu(-)/Arg(+), Asp(-)/Arg(+), or Glu(-)/Lys(+)). Within each set of peptides, the distance and order of the oppositely charged residues in the peptide sequence differ, such that they have different capabilities of forming salt bridges. Our results indicate that stabilizing salt bridges (in which the interacting residues are spaced and ordered such that they favor helix formation) speed up α-helix formation by up to 50% and slow down the unfolding of the α-helix, whereas salt bridges with an unfavorable geometry have the opposite effect. Comparing the peptides with different types of charge pairs, we observe that salt bridges between side chains of Glu(-) and Arg(+) are most favorable for the speed of folding, probably because of the larger conformational space of the salt-bridging Glu(-)/Arg(+) rotamer pairs compared to Asp(-)/Arg(+) and Glu(-)/Lys(+). We speculate that the observed impact of salt bridges on the folding kinetics might explain why some proteins contain salt bridges that do not stabilize the final, folded conformation.

Free energy landscape of a minimalist salt bridge model

Salt bridges are essential to protein stability and dynamics. Despite the importance, there has been scarce of detailed discussion on how salt bridge partners interact with each other in distinct solvent exposed environments. In this study, employing a recent generalized orthogonal space tempering (gOST) method, we enabled efficient molecular dynamics simulation of repetitive breaking and reforming of salt bridge structures within a minimalist salt-bridge model, the Asp-Arg dipeptide and thereby were able to map its detailed free energy landscape in aqueous solution. Free energy surface analysis shows that although individually-solvated states are more favorable, salt-bridge states still occupy a noticeable portion of the overall population. Notably, the competing forces, e.g. intercharge attractions that drive the formation of salt bridges and solvation forces that pull the charged groups away from each other, are energetically comparable. As the result, the salt bridge stability is highly tunable by local environments; for instance when local water molecules are perturbed to interact more strongly with each other, the population of the salt-bridge states is likely to increase. Our results reveal the critical role of local solvent structures in modulating salt-bridge partner interactions and imply the importance of water fluctuations on conformational dynamics that involves solvent accessible salt bridge formations.

Solvent-Exposed Salt Bridges Influence the Kinetics of α-Helix Folding and Unfolding

Salt bridges are known to play an essential role in the thermodynamic stability of the folded conformation of many proteins, but their influence on the kinetics of folding remains largely unknown. Here, we investigate the effect of Glu-Arg salt bridges on the kinetics of α-helix folding using temperature-jump transient-infrared spectroscopy and steady-state UV circular dichroism. We find that geometrically optimized salt bridges (Glu^- and Arg⁺ are spaced four peptide units apart, and the Glu/Arg order is such that the side-chain rotameric preferences favor salt-bridge formation) significantly speed up folding and slow down unfolding, whereas salt bridges with unfavorable geometry slow down folding and slightly speed up unfolding. Our observations suggest a possible explanation for the surprising fact that many biologically active proteins contain salt bridges that do not stabilize the native conformation: these salt bridges might have a kinetic rather than a thermodynamic function.

Stability mechanisms of a thermophilic laccase probed by molecular dynamics

Laccases are highly stable, industrially important enzymes capable of oxidizing a large range of substrates. Causes for their stability are, as for other proteins, poorly understood. In this work, multiple-seed molecular dynamics (MD) was applied to a Trametes versicolor laccase in response to variable ionic strengths, temperatures, and glycosylation status. Near-physiological conditions provided excellent agreement with the crystal structure (average RMSD ∼0.92 Å) and residual agreement with experimental B-factors. The persistence of backbone hydrogen bonds was identified as a key descriptor of structural response to environment, whereas solvent-accessibility, radius of gyration, and fluctuations were only locally relevant. Backbone hydrogen bonds decreased systematically with temperature in all simulations (∼9 per 50 K), probing structural changes associated with enthalpy-entropy compensation. Approaching T_opt (∼350 K) from 300 K, this change correlated with a beginning “unzipping” of critical β-sheets. 0 M ionic strength triggered partial denucleation of the C-terminal (known experimentally to be sensitive) at 400 K, suggesting a general salt stabilization effect. In contrast, F⁻ (but not Cl⁻) specifically impaired secondary structure by formation of strong hydrogen bonds with backbone NH, providing a mechanism for experimentally observed small anion destabilization, potentially remedied by site-directed mutagenesis at critical intrusion sites. N-glycosylation was found to support structural integrity by increasing persistent backbone hydrogen bonds by ∼4 across simulations, mainly via prevention of F⁻ intrusion. Hydrogen-bond loss in distinct loop regions and ends of critical β-sheets suggest potential strategies for laboratory optimization of these industrially important enzymes.

Structural study of carboxylesterase from hyperthermophilic bacteria Geobacillus stearothermophilus by molecular dynamics simulation

Carboxylesterases are ubiquitous enzymes with important physiological, industrial and medical applications such as synthesis and hydrolysis of stereo specific compounds, including the metabolic processing of drugs, and antimicrobial agents. Here, we have performed molecular dynamics simulations of carboxylesterase from hyperthermophilic bacterium Geobacillus stearothermophilus (GsEst) for 10ns each at five different temperatures namely at 300K, 343K, 373K, 473K and 500K. Profiles of root mean square fluctuation (RMSF) identify thermostable and thermosensitive regions of GsEst. Unfolding of GsEst initiates at the thermosensitive alpha-helices and proceeds to the thermostable beta-sheets. Five ion-pairs have been identified as critical ion-pairs for thermostability and are maintained stably throughout the higher temperature simulations. A detailed investigation of the active site residues of this enzyme suggests that the geometry of this site is well preserved up to 373K. Furthermore, the hydrogen bonds between Asp188 and His218 of the active site are stably maintained at higher temperatures imparting stability of this site. Radial distribution functions (RDFs) show similar pattern of solvent ordering and water penetration around active site residues up to 373K. Principal component analysis suggests that the motion of the entire protein as well as the active site is similar at 300K, 343K and 373K. Our study may help to identify the factors responsible for thermostability of GsEst that may endeavor to design enzymes with enhanced thermostability.

Insights into thermal stability of thermophilic nitrile hydratases by molecular dynamics simulation

Thermal stability is of great importance for industrial enzymes. Here we explored the thermal-stable mechanism of thermophilic nitrile hydratases (NHases) utilizing a molecular dynamic simulation. At a nanosecond timescale, profiles of root mean square fluctuation (RMSF) of two thermophilic NHases, 1UGQ and 1V29, under enhancing thermal stress were carried out at 300 K, 320 K, 350 K and 370 K, respectively. Results showed that the region A1 (211-231 aa) and A2 (305-316 aa) in 1UGQ, region B1 (186-192 aa) in 1V29, and most of terminal ends in both enzymes are hyper-sensitive. Salt-bridge analyses revealed that in one hand, salt-bridges contributed to maintaining the rigid structure and stable performance of the thermophilic 1UGQ and 1V29; in the other hand, salt-bridges involved in thermal sensitive regions are relatively weak and prone to be broken at elevated temperature, thereby cannot hold the stable conformation of the spatial neighborhood. In 1V29, region A1 was stabilized by a well-organized hook-hook like cluster with multiple salt-bridge interactions, region A2 was stabilized by two strong salt-bridge interactions of GLU52-ARG332 and GLU334-ARG332. In 1UGQ, the absence of a charged residue decreased its thermal sensitivity of region B1, and the formation of a small beta-sheet containing a stable salt-bridge in C-beta-terminal significantly enhanced its thermal stability. By radius of gyration calculation containing or eliminating the thermal sensitive regions, we quantified the contribution of thermal sensitive regions for thermal sensitivity of 1UGQ and 1V29. Consequently, we presented strategies to improve thermal stability of the industrialized mesophilic NHase by introducing stable salt-bridge interactions into its thermal sensitive regions.

Analyzing large-scale structural change in proteins: comparison of principal component projection and Sammon mapping

Effective analysis of large-scale conformational transitions in macromolecules requires transforming them into a lower dimensional representation that captures the dominant motions. Herein, we apply and compare two different dimensionality reduction techniques, namely, principal component analysis (PCA), a linear method, and Sammon mapping, which is nonlinear. The two methods are used to analyze four different protein transition pathways of varying complexity, obtained by using either the conjugate peak refinement method or constrained molecular dynamics. For the return-stroke in myosin, both Sammon mapping and PCA show that the conformational change is dominated by a simple rotation of a rigid body. Also, in the case of the T-->R transition in hemoglobin, both methods are able to identify the two main quaternary transition events. In contrast, in the cases of the unfolding transition of staphylococcal nuclease or the signaling switch of Ras p21, which are both more complex conformational transitions, only Sammon mapping is able to identify the distinct phases of motion.

Molecular dynamics simulations of the periplasmic ferric-hydroxamate binding protein FhuD

FhuD is a periplasmic binding protein (PBP) that, under iron-limiting conditions, transports various hydroxamate-type siderophores from the outer membrane receptor (FhuA) to the inner membrane ATP-binding cassette transporter (FhuBC). Unlike many other PBPs, FhuD possesses two independently folded domains that are connected by an alpha-helix rather than two or three central beta-strands. Crystal structures of FhuD with and without bound gallichrome have provided some insight into the mechanism of siderophore binding as well as suggested a potential mechanism for FhuD binding to FhuB. Since the alpha-helix connecting the two domains imposes greater rigidity on the structure relative to the beta-strands in other 'classical' PBPs, these structures reveal no large conformational change upon binding a hydroxamate-type siderophore. Therefore, it is difficult to explain how the inner membrane transporter FhuB can distinguish between ferrichrome-bound and ferrichrome-free FhuD. In the current study, we have employed a 30 ns molecular dynamics simulation of FhuD with its bound siderophore removed to explore the dynamic behavior of FhuD in the substrate-free state. The MD simulation suggests that FhuD is somewhat dynamic with a C-terminal domain closure of 6 degrees upon release of its siderophore. This relatively large motion suggests differences that would allow FhuB to distinguish between ferrichrome-bound and ferrichrome-free FhuD.

Charged residues in the alpha1 and beta2 pre-M1 regions involved in GABAA receptor activation

For Cys-loop ligand-gated ion channels (LGIC), the protein movements that couple neurotransmitter binding to channel gating are not well known. The pre-M1 region, which links the extracellular agonist-binding domain to the channel-containing transmembrane domain, is in an ideal position to transduce binding site movements to gating movements. A cluster of cationic residues in this region is observed in all LGIC subunits, and in particular, an arginine residue is absolutely conserved. We mutated charged pre-M1 residues in the GABAA receptor alpha1 (K219, R220, K221) and beta2 (K213, K215, R216) subunits to cysteine and expressed the mutant subunits with wild-type beta2 or alpha1 in Xenopus oocytes. Cysteine substitution of beta2R216 abolished channel gating by GABA without altering the binding of the GABA agonist [3H]muscimol, indicating that this residue plays a key role in coupling GABA binding to gating. Tethering thiol-reactive methanethiosulfonate (MTS) reagents onto alpha1K219C, beta2K213C, and beta2K215C increased maximal GABA-activated currents, suggesting that structural perturbations of the pre-M1 regions affect channel gating. GABA altered the rates of sulfhydryl modification of alpha1K219C, beta2K213C, and beta2K215C, indicating that the pre-M1 regions move in response to channel activation. A positively charged MTS reagent modified beta2K213C and beta2K215C significantly faster than a negatively charged reagent, and GABA activation eliminated modification of beta2K215C by the negatively charged reagent. Overall, the data indicate that the pre-M1 region is part of the structural machinery coupling GABA binding to gating and that the transduction of binding site movements to channel movements is mediated, in part, by electrostatic interactions.

Structural mechanism of the recovery stroke in the myosin molecular motor

The power stroke pulling myosin along actin filaments during muscle contraction is achieved by a large rotation ( approximately 60 degrees ) of the myosin lever arm after ATP hydrolysis. Upon binding the next ATP, myosin dissociates from actin, but its ATPase site is still partially open and catalytically off. Myosin must then close and activate its ATPase site while returning the lever arm for the next power stroke. A mechanism for this coupling between the ATPase site and the distant lever arm is determined here by generating a continuous series of optimized intermediates between the crystallographic end-states of the recovery stroke. This yields a detailed structural model for communication between the catalytic and the force-generating regions that is consistent with experimental observations. The coupling is achieved by an amplifying cascade of conformational changes along the relay helix lying between the ATPase and the domain carrying the lever arm.

Nonuniform charge scaling (NUCS): A practical approximation of solvent electrostatic screening in proteins

In molecular mechanics calculations, electrostatic interactions between chemical groups are usually represented by a Coulomb potential between the partial atomic charges of the groups. In aqueous solution these interactions are modified by the polarizable solvent. Although the electrostatic effects of the polarized solvent on the protein are well described by the Poisson--Boltzmann equation, its numerical solution is computationally expensive for large molecules such as proteins. The procedure of nonuniform charge scaling (NUCS) is a pragmatic approach to implicit solvation that approximates the solvent screening effect by individually scaling the partial charges on the explicit atoms of the macromolecule so as to reproduce electrostatic interaction energies obtained from an initial Poisson--Boltzmann analysis. Once the screening factors have been determined for a protein the scaled charges can be easily used in any molecular mechanics program that implements a Coulomb term. The approach is particularly suitable for minimization-based simulations, such as normal mode analysis, certain conformational reaction path or ligand binding techniques for which bulk solvent cannot be included explicitly, and for combined quantum mechanical/molecular mechanical calculations when the interface to more elaborate continuum solvent models is lacking. The method is illustrated using reaction path calculations of the Tyr 35 ring flip in the bovine pancreatic trypsin inhibitor.