Model of the 3-D structure of the GLUT3 glucose transporter and molecular dynamics simulation of glucose transport

A molecular model of the three-dimensional (3-D) structure of the glucose transport protein, GLUT3, has been derived by homology modeling. The model was built on the basis of structural data from the MscL protein, which is a mechanosensitive ion channel, and general insights from aquaporin (a water permeation pore). Structurally conserved regions were defined by amino acid sequence comparisons, optimum interconnecting loops were selected from the protein databank, and amino (N)- and carboxy (C)-terminal ends of the protein were generated as random coil structures. The model was then subjected to energy minimization and molecular dynamics simulations in the presence of bound substrate (D-glucose). In the proposed structure of GLUT3, the 12 transmembrane (TM) helices form a right-hand barrel with a central hydrophilic pore. The pore is shaped like a funnel with dimensions of approximately 5-6 A by 8 A at its narrowest point. A network of polar and aromatic amino acids line the pore region and may facilitate the movement of glucose along the channel. A putative binding site for inhibitory ligands, such as forskolin and cytochalasin B, was identified on an intracellular aspect of the protein. Molecular dynamics studies showed that changes in the tilt and flexibility of key TM helices may modulate the opening of the pore to effect glucose transport. The proposed structure of GLUT3 may prove useful in guiding future experiments aimed at more precisely defining various functional regions of the transporter and may encourage efforts to develop models of other complex membrane proteins.     

A Hydrophobic Gold Surface Triggers Misfolding and Aggregation of the Amyloidogenic Josephin Domain in Monomeric Form,While Leaving the Oligomers Unaffected

Protein misfolding and aggregation in intracellular and extracellular spaces is regarded as a main marker of the presence of degenerative disorders such as amyloidoses. To elucidate the mechanisms of protein misfolding, the interaction of proteins with inorganic surfaces is of particular relevance, since surfaces displaying different wettability properties may represent model systems of the cell membrane. Here, we unveil the role of surface hydrophobicity/hydrophilicity in the misfolding of the Josephin domain (JD), a globular-shaped domain of ataxin-3, the protein responsible for the spinocerebellar ataxia type 3. By means of a combined experimental and theoretical approach based on atomic force microscopy, Fourier transform infrared spectroscopy and molecular dynamics simulations, we reveal changes in JD morphology and secondary structure elicited by the interaction with the hydrophobic gold substrate, but not by the hydrophilic mica. Our results demonstrate that the interaction with the gold surface triggers misfolding of the JD when it is in native-like configuration, while no structural modification is observed after the protein has undergone oligomerization. This raises the possibility that biological membranes would be unable to affect amyloid oligomeric structures and toxicity.     

Disulfide bonds convert small heat shock protein Hsp16.3 from a chaperone to a non-chaperone: implications for the evolution of cysteine in molecular chaperones

Molecular chaperones mainly function in assisting newly synthesized polypeptide folding and protect non-native proteins from aggregation, with known structural features such as the ability of spontaneous folding/refolding and high conformational flexibility. In this report, we verified the assumption that the lack of disulfide bonds in molecular chaperones is a prerequisite for such unique structural features. Using small heat shock protein (one sub-class of chaperones) Hsp16.3 as a model system, our results show the following: (1) Cysteine-free Hsp16.3 wild type protein can efficiently exhibit chaperone activity and spontaneously refold/reassemble with high conformational flexibility. (2) Whereas Hsp16.3 G89C mutant with inter-subunit disulfide bonds formed seems to lose the nature of chaperone proteins, i.e., under stress conditions, it neither acts as molecular chaperone nor spontaneously refolds/reassembles. Structural analysis indicated that the mutant exists as an unstable molten globule-like state, which incorrectly exposes hydrophobic surfaces and irreversibly tends to form aggregates that can be suppressed by the other molecular chaperone (alpha-crystallin). By contrast, reduction of disulfide bond in the Hsp16.3 G89C mutant can significantly recover its character as a molecular chaperone. In light of these results, we propose that disulfide bonds could severely disturb the structure/function of molecular chaperones like Hsp16.3. Our results might not only provide insights into understanding the structural basis of chaperone upon binding substrates, but also explain the observation that the occurrence of cysteine in molecular chaperones is much lower than that in other protein families, subsequently being helpful to understand the evolution of protein family.     

Structure, Dynamics, and Ion Conductance of the Phospholamban Pentamer

A 52-residue membrane protein, phospholamban (PLN) is an inhibitor of an adenosine-5′-triphosphate-driven calcium pump, the Ca²⁺-ATPase. Although the inhibition of Ca²⁺-ATPase involves PLN monomers, in a lipid bilayer membrane, PLN monomers form stable pentamers of unknown biological function. The recent NMR structure of a PLN pentamer depicts cytoplasmic helices extending normal to the bilayer in what is known as the bellflower conformation. The structure shows transmembrane helices forming a hydrophobic pore 4 Å in diameter, which is reminiscent of earlier reports of possible ion conductance through PLN pentamers. However, recent FRET measurements suggested an alternative structure for the PLN pentamer, known as the pinwheel model, which features a narrower transmembrane pore and cytoplasmic helices that lie against the bilayer. Here, we report on structural dynamics and conductance properties of the PLN pentamers from all-atom (AA) and coarse-grained (CG) molecular dynamics simulations. Our AA simulations of the bellflower model demonstrate that in a lipid bilayer membrane or a detergent micelle, the cytoplasmic helices undergo large structural fluctuations, whereas the transmembrane pore shrinks and becomes asymmetric. Similar asymmetry of the transmembrane region was observed in the AA simulations of the pinwheel model; the cytoplasmic helices remained in contact with the bilayer. Using the CG approach, structural dynamics of both models were investigated on a microsecond timescale. The cytoplasmic helices of the CG bellflower model were observed to fall against the bilayer, whereas in the CG pinwheel model the conformation of the cytoplasmic helices remained stable. Using steered molecular dynamics simulations, we investigated the feasibility of ion conductance through the pore of the bellflower model. The resulting approximate potentials of mean force indicate that the PLN pentamer is unlikely to function as an ion channel.     

Structural features that predict real-value fluctuations of globular proteins

It is crucial to consider dynamics for understanding the biological function of proteins. We used a large number of molecular dynamics (MD) trajectories of nonhomologous proteins as references and examined static structural features of proteins that are most relevant to fluctuations. We examined correlation of individual structural features with fluctuations and further investigated effective combinations of features for predicting the real value of residue fluctuations using the support vector regression (SVR). It was found that some structural features have higher correlation than crystallographic B‐factors with fluctuations observed in MD trajectories. Moreover, SVR that uses combinations of static structural features showed accurate prediction of fluctuations with an average Pearson's correlation coefficient of 0.669 and a root mean square error of 1.04 Å. This correlation coefficient is higher than the one observed in predictions by the Gaussian network model (GNM). An advantage of the developed method over the GNMs is that the former predicts the real value of fluctuation. The results help improve our understanding of relationships between protein structure and fluctuation. Furthermore, the developed method provides a convienient practial way to predict fluctuations of proteins using easily computed static structural features of proteins. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

Homology Modeling of Adenylosuccinate Synthetase from Saccharomyces Cerevisiae Reveals a Possible Binding Region for Single-Stranded ARS Sequences

Abstract

Adenylosuccinate synthetase from Saccharomyces cerevisiae was investigated in order to find a structural explanation for its ability to bind specifically to single-stranded ARS elements (autonomously replicating sequences). Using the E. coli enzyme as template, a model for the structure of adenylosuccinate synthetase from S. cerevisiae was generated and subsequently refined by molecular dynamics techniques.

The resulting three-dimensional structure offers an explanation for the DNA binding activity of the yeast enzyme by revealing a distinct basic region that is not present in the homologous enzymes from other organisms.

The model is also in good agreement with biochemical data available for a mutant protein in which Glycine 252 is replaced by Aspartate. On the basis of the model a significant structural distortion near the catalytic center was predicted for this mutant, corresponding well to the enzymatic inactivity observed. The mutant enzyme shows larger structural fluctuations than the wild-type protein according to the results of two independent molecular dynamics simulations.     

Cell adhesion and proliferation on poly(N-vinylacetamide) hydrogels and double network approaches for changing cellular affinities

Poly(N-vinylacetamide) hydrogels (PNVA gels) were synthesized to investigate their basic characteristics for biomedical applications such as water contact angles, protein uptake, and mouse fibroblasts (L-929) cell adhesion. Because PNVA gels show hydrophilic features, double network (DN) hydrogels were prepared by the secondary polymerization of N-vinylacetamide (NVA) or acrylamide (AAm) in PNVA gels (NVA/NVA DN gels and NVA/AAm DN gels, respectively), in order to vary PNVA gel features for biocompatibility. Contact angles for both DN gels decreased to around 20 degrees, whereas both PNVA and PAAm gels were over 30 degrees. On the other hand, more protein tended to adsorb to DN gels than single network hydrogels. Compared to PNVA gel, cell adhesion and proliferation on NVA/NVA DN gel were improved with less swelling ratio and much protein uptake, while no significant difference was observed on NVA/AAm DN gel, probably due to more hydrophilic character, supported by lowest water contact angle. These complicated structure change in DN gels would provide a new methodology for tuning the biocompatibility of hydrogels and for controlling surface hydrophilic characteristics and network structures.     

Does there exist an intrinsic relationship between the flexibility and self-assembly of pepfactants?

Short peptide surfactants (pepfactants) bearing the amphiphilic molecular architecture of hydrophobic and hydrophilic moieties are attractive for a wide range of biological and medical applications. Understanding of structural basis and molecular mechanism underlying the self-assembling behaviour of pepfactants is fundamentally important for rationally designing surfactant-like peptides to function as diverse biomaterials. To date, however, the relationship between the physicochemical properties and self-assembly of pepfactants still remains largely unexplored. In this study, we attempt to elucidate the role of structural flexibility and dynamics in peptide self-assembly. Two fast and reliable quantitative structure–property relationship predictors are carefully developed with the sophisticated genetic algorithm/partial least squares method; the predictors are then applied to estimate molecular flexibility and self-assembling ability for 10,000 randomly generated surfactant-like peptides. As a result, a significant negative correlation between the flexibility and self-assembly is observed, which can be fitted fairly well using an exponential curve. Furthermore, atomistic molecular dynamics simulations also reveal a noticeable difference of flexibility profile between strong and weak self-assembling peptides. All these come together to suggest that the self-assembling behaviour of pepfactants is highly dependent on their structural flexibility.     

Molecular modeling study of the norfloxacin-DNA complex

Molecular modeling and molecular dynamics were performed to investigate the interaction of norfloxacin with the DNA oligonucleotide 5'-d(ATACGTAT)(2). Eight quinolone-DNA binding structures were built by molecular modeling on the basis of experimental results. A 100ps molecular dynamics calculation was carried out on two groove binding models and six partially intercalating models. The resulting average structures were compared with each other and to free DNA structure as a reference. The favorable binding mode of norfloxacin to a DNA substrate was pursued by structural assess including steric hindrance, presence of hydrogen-bonding, non-bonding energies of the complex and presence of abnormal structural distortion. Although two of the intercalative models showed the highest binding energy and the lowest non-bonding interaction energy, they presented structural features which contrast with experimental results. On the other hand, one groove binding model demonstrated the most acceptable structure when the experimental observation was accounted. In this model, hydrogen bonding of the carbonyl and carboxyl group of the norfloxacin rings with the DNA bases was present, and norfloxacin binds to the amine group of the guanine base which protrudes toward the minor groove of B-DNA.     

Molecular modelling of membrane sterols with the use of the GROMOS 96 force field

Membrane located sterols determine the structure and function of eucariotic cell membranes. Moreover, they are targets for important antifungal antibiotic amphotericin B. Knowledge about the geometry and dynamics of sterols in the environment of lipidic membranes is necessary to understand their functions. However, due to the dynamic character of the membrane, no experimental data about sterol behaviour on the molecular level is available. Hence molecular modelling simulations could be a source of useful information. The main goal of this paper is to prove the adequacy of the GROMOS 96 force field for molecular simulations of membrane sterols. We focused our attention on the reproduction of characteristic geometrical features observed in the crystal of cholesterol hemiethanolate by molecular dynamics simulations. The results presented clearly indicate that the GROMOS 96 force field can be a useful tool to simulate the highly lipophilic systems. Moreover, interactions responsible for the stability of such systems can also be recognised.     

Membrane protein dynamics in different environments: simulation study of the outer membrane protein X in a lipid bilayer and in a micelle

The bacterial outer membrane protein OmpX from Escherichia coli has been investigated by molecular dynamics simulations when embedded in a phospholipid bilayer and as a protein-micelle aggregate. The resulting simulation trajectories were analysed in terms of structural and dynamic properties of the membrane protein. In agreement with experimental observations, highest relative stability was found for the β-barrel region that is embedded in the lipophilic phase, whereas an extracellular protruding β-sheet, which is a unique structural feature of OmpX that supposedly plays an important role in cell adhesion and invasion, shows larger structure fluctuations. Additionally, we investigated water permeation into the core of the β-barrel protein, which contains a tight salt-bridge and hydrogen-bond network, so that extensive water flux is unlikely. Differences between the bilayer and the micellar system were observed in the length of the barrel and its position inside the lipid environment, and in the protein interactions with the hydrophilic part of the lipids near the lipid/water interface. Those variations suggest that micelles and other detergent environments might not offer a wholly membrane-like milieu to promote adoption of the physiological conformational state by OmpX.     

A strategic fluorescence labeling of D-galactose/D-glucose-binding protein from Escherichia coli helps to shed light on the protein structural stability and dynamics

The D-glucose/D-galactose-binding protein (GGBP) of Escherichia coli serves as an initial component for both chemotaxis toward D-galactose and D-glucose and high-affinity active transport of the two sugars. GGBP is a monomer with a molecular weight of about 32 kDa that binds glucose with micromolar affinity. The sugar-binding site is located in the cleft between the two lobes of the bilobate protein. In this work, the local and global structural features of GGBP were investigated by a strategic fluorescence labeling procedure and spectroscopic methodologies. A mutant form of GGBP containing the amino acid substitution Met to Cys at position 182 was realized and fluorescently labeled to probe the effect of glucose binding on the local and overall structural organization of the protein. The labeling of the N-terminus with a fluorescence probe as well as the protein intrinsic fluorescence were also used to obtain a complete picture of the GGBP structure and dynamics. Our results showed that the binding of glucose to GGBP resulted in no stabilizing effect on the N-terminus portion of GGBP and in a moderate stabilization of the protein matrix in the vicinity of the ligand-binding site. On the contrary, it was observed that the binding of glucose has a strong stabilization effect on the C-terminal domain of the GGBP structure.     

Hydrophilic Linkers and Polar Contacts Affect Aggregation of FG Repeat Peptides

Transport of large proteins into the nucleus involves two events, binding of the cargo protein to a transport receptor in the cytoplasm and passage of the cargo-transporter complex through the selective permeability barrier of the nuclear pore complex. The permeability barrier is formed by largely disordered polypeptides, each containing a number of conserved hydrophobic phenylalanine-glycine (FG) sequence motifs, connected by hydrophilic linkers of varying sequence (FG nups). How the motifs interact to form the permeability barrier, however, is not yet known. We have, therefore, carried out molecular dynamics simulations on various model FG repeat peptides to study the aggregation propensity of FG nups and the specific roles of the hydrophobic FG motifs and the hydrophilic linkers. Our simulations show spontaneous aggregation of the model nups into hydrated aggregates, which exhibit structural features assumed to be part of the permeability barrier. Our simulations suggest that short β-sheets are an important structural feature of the aggregates and that Phe residues are sufficiently exposed to allow rapid binding of transport receptors. A surprisingly large influence of the amino acid composition of the hydrophilic linkers on aggregation is seen, as well as a major contribution of hydrogen-bonding patterns.     

Molecular Modeling Study of the Norfloxacin-DNA Complex

Abstract

Molecular modeling and molecular dynamics were performed to investigate the interaction of norfloxacin with the DNA oligonucleotide 5′-d(ATACGTAT)₂. Eight quinolone-DNA binding structures were built by molecular modeling on the basis of experimental results. A 100ps molecular dynamics calculation was carried out on two groove binding models and six partially intercalating models. The resulting average structures were compared with each other and to free DNA structure as a reference. The favorable binding mode of norfloxacin to a DNA substrate was pursued by structural assess including steric hindrance, presence of hydrogen-bonding, non-bonding energies of the complex and presence of abnormal structural distortion. Although two of the intercalative models showed the highest binding energy and the lowest non-bonding interaction energy, they presented structural features which contrast with experimental results. On the other hand, one groove binding model demonstrated the most acceptable structure when the experimental observation was accounted. In this model, hydrogen bonding of the carbonyl and carboxyl group of the norfloxacin rings with the DNA bases was present, and norfloxacin binds to the amine group of the guanine base which protrudes toward the minor groove of B-DNA.     

Homology modeling of a heme protein, lignin peroxidase, from the crystal structure of cytochrome c peroxidase.

A 3-dimensional model of lignin peroxidase (LiP) was constructed based on its sequence homology with other peroxidases, particularly cytochrome c peroxidase, the only protein with a known crystal structure in the peroxidase family. The construction of initial conformations of insertions and deletions was assisted by secondary structure predictions, amphipathic helix predictions, and consideration of the specific protein environment. A succession of molecular dynamics simulations of these regions with surrounding residues as constraints were carried out to relax the bond lengths and angles. Full protein molecular dynamics simulations with explicit consideration of bound waters were performed to relax the geometry and to identify dynamically flexible regions of the successive models for further refinement. Among the important functionally relevant structural features predicted are: (i) four disulfide bonds are predicted to be formed between Cys3 and Cys15, Cys14 and Cys285, Cys34 and Cys120 and Cys249 and Cys317; (ii) a glycosylation site, Asn257, was located on the surface; (iii) Glu40 was predicted to form a salt bridge with Arg43 on the distal side of the heme and was considered as a possible origin for the pH dependence of compound I formation; and (iv) two candidate substrate binding sites with a cluster of surface aromatic residues and flexible backbones were found in the refined model, consistent with the nature of known substrates of LiP. Based on these predicted structural features of the model, further theoretical and experimental studies are proposed to continue to elucidate the structure and function of LiP.     

Modeling of the metallo-beta-lactamase from B. fragilis: structural and dynamic effects of inhibitor binding

The structure and dynamics of an inhibitor-bound complex of the metallo-beta-lactamase from Bacteroides fragilis are studied by using molecular dynamics. A search of the conformational space was performed to obtain three distinct models of the complex, which were then subjected to solvated molecular dynamics. A solvated molecular dynamics study of the apo protein was performed to serve as a baseline for comparison with the bound simulations. We find loop conformation changes due to binding as well as a decrease in flexibility of the protein as a whole and especially in the major loop of the beta-lactamase. We report the structural and dynamical features of the inhibitor-bound and apo models, as well as experimentally measurable quantities, which should be capable of distinguishing the two binding modes we have determined.     

Monte Carlo studies on equilibrium globular protein folding. II. Beta-barrel globular protein models

In the context of dynamic Monte Carlo simulations on a model protein confined to a tetrahedral lattice, the interplay of protein size and tertiary structure, and the requirements for an all-or-none transition to a unique native state, are investigated. Small model proteins having a primary sequence consisting of a central bend neutral region flanked by two tails having an alternating hydrophobic/hydrophilic pattern of residues are seen to undergo a continuous transition to a beta-hairpin collapsed state. On increasing the length of the tails, the beta-hairpin structural motif is found to be in equilibrium with a four-member beta-barrel. Further increase of the tail length results in the shift of the structural equilibrium to the four-member beta-barrel. The random coil to beta-barrel transition is of an all-or-none character, but while the central turn is always the desired native bend, the location of the turns involving the two external strands is variable. That is, beta-barrels having the external stands that are two residues out of register are also observed in the transition region. Introduction into the primary sequence of two additional regions that are at the very least neutral toward turn formation produces an all-or-none transition to the unique, native, four-member beta-barrel. Various factors that can augment the stability of the native conformation are explored. Overall, these folding simulations strongly indicate that the general rules of globular protein folding are rather robust--namely, one requires a general pattern of hydrophobic/hydrophilic residues that allow the protein to have a well-defined interior and exterior and the presence of regions in the amino acid sequence that at the very least are locally indifferent to turn formation. Since no site-specific interactions between hydrophobic and hydrophilic residues are required to produce a unique four-member beta-barrel, these simulations strongly suggest that site specificity is involved in structural fine-tuning.     

Structural insights into human GPCR protein OA1: a computational perspective

Human ocular albinism type 1 protein (OA1)—a member of the G-protein coupled receptor (GPCR) superfamily—is an integral membrane glycoprotein expressed exclusively by intracellular organelles known as melanocytes, and is responsible for the proper biogenesis of melanosomes. Mutations in the Oa1 gene are responsible for the disease ocular albinism. Despite its clinical importance, there is a lack of in-depth understanding of its structure and mechanism of activation due to the absence of a crystal structure. In the present study, homology modeling was applied to predicting OA1 structure following thorough sequence analysis and secondary structure predictions. The predicted model had the signature residues and motifs expected of GPCRs, and was used for carrying out molecular docking studies with an endogenous ligand, l-DOPA and an antagonist, dopamine; the results agreed quite well with the available experimental data. Finally, three sets of explicit molecular dynamics simulations were carried out in lipid bilayer, the results of which not only confirmed the stability of the predicted model, but also helped witness some differences in structural features such as rotamer toggle switch, helical tilts and hydrogen bonding pattern that helped distinguish between the agonist- and antagonist-bound receptor forms. In place of the typical “D/ERY”-motif-mediated “ionic lock”, a hydrogen bond mediated by the “DAY” motif was observed that could be used to distinguish the agonist and antagonist bound forms of OA1. In the absence of a crystal structure, this study helped to shed some light on the structural features of OA1, and its behavior in the presence of an agonist and an antagonist, which might be helpful in the future drug discovery process for ocular albinism.     

Interpretation of HDX Data by Maximum-Entropy Reweighting of Simulated Structural Ensembles

Hydrogen-deuterium exchange combined with mass spectrometry (HDX-MS) is a widely applied biophysical technique that probes the structure and dynamics of biomolecules without the need for site-directed modifications or bio-orthogonal labels. The mechanistic interpretation of HDX data, however, is often qualitative and subjective, owing to a lack of quantitative methods to rigorously translate observed deuteration levels into atomistic structural information. To help address this problem, we have developed a methodology to generate structural ensembles that faithfully reproduce HDX-MS measurements. In this approach, an ensemble of protein conformations is first generated, typically using molecular dynamics simulations. A maximum-entropy bias is then applied post hoc to the resulting ensemble such that averaged peptide-deuteration levels, as predicted by an empirical model, agree with target values within a given level of uncertainty. We evaluate this approach, referred to as HDX ensemble reweighting (HDXer), for artificial target data reflecting the two major conformational states of a binding protein. We demonstrate that the information provided by HDX-MS experiments and by the model of exchange are sufficient to recover correctly weighted structural ensembles from simulations, even when the relevant conformations are rarely observed. Degrading the information content of the target data—e.g., by reducing sequence coverage, by averaging exchange levels over longer peptide segments, or by incorporating different sources of uncertainty—reduces the structural accuracy of the reweighted ensemble but still allows for useful insights into the distinctive structural features reflected by the target data. Finally, we describe a quantitative metric to rank candidate structural ensembles according to their correspondence with target data and illustrate the use of HDXer to describe changes in the conformational ensemble of the membrane protein LeuT. In summary, HDXer is designed to facilitate objective structural interpretations of HDX-MS data and to inform experimental approaches and further developments of theoretical exchange models.