Simulating proteins at constant pH: An approach combining molecular dynamics and Monte Carlo simulation

For the structure and function of proteins, the pH of the solution is one of the determining parameters. Current molecular dynamics (MD) simulations account for the solution pH only in a limited way by keeping each titratable site in a chosen protonation state. We present an algorithm that generates trajectories at a Boltzmann distributed ensemble of protonation states by a combination of MD and Monte Carlo (MC) simulation. The algorithm is useful for pH-dependent structural studies and to investigate in detail the titration behavior of proteins. The method is tested on the acidic residues of the protein hen egg white lysozyme. It is shown that small structural changes may have a big effect on the pK(A) values of titratable residues.     

The pH dependence of flavivirus envelope protein structure: insights from molecular dynamics simulations

The flavivirus membrane fusion is triggered by the acid pH of the endosomes after virus endocytosis. The proposed mechanism involves changes in the protonation state of conserved histidine residues of the E protein present in the viral surface that undergoes a series of structural rearrangements that result in the fusion between the endosome and viral bilayers. We studied the pH dependence of E protein rearrangements of dengue virus type 2, used as a model, in the pH range experimented by the virus along the fusion process. We employed a low computational cost scheme to explore the behavior of the E protein by molecular dynamics (MD) simulations of complete systems that include the protein, the solvent, and ions. The procedure alternates cyclically the update of the ionization states of the protein residues with common MD steps applied to the new ionization configuration. Important pH-dependent protein structure rearrangements consistent with the changes of the protonation states of conserved histidine residues were observed. The involvement of other conserved residues in the flavivirus in the rearrangements was also identified. The results show interesting correlations with a proposed model for the fusion mechanism, as well as the experimentally identified key residues, contributing to a better understanding of the structural changes in protein E that lead to the fusion process.     

Characterizing a Histidine Switch Controlling pH-Dependent Conformational Changes of the Influenza Virus Hemagglutinin

During the fusion of the influenza virus to the host cell, bending of the HA2 chain of hemagglutinin into a hairpin-shaped structure in a pH-dependent manner facilitates the fusion of the viral envelope and the endosomal membrane. To characterize the structural and dynamical responses of the hinge region of HA2 to pH changes and examine the role of a conserved histidine in this region (the hinge histidine), we have performed an extensive set of molecular dynamics (MD) simulations of 26-residue peptides encompassing the hinge regions of several hemagglutinin subtypes under both neutral and low pH conditions, modeled by the change of the protonation state of the hinge histidine. More than 70 sets of MD simulations (collectively amounting to 25.1 μs) were performed in both implicit and explicit solvents to study the effect of histidine protonation on structural dynamics of the hinge region. In both explicit and implicit solvent simulations, hinge bending was consistently observed upon the protonation of the histidine in all the simulations starting with an initial straight helical conformation, whereas the systems with a neutral histidine retained their primarily straight conformation throughout the simulations. Conversely, the MD simulations starting from an initially bent conformation resulted in the formation of a straight helical structure upon the neutralization of the hinge histidine, whereas the bent structure was maintained when the hinge histidine remained protonated. Finally, mutation of the hinge histidine to alanine abolishes the bending response of the peptide altogether. A molecular mechanism based on the interaction of the hinge histidine with neighboring acidic residues is proposed to be responsible for its role in controlling the conformation of the hinge. We propose that this might present a common mechanism for pH-controlled structural changes in helical structures when histidines act as the pH sensor.     

Characterizing a Histidine Switch Controlling pH-Dependent Conformational Changes of the Influenza Virus Hemagglutinin

During the fusion of the influenza virus to the host cell, bending of the HA2 chain of hemagglutinin into a hairpin-shaped structure in a pH-dependent manner facilitates the fusion of the viral envelope and the endosomal membrane. To characterize the structural and dynamical responses of the hinge region of HA2 to pH changes and examine the role of a conserved histidine in this region (the hinge histidine), we have performed an extensive set of molecular dynamics (MD) simulations of 26-residue peptides encompassing the hinge regions of several hemagglutinin subtypes under both neutral and low pH conditions, modeled by the change of the protonation state of the hinge histidine. More than 70 sets of MD simulations (collectively amounting to 25.1 μs) were performed in both implicit and explicit solvents to study the effect of histidine protonation on structural dynamics of the hinge region. In both explicit and implicit solvent simulations, hinge bending was consistently observed upon the protonation of the histidine in all the simulations starting with an initial straight helical conformation, whereas the systems with a neutral histidine retained their primarily straight conformation throughout the simulations. Conversely, the MD simulations starting from an initially bent conformation resulted in the formation of a straight helical structure upon the neutralization of the hinge histidine, whereas the bent structure was maintained when the hinge histidine remained protonated. Finally, mutation of the hinge histidine to alanine abolishes the bending response of the peptide altogether. A molecular mechanism based on the interaction of the hinge histidine with neighboring acidic residues is proposed to be responsible for its role in controlling the conformation of the hinge. We propose that this might present a common mechanism for pH-controlled structural changes in helical structures when histidines act as the pH sensor.     

Protonation states of important acidic residues in the central Ca²⁺ ion binding sites of the Ca²⁺-ATPase: a molecular modeling study

The P-type ATPases are responsible for the transport of cations across cell membranes. The sarco(endo)plasmic reticulum Ca2?-ATPase (SERCA) transports two Ca2? ions from the cytoplasm to the lumen of the sarco(endo)plasmic reticulum and countertransports two or three protons per catalytic cycle. Two binding sites for Ca2? ions have been located via protein crystallography, including four acidic amino acid residues that are essential to the ion coordination. In this study, we present molecular dynamics (MD) simulations examining the protonation states of these amino acid residues in a Ca2?-free conformation of SERCA. Such knowledge will be important for an improved understanding of atomistic details of the transport mechanism of protons and Ca2? ions. Eight combinations of the protonation of four central acidic residues, Glu309, Glu771, Asp800, and Glu908, are tested from 10 ns MD simulations with respect to protein stability and ability to maintain a structure similar to the crystal structure. The trajectories for the most prospective combinations of protonation states were elongated to 50 ns and subjected to more detailed analysis, including prediction of pK(a) values of the four acidic residues over the trajectories. From the simulations we find that the combination leaving only Asp800 as charged is most likely. The results are compared to available experimental data and explain the observed destabilization upon full deprotonation, resulting in the entry of cytoplasmic K? ions into the Ca2? binding sites during the simulation in which Ca2? ions are absent. Furthermore, a hypothesis for the exchange of protons from the central binding cavity is proposed.     

Insight into the structural requirements of aminopyrimidine derivatives for good potency against both purified enzyme and whole cells of M. tuberculosis: combination of HQSAR,CoMSIA, and MD simulation studies

The Mycobacterium tuberculosis protein kinase B (PknB) is critical for growth and survival of M. tuberculosis within the host. The series of aminopyrimidine derivatives show impressive activity against PknB (IC₅₀ < .5 μM). However, most of them show weak or no cellular activity against M. tuberculosis (MIC > 63 μM). Consequently, the key structural features related to activity against of both PknB and M. tuberculosis need to be investigated. Here, two- and three-dimensional quantitative structure–activity relationship (2D and 3D QSAR) analyses combined with molecular dynamics (MD) simulations were employed with the aim to evaluate these key structural features of aminopyrimidine derivatives. Hologram quantitative structure–activity relationship (HQSAR) and CoMSIA models constructed from IC₅₀ and MIC values of aminopyrimidine compounds could establish the structural requirements for better activity against of both PknB and M. tuberculosis. The NH linker and the R₁ substituent of the template compound are not only crucial for the biological activity against PknB but also for the biological activity against M. tuberculosis. Moreover, the results obtained from MD simulations show that these moieties are the key fragments for binding of aminopyrimidine compounds in PknB. The combination of QSAR analysis and MD simulations helps us to provide a structural concept that could guide future design of PknB inhibitors with improved potency against both the purified enzyme and whole M. tuberculosis cells.     

Resuscitation Promoting Factors: a Family of Microbial Proteins in Survival and Resuscitation of Dormant Mycobacteria

Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), is an extraordinarily successful pathogen of humankind. It has been estimated that up to one-third of the world’s population is infected with M. tuberculosis, and this population is an important reservoir for disease reactivation. Resuscitation promoting factor (Rpf) is a secretory protein, which was first reported in Micrococcus luteus. There are five functionally redundant Rpf-like proteins found in M. tuberculosis. Rpf promotes the resuscitation of dormant bacilli to yield normal, viable colony forming bacteria. All Rpfs share a conserved domain of about 70 amino acids and possess a lysozyme-like activity. The structural studies of the conserved domain suggest that Rpfs could be considered as a c-type lysozyme and lytic transglycosylases. Recently a novel class of nitrophenylthiocyanates (NPT) inhibitors of the muralytic activity of Rpf were reported which opens a new approach in the study of cell-wall hydrolyzing enzymes. This review describes molecular and structural studies conducted on Rpf proteins, their role in the resuscitation of dormant bacteria, in the reactivation of latent infection and identification of low molecular weight inhibitors of resuscitation promoting factors.     

Prediction of titration properties of structures of a protein derived from molecular dynamics trajectories.

This paper explores the dependence of the molecular dynamics (MD) trajectory of a protein molecule on the titration state assigned to the molecule. Four 100-ps MD trajectories of bovine pancreatic trypsin inhibitor (BPTI) were generated, starting from two different structures, each of which was held in two different charge states. The two starting structures were the X-ray crystal structure and one of the solution structures determined by NMR, and the charge states differed only in the ionization state of N terminus. Although it is evident that the MD simulations were too short to sample fully the equilibrium distribution of structures in each case, standard Poisson-Boltzmann titration state analysis of the resulting configurations shows general agreement between the overall titration behavior of the protein and the charge state assumed during MD simulation: at pH 7, the total net charge of the protein resulting from the titration analysis is consistently lower for the protein with the N terminus assumed to be neutral than for the protein with the N terminus assumed to be charged. For most of the ionizable residues, the differences in the calculated pKaS among the four trajectories are statistically negligible and remain in good agreement with the data obtained by crystal structure titration and by experiment. The exceptions include the N terminus, which responds directly to the change of its imposed charge; the C terminus, which in the NMR structure interacts strongly with the former; and a few other residues (Arg 1, Glu 7, Tyr 35, and Arg 42) whose pKaS reflect the initial structure and the limited trajectory lengths. This study illustrates the importance of the careful assignment of protonation states at the start of MD simulations and points to the need for simulation methods that allow for the variation of the protonation state in the calculation of equilibrium properties.     

Identification of a conserved 8 aa insert in the PIP5K protein in the Saccharomycetaceae family of fungi and the molecular dynamics simulations and structural analysis to investigate its potential functional role

Homologs of the phosphatidylinositol‐4‐phosphate‐5‐kinase (PIP5K), which controls a multitude of essential cellular functions, contain a 8 aa insert in a conserved region that is specific for the Saccharomycetaceae family of fungi. Using structures of human PIP4K proteins as templates, structural models were generated of the Saccharomyces cerevisiae and human PIP5K proteins. In the modeled S. cerevisiae PIP5K, the 8 aa insert forms a surface exposed loop, present on the same face of the protein as the activation loop of the kinase domain. Electrostatic potential analysis indicates that the residues from 8 aa conserved loop form a highly positively charged surface patch, which through electrostatic interaction with the anionic portions of phospholipid head groups, is expected to play a role in the membrane interaction of the yeast PIP5K. To unravel this prediction, molecular dynamics (MD) simulations were carried out to examine the binding interaction of PIP5K, either containing or lacking the conserved signature insert, with two different membrane lipid bilayers. The results from MD studies provide insights concerning the mechanistic of interaction of PIP5K with lipid bilayer, and support the contention that the identified 8 aa conserved insert in fungal PIP5K plays an important role in the binding of this protein with membrane surface. Proteins 2017; 85:1454–1467. © 2017 Wiley Periodicals, Inc.     

Molecular dynamics simulations reveal structural differences among wild-type NPC1 protein and its mutant forms

Abstract

NPC1 is a 25-exon gene located on the long arm of chromosome 18q11.2 and encodes NPC1, a transmembrane protein comprising 1278 amino acid residues. Mutations in the NPC1 gene can cause Niemann-Pick disease type C (NP-C), a rare autosomal-recessive neurovisceral disease. We assessed mutant protein folding using computer-based molecular dynamics (MD) simulations and molecular docking of the three most common NPC1 mutations, all of which result in changes in a cysteine-rich luminal loop region of the protein: a) I1061T is the most commonly detected variant in patients with NP-C worldwide; b) P1007A is the second most common variant, frequently detected in Portuguese, British and German patients; c) G992W occurs most often in patients of Acadian descent. Analyses of molecular structural information and related cellular physiological processes revealed that mutant NPC1 proteins exhibited altered function despite being far from the N-terminal domain cholesterol binding. MD simulations revealed that mutant I1061T protein shows remarkable instability in comparison the WT and also de other mutants, and interestingly this mutant has been identified as the most common variant. In the case of the mutant P1007A, it is presumed that this substitution promotes larger structural changes than proline due to their greater hydrophobic properties.

Structural changes related to the G992W mutation may affect the physicochemical space of G992W variant protein because tryptophan induces hydrophobic interactions. Cholesterol docking studies focused on binding recognition showed differences in the binding positions of variants versus the wild-type protein that go some way to explaining the molecular pathogenesis.

Communicated by Ramaswamy H. Sarma     

Protonation state of the selectivity filter of bacterial voltage-gated sodium channels is modulated by ions

The selectivity filter (SF) of bacterial voltage-gated sodium channels consists of four glutamate residues arranged in a C₄ symmetry. The protonation state population of this tetrad is unclear. To address this question, we simulate the pore domain of bacterial voltage-gated sodium channel of Magnetococcus sp. (Na_vMs) through constant pH methodology in explicit solvent and free energy perturbation calculations. We find that at physiological pH the fully deprotonated as well as singly and doubly protonated states of the SF appear feasible, and that the calculated pKa decreases with each additional bound ion, suggesting that a decrease in the number of ions in the pore can lead to protonation of the SF. Previous molecular dynamics simulations have suggested that protonation can lead to a decrease in the conductance, but no pKa calculations were performed. We confirm a decreased ionic population of the pore with protonation, and also observe structural symmetry breaking triggered by protonation; the SF of the deprotonated channel is closest to the C₄ symmetry observed in crystal structures of the open state, while the SF of protonated states display greater levels of asymmetry which could lead to transition to the inactivated state which possesses a C₂ symmetry in the crystal structure. We speculate that the decrease in the number of ions near the mouth of the channel, due to either random fluctuations or ion depletion due to conduction, could be a self-regulatory mechanism resulting in a nonconducting state that functionally resembles inactivated states.     

Structural plasticity and thermal stability of the histone-like protein from Spiroplasma melliferum are due to phenylalanine insertions into the conservative scaffold

The histone-like (HU) protein is one of the major nucleoid-associated proteins of the bacterial nucleoid, which shares high sequence and structural similarity with IHF but differs from the latter in DNA-specificity. Here, we perform an analysis of structural-dynamic properties of HU protein from Spiroplasma melliferum and compare its behavior in solution to that of another mycoplasmal HU from Mycoplasma gallisepticum. The high-resolution heteronuclear NMR spectroscopy was coupled with molecular-dynamics study and comparative analysis of thermal denaturation of both mycoplasmal HU proteins. We suggest that stacking interactions in two aromatic clusters in the HUSpm dimeric interface determine not only high thermal stability of the protein, but also its structural plasticity experimentally observed as slow conformational exchange. One of these two centers of stacking interactions is highly conserved among the known HU and IHF proteins. Second aromatic core described recently in IHFs and IHF-like proteins is considered as a discriminating feature of IHFs. We performed an electromobility shift assay to confirm high affinities of HUSpm to both normal and distorted dsDNA, which are the characteristics of HU protein. MD simulations of HUSpm with alanine mutations of the residues forming the non-conserved aromatic cluster demonstrate its role in dimer stabilization, as both partial and complete distortion of the cluster enhances local flexibility of HUSpm.     

Distance restraints from crosslinking mass spectrometry: Mining a molecular dynamics simulation database to evaluate lysine–lysine distances

Integrative structural biology attempts to model the structures of protein complexes that are challenging or intractable by classical structural methods (due to size, dynamics, or heterogeneity) by combining computational structural modeling with data from experimental methods. One such experimental method is chemical crosslinking mass spectrometry (XL‐MS), in which protein complexes are crosslinked and characterized using liquid chromatography‐mass spectrometry to pinpoint specific amino acid residues in close structural proximity. The commonly used lysine‐reactive N‐hydroxysuccinimide ester reagents disuccinimidylsuberate (DSS) and bis(sulfosuccinimidyl)suberate (BS³) have a linker arm that is 11.4 Å long when fully extended, allowing C_α (alpha carbon of protein backbone) atoms of crosslinked lysine residues to be up to ～24 Å apart. However, XL‐MS studies on proteins of known structure frequently report crosslinks that exceed this distance. Typically, a tolerance of ～3 Å is added to the theoretical maximum to account for this observation, with limited justification for the chosen value. We used the Dynameomics database, a repository of high‐quality molecular dynamics simulations of 807 proteins representative of diverse protein folds, to investigate the relationship between lysine–lysine distances in experimental starting structures and in simulation ensembles. We conclude that for DSS/BS³, a distance constraint of 26–30 Å between C_α atoms is appropriate. This analysis provides a theoretical basis for the widespread practice of adding a tolerance to the crosslinker length when comparing XL‐MS results to structures or in modeling. We also discuss the comparison of XL‐MS results to MD simulations and known structures as a means to test and validate experimental XL‐MS methods.     

Monitoring anthrax toxin receptor dissociation from the protective antigen by NMR

The binding of the Bacillus anthracis protective antigen (PA) to the host cell receptor is the first step toward the formation of the anthrax toxin, a tripartite set of proteins that include the enzymatic moieties edema factor (EF), and lethal factor (LF). PA is cleaved by a furin‐like protease on the cell surface followed by the formation of a donut‐shaped heptameric prepore. The prepore undergoes a major structural transition at acidic pH that results in the formation of a membrane spanning pore, an event which is dictated by interactions with the receptor and necessary for entry of EF and LF into the cell. We provide direct evidence using 1‐dimensional ¹³C‐edited ¹H NMR that low pH induces dissociation of the Von‐Willebrand factor A domain of the receptor capillary morphogenesis protein 2 (CMG2) from the prepore, but not the monomeric full length PA. Receptor dissociation is also observed using a carbon‐13 labeled, 2‐fluorohistidine labeled CMG2, consistent with studies showing that protonation of His‐121 in CMG2 is not a mechanism for receptor release. Dissociation is likely caused by the structural transition upon formation of a pore from the prepore state rather than protonation of residues at the receptor PA or prepore interface.     

Theoretical and Experimental Study of DNA Helix-Coil Transition in Acidic and Alkaline Medium

Abstract

The theoretical approach to the calculation of the influence of selective binding of small ligands on DNA helix-coil transition has been described in the previous paper (Lando D.Yu., J. Biomol. Struct. Dyrt., (1994)). In the present paper that method is used for the study of DNA protonation and deprotonation in acidic and alkaline medium by theoretical analysis of pH effect on DNA heat denaturation.

The mechanism of DNA protonation in acidic medium and pK values of nucleotides are well known. It gave us an opportunity to check the theory without any fitting of pK values. A good agreement between experimental and calculated functions T_m(pH) and ΔT(pH) (melting temperature and melting range width) obtained for acidic medium proved the validity of the theory. However, for alkaline medium there was not even qualitative agreement when the agreed-upon mechanism of deprotonation was considered. Looking into the cause of the discrepancy, we have studied the DNA melting for different mechanisms of deprotonation by calculation of T_m(pH) and ΔT(pH). As a result, it has been established that the discrepancy is due to deprotonation of bonded GC base pairs of helical DNA regions (pK= 11).

It was shown that the early known protonation and newly found deprotonation of helical DNA essentially stabilised double helix in alkaline and acidic medium.     

All-atom molecular dynamics simulations reveal significant differences in interaction between antimycin and conserved amino acid residues in bovine and bacterial bc1 complexes

Antimycin A is the most frequently used specific and powerful inhibitor of the mitochondrial respiratory chain. We used all-atom molecular dynamics (MD) simulations to study the dynamic aspects of the interaction of antimycin A with the Q_i site of the bacterial and bovine bc₁ complexes embedded in a membrane. The MD simulations revealed considerable conformational flexibility of antimycin and significant mobility of antimycin, as a whole, inside the Q_i pocket. We conclude that many of the differences in antimycin binding observed in high-resolution x-ray structures may have a dynamic origin and result from fluctuations of protein and antimycin between multiple conformational states of similar energy separated by low activation barriers, as well as from the mobility of antimycin within the Q_i pocket. The MD simulations also revealed a significant difference in interaction between antimycin and conserved amino acid residues in bovine and bacterial bc₁ complexes. The strong hydrogen bond between antimycin and conserved Asp-228 (bovine numeration) was observed to be frequently broken in the bacterial bc₁ complex and only rarely in the bovine bc₁ complex. In addition, the distances between antimycin and conserved His-201 and Lys-227 were consistently larger in the bacterial bc₁ complex. The observed differences could be responsible for a weaker interaction of antimycin with the bacterial bc₁ complex.     

Interaction of Melanin with Proteins – The Importance of an Acidic Intramelanosomal pH

Melanin is a highly irregular heteropolymer consisting of monomeric units derived from the enzymatic oxidation of the amino acid tyrosine. The process of melanin formation takes place in specialized acidic organelles (melanosomes) in melanocytes. The process of melanin polymerization requires an alkaline pH in vitro, and therefore, the purpose of an acidic environment in vivo remains a mystery. It is known that melanin is always bound to protein in vivo. It is also seen that polymerization in vitro at an acidic pH necessarily requires the presence of proteins. The effect of various model proteins on melanin synthesis and their interaction with melanin was studied. It was seen that many proteins could increase melanin synthesis at an acidic pH, and that different proteins resulted in the formation of different states of melanin, i.e., a precipitate or a soluble, protein‐bound form. We also present evidence to show that soluble protein‐bound melanin is present in vivo (in B16 cells as well as in B16 melanoma tissue). An acidic pH appeared to be necessary to ensure the formation of a uniform, very high molecular weight melano–protein complex. The interaction between melanin and proteins appears to be largely charge‐dependent as evidenced by zeta potential measurements, and this interaction is also increased in an acidic pH. Thus, it appears that an acidic intramelanosomal pH is essential to ensure maximum interaction between protein and melanin, and also to ensure that all the melanin formed is protein‐bound.     

Constant pH molecular dynamics of proteins in explicit solvent with proton tautomerism

pH is a ubiquitous regulator of biological activity, including protein‐folding, protein‐protein interactions, and enzymatic activity. Existing constant pH molecular dynamics (CPHMD) models that were developed to address questions related to the pH‐dependent properties of proteins are largely based on implicit solvent models. However, implicit solvent models are known to underestimate the desolvation energy of buried charged residues, increasing the error associated with predictions that involve internal ionizable residue that are important in processes like hydrogen transport and electron transfer. Furthermore, discrete water and ions cannot be modeled in implicit solvent, which are important in systems like membrane proteins and ion channels. We report on an explicit solvent constant pH molecular dynamics framework based on multi‐site λ‐dynamics (CPHMD^MSλD). In the CPHMD^MSλD framework, we performed seamless alchemical transitions between protonation and tautomeric states using multi‐site λ‐dynamics, and designed novel biasing potentials to ensure that the physical end‐states are predominantly sampled. We show that explicit solvent CPHMD^MSλD simulations model realistic pH‐dependent properties of proteins such as the Hen‐Egg White Lysozyme (HEWL), binding domain of 2‐oxoglutarate dehydrogenase (BBL) and N‐terminal domain of ribosomal protein L9 (NTL9), and the pK_a predictions are in excellent agreement with experimental values, with a RMSE ranging from 0.72 to 0.84 pK_a units. With the recent development of the explicit solvent CPHMD^MSλD framework for nucleic acids, accurate modeling of pH‐dependent properties of both major class of biomolecules—proteins and nucleic acids is now possible. Proteins 2014; 82:1319–1331. © 2013 Wiley Periodicals, Inc.     

Enhanced conformational flexibility of the histone-like (HU) protein from Mycoplasma gallisepticum

The histone-like (HU) protein is one of the major nucleoid-associated proteins involved in DNA supercoiling and compaction into bacterial nucleoid as well as in all DNA-dependent transactions. This small positively charged dimeric protein binds DNA in a non-sequence specific manner promoting DNA super-structures. The majority of HU proteins are highly conserved among bacteria; however, HU protein from Mycoplasma gallisepticum (HUMgal) has multiple amino acid substitutions in the most conserved regions, which are believed to contribute to its specificity to DNA targets unusual for canonical HU proteins. In this work, we studied the structural dynamic properties of the HUMgal dimer by NMR spectroscopy and MD simulations. The obtained all-atom model displays compliance with the NMR data and confirms the heterogeneous backbone flexibility of HUMgal. We found that HUMgal, being folded into a dimeric conformation typical for HU proteins, has a labile α-helical body with protruded β-stranded arms forming DNA-binding domain that are highly flexible in the absence of DNA. The amino acid substitutions in conserved regions of the protein are likely to affect the conformational lability of the HUMgal dimer that can be responsible for complex functional behavior of HUMgal in vivo, e.g. facilitating its spatial adaptation to non-canonical DNA-targets.     

pH‐dependent molecular dynamics of vesicular stomatitis virus glycoprotein G

Vesicular stomatitis virus glycoprotein G (VSV‐G) belongs to a new class of viral fusion proteins (Class III). The structure of VSV‐G has been solved in two different conformations and fusion is known to be triggered by low pH. To investigate Class III fusion mechanisms, molecular dynamics simulations were performed on the VSV‐G prefusion structure in two different protonation states: at physiological pH (pH 7) and low pH present in the endosome (pH 5). Domain IV containing the fusion loops, which need to interact with the target membrane, exhibits the highest mobility. Energetic analyses revealed weakened interaction between Domain IV and the protein core at pH 5, which can be attributed to two pairs of structurally neighboring conserved and differentially protonated residues in the Domain IV–core interface. Energetic calculations also demonstrated that the interaction between the subunits in the core of the trimeric VSV‐G is strengthened at pH 5, mainly due to newly formed interactions between the C‐terminal loop of Domain II and the N‐terminus of the adjacent subunit. A pair of interacting residues in this interface that is affected by differential protonation was shown to be the main effectors of this phenomenon. The results of this study thus enhance the mechanistic understanding of the effects of protonation changes in VSV‐G. Proteins 2012. © 2012 Wiley Periodicals, Inc.