Calculating pKa values in the cAMP-dependent protein kinase: the effect of conformational change and ligand binding

The conformational change observed upon ligand binding and phosphorylation for the cAMP-dependent protein kinase (protein kinase A-PKA) is of high importance for the regulation of its activity. We calculate pKa values and net charges for 18 3D structures of PKA in various conformations and liganded states to examine the role of electrostatics in ligand binding and activation. We find that the conformational change of PKA takes place without any significant net proton uptake/release at all pH values, thus indicating that PKA has evolved to reduce any pH-dependent barriers to the conformational motion. We furthermore find that the binding of ligands induces large changes in the net charge of PKA at most pH values, but significantly, we find that the net charge difference at physiological pH is close to zero, thus indicating that the active-site pKa values have been preorganized for substrate binding. We are unable to unequivocally resolve the identity of the groups responsible for determining the pH-activity profile of PKA but speculate that the titration of Lys 168 or the titration of ATP itself could be responsible for the loss of activity at high pH values. Finally, we examine the effect of point mutations on the pKa values of the PKA catalytic residues and find these to be relatively insensitive to both noncharge-altering and charge-altering mutations.     

Microscopic mechanisms that govern the titration response and pKa values of buried residues in staphylococcal nuclease mutants

To probe the microscopic mechanisms that govern the titration behavior of buried ionizable groups, microsecond explicit solvent molecular dynamics simulations are carried out for several mutants of Staphylococcal nuclease using both fixed charge and polarizable force fields. While the ionization of Asp 66, Glu 66, and Lys 125 lead to enhanced structural fluctuations and partial unfolding of adjacent α‐helical regions, the ionization of Lys 25 causes local unfolding of adjacent β sheets. Using the sampled conformational ensembles, good agreement with experimental pK_a values is obtained with Poisson–Boltzmann calculations using a protein dielectric constant of 2–4 for V66D/E; slightly larger dielectric constants are needed for Lys mutants especially L25K, suggesting that structural responses beyond microseconds are involved in ionization of Lys 25. Overall, the set of unbiased simulations provides insights into the spatial and temporal scales of protein and solvent motions that dictate the diverse titration behaviors of buried protein residues. Proteins 2017; 85:268–281. © 2016 Wiley Periodicals, Inc.     

Intrinsically disordered inhibitor of glutamine synthetase is a functional protein with random‐coil‐like pKa values

The sequential action of glutamine synthetase (GS) and glutamate synthase (GOGAT) in cyanobacteria allows the incorporation of ammonium into carbon skeletons. In the cyanobacterium Synechocystis sp. PCC 6803, the activity of GS is modulated by the interaction with proteins, which include a 65‐residue‐long intrinsically disordered protein (IDP), the inactivating factor IF7. This interaction is regulated by the presence of charged residues in both IF7 and GS. To understand how charged amino acids can affect the binding of an IDP with its target and to provide clues on electrostatic interactions in disordered states of proteins, we measured the pK_a values of all IF7 acidic groups (Glu32, Glu36, Glu38, Asp40, Asp58, and Ser65, the backbone C‐terminus) at 100 mM NaCl concentration, by using NMR spectroscopy. We also obtained solution structures of IF7 through molecular dynamics simulation, validated them on the basis of previous experiments, and used them to obtain theoretical estimates of the pK_a values. Titration values for the two Asp and three Glu residues of IF7 were similar to those reported for random‐coil models, suggesting the lack of electrostatic interactions around these residues. Furthermore, our results suggest the presence of helical structure at the N‐terminus of the protein and of conformational changes at acidic pH values. The overall experimental and in silico findings suggest that local interactions and conformational equilibria do not play a role in determining the electrostatic features of the acidic residues of IF7.     

The 3A6-TCR/superagonist/HLA-DR2a complex shows similar interface and reduced flexibility compared to the complex with self-peptide

T-cell receptor (TCR) recognition of the myelin basic protein (MBP) peptide presented by major histocompatibility complex (MHC) protein HLA-DR2a, one of the MHC class II alleles associated with multiple sclerosis, is highly variable. Interactions in the trimolecular complex between the TCR of the MBP83-99-specific T cell clone 3A6 with the MBP-peptide/HLA-DR2a (abbreviated TCR/pMHC) lead to substantially different proliferative responses when comparing the wild-type decapeptide MBP90-99 and a superagonist peptide, which differs mainly in the residues that point toward the TCR. Here, we investigate the influence of the peptide sequence on the interface and intrinsic plasticity of the TCR/pMHC trimolecular and pMHC bimolecular complexes by molecular dynamics simulations. The intermolecular contacts at the TCR/pMHC interface are similar for the complexes with the superagonist and the MBP self-peptide. The orientation angle between TCR and pMHC fluctuates less in the complex with the superagonist peptide. Thus, the higher structural stability of the TCR/pMHC tripartite complex with the superagonist peptide, rather than a major difference in binding mode with respect to the self-peptide, seems to be responsible for the stronger proliferative response.     

The role of flexibility and hydration on the sequence-specific DNA recognition by the Tn916 integrase protein: a molecular dynamics analysis

The N-terminal domain of the Tn916 integrase protein (INT-DBD) is responsible for DNA binding in the process of strand cleavage and joining reactions required for transposition of the Tn916 conjugative transposon. Site-specific association is facilitated by numerous protein-DNA contacts from the face of a three-stranded beta-sheet inserted into the major groove. The protein undergoes a subtle conformational transition and is slightly unfolded in the protein-DNA complex. The conformation of many charged residues is poorly defined by NMR data but mutational studies have indicated that removal of polar side chains decreases binding affinity, while non-polar contacts are malleable. Based on analysis of the binding enthalpy and binding heat capacity, we have reasoned that dehydration of the protein-DNA interface is incomplete. This study presents results from a molecular dynamics investigation of the INT-DBD-DNA complex aimed at a more detailed understanding of the role of conformational dynamics and hydration in site-specific binding. Comparison of simulations (total of 13 ns) of the free protein and of the bound protein conformation (in isolation or DNA-bound) reveals intrinsic flexibility in certain parts of the molecule. Conformational adaptation linked to partial unfolding appears to be induced by protein-DNA contacts. The protein-DNA hydrogen-bonding network is highly dynamic. The simulation identifies protein-DNA interactions that are poorly resolved or only surmised from the NMR ensemble. Single water molecules and water clusters dynamically optimize the complementarity of polar interactions at the 'wet' protein-DNA interface. The simulation results are useful to establish a qualitative link between experimental data on individual residue's contribution to binding affinity and thermodynamic properties of INT-DBD alone and in complex with DNA.     

Peptide dynamic fingerprints: a tool for investigating the role of conformational flexibility for GLP-1 analogs affinity.

Glucagon-like peptide-1 (GLP-1) is a 30-residue peptide implicated in short-term appetite regulation. Its analogs are presumed to be potential drugs against obesity and non-insulin dependent diabetes mellitus (NIDDM or type 2 diabetes). This study examined how the dynamic fingerprints can be used for establishing dynamics-activity relationships in a series of peptides for which the mechanism of action is unknown and in which mutations can cause an increase or decrease in biological activity. The 3D autocorrelation method was used to generate maps of both active and inactive analogs. As the active conformation of GLP-1 is not yet clearly defined, the dynamic fingerprints of peptides in an aqueous environment were compared to explain the high affinity of the peptide for its receptor. The suggestion that the peptide could bind to the receptor in a folded conformation has been examined. In the case of the GLP-1 analogs, it was shown that the folding tendency cannot be directly related to affinity values and the results do not favor a folded active conformation model of GLP-1.     

A mechanistic model for understanding invasions: using the environment as a predictor of population success

Aim We set out to develop a temperature‐ and salinity‐dependent mechanistic population model for copepods that can be used to understand the role of environmental parameters in population growth or decline. Models are an important tool for understanding the dynamics of invasive species; our model can be used to determine an organism’s niche and explore the potential for invasion of a new habitat. Location Strait of Georgia, British Columbia, Canada. Methods We developed a birth rate model to determine the environmental niche for an estuarine copepod. We conducted laboratory experiments to estimate demographic parameters over a range of temperatures and salinities for Eurytemora affinis collected from the Nanaimo Estuary, British Columbia (BC). The parameterized model was then used to explore what environmental conditions resulted in population growth vs. decline. We then re‐parameterized our model using previously published data for E. affinis collected in the Seine Estuary, France (SE), and compared the dynamics of the two populations. Results We established regions in temperature–salinity space where E. affinis populations from BC would likely grow vs. decline. In general, the population from BC exhibited positive and higher intrinsic growth rates at higher temperatures and salinities. The population from SE exhibited positive and higher growth rates with increasing temperature and decreasing salinity. These different relationships with environmental parameters resulted in predictions of complex interactions among temperature, salinity and growth rates if the two subspecies inhabited the same estuary. Main conclusions We developed a new mechanistic model that describes population dynamics in terms of temperature and salinity. This model may prove especially useful in predicting the potential for invasion by copepods transported to Pacific north‐west estuaries via ballast water, or in any system where an ecosystem is subject to invasion by a species that shares demographic characteristics with an established (sub)species.     

A bacterial virulence factor with a dual role as an adhesin and a solute-binding protein: the crystal structure at 1.5 A resolution of the PEB1a protein from the food-borne human pathogen Campylobacter jejuni

The PEB1a protein is an antigenic factor exposed on the surface of the food-borne human pathogen Campylobacter jejuni, which has a major role in adherence and host colonisation. PEB1a is also the periplasmic binding protein component of an aspartate/glutamate ABC transporter essential for optimal microaerobic growth on these dicarboxylic amino acids. Here, we report the crystal structure of PEB1a at 1.5 A resolution. The protein has a typical two-domain alpha/beta structure, characteristic of periplasmic extracytoplasmic solute receptors and a chain topology related to the type II subfamily. An aspartate ligand, clearly defined by electron density in the interdomain cleft, forms extensive polar interactions with the protein, the majority of which are made with the larger domain. Arg89 and Asp174 form ion-pairing interactions with the main chain alpha-carboxyl and alpha-amino-groups, respectively, of the ligand, while Arg67, Thr82, Lys19 and Tyr156 co-ordinate the ligand side-chain carboxyl group. Lys19 and Arg67 line a positively charged groove, which favours binding of Asp over the neutral Asn. The ligand-binding cleft is of sufficient depth to accommodate a glutamate. This is the first structure of an ABC-type aspartate-binding protein, and explains the high affinity of the protein for aspartate and glutamate, and its much weaker binding of asparagine and glutamine. Stopped-flow fluorescence spectroscopy indicates a simple bimolecular mechanism of ligand binding, with high association rate constants. Sequence alignments and phylogenetic analyses revealed PEB1a homologues in some Gram-positive bacteria. The alignments suggest a more distant homology with GltI from Escherichia coli, a known glutamate and aspartate-binding protein, but Lys19 and Tyr156 are not conserved in GltI. Our results provide a structural basis for understanding both the solute transport and adhesin/virulence functions of PEB1a.     

The structure of a glycogen phosphorylase glucopyranose spirohydantoin complex at 1.8 A resolution and 100 K: the role of the water structure and its contribution to binding.

A glucopyranose spirohydantoin (a pyranose analogue of the potent herbicide, hydantocidin) has been identified as the highest affinity glucose analogue inhibitor of glycogen phosphorylase b (GPb). In order to elucidate the structural features that contribute to the binding, the structures of GPb in the native T state conformation and in complex with glucopyranose spirohydantoin have been determined at 100 K to 2.0 A and 1.8 A resolution, respectively, and refined to crystallographic R values of 0.197 (R[free] 0.248) and 0.182 (R[free] 0.229), respectively. The low temperature structure of GPb is almost identical to that of the previously determined room temperature structure, apart from a decrease in overall atomic temperature factors ((B) room temperature GPb = 34.9 A2; (B) 100 K GPb = 23.4 A2). The glucopyranose spirohydantoin inhibitor (Ki = 3.0 microM) binds at the catalytic site and induces small changes in two key regions of the protein: the 280s loop (residues 281-286) that results in a decrease in mobility of this region, and the 380s loop (residues 377-385) that undergoes more significant shifts in order to optimize contact to the ligand. The hydantoin group, that is responsible for increasing the affinity of the glucose compound by a factor of 10(3), makes only one hydrogen bond to the protein, from one of its NH groups to the main chain oxygen of His377. The other polar groups of the hydantoin group form hydrogen bonds to five water molecules. These waters are involved in extensive networks of hydrogen bonds and appear to be an integral part of the protein structure. Analysis of the water structure at the catalytic site of the native enzyme, shows that five waters are displaced by ligand binding and that there is a significant decrease in mobility of the remaining waters on formation of the GPb-hydantoin complex. The ability of the inhibitor to exploit existing waters, to displace waters and to recruit new waters appears to be important for the high affinity of the inhibitor.     

基于生态需水保障的农业用水安全评价——以山东省引黄灌区为例

保障农业用水安全和生态安全是流域水资源管理的重点,针对黄河口和山东引黄灌区的用水矛盾,采用阈值分析与地统计学方法,考虑作物蒸散发和有效降雨计算山东省引黄灌区灌溉需水量,在优先保障黄河口不同等级生态需水条件下分析灌区（划分为不同调控区）可用水量的响应特征,以地理信息系统（GIS）为平台计算具有时空差异的山东省引黄灌区农业用水安全压力指数,进而评价多时空尺度下的山东省引黄灌区农业用水安全。结果表明,在平水年,保障适宜等级的生态需水后大概有33%的年份农业用水安全存在压力,保障最低等级的生态需水约有27%的农业用水面临短缺,但是大部分情况下农业用水安全压力指数都在30%之下,然而保障最高等级生态需水后,有50%以上的年份存在农业用水安全压力,这种压力无论是从出现频次还是在指数强度上都有明显增加。在空间尺度上,以打渔张、刘春家、麻湾和簸箕李等灌区为代表的调控区2、11、13、14和17的农业用水压力显著,代表年内保障最低等级的生态需水后调控区的农业用水安全压力指数超过了20%,保障适宜等级的生态需水后,大部分调控区的农业用水安全压力指数超过了20%,调控区2和13的压力指数超过了30%,保障最高等级的生态需水后大部分调控区的农业用水安全压力指数超过了60%。在引黄水量调配过程中,应该按照水文年的不同保障合理的生态水量,同时充分考虑水文气象因子的空间差异性,在不同调控区细化分配方案,平水年份保障适宜等级生态需水后农业用水短缺由55.28降低到18.25亿m³。该评价方法反映了优先保障生态需水后灌区农业水资源保障情况,并能有效降低农业和生态用水之间的矛盾,为管理部门进行"精细配水"提供依据。