Coenzyme reorientations in the active sites of aspartate aminotransferase isoenzymes studied by linear dichroism method

Cytosolic and mitochondrial pig heart aspartate aminotransferases (cAspAT and mAspAT) and chicken heart cAspAT have been oriented in a compressed slab of polyacrylamide gel and their linear dichroism LD spectra have been recorded. The coenzyme's tilt angles in the active sites of chicken cAspAT and pig mAspAT and their quasisubstrate complexes imitating catalytic intermediates have been computed. The computations are based on reduced linear dichroism values (delta A/A), the known directions of the transition dipole moments in the coenzyme ring and atomic coordinates of the coenzyme obtained by X-ray crystallography. It has been found that formation of the enzyme complex with glutarate and protonation of the internal pyridoxal-lysine aldimine induce reorientations of the coenzyme. As a result of protonation, the coenzyme ring tilts by 27 degrees in cAspAT and 13 degrees in mAspAT. Formation of the external aldimine with 2-methylaspartate is accompanied by tilting of the coenzyme ring by 44 degrees in cAspAT and 39 degrees in mAspAT. For the quinonoid complex with erythro-3-hydroxyaspartate, the tilt angles were found to be 63 degrees in cAspAT and 53 degrees in mAspAT. It is inferred that the basic features of the active site dynamics are similar in the three AspAT's studied. The differences in the coenzyme tilt angles between cAspAT and mAspAT may be linked to catalytic and structural peculiarities of the isoenzymes.     

Determination of the protonation state of the Asp dyad: conventional molecular dynamics versus thermodynamic integration

The protonation state of the Asp dyad is important as it can reveal enzymatic mechanisms, and the information this provides can be used in the development of drugs for proteins such as memapsin 2 (BACE-1), HIV-1 protease, and rennin. Conventional molecular dynamics (MD) simulations have been successfully used to determine the preferred protonation state of the Asp dyad. In the present work, we demonstrate that the results obtained from conventional MD simulations can be greatly influenced by the particular force field applied or the values used for control parameters. In principle, free-energy changes between possible protonation states can be used to determine the protonation state. We show that protonation state prediction by the thermodynamic integration (TI) method is insensitive to force field version or to the cutoff for calculating nonbonded interactions (a control parameter). In the present study, the protonation state of the Asp dyad predicted by TI calculations was the same regardless of the force field and cutoff value applied. Contrary to the intuition that conventional MD is more efficient, our results clearly show that the TI method is actually more efficient and more reliable for determining the protonation state of the Asp dyad.     

Long time dynamics of Met-enkephalin: comparison of explicit and implicit solvent models

Met-enkephalin is one of the smallest opiate peptides. Yet, its dynamical structure and receptor docking mechanism are still not well understood. The conformational dynamics of this neuron peptide in liquid water are studied here by using all-atom molecular dynamics (MD) and implicit water Langevin dynamics (LD) simulations with AMBER potential functions and the three-site transferable intermolecular potential (TIP3P) model for water. To achieve the same simulation length in physical time, the full MD simulations require 200 times as much CPU time as the implicit water LD simulations. The solvent hydrophobicity and dielectric behavior are treated in the implicit solvent LD simulations by using a macroscopic solvation potential, a single dielectric constant, and atomic friction coefficients computed using the accessible surface area method with the TIP3P model water viscosity as determined here from MD simulations for pure TIP3P water. Both the local and the global dynamics obtained from the implicit solvent LD simulations agree very well with those from the explicit solvent MD simulations. The simulations provide insights into the conformational restrictions that are associated with the bioactivity of the opiate peptide dermorphin for the delta-receptor.     

Coenzyme Reorientations in the Active Sites of Aspartate Aminotransferase Isoenzymes Studied by Linear Dichroism Method

Abstract

Cytosolic and mitochondrial pig heart aspartate aminotransferases (cAspAT and mAspAT) and chicken heart cAspAT have been oriented in a compressed slab of polyacrylamide gel and their linear dichroism LD spectra have been recorded. The coenzyme's tilt angles in the active sites of chicken cAspAT and pig mAspAT and their quasisubstrate complexes imitating catalytic intermediates have been computed. The computations are based on reduced linear dichroism values (ΔA/A), the known directions of the transition dipole moments in the coenzyme ring and atomic coordinates of the coenzyme obtained by X-ray crystallography. It has been found that formation of the enzyme complex with glutarate and protonation of the internal pyridoxal-lysine aldimine induce reorientations of the coenzyme. As a result of protonation, the coenzyme ring tilts by 27° in cAspAT and 13° in mAspAT. Formation of the external aldimine with 2-mehtylaspartate is accompanied by tilting of the coenzyme ring by 44° in cAspAT and 39° in mAspAT. For the quinonoid complex with erythro-3-hydroxyaspartate, the tilt angles were found to be 63° in cAspAT and 53° in mAspAT. It is inferred that the basic features of the active site dynamics are similar in the three AspATs studied. The diiferences in the coenzyme tilt angles between cAspAT and mAspAT may be linked to catalytic and structural peculiarities of the isoenzymes.     

Evidence from 13C NMR for protonation of carbamyl-P and N-(phosphonacetyl)-L-aspartate in the active site of aspartate transcarbamylase.

Nuclear magnetic resonance has been used to study the binding of [13C]carbamyl-P (90% enriched) to the catalytic subunit of Escherichia coli aspartate transcarbamylase. Upon forming a binary complex, there is a small change in the chemical shift of the carbonyl carbon resonance, 2 Hz upfield at pH 7.0, indicating that the environments of the carbonyl group in the active site and in water are similar. When succinate, an analog of L-aspartate, is added to form a ternary complex, there is a large downfield change in the chemical shift for carbamyl-P, consistent with interaction between the carbonyl group and a proton donor of the enzyme. The change might also be caused by a ring current froma nearby aromatic amino acid residue. From the pH dependence of this downfield change and from the effects of L-aspartate analogs other than succinate, the form of the enzyme involved is proposed to be an isomerized ternary complex, previously observed in temperature jump and proton NMR studies. The downfield change to chemical shift for carbamyl-P bound to the isomerized complex is 17.7 +/- 1.0 Hz. Using this value, the relative ability of other four-carbon dicarboxylic acids to form isomerized ternary complexes with the enzyme and carbamyl-P has been evaluated quantitatively. The 13C peak for the transition state analog N-(phosphonacetyl)-L-aspartate (PALA), 90% enriched specifically at the amide carbonyl group, is shifted 20 Hz downfield of the peak for free PALA upon binding to the catalytic subunit at pH 7.0. In contrast, the peak for [1-13C] phosphonaceatmide shifts upfield by about 6 Hz upon binding. Since PALA induces isomerization of the enzyme and phosphonacetamide does not, these data provide further evidence consistent with protonation of the carbonyl group only upon isomerization. The degrees of protonation is strong acids of the carbonyl groups of PALA, phosphonacetamide and urethan (a model for the labile carbamyl-P) have been determined, as have the chemical shifts for these compounds upon full protonation. From these data it is calculated that the amide carbonyl groups of carbamyl-P and PALA might be protonated to a maximum of about 20% in the isomerized complexes at pH 7.0. The change in conformation of the enzyme-carbamyl-P complex upon binding L-aspartate, previously proposed to aid catalysis by compressing the two substrates together in the active site, may be accompanied by polarization of the C=O bond, making this ordinarily unreactive group a much better electrophile. A keto analog of PALA, 4,5-dicarboxy-2-ketopentyl phosphonate, also binds tightly to the catalytic subunit and induces a very similar conformational change, whereas an alcohol analog, 4,5-dicarboxy-2-hydroxypentyl phosphonate, does not bind tightly, indicating the critical importance of an unhindered carbonyl group with trigonal geometry.     

Protonation and deprotonation enthalpies of alloxan and implications for the structure and energy of its complexes with water: a computational study

The optimized geometries, harmonic vibrational frequencies, and energies of the structures of monohydrated alloxan were computed at the DFT/ωB97X-D and B3LYP/6–311++G** level of theory. Results confirm that the monohydrate exists as a dipolar alloxan–water complex which represents a global minimum on the potential energy surface (PES). Trajectory dynamics simulations show that attempt to reorient this monohydrate, to a more favorable orientation for H-bonding, is opposed by an energy barrier of 25.07?kJ/mol. Alloxan seems to prefer acting as proton donor than proton acceptor. A marked stabilization due to the formation of N–H–OH₂ bond is observed. The concerted proton donor–acceptor interaction of alloxan with one H₂O molecule does not increase the stability of the alloxan–water complex. The proton affinity of the O and N atoms and the deprotonation enthalpy of the NH bond of alloxan are computed at the same level of theory. Results are compared with recent data on uracil, thymine, and cytosine. The intrinsic acidities and basicities of the four pyrimidines were discussed. Results of the present study reveal that alloxan is capable of forming stronger H-bonds and more stable cyclic complex with water; yet it is of much lower basicity than other pyrimidines.     

pH‐dependent conformational dynamics of beta‐secretase 1: A molecular dynamics study

Beta‐secretase 1 (BACE‐1) is an aspartyl protease implicated in the overproduction of β‐amyloid fibrils responsible for Alzheimer disease. The process of β‐amyloid genesis is known to be pH dependent, with an activity peak between solution pH of 3.5 and 5.5. We have studied the pH‐dependent dynamics of BACE‐1 to better understand the pH dependent mechanism. We have implemented support for graphics processor unit (GPU) accelerated constant pH molecular dynamics within the AMBER molecular dynamics software package and employed this to determine the relative population of different aspartyl dyad protonation states in the pH range of greatest β‐amyloid production, followed by conventional molecular dynamics to explore the differences among the various aspartyl dyad protonation states. We observed a difference in dynamics between double‐protonated, mono‐protonated, and double‐deprotonated states over the known pH range of higher activity. These differences include Tyr 71‐aspartyl dyad proximity and active water lifetime. This work indicates that Tyr 71 stabilizes catalytic water in the aspartyl dyad active site, enabling BACE‐1 activity.     

Conformational studies of parathyroid hormone (PTH)/PTH-related protein (PTHrp) chimeric peptides

The N-terminal 1-34 segments of both parathyroid hormone (PTH) and parathyroid hormone-related protein (PTHrP) bind and activate the same membrane-embedded G protein-coupled receptor (PTH1 Rc) present on the surface of cells in target tissues such as bone and kidney. This binding occurs in spite of major differences between the two hormones in their amino acid sequence. Recently, it was shown that in (1-34) PTH/PTHrP hybrid peptides, the N-terminal 1-14 segment of PTHrP is incompatible with the C-terminal 15-34 region of PTH in terms of bioactivity. The sites of incompatibility were identified at positions 5 in PTHrP and 19 in PTH. In the present paper we describe the synthesis, biological evaluation, and conformational characterization of two segmental hybrids: PTHrP(1-27)-[Tyr(34)]bPTH(28-34)-NH(2) (hybrid I) and PTHrP(1-18)-[Nal(23), Tyr(34)]bPTH(19-34)-NH(2) (hybrid II). Hybrid I is as active as PTH(1-34)NH(2) and more than two orders of magnitude more active than hybrid II. The conformational properties of the hybrids were studied in water/trifluoroethanol (TFE) mixtures and in aqueous solutions containing dodecylphosphocholine (DPC) micelles by CD, two-dimensional nmr and computer simulations. Upon addition of TFE to the aqueous solution, both hybrids undergo a coil-helix transition. The helix content in 1:1 water/TFE obtained by CD data is about 75% for both hybrids. In the presence of DPC, helix formation is observed at detergent concentrations above critical micellar concentration and the maximum helix content is of approximately 35 and approximately 30% for hybrid I and II, respectively. Combined nmr analysis, distance geometry, and molecular dynamics calculations suggest that, in both solvent systems, the biologically active hybrid I exhibits two flexible sites, centered at residues 12 and 19, connecting helical segments. The flexibility point at position 19 is not present in the poorly active hybrid II. Our findings support the hypothesis, proposed in our previous work, that in bioactive PTH analogues the presence and location of flexibility points between helical segments are essential for enabling them to fold into the bioactive conformation upon interaction with the PTH1 receptor.     

Restrained least squares refinement of native (calcium) and cadmium-substituted carp parvalbumin using X-ray crystallographic data at 1.6-A resolution

Carp parvalbumin coordinates calcium through one carbonyl oxygen atom and the oxygen-containing side chains of 5 amino acid residues, or 4 residues and a water molecule, in a helix-loop-helix structural motif. Other calcium-binding proteins, including calmodulin and troponin C, also possess this unique calcium-binding design, which is designated EF-hand or calmodulin fold. Parvalbumin has two such sites, labeled CD and EF. Each of the calcium-binding sites of refined structures of proteins belonging to this group has a 7-oxygen coordination sphere except those of the structure of parvalbumin as it was reported in 1975. This structure had been refined at 1.9 A using difference Fourier techniques on film data. The CD site appeared to be 6-coordinate and the EF site 8-coordinate. Results of NMR experiments using 113Cd-substituted parvalbumin, however, indicate that the sites are similar to one another with coordination number greater than 6. To resolve the inconsistency between crystallographic and NMR results, 1.6 A area detector data was collected for native and cadmium-substituted parvalbumin; the structures have been refined to R factors of 18.7% and 16.4%, respectively, with acceptable geometry and low errors in atomic coordinates. Differences between the parvalbumin structure described in 1975 and the present structure are addressed, including the discovery of 7-coordination for both the CD and EF sites.     

CD of gels and suspensions: apparent CD in the soret region of sickle hemoglobin gels

CD spectra in the soret region of sickle-cell deoxyhemoglobin (deoxy-HbS) fiber gels are radically different from the CD of deoxy-HbS in solution. An explanation is found using the Stokes–Mueller representation of the interaction of a polarized beam with the instrument optical train and sample to derive expressions for the apparent CD of gels and suspensions of optically active molecules that consist of randomly oriented domains or particles that are linearly dichroic and linearly birefringent. These theoretical considerations show that the apparent CD spectra from such systems have contributions from the LD and birefringence of each domain even if no net linear birefringence and dichroism is apparent in the sample. Thus, the interpretation of the CD from gels and suspensions is problematic, unless it can be demonstrated that each domain or particle has extremely small absorbance or that the LD and birefringence of each is a very small fraction of the total absorbance. As a result, we conclude that the spectra of HbS gels are not due to the CD of the heme per se; rather, they also reflect the randomly oriented domain structure of the gels and the LD and linear birefringence associated with each domain.     

Influence of the acetylcholinesterase active site protonation on omega loop and active site dynamics

Existence of alternative entrances in acetylcholinesterase (AChE) could explain the contrast between the very high AChE catalytic efficiency and the narrow and long access path to the active site revealed by X-ray crystallography. Alternative entrances could facilitate diffusion of the reaction products or at least water and ions from the active site. Previous molecular dynamics simulations identified side door and back door as the most probable alternative entrances. The simulations of non-inhibited AChE suggested that the back door opening events occur only rarely (0.8% of the time in the 10ns trajectory). Here we present a molecular dynamics simulation of non-inhibited AChE, where the back door opening appears much more often (14% of the time in the 12ns trajectory) and where the side door opening was observed quite frequently (78% of trajectory time). We also present molecular dynamics, where the back door does not open at all, or where large conformational changes of the AChE omega loop occur together with alternative passage opening events. All these differences in AChE dynamical behavior are caused by different protonation states of two glutamate residues located on bottom of the active site gorge (Glu202 and G450 in Mus musculus AChE). Our results confirm the results of previous molecular dynamics simulations, expand the view and suggest the probable reasons for the overall conformational behavior of AChE omega loop.     

Protonation of crotonyl-CoA dienolate by human glutaryl-CoA dehydrogenase occurs by solvent-derived protons

The protonation of crotonyl-CoA dienolate following decarboxylation of glutaconyl-CoA by glutaryl-CoA dehydrogenase was investigated. Although it is generally held that the active sites of acyl-CoA dehydrogenases are desolvated when substrate binds, recent evidence has established that water has access to the active site in these binary complexes of glutaryl-CoA dehydrogenase. The present investigation shows that the dehydrogenase catalyzes (a) a rapid exchange of C-4 methyl protons of crotonyl-CoA with bulk solvent and (b) protonation of crotonyl-CoA dienolate by solvent-derived protons under single turnover conditions. Both of the reactions require the catalytic base, Glu370. These findings indicate that decarboxylation proceeds via a dienolate intermediate. The involvement of water in catalysis by glutaryl-CoA dehydrogenase was previously unrecognized and is in conflict with a classically held intramolecular 1,3-prototropic shift for protonation of crotonyl-CoA dienolate.     

Elucidation of the protonation states of the catalytic residues in mtKasA: implications for inhibitor design

KasA (β-ketoacyl ACP synthase I) is involved in the biosynthetic pathway of mycolic acids, an essential component of the cell wall in Mycobacterium tuberculosis. It was shown that KasA is essential for the survival of the pathogen and thus could serve as a new drug target for the treatment of tuberculosis. The active site of KasA was previously characterized by X-ray crystallography. However, questions regarding the protonation state of specific amino acids, the orientation of the histidine groups within the active site, and additional conformers being accessible at ambient temperatures remain open and have to be addressed prior to the design of new inhibitors. We investigate the active site of KasA in this work by means of structural motifs and relative energies. Molecular dynamics (MD) simulations, free energy perturbation computations, and calculations employing the hybrid quantum mechanics/molecular mechanics (QM/MM) method made it possible to determine the protonation status and reveal important details about the catalytic mechanism of KasA. Additionally, we can rationalize the molecular basis for the acyl-transfer activity in the H311A mutant. Our data strongly suggest that inhibitors should be able to inhibit different protonation states because the enzyme can switch easily between a zwitterionic and neutral state.     

Hydration energy landscape of the active site cavity in cytochrome P450cam

Hydration of protein cavities influences protein stability, dynamics, and function. Protein active sites usually contain water molecules that, upon ligand binding, are either displaced into bulk solvent or retained to mediate protein–ligand interactions. The contribution of water molecules to ligand binding must be accounted for to compute accurate values of binding affinities. This requires estimation of the extent of hydration of the binding site. However, it is often difficult to identify the water molecules involved in the binding process when ligands bind on the surface of a protein. Cytochrome P450cam is, therefore, an ideal model system because its substrate binds in a buried active site, displacing partially disordered solvent, and the protein is well characterized experimentally. We calculated the free energy differences for having five to eight water molecules in the active site cavity of the unliganded enzyme from molecular dynamics simulations by thermodynamic integration employing a three-stage perturbation scheme. The computed free energy differences between the hydration states are small (within 12 kJ mol⁻¹) but distinct. Consistent with the crystallographic determination and studies employing hydrostatic pressure, we calculated that, although ten water molecules could in principle occupy the volume of the active site, occupation by five to six water molecules is thermodynamically most favorable. Proteins 32:381–396, 1998. © 1998 Wiley-Liss, Inc.     

Effect of frictional force on urea and thiourea during an intramembrane transport: an explicit solvent Langevin dynamics simulation method

An explicit solvent Langevin dynamics (LD) simulation method is employed to investigate the effect of frictional force on the solute–solvent interactions in urea and thiourea, which are subjected to varied degree of frictional forces akin to diffusional transport across the membrane system. Utilizing the frictional force concept in LD method, the effect of increasing frictional coefficient is studied on the microscopic properties such as single-point energy and N–C–N bond angle of urea and thiourea in aqueous solutions. The present work follows a self-conscious approach by explicitly including the solvent molecules in the LD simulations. The results are further compared with the standard Langevin approach without explicitly specifying solvent molecules. An analysis of the results reveals that the two cases differ strongly in the intermediate (Brownian) regime of frictional coefficient.     

Calculating pH-dependent free energy of proteins by using Monte Carlo protonation probabilities of ionizable residues

Protein folding, stability, and function are usually influenced by pH. And free energy plays a fundamental role in analysis of such pH-dependent properties. Electrostatics-based theoretical framework using dielectric solvent continuum model and solving Poisson-Boltzmann equation numerically has been shown to be very successful in understanding the pH-dependent properties. However, in this approach the exact computation of pH-dependent free energy becomes impractical for proteins possessing more than several tens of ionizable sites (e.g.>30), because exact evaluation of the partition function requires a summation over a vast number of possible protonation microstates. Here we present a method which computes the free energy using the average energy and the protonation probabilities of ionizable sites obtained by the well-established Monte Carlo sampling procedure. The key feature is to calculate the entropy by using the protonation probabilities. We used this method to examine a well-studied protein (lysozyme) and produced results which agree very well with the exact calculations. Applications to the optimum pH of maximal stability of proteins and protein–DNA interactions have also resulted in good agreement with experimental data. These examples recommend our method for application to the elucidation of the pH-dependent properties of proteins.     

Coupling dynamics and evolutionary information with structure to identify protein regulatory and functional binding sites

Binding sites in proteins can be either specifically functional binding sites (active sites) that bind specific substrates with high affinity or regulatory binding sites (allosteric sites), that modulate the activity of functional binding sites through effector molecules. Owing to their significance in determining protein function, the identification of protein functional and regulatory binding sites is widely acknowledged as an important biological problem. In this work, we present a novel binding site prediction method, Active and Regulatory site Prediction (AR-Pred), which supplements protein geometry, evolutionary, and physicochemical features with information about protein dynamics to predict putative active and allosteric site residues. As the intrinsic dynamics of globular proteins plays an essential role in controlling binding events, we find it to be an important feature for the identification of protein binding sites. We train and validate our predictive models on multiple balanced training and validation sets with random forest machine learning and obtain an ensemble of discrete models for each prediction type. Our models for active site prediction yield a median area under the curve (AUC) of 91% and Matthews correlation coefficient (MCC) of 0.68, whereas the less well-defined allosteric sites are predicted at a lower level with a median AUC of 80% and MCC of 0.48. When tested on an independent set of proteins, our models for active site prediction show comparable performance to two existing methods and gains compared to two others, while the allosteric site models show gains when tested against three existing prediction methods. AR-Pred is available as a free downloadable package at https://github.com/sambitmishra0628/AR-PRED_source .     

OmpT: molecular dynamics simulations of an outer membrane enzyme

Five molecular dynamics simulations (total duration >25 ns) have been performed on the Escherichia coli outer membrane protease OmpT embedded in a dimyristoylphosphatidylcholine lipid bilayer. Globally the protein is conformationally stable. Some degree of tilt of the beta-barrel is observed relative to the bilayer plane. The greatest degree of conformational flexibility is seen in the extracellular loops. A complex network of fluctuating H-bonds is formed between the active site residues, such that the Asp210-His212 interaction is maintained throughout, whereas His212 and Asp83 are often bridged by a water molecule. This supports a catalytic mechanism whereby Asp83 and His212 bind a water molecule that attacks the peptide carbonyl. A configuration yielded by docking calculations of OmpT simulation snapshots and a model substrate peptide Ala-Arg-Arg-Ala was used as the starting point for an extended Huckel calculation on the docked peptide. These placed the lowest unoccupied molecular orbital mainly on the carbon atom of the central C=O in the scissile peptide bond, thus favoring attack on the central peptide by the water held by residues Asp83 and His212. The trajectories of water molecules reveal exchange of waters between the intracellular face of the membrane and the interior of the barrel but no exchange at the extracellular mouth. This suggests that the pore-like region in the center of OmpT may enable access of water to the active site from below. The simulations appear to reveal the presence of specific lipid interaction sites on the surface of the OmpT barrel. This reveals the ability of extended MD simulations to provide meaningful information on protein-lipid interactions.     

Distinguishing two groups of flavin reductases by analyzing the protonation state of an active site carboxylic acid

Flavin-containing reductases are involved in a wide variety of physiological reactions such as photosynthesis, nitric oxide synthesis, and detoxification of foreign compounds, including therapeutic drugs. Ferredoxin-NADP(H)-reductase (FNR) is the prototypical enzyme of this family. The fold of this protein is highly conserved and occurs as one domain of several multidomain enzymes such as the members of the diflavin reductase family. The enzymes of this family have emerged as fusion of a FNR and a flavodoxin. Although the active sites of these enzymes are very similar, different enzymes function in opposite directions, that is, some reduce oxidized nicotinamide adenine dinucleotide phosphate (NADP(+)) and some oxidize reduced nicotinamide adenine dinucleotide phosphate (NADPH). In this work, we analyze the protonation behavior of titratable residues of these enzymes through electrostatic calculations. We find that a highly conserved carboxylic acid in the active site shows a different titration behavior in different flavin reductases. This residue is deprotonated in flavin reductases present in plastids, but protonated in bacterial counterparts and in diflavin reductases. The protonation state of the carboxylic acid may also influence substrate binding. The physiological substrate for plastidic enzymes is NADP(+), but it is NADPH for the other mentioned reductases. In this article, we discuss the relevance of the environment of this residue for its protonation and its importance in catalysis. Our results allow to reinterpret and explain experimental data.