Exploration of conformational changes in lactose permease upon sugar binding and proton transfer through coarse‐grained simulations

Escherichia coli lactose permease (LacY) actively transports lactose and other galactosides across cell membranes through lactose/H⁺ symport process. Lactose/H⁺ symport is a highly complex process that involves sugar translocation, H⁺ transfer, and large‐scale protein conformational changes. The complete picture of lactose/H⁺ symport is largely unclear due to the complexity and multiscale nature of the process. In this work, we develop the force field for sugar molecules compatible with PACE, a hybrid and coarse‐grained force field that couples the united‐atom protein models with the coarse‐grained MARTINI water/lipid. After validation, we implement the new force field to investigate the binding of a ‐d ‐galactopyranosyl‐1‐thio‐ ‐d ‐galactopyranoside (TDG) molecule to a wild‐type LacY. Results show that the local interactions between TDG and LacY at the binding pocket are consistent with the X‐ray experiment. Transitions from inward‐facing to outward‐facing conformations upon TDG binding and protonation of Glu269 have been achieved from ～5.5 µs simulations. Both the opening of the periplasmic side and closure of the cytoplasmic side of LacY are consistent with double electron–electron resonance and thiol cross‐linking experiments. Our analysis suggests that the conformational changes of LacY are a cumulative consequence of interdomain H‐bonds breaking at the periplasmic side, interdomain salt‐bridge formation at the cytoplasmic side, and the TDG orientational changes during the transition.     

Coarse‐grained and all‐atom modeling of structural states and transitions in hemoglobin

Hemoglobin (Hb), an oxygen‐binding protein composed of four subunits (α1, α2, β1, and β2), is a well‐known example of allosteric proteins that are capable of cooperative ligand binding. Despite decades of studies, the structural basis of its cooperativity remains controversial. In this study, we have integrated coarse‐grained (CG) modeling, all‐atom simulation, and structural data from X‐ray crystallography and wide‐angle X‐ray scattering (WAXS), aiming to probe dynamic properties of the two structural states of Hb (T and R state) and the transitions between them. First, by analyzing the WAXS data of unliganded and liganded Hb, we have found that the structural ensemble of T or R state is dominated by one crystal structure of Hb with small contributions from other crystal structures of Hb. Second, we have used normal mode analysis to identify two distinct quaternary rotations between the α1β1 and α2β2 dimer, which drive the transitions between T and R state. We have also identified the hot‐spot residues whose mutations are predicted to greatly change these quaternary motions. Third, we have generated a CG transition pathway between T and R state, which predicts a clear order of quaternary and tertiary changes involving α and β subunits in Hb. Fourth, we have used the accelerated molecular dynamics to perform an all‐atom simulation starting from the T state of Hb, and we have observed a transition toward the R state of Hb. Further analysis of crystal structural data and the all‐atom simulation trajectory has corroborated the order of quaternary and tertiary changes predicted by CG modeling. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Cooperativity network of Trp‐cage miniproteins: probing salt‐bridges

Trp‐cage miniprotein was used to investigate the role of a salt‐bridge (Asp⁹–Arg¹⁶) in protein formation, by mutating residues at both sides, we mapped its contribution to overall stability and its role in folding mechanism. We found that both of the above side‐chains are also part of a dense interaction network composed of electrostatic, H‐bonding, hydrophobic, etc. components. To elucidate the fold stabilizing effects, we compared and contrasted electronic circular dichroism and NMR data of miniproteins equipped with a salt‐bridge with those of the salt‐bridge deleted mutants. Data were acquired both in neutral and in acidic aqueous solutions to decipher the pH dependency of both fully and partially charged partners. Our results indicate that the folding of Trp‐cage miniproteins is more complex than a simple two‐state process as we detected an intermediate state that differs significantly from the native fold. The intermediate formation is related to the salt‐bridge stabilization; in the miniprotein variants equipped with salt‐bridge the population of the intermediate state at acidic pH is significantly higher than it is for the salt‐bridge deleted mutants. In this molecular framework Arg¹⁶ stabilizes more than Asp⁹ does, because of its higher degree of 3D‐fold cooperation. In conclusion, the Xxx$^{9} \leftrightarrow$ Yyy¹⁶ salt‐bridge is not an isolated entity of this fold; rather it is an integrated part of a complex interaction network. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Coarse‐grained modeling of the intrinsically disordered protein Histatin 5 in solution: Monte Carlo simulations in combination with SAXS

Monte Carlo simulations and coarse‐grained modeling have been used to analyze Histatin 5, an unstructured short cationic salivary peptide known to have anticandidical properties. The calculated scattering functions have been compared with intensity curves and the distance distribution function P(r) obtained from small angle X‐ray scattering (SAXS), at both high and low salt concentrations. The aim was to achieve a molecular understanding and a physico‐chemical insight of the obtained SAXS results and to gain information of the conformational changes of Histatin 5 due to altering salt content, charge distribution, and net charge. From a modeling perspective, the accuracy of the electrostatic interactions are of special interest. The used coarse‐grained model was based on the primitive model in which charged hard spheres differing in charge and in size represent the ionic particles, and the solvent only enters the model through its relative permittivity. The Hamiltonian of the model comprises three different contributions: (i) excluded volumes, (ii) electrostatic, and (iii) van der Waals interactions. Even though the model can be considered as gross omitting all atomistic details, a great correspondence is obtained with the experimental results. Proteins 2016; 84:777–791. © 2016 Wiley Periodicals, Inc.     

High salt‐treatment‐induced Na+ extrusion and low salt‐treatment‐induced Na+ accumulation in suspension‐cultured cells of the mangrove plant,Bruguiera sexangula

A suspension‐cultured cell strain of the mangrove plant (Bruguiera sexangula) was established from a callus culture and maintained in an amino acid medium in the absence of NaCl. NaCl non‐adapted cells were transferred to media containing 0–200 mm NaCl. The initial growth rate decreased gradually with increasing salt concentrations. However, at up to 150 mm NaCl, cell number growth at the highest point was almost the same as that at lower salt concentrations. Cells even continued to grow in the presence of 200 mm NaCl. Cells incubated in a medium containing 50 mm NaCl for 3 weeks accumulated Na⁺, while those incubated in 150 mm NaCl for 2 d showed only a transient increase in Na⁺ and Cl^– concentrations. In the latter treatment, the intracellular concentration of Na⁺ returned to the original low level within 2 weeks. It took a longer time for Cl^– to return to its original level. As a result, the Na⁺ and Cl^– concentrations in cells cultured with 50 mm NaCl were much larger than those in cells cultured with 150 mm NaCl. The intracellular distribution of ions after transfer to the medium containing 150 mm NaCl was analysed by isolating the vacuoles. Treatment with amiloride, an inhibitor of the Na⁺/H⁺ antiporter, suppressed the recovery of Na⁺ to the original level in the cells. Treatment with 150 mm NaCl for 3 d stimulated the activities of both the vanadate‐dependent H⁺‐ATPase and the Na⁺/H⁺ antiporter in the plasma membrane fraction.     

A coarse‐grained potential for fold recognition and molecular dynamics simulations of proteins

A coarse‐grained potential for protein simulations and fold ranking is presented. The potential is based on a two‐point model of individual amino acids and a specific implementation of hydrogen bonding. Parameters are determined for distance dependent pair interactions, pseudo bonds, angles, and torsions. A scaling factor for a hydrogen bonding term is also determined. Iterative sampling for 4867 proteins reproduces distributions of internal coordinates and distances observed in the Protein Data Bank. The adjustment of the potential and resampling are in the spirit of the generalized ensemble approach. No native structure information (e.g., secondary structure) is used in the calculation of the potential or in the simulation of a particular protein. The potential is subject to two tests as follows: (i) simulations of 956 globular proteins in the neighborhood of their native folds (these proteins were not used in the training set) and (ii) discrimination between native and decoy structures for 2470 proteins with 305,000 decoys and the “Decoys ‘R’ Us” dataset. In the first test, 58% of tested proteins stay within 5 Å from the native fold in Molecular Dynamics simulations of more than 20 nanoseconds using the new potential. The potential is also useful in differentiating between correct and approximate folds providing significant signal for structure prediction algorithms. Sampling with the potential consistently regenerates the distribution of distances and internal coordinates it learned. Nevertheless, during Molecular Dynamics simulations structures are found that reproduce the learned distributions but are far from the native fold. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Ligand-induced conformational changes in the lactose permease of Escherichia coli: evidence for two binding sites.

By using a lactose permease mutant containing a single Cys residue in place of Val 331 (helix X), conformational changes induced by ligand binding were studied. With right-side-out membrane vesicles containing Val 331-->Cys permease, lactose transport is inactivated by either N-ethylmaleimide (NEM) or 7-diethylamino-3-(4'-maleimidylphenyl)-4-methylcoumarin (CPM). Remarkably, beta,D-galactopyranosyl 1-thio-beta,D-galactopyranoside (TDG) enhances the rate of inactivation by CPM, a hydrophobic sulfhydryl reagent, whereas NEM inactivation is attenuated by the ligand. Val 331-->Cys permease was then purified and studied in dodecyl-beta,D-maltoside by site-directed fluorescence spectroscopy. The reactivity of Val 331-->Cys permease with 2-(4'-maleimidylanilino)-naphthalene-6-sulfonic acid (MIANS) is not changed over a low range of TDG concentrations (< 0.8 mM), but the fluorescence of the MIANS-labeled protein is quenched in a saturable manner (apparent Kd approximately equal to 0.12 mM) without a change in emission maximum. In contrast, over a higher range of TDG concentrations (1-10 mM), the reactivity of Val 331-->Cys permease with MIANS is enhanced and the emission maximum of MIANS-labeled permease is blue shifted by 3-7 nm. Furthermore, the fluorescence of MIANS-labeled Val 331 -->Cys permease is quenched by both acrylamide and iodide, but the former is considerably more effective. A low concentration of TDG (0.2 mM) does not alter quenching by either compound, whereas a higher concentration of ligand (10 mM) decreases the quenching constant for iodide by about 50% and for acrylamide by about 20%.(ABSTRACT TRUNCATED AT 250 WORDS)     

Hierarchical domain‐motion analysis of conformational changes in sarcoplasmic reticulum Ca2+‐ATPase

Sarco(endo)plasmic reticulum Ca²⁺‐ATPase transports two Ca²⁺ per ATP‐hydrolyzed across biological membranes against a large concentration gradient by undergoing large conformational changes. Structural studies with X‐ray crystallography revealed functional roles of coupled motions between the cytoplasmic domains and the transmembrane helices in individual reaction steps. Here, we employed “Motion Tree (MT),” a tree diagram that describes a conformational change between two structures, and applied it to representative Ca²⁺‐ATPase structures. MT provides information of coupled rigid‐body motions of the ATPase in individual reaction steps. Fourteen rigid structural units, “common rigid domains (CRDs)” are identified from seven MTs throughout the whole enzymatic reaction cycle. CRDs likely act as not only the structural units, but also the functional units. Some of the functional importance has been newly revealed by the analysis. Stability of each CRD is examined on the morphing trajectories that cover seven conformational transitions. We confirmed that the large conformational changes are realized by the motions only in the flexible regions that connect CRDs. The Ca²⁺‐ATPase efficiently utilizes its intrinsic flexibility and rigidity to response different switches like ligand binding/dissociation or ATP hydrolysis. The analysis detects functional motions without extensive biological knowledge of experts, suggesting its general applicability to domain movements in other membrane proteins to deepen the understanding of protein structure and function. Proteins 2015; 83:746–756. © 2015 Wiley Periodicals, Inc.     

Anchorage‐dependent binding of integrin I‐domain to adhesion ligands

The dynamic interactions between leukocyte integrin receptors and ligands in the vascular endothelium, extracellular matrix, or invading pathogens result in leukocyte adhesion, extravasation, and phagocytosis. This work examined the mechanical strength of the connection between iC3b, a complement component that stimulates phagocytosis, and the ligand‐binding domain, the I‐domain, of integrin α_Mβ₂. Single‐molecule force measurements of α_M I‐domain–iC3b complexes were conducted by atomic force microscope. Strikingly, depending on loading rates, immobilization of the I‐domain via its C‐terminus resulted in a 1.3‐fold to 1.5‐fold increase in unbinding force compared with I‐domains immobilized via the N‐terminus. The force spectra (unbinding force versus loading rate) of the I‐domain–iC3b complexes revealed that the enhanced mechanical strength is due to a 2.4‐fold increase in the lifetime of the I‐domain–iC3b bond. Given the structural and functional similarity of all integrin I‐domains, our result supports the existing allosteric regulatory model by which the ligand binding strength of integrin can be increased rapidly when a force is allowed to stretch the C‐terminus of the I‐domain. This type of mechanism may account for the rapid ligand affinity adjustment during leukocyte migration. Copyright © 2015 John Wiley & Sons, Ltd.     

Structure of HIV‐1 reverse transcriptase/d4TTP complex: Novel DNA cross‐linking site and pH‐dependent conformational changes

Stavudine (d4T, 2′,3′‐didehydro‐2′,3′‐dideoxythymidine) was one of the first chain‐terminating nucleoside analogs used to treat HIV infection. We present the first structure of the active, triphosphate form of d4T (d4TTP) bound to a catalytic complex of HIV‐1 RT/dsDNA template‐primer. We also present a new strategy for disulfide (S–S) chemical cross‐linking between N⁶ of a modified adenine at the second overhang base to I63C in the fingers subdomain of RT. The cross‐link site is upstream of the duplex‐binding region of RT, however, the structure is very similar to published RT structures with cross‐linking to Q258C in the thumb, which suggests that cross‐linking at either site does not appreciably perturb the RT/DNA structures. RT has a catalytic maximum at pH 7.5. We determined the X‐ray structures of the I63C‐RT/dsDNA/d4TTP cross‐linked complexes at pH 7, 7.5, 8, 8.5, 9, and 9.5. We found small (~0.5 Å), pH‐dependent motions of the fingers subdomain that folds in to form the dNTP‐binding pocket. We propose that the pH‐activity profile of RT relates to this motion of the fingers. Due to side effects of neuropathy and lipodystrophy, use of d4T has been stopped in most countries, however, chemical modification of d4T might lead to the development of a new class of nucleoside analogs targeting RNA and DNA polymerases.     

Heparin‐induced amyloid fibrillation of β2‐microglobulin explained by solubility and a supersaturation‐dependent conformational phase diagram

Amyloid fibrils are fibrillar deposits of denatured proteins associated with amyloidosis and are formed by a nucleation and growth mechanism. We revisited an alternative and classical view of amyloid fibrillation: amyloid fibrils are crystal‐like precipitates of denatured proteins formed above solubility upon breaking supersaturation. Various additives accelerate and then inhibit amyloid fibrillation in a concentration‐dependent manner, suggesting that the combined effects of stabilizing and destabilizing forces affect fibrillation. Heparin, a glycosaminoglycan and anticoagulant, is an accelerator of fibrillation for various amyloidogenic proteins. By using β₂‐microglobulin, a protein responsible for dialysis‐related amyloidosis, we herein examined the effects of various concentrations of heparin on fibrillation at pH 2. In contrast to previous studies that focused on accelerating effects, higher concentrations of heparin inhibited fibrillation, and this was accompanied by amorphous aggregation. The two‐step effects of acceleration and inhibition were similar to those observed for various salts. The results indicate that the anion effects caused by sulfate groups are one of the dominant factors influencing heparin‐dependent fibrillation, although the exact structures of fibrils and amorphous aggregates might differ between those formed by simple salts and matrix‐forming heparin. We propose that a conformational phase diagram, accommodating crystal‐like amyloid fibrils and glass‐like amorphous aggregates, is important for understanding the effects of various additives.     

Norovirus RNA‐dependent RNA polymerase: A computational study of metal‐binding preferences

Norovirus (NV) RNA‐dependent RNA polymerase (RdRP) is essential for replicating the genome of the virus, which makes this enzyme a key target for the development of antiviral agents against NV gastroenteritis. In this work, a complex of NV RdRP bound to manganese ions and an RNA primer‐template duplex was investigated using X‐ray crystallography and hybrid quantum chemical/molecular mechanical simulations. Experimentally, the complex crystallized in a tetragonal crystal form. The nature of the primer/template duplex binding in the resulting structure indicates that the complex is a closed back‐tracked state of the enzyme, in which the ‐end of the primer occupies the position expected for the post‐incorporated nucleotide before translocation. Computationally, it is found that the complex can accept a range of divalent metal cations without marked distortions in the active site structure. The highest binding energy is for copper, followed closely by manganese and iron, and then by zinc, nickel, and cobalt. Proteins 2017; 85:1435–1445. © 2017 Wiley Periodicals, Inc.     

Helix stability of oligoglycine,oligoalanine, and oligo‐β‐alanine dodecamers reflected by hydrogen‐bond persistence

Helices are important structural/recognition elements in proteins and peptides. Stability and conformational differences between helices composed of α‐ and β‐amino acids as scaffolds for mimicry of helix recognition has become a theme in medicinal chemistry. Furthermore, helices formed by β‐amino acids are experimentally more stable than those formed by α‐amino acids. This is paradoxical because the larger sizes of the hydrogen‐bonding rings required by the extra methylene groups should lead to entropic destabilization. In this study, molecular dynamics simulations using the second‐generation force field, AMOEBA (Ponder, J.W., et al., Current status of the AMOEBA polarizable force field. J Phys Chem B, 2010. 114 (8): p. 2549–64.) explored the stability and hydrogen‐bonding patterns of capped oligo‐β‐alanine, oligoalanine, and oligoglycine dodecamers in water. The MD simulations showed that oligo‐β‐alanine has strong acceptor+2 hydrogen bonds, but surprisingly did not contain a large content of 3₁₂‐helical structures, possibly due to the sparse distribution of the 3₁₂‐helical structure and other structures with acceptor+2 hydrogen bonds. On the other hand, despite its backbone flexibility, the β‐alanine dodecamer had more stable and persistent <3.0 Å hydrogen bonds. Its structure was dominated more by multicentered hydrogen bonds than either oligoglycine or oligoalanine helices. The 3₁ (PII) helical structure, prevalent in oligoglycine and oligoalanine, does not appear to be stable in oligo‐β‐alanine indicating its competition with other structures (stacking structure as indicated by MD analyses). These differences are among the factors that shape helical structural preferences and the relative stabilities of these three oligopeptides. Proteins 2014; 82:3043–3061. © 2014 Wiley Periodicals, Inc.     

Molecular modeling of the binding modes of the iron‐sulfur protein to the Jac1 co‐chaperone from Saccharomyces cerevisiae by all‐atom and coarse‐grained approaches

The iron‐sulfur protein 1 (Isu1) and the J‐type co‐chaperone Jac1 from yeast are part of a huge ATP‐dependent system, and both interact with Hsp70 chaperones. Interaction of Isu1 and Jac1 is a part of the iron‐sulfur cluster biogenesis system in mitochondria. In this study, the structure and dynamics of the yeast Isu1–Jac1 complex has been modeled. First, the complete structure of Isu1 was obtained by homology modeling using the I‐TASSER server and YASARA software and thereafter tested for stability in the all‐atom force field AMBER. Then, the known experimental structure of Jac1 was adopted to obtain initial models of the Isu1–Jac1 complex by using the ZDOCK server for global and local docking and the AutoDock software for local docking. Three most probable models were subsequently subjected to the coarse‐grained molecular dynamics simulations with the UNRES force field to obtain the final structures of the complex. In the most probable model, Isu1 binds to the left face of the Γ‐shaped Jac1 molecule by the β‐sheet section of Isu1. Residues L₁₀₅, L₁₀₉, and Y₁₆₃ of Jac1 have been assessed by mutation studies to be essential for binding (Ciesielski et al., J Mol Biol 2012; 417:1–12). These residues were also found, by UNRES/molecular dynamics simulations, to be involved in strong interactions between Isu1 and Jac1 in the complex. Moreover, N₉₅, T₉₈, P₁₀₂, H₁₁₂, V₁₅₉, L₁₆₇, and A₁₇₀ of Jac1, not yet tested experimentally, were also found to be important in binding. Proteins 2015; 83:1414–1426. © 2015 Wiley Periodicals, Inc.     

Spike‐dependent calcium influx in dendrites of the cricket giant interneuron

Identified wind‐sensitive giant interneurons in the cricket's cercal sensory system integrate cercal afferent signals and release an avoidance behavior. A calcium‐imaging technique was applied to the giant interneurons to examine the presence of the voltage‐dependent Ca²⁺ channels (VDCCs) in their dendrites. We found that presynaptic stimuli to the cercal sensory nerve cords elevated the cytosolic Ca²⁺ concentration ([Ca²⁺]_i) in the dendrites of the giant interneurons. The dendritic Ca²⁺ rise coincided with the spike burst of the giant interneurons, and the rate of Ca²⁺ rise depended on the frequency of the action potentials. These results suggest that the action potentials directly caused [Ca²⁺]_i increase. Observation of the [Ca²⁺]_i elevation induced by depolarizing current injection demonstrates the presence of the VDCCs in the dendrites. Although hyperpolarizing current injection into the giant interneuron suppressed action potential generation, EPSPs could induce no [Ca²⁺]_i increase. This result means that ligand‐gated channels do not contribute to the synaptically stimulated Ca²⁺ elevation. On the other hand, antidromically stimulated spikes also increased [Ca²⁺]_i in all cellular regions including the dendrites. And bath application of a mixture of Ni²⁺, Co²⁺, and Cd²⁺ or tetrodotoxin inhibited the [Ca²⁺]_i elevation induced by the antidromic stimulation. From these findings, we suppose that the axonal spikes antidromically propagate and induce the Ca²⁺ influx via VDCCs in the dendrites. The spike‐dependent Ca²⁺ elevation may regulate the sensory signals processing via second‐messenger cascades in the giant interneurons. © 2000 John Wiley & Sons, Inc. J Neurobiol 44: 45–56, 2000     

Maize ZmVPP5 is a truncated Vacuole H+‐PPase that confers hypersensitivity to salt stress

In plants, Vacuole H⁺‐PPases (VPPs) are important proton pumps and encoded by multiple genes. In addition to full‐length VPPs, several truncated forms are expressed, but their biological functions are unknown. In this study, we functionally characterized maize vacuole H⁺‐PPase 5 (ZmVPP5), a truncated VPP in the maize genome. Although ZmVPP5 shares high sequence similarity with ZmVPP1, ZmVPP5 lacks the complete structure of the conserved proton transport and the inorganic pyrophosphatase‐related domain. Phylogenetic analysis suggests that ZmVPP5 might be derived from an incomplete gene duplication event. ZmVPP5 is expressed in multiple tissues, and ZmVPP5 was detected in the plasma membrane, vacuole membrane and nuclei of maize cells. The overexpression of ZmVPP5 in yeast cells caused a hypersensitivity to salt stress. Transgenic maize lines with overexpressed ZmVPP5 also exhibited the salt hypersensitivity phenotype. A yeast two‐hybrid analysis identified the ZmBag6‐like protein as a putative ZmVPP5‐interacting protein. The results of bimolecular luminescence complementation (BiLC) assay suggest an interaction between ZmBag6‐like protein and ZmVPP5 in vivo. Overall, this study suggests that ZmVPP5 might act as a VPP antagonist and participate in the cellular response to salt stress. Our study of ZmVPP5 has expanded the understanding of the origin and functions of truncated forms of plant VPPs.     

Salt effects on hydrophobic‐core formation in folding of a helical miniprotein studied by molecular dynamics simulations

We have investigated effects of salt ions on folding events of a helical miniprotein chicken villin headpiece subdomain HP36. Low concentrations of ions alter electrostatic interactions between charged groups of a protein and can change the populations of conformers. Here, we compare two data sets of folding simulations of HP36 in explicit water solvent with or without ions. For efficient sampling of the conformational space of HP36, the multicanonical replica‐exchange molecular dynamics method was employed. Our analyses suggest that salt alters salt‐bridging nature of the protein at later stages of folding at room temperature. Especially, more nonnative, nonlocal salt bridges are formed at near‐native conformations in pure water. Our analyses also show that such salt‐bridge formation hinders the fully native hydrophobic‐core packing at the final stages of folding. Proteins 2014; 82:933–943. © 2013 Wiley Periodicals, Inc.