Knut Teigen

Lipid14: The Amber Lipid Force Field

The AMBER lipid force field has been updated to create Lipid14, allowing tensionless simulation of a number of lipid types with the AMBER MD package. The modular nature of this force field allows numerous combinations of head and tail groups to create different lipid types, enabling the easy insertion of new lipid species. The Lennard-Jones and torsion parameters of both the head and tail groups have been revised and updated partial charges calculated. The force field has been validated by simulating bilayers of six different lipid types for a total of 0.5 μs each without applying a surface tension; with favorable comparison to experiment for properties such as area per lipid, volume per lipid, bilayer thickness, NMR order parameters, scattering data, and lipid lateral diffusion. As the derivation of this force field is consistent with the AMBER development philosophy, Lipid14 is compatible with the AMBER protein, nucleic acid, carbohydrate, and small molecule force fields.

Epac1 and cAMP-dependent Protein Kinase Holoenzyme Have Similar cAMP Affinity, but Their cAMP Domains Have Distinct Structural Features and Cyclic Nucleotide Recognition

The cAMP-dependent protein kinase (PKA I and II) and the cAMP-stimulated GDP exchange factors (Epac1 and -2) are major cAMP effectors. The cAMP affinity of the PKA holoenzyme has not been determined previously. We found that cAMP bound to PKA I with a Kd value (2.9 μm) similar to that of Epac1. In contrast, the free regulatory subunit of PKA type I (RI) had Kd values in the low nanomolar range. The cAMP sites of RI therefore appear engineered to respond to physiological cAMP concentrations only when in the holoenzyme form, whereas Epac can respond in its free form. Epac is phylogenetically younger than PKA, and its functional cAMP site has presumably evolved from site B of PKA. A striking feature is the replacement of a conserved Glu in PKA by Gln (Epac1) or Lys (Epac2). We found that such a switch (E326Q) in site B of human RIα led to a 280-fold decreased cAMP affinity. A similar single switch early in Epac evolution could therefore have decreased the high cAMP affinity of the free regulatory subunit sufficiently to allow Epac to respond to physiologically relevant cAMP levels. Molecular dynamics simulations and cAMP analog mapping indicated that the E326Q switch led to flipping of Tyr-373, which normally stacks with the adenine ring of cAMP. Combined molecular dynamics simulation, GRID analysis, and cAMP analog mapping of wild-type and mutated BI and Epac1 revealed additional differences, independent of the Glu/Gln switch, between the binding sites, regarding space (roominess), hydrophobicity/polarity, and side chain flexibility. This helped explain the specificity of current cAMP analogs and, more importantly, lays a foundation for the generation of even more discriminative analogs.

LIPID11: A Modular Framework for Lipid Simulations using Amber

Accurate simulation of complex lipid bilayers has long been a goal in condensed phase molecular dynamics (MD). Structure and function of membrane-bound proteins are highly dependent on the lipid bilayer environment and are challenging to study through experimental methods. Within Amber, there has been limited focus on lipid simulations, although some success has been seen with the use of the General Amber Force Field (GAFF). However, to date there are no dedicated Amber lipid force fields. In this paper we describe a new charge derivation strategy for lipids consistent with the Amber RESP approach and a new atom and residue naming and type convention. In the first instance, we have combined this approach with GAFF parameters. The result is LIPID11, a flexible, modular framework for the simulation of lipids that is fully compatible with the existing Amber force fields. The charge derivation procedure, capping strategy, and nomenclature for LIPID11, along with preliminary simulation results and a discussion of the planned long-term parameter development are presented here. Our findings suggest that LIPID11 is a modular framework feasible for phospholipids and a flexible starting point for the development of a comprehensive, Amber-compatible lipid force field.

Selectivity and Affinity Determinants for Ligand Binding to the Aromatic Amino Acid Hydroxylases

Hydroxylation of the aromatic amino acids phenylalanine, tyrosine and tryptophan is carried out by a family of non-heme iron and tetrahydrobiopterin (BH4) dependent enzymes, i.e. the aromatic amino acid hydroxylases (AAHs). The reactions catalyzed by these enzymes are important for biomedicine and their mutant forms in humans are associated with phenylketonuria (phenylalanine hydroxylase), Parkinsons disease and DOPA-responsive dystonia (tyrosine hydroxylase), and possibly neuropsychiatric and gastrointestinal disorders (tryptophan hydroxylase 1 and 2). We attempt to rationalize current knowledge about substrate and inhibitor specificity based on the three-dimensional structures of the enzymes and their complexes with substrates, cofactors and inhibitors. In addition, further insights on the selectivity and affinity determinants for ligand binding in the AAHs were obtained from molecular interaction field (MIF) analysis. We applied this computational structural approach to a rational analysis of structural differences at the active sites of the enzymes, a strategy that can help in the design of novel selective ligands for each AAH.

Phosphorylation and Mutations of Ser16 in Human Phenylalanine Hydroxylase KINETIC AND STRUCTURAL EFFECTS

Phosphorylation of phenylalanine hydroxylase (PAH) at Ser16 by cyclic AMP-dependent protein kinase is a post-translational modification that increases its basal activity and facilitates its activation by the substrate l-Phe. So far there is no structural information on the flexible N-terminal tail (residues 1–18), including the phosphorylation site. To get further insight into the molecular basis for the effects of phosphorylation on the catalytic efficiency and enzyme stability, molecular modeling was performed using the crystal structure of the recombinant rat enzyme. The most probable conformation and orientation of the N-terminal tail thus obtained indicates that phosphorylation of Ser16 induces a local conformational change as a result of an electrostatic interaction between the phosphate group and Arg13 as well as a repulsion by Glu280 in the loop at the entrance of the active site crevice structure. The modeled reorientation of the N-terminal tail residues (Met1–Leu15) on phosphorylation is in agreement with the observed conformational change and increased accessibility of the substrate to the active site, as indicated by circular dichroism spectroscopy and the enzyme kinetic data for the full-length phosphorylated and nonphosphorylated human PAH. To further validate the model we have prepared and characterized mutants substituting Ser16 with a negatively charged residue and found that S16E largely mimics the effects of phosphorylation of human PAH. Both the phosphorylated enzyme and the mutants with acidic side chains instead of Ser16 revealed an increased resistance toward limited tryptic proteolysis and, as indicated by circular dichroism spectroscopy, an increased content of α-helical structure. In agreement with the modeled structure, the formation of an Arg13 to Ser16 phosphate salt bridge and the conformational change of the N-terminal tail also explain the higher stability toward limited tryptic proteolysis of the phosphorylated enzyme. The results obtained with the mutant R13A and E381A further support the model proposed for the molecular mechanism for the activation of the enzyme by phosphorylation.

Specific interaction of the diastereomers 7(R)- and 7(S)-tetrahydrobiopterin with phenylalanine hydroxylase: implications for understanding primapterinuria and vitiligo

Pterin‐4a‐carbinolamine dehydratase (PCD) is an essential component of the phenylalanine hydroxylase (PAH) system, catalyzing the regeneration of the essential cofactor 6(R)‐L‐erythro‐5,6,7,8‐tetrahydrobiopterin [6(R)BH4]. Mutations in PCD or its deactivation by hydrogen peroxide result in the generation of 7(R,S)BH4, which is a potent inhibitor of PAH that has been implicated in primapterinuria, a variant form of phenylketonuria, and in the skin depigmentation disorder vitiligo. We have synthesized and separated the 7(R) and 7(S) diastereomers confirming their structure by NMR. Both 7(R)‐ and 7(S)BH4 function as poor cofactors for PAH, whereas only 7(S)BH4 acts as a potent competitive inhibitor vs. 6(R)BH4 (Ki2.3–4.9 µM). Kinetic and binding studies, as well as characterization of the pterin‐enzyme complexes by fluorescence spectroscopy, revealed that the inhibitory effects of 7(R,S)BH4 on PAH are in fact specifically based on 7(S)BH4 binding. The molecular dynamics simulated structures of the pterin‐PAH complexes indicate that 7(S)BH4 inhibition is due to its interaction with the polar region at the pterin binding site close to Ser‐251, whereas its low efficiency as cofactor is related to a suboptimal positioning toward the catalytic iron. 7(S)BH4 is not an inhibitor for tyrosine hydroxylase (TH) in the physiological range, presumably due to the replacement of Ser‐251 by the corresponding Ala297. Taken together, our results identified structural determinants for the specific regulation of PAH and TH by 7(S)BH4, which in turn aid in the understanding of primapterinuria and acute vitiligo. —Pey, A. L., Martinez, A., Charubala, R., Maitland, D. J., Teigen, K., Calvo, A., Pfleiderer, W., Wood, J. M., Schallreuter, K. U. Specific interaction of the diastereomers 7(R)‐ and 7(S)‐tetrahydrobiopterin with phenylalanine hydroxylase: implications for understanding primapterinuria and vitiligo FASEB J. 20, E1451–E1464 (2006)

Rescuing Proteins of Low Kinetic Stability by Chaperones and Natural Ligands: Phenylketonuria, a Case Study

Publisher Summary Phenylketonuria (PKU) is a disease caused by deleterious mutations in phenylalanine hydroxylase (PAH) and constitutes a paradigm for misfolding diseases. Folding is the process by which a protein reaches a functional and stable native structure, while misfolding can be seen as the failure to attain this fully functional conformation. Natural substrates, cofactors, and inhibitors have effects on protein stability beyond their functional role in enzyme function by the same arguments as for other specific ligands and can be considered as natural chaperone ligands. To avoid pathogenic misfolding, the cell is equipped with protein quality control systems (QCS) mainly including chaperones, the ubiquitin proteasome pathway (UPP) and, in some instances, the aggresome. Binding of a ligand to a specific binding site on the native state of a protein will influence the unfolding equilibrium which will be shifted towards the natively folded state, resulting in an increase in protein stability.

Intramolecular hydrogen bonding in articaine can be related to superior bone tissue penetration: A molecular dynamics study

Local anesthetics (LAs) are drugs that cause reversible loss of nociception during surgical procedures. Articaine is a commonly used LA in dentistry that has proven to be exceptionally effective in penetrating bone tissue and induce anesthesia on posterior teeth in maxilla and mandibula. In the present study, our aim was to gain a deeper understanding of the penetration of articaine through biological membranes by studying the interactions of articaine with a phospholipid membrane. Our approach involves Langmuir monolayer experiments combined with molecular dynamics simulations. Membrane permeability of LAs can be modulated by pH due to a titratable amine group with a pKa value close to physiological pH. A change in protonation state is thus known to act as a lipophilicity switch in LAs. Our study shows that articaine has an additional unique lipophilicity switch in its ability to form an intramolecular hydrogen bond. We suggest this intramolecular hydrogen bond as a novel and additional solvent-dependent mechanism for modulation of lipophilicity of articaine which may enhance its diffusion through membranes and connective tissue.

Probing cofactor specificity in phenylalanine hydroxylase by molecular dynamics simulations.

Abstract Phenylalanine hydroxylase (PAH) is a tetrahydrobiopterin-dependent enzyme that catalyzes the hydroxylation of L-phenylalanine (L-Phe) to L-tyrosine using dioxygen as an additional substrate. The requirement of PAH for a cofactor is absolute, but several cofactor analogs are able to substitute the natural cofactor in catalysis. However, it is only the natural cofactor 6R-tetrahydrobiopterin (6R-BH4) that induces a negative regulatory effect on the enzyme. In order to get further insights on the molecular basis for this specificity, we studied the structure of the cofactor-enzyme complex and the conformational changes induced by cofactor binding by molecular dynamics simulations. Simulations were carried out on the enzyme alone and complexed with 6R-BH4 and with two cofactor analogs, 6S-BH4 and 6-methyl-tetrahydropterin (6M-PH4). In the resting unbound enzyme Tyr377 in the catalytic domain is hydrogen bonded to both Ser23 and Glu21 of the autoregulatory N-terminal sequence. This hydrogen bonding network is disturbed by the binding of BH4, which interacts with Ser23. By doing so, 6R-BH4 facilitates an interaction between Glu21 and the active site iron, further pulling the N-terminal into the active site of PAH and blocking the L-Phe binding site. Thus, in the 6R-BH4 complexed enzyme, the N-terminal functions as an intrinsic amino acid regulatory sequence (IARS). Neither 6M-PH4 nor 6S-BH4 can interact favorably with Ser23, and do not induce an inhibitory effect on PAH. These simulations thus explain the previous findings that the two hydroxyl groups in the side chain of the 6R epimer of BH4 are essential for the inhibitory regulatory effect on PAH.

Structural and stability effects of phosphorylation: Localized structural changes in phenylalanine hydroxylase

Phosphorylation of phenylalanine hydroxylase (PAH) at Ser16 by cAMP‐dependent protein kinase increases the basal activity of the enzyme and its resistance to tryptic proteolysis. The modeled structures of the full‐length phosphorylated and unphosphorylated enzyme were subjected to molecular dynamics simulations, and we analyzed the energy of charge–charge interactions for individual ionizable residues in the final structures. These calculations showed that the conformational changes induced by incorporation of phosphate were localized and limited mostly to the region around the phosphoserine (Arg13–Asp17) and a region around the active site in the catalytic domain that includes residues involved in the binding of the iron and the substrate L‐Phe (Arg270 and His285). The absence of a generalized conformational change was confirmed by differential scanning calorimetry, thermal‐dependent circular dichroism, fluorescence spectroscopy, and limited chymotryptic proteolysis of the phosphorylated and unphosphorylated PAH. Our results explain the effect of phosphorylation of PAH on both the resistance to proteolysis specifically by trypsin‐like enzymes and on the increase in catalytic efficiency.