Theoretical study of the conformation of the lipoamide arm in a mutant H protein

The lipoamide arm of the H protein plays a pivotal role in the catalytic cycle of the glycine decarboxylase complex (GDC) by being successively methylamine loaded (Hmet), reduced (Hred), and oxidized (Hox). In a previous study, we calculated free-energy surfaces as a function of the lipoamide arm position of the three forms of the wild-type protein and found close agreement with the available experimental data. Our simulations, together with crystallographic and NMR data, showed that the methylamine-loaded arm is locked in a cavity by interaction with Ser12, Glu14, and Asp67. In this work, we investigate the behavior of the methylamine-loaded form of a mutant H protein (HEA) where Glu14 has been replaced by Ala. We find that the arm can still be held in the cavity but that the energy barrier to release of the arm is halved from approximately 40 kcal mol(-1) for Hmet to approximately 12 kcal mol(-1) for HEA. To compensate for the loss of Glu14, the methylamine group shifts toward Ser66 in the mutant form. These results provide a structural basis for the equilibrium between the loaded and the unloaded forms of the arm observed by Gueguen et al. (Gueguen et al., J Biol Chem 1999;274:26344-26352) in HEA.     

Combined structural and biochemical analysis of the H-T complex in the glycine decarboxylase cycle: evidence for a destabilization mechanism of the H-protein

The lipoate containing H-protein plays a pivotal role in the catalytic cycle of the glycine decarboxylase complex (GDC), undergoing reducing methylamination, methylene transfer, and oxidation. The transfer of the CH(2) group is catalyzed by the T-protein, which forms a 1:1 complex with the methylamine-loaded H-protein (Hmet). The methylamine group is then deaminated and transferred to the tetrahydrofolate-polyglutamate (H(4)FGlu(n)) cofactor of T-protein, forming methylenetetrahydrofolate-polyglutamate. The methylamine group is buried inside the protein structure and highly stable. Experimental data show that the H(4)FGlu(n) alone does not induce transfer of the methylene group, and molecular modeling also indicates that the reaction cannot take place without significant structural perturbations of the H-protein. We have, therefore, investigated the effect of the presence of the T-protein on the stability of Hmet. Addition of T-protein without H(4)FGlu(n) greatly increases the rate of the unloading reaction of Hmet, reducing the activation energy by about 20 kcal mol(-1). Differences of the (1)H and (15)N chemical shifts of the H-protein in its isolated form and in the complex with the T-protein show that the interaction surface for the H-protein is localized on one side of the cleft where the lipoate arm is positioned. This suggests that the role of the T-protein is not only to locate the tetrahydrofolate cofactor in a position favorable for a nucleophilic attack on the methylene carbon but also to destabilize the H-protein in order to facilitate the unlocking of the arm and initiate the reaction.     

Interaction between the lipoamide-containing H-protein and the lipoamide dehydrogenase (L-protein) of the glycine decarboxylase multienzyme system 2. Crystal structures of H- and L-proteins.

The glycine decarboxylase complex consists of four different component enzymes (P-, H-, T- and L-proteins). The 14-kDa lipoamide-containing H-protein plays a pivotal role in the complete sequence of reactions as its prosthetic group (lipoic acid) interacts successively with the three other components of the complex and undergoes a cycle of reductive methylamination, methylamine transfer and electron transfer. With the aim to understand the interaction between the H-protein and its different partners, we have previously determined the crystal structure of the oxidized and methylaminated forms of the H-protein. In the present study, we have crystallized the H-protein in its reduced state and the L-protein (lipoamide dehydrogenase or dihydrolipoamide dehydrogenase). The L-protein has been overexpressed in Escherichia coli and refolded from inclusion bodies in an active form. Crystals were obtained from the refolded L-protein and the structure has been determined by X-ray crystallography. This first crystal structure of a plant dihydrolipoamide dehydrogenase is similar to other known dihydrolipoamide dehydrogenase structures. The crystal structure of the H-protein in its reduced form has been determined and compared to the structure of the other forms of the protein. It is isomorphous to the structure of the oxidized form. In contrast with methylaminated H-protein where the loaded lipoamide arm was locked into a cavity of the protein, the reduced lipoamide arm appeared freely exposed to the solvent. Such a freedom is required to allow its targeting inside the hollow active site of L-protein. Our results strongly suggest that a direct interaction between the H- and L-proteins is not necessary for the reoxidation of the reduced lipoamide arm bound to the H-protein. This hypothesis is supported by biochemical data [Neuburger, M., Polidori, A.M., Piètre, E., Faure, M., Jourdain, A., Bourguignon, J., Pucci, B. & Douce, R. (2000) Eur. J. Biochem. 267, 2882-2889] and by small angle X-ray scattering experiments reported herein.     

Investigation of the local structure and dynamics of the H subunit of the mitochondrial glycine decarboxylase using heteronuclear NMR spectroscopy.

The lipoate-dependent H protein plays a pivotal role in the catalytic cycle of the glycine decarboxylase complex (GDC), undergoing reducing methylamination, methylene transfer, and oxidation. The local structure and backbone dynamics of the methylamine-loaded H (Hmet), oxidized H (Hox), and H apoprotein (Hapo) have been investigated in solution. Filtered NOESY experiments using a [13C]Hmet as well as comparison of the heteronuclear shifts between the Hox and Hmet proteins demonstrate that the methylamine group is located inside a cleft of the protein. Furthermore, this group appears to be locked in this configuration as indicated by the high value of the activation energy (37 kcal/mol) of the global unloading reaction and by its restricted mobility, deduced from 13C relaxation measurements. Comparisons of the 1H and 15N chemical shifts and 15N relaxation in the three forms suggest that part of the lipoyl-lysine arm interacts with the protein polypeptide in the Hox and Hmet. The major change induced by the loading of the methylamine group concerns the C-terminal helix whose mobility becomes completely restricted compared to those of the Hox and Hapo. This C-terminal helix exhibits different reorientational characteristics in the three forms, which can be explained in the Hapo by a model consisting of a twisting motion about an axis passing through the helix. Our results indicate that the model of a freely swinging arm proposed for other lipoate-containing proteins is not acceptable in solution for the GDC. The implication of this observation in terms of the mechanism of the interaction of the H protein with the T protein, its physiological partner during the catalytic cycle, is discussed.     

Crystal Structure of Aminomethyltransferase in Complex with Dihydrolipoyl-H-Protein of the Glycine Cleavage System: IMPLICATIONS FOR RECOGNITION OF LIPOYL PROTEIN SUBSTRATE,DISEASE-RELATED MUTATIONS,AND REACTION MECHANISM*

Aminomethyltransferase, a component of the glycine cleavage system termed T-protein, reversibly catalyzes the degradation of the aminomethyl moiety of glycine attached to the lipoate cofactor of H-protein, resulting in the production of ammonia, 5,10-methylenetetrahydrofolate, and dihydrolipoate-bearing H-protein in the presence of tetrahydrofolate. Several mutations in the human T-protein gene are known to cause nonketotic hyperglycinemia. Here, we report the crystal structure of Escherichia coli T-protein in complex with dihydrolipoate-bearing H-protein and 5-methyltetrahydrofolate, a complex mimicking the ternary complex in the reverse reaction. The structure of the complex shows a highly interacting intermolecular interface limited to a small area and the protein-bound dihydrolipoyllysine arm inserted into the active site cavity of the T-protein. Invariant Arg²⁹² of the T-protein is essential for complex assembly. The structure also provides novel insights in understanding the disease-causing mutations, in addition to the disease-related impairment in the cofactor-enzyme interactions reported previously. Furthermore, structural and mutational analyses suggest that the reversible transfer of the methylene group between the lipoate and tetrahydrofolate should proceed through the electron relay-assisted iminium intermediate formation.     

Computational detection of the binding-site hot spot at the remodeled human growth hormone-receptor interface

A hierarchical computational approach is used to identify the engineered binding-site cavity at the remodeled intermolecular interface between the mutants of human growth hormone (hGH) and the extracellular domain of its receptor (hGHbp). Multiple docking simulations are conducted with the remodeled hGH-hGHbp complex for a panel of potent benzimidazole-containing inhibitors that can restore the binding affinity of the wild-type complex, and for a set of known nonactive small molecules that contain different heterocyclic motifs. Structural clustering of ligand-bound conformations and binding free-energy calculations, using the AMBER force field and a continuum solvation model, can rapidly locate and screen numerous ligand-binding modes on the protein surface and detect the binding-site hot spot at the intermolecular interface. Structural orientation of the benzimidazole motif in the binding-site cavity closely mimics the position of the hot spot residue W104 in the crystal structure of the wild-type complex, which is recognized as an important structural requirement for restoring binding affinity. Despite numerous pockets on the protein surface of the mutant hGH-hGHbp complex, the binding-site cavity presents the energetically favorable hot spot for the benzimidazole-containing inhibitors, whereas for a set of nonactive molecules, the lowest energy ligand conformations do not necessarily bind in the engineered cavity. The results reveal a dominant role of the intermolecular van der Waals interactions in providing favorable ligand-protein energetics in the redesigned interface, in agreement with the experimental and computational alanine scanning of the hGH-hGHbp complex.     

Effects of environment on the structure of Pyrococcus furiosus rubredoxin: a molecular dynamics study

Molecular dynamics simulations based on a 0.95-A resolution crystal structure of Pyrococcus furiosus have been performed to elucidate the effects of the environment on the structure of rubredoxin, and proteins in general. Three 1-ns simulations are reported here: two crystalline state simulations at 123 and 300 K, and a solution state simulation at 300 K. These simulations show that temperature has a greater impact on the protein structure than the close molecular contacts of the crystal matrix in rubredoxin, although both have an effect on its dynamic properties. These results indicate that differences between NMR solution structures and X-ray crystal structures will be relatively minor if they are done at similar temperatures. In addition, the crystal simulations appears to mimic previous crystallographic experiments on the effects of cryo-temperature on temperature factors, and might provide a useful tool in the structural analysis of protein structures solved at cryo-temperatures.     

On the reaction mechanism of class Pi glutathione S-transferase

Theoretical calculations were performed to examine the ionization of the phenolic group of Tyr7 and the thiol group of glutathione in aqueous solution and in the protein class-pi glutathione S-transferase (GST-Pi). Three model systems were considered for simulations in the protein environment: the free enzyme, the complex between glutathione and the enzyme, and the complex between 1-chloro-2.4-dinitrobenzene, glutathione, and the enzyme. The structures derived from Molecular Dynamics simulations were compared with the crystallographic data available for the complex between the inhibitor S-(p-nitrobenzyl)glutathione and GST-Pi, the glutathione-bound form of GST-Pi, and the free enzyme carboxymethylated in Cys47. Free-energy perturbation techniques were used to determine the thermodynamics quantities for ionization of the phenol and thiol groups. The functional implications of Tyr7 in the activation of the glutathione thiol group are discussed in the light of present results, which in agreement with previous studies suggest that Tyr7 in un-ionized form contributes to the catalytic process of glutathione S-transferase, the thiolate anion being stabilized by hydrogen bond with Tyr7 and by interactions with hydrating water molecules. Proteins 28:530–542, 1997 © 1997 Wiley-Liss, Inc.     

Crystal structure of lipoate-protein ligase A from Escherichia coli. Determination of the lipoic acid-binding site

Lipoate-protein ligase A (LplA) catalyzes the formation of lipoyl-AMP from lipoate and ATP and then transfers the lipoyl moiety to a specific lysine residue on the acyltransferase subunit of alpha-ketoacid dehydrogenase complexes and on H-protein of the glycine cleavage system. The lypoyllysine arm plays a pivotal role in the complexes by shuttling the reaction intermediate and reducing equivalents between the active sites of the components of the complexes. We have determined the X-ray crystal structures of Escherichia coli LplA alone and in a complex with lipoic acid at 2.4 and 2.9 angstroms resolution, respectively. The structure of LplA consists of a large N-terminal domain and a small C-terminal domain. The structure identifies the substrate binding pocket at the interface between the two domains. Lipoic acid is bound in a hydrophobic cavity in the N-terminal domain through hydrophobic interactions and a weak hydrogen bond between carboxyl group of lipoic acid and the Ser-72 or Arg-140 residue of LplA. No large conformational change was observed in the main chain structure upon the binding of lipoic acid.     

Influence of RNA Binding on the Structure and Dynamics of the Lassa Virus Nucleoprotein

Lassa virus protects its viral genome through the formation of a ribonucleoprotein complex in which the nucleoprotein (NP) encapsidates the single-stranded RNA genome. Crystal structures provide evidence that a conformational change must occur to allow for RNA binding. In this study, the mechanism by which NP binds to RNA and how the conformational changes in NP are achieved was investigated with molecular-dynamics simulations. NP was structurally characterized in an open configuration when bound to RNA and in a closed form in the absence of RNA. Our results show that when NP is bound to RNA, the protein is highly dynamic and the system undergoes spontaneous deviations away from the open-state configuration. The equilibrium simulations are supported by free-energy calculations that quantify the influence of RNA on the free-energy surface, which governs NP dynamics. We predict that the globally stable states are qualitatively in agreement with the observed crystal structures, but that both open and closed conformations are thermally accessible in the presence of RNA. The free-energy calculations also provide a prediction of the location of the transition state for RNA binding and identify an intermediate metastable state that exhibits correlated motions that could promote RNA binding.     

Exploration of the conformational space of myosin recovery stroke via molecular dynamics

Muscle contractions are driven by cyclic conformational changes of myosin, whose molecular mechanisms of operation are being elucidated by recent advances in crystallographic studies and single molecule experiments. To complement such structural studies and consider the energetics of the conformational changes of myosin head, umbrella sampling molecular dynamics (MD) simulations were performed with the all-atom model of the scallop myosin sub-fragment 1 (S1) with a bound ATP in solution in explicit water using the crystallographic near-rigor and transition state conformations as two references. The constraints on RMSD reaction coordinates used for the umbrella sampling were found to steer the conformational changes efficiently, and relatively close correlations have been observed between the set of characteristic structural changes including the lever arm rotation and the closing of the nucleotide binding pocket. The lever arm angle and key residue interaction distances in the nucleotide binding pocket and the relay helix show gradual changes along the recovery stroke reaction coordinate, consistent with previous crystallographic and computational minimum energy studies. Thermal fluctuations, however, appear to make the switch-2 coordination of ATP more flexible than suggested by crystal structures. The local solvation environment of the fluorescence probe, Trp 507 (scallop numbering), also appears highly mobile in the presence of thermal fluctuations.     

Molecular dynamics simulations of papilloma virus E2 DNA sequences: dynamical models for oligonucleotide structures in solution

The specificity of papilloma virus E2 protein-DNA binding depends critically upon the sequence of a region of the DNA not in direct contact with the protein, and represents one of the simplest known examples of indirect readout. A detailed characterization of this system in solution is important to the further investigation hypothesis of a structural code for DNA recognition by regulatory proteins. In the crystalline state, the E2 DNA oligonucleotide sequence, d(ACCGAATTCGGT), exhibits three different structural forms. We report herein studies of the structure of E2 DNA in solution based on a series of molecular dynamics (MD) simulations including counterions and water, utilizing both the canonical and various crystallographic structures as initial points of departure. All MDs converged on a single dynamical structure of d(ACCGAATTCGGT) in solution. The predicted structure is in close accord with two of the three crystal structures, and indicates that a significant kink in the double helix at the central ApT step in the other crystal molecule may be a packing effect. The dynamical fine structure was analyzed on the basis of helicoidal parameters. The calculated curvature in the sequence was found to originate primarily from YPR steps in the regions flanking the central AATT tract. In order to study the role of structural adaptation of the DNA in the binding process, a subsequent simulation on the 16-mer cognate sequence d(CAACCGAATTCGGTTG) was initiated from the crystallographic coordinates of the bound DNA in the crystal structure of the protein DNA complex. MD simulations starting with the protein-bound form relaxed rapidly back to the dynamical structure predicted from the previous simulations on the uncomplexed DNA. The MD results show that the bound form E2 DNA is a dynamically unstable structure in the absence of protein, and arises as a consequence of both structural changes intrinsic to the sequence and induced by the interaction with protein.     

Catalytic mechanism of cyclophilin as observed in molecular dynamics simulations: pathway prediction and reconciliation of X-ray crystallographic and NMR solution data

Cyclophilins are proteins that catalyze X-proline cis-trans interconversion, where X represents any amino acid. Its mechanism of action has been investigated over the past years but still generates discussion, especially because until recently structures of the ligand in the cis and trans conformations for the same system were lacking. X-ray crystallographic structures for the complex cyclophilin A and HIV-1 capsid mutants with ligands in the cis and trans conformations suggest a mechanism where the N-terminal portion of the ligand rotates during the cis-trans isomerization. However, a few years before, a C-terminal rotating ligand was proposed to explain NMR solution data. In the present study we use molecular dynamics (MD) simulations to generate a trans structure starting from the cis structure. From simulations starting from the cis and trans structures obtained through the rotational pathways, the seeming contradiction between the two sets of experimental data could be resolved. The simulated N-terminal rotated trans structure shows good agreement with the equivalent crystal structure and, moreover, is consistent with the NMR data. These results illustrate the use of MD simulation at atomic resolution to model structural transitions and to interpret experimental data.     

Molecular dynamics simulations of the acyl-enzyme and the tetrahedral intermediate in the deacylation step of serine proteases

Despite the availability of many experimental data and some modeling studies, questions remain as to the precise mechanism of the serine proteases. Here we report molecular dynamics simulations on the acyl-enzyme complex and the tetrahedral intermediate during the deacylation step in elastase catalyzed hydrolysis of a simple peptide. The models are based on recent crystallographic data for an acyl-enzyme intermediate at pH 5 and a time-resolved study on the deacylation step. Simulations were carried out on the acyl enzyme complex with His-57 in protonated (as for the pH 5 crystallographic work) and deprotonated forms. In both cases, a water molecule that could provide the nucleophilic hydroxide ion to attack the ester carbonyl was located between the imidazole ring of His-57 and the carbonyl carbon, close to the hydrolytic position assigned in the crystal structure. In the "neutral pH" simulations of the acyl-enzyme complex, the hydrolytic water oxygen was hydrogen bonded to the imidazole ring and the side chain of Arg-61. Alternative stable locations for water in the active site were also observed. Movement of the His-57 side-chain from that observed in the crystal structure allowed more solvent waters to enter the active site, suggesting that an alternative hydrolytic process directly involving two water molecules may be possible. At the acyl-enzyme stage, the ester carbonyl was found to flip easily in and out of the oxyanion hole. In contrast, simulations on the tetrahedral intermediate showed no significant movement of His-57 and the ester carbonyl was constantly located in the oxyanion hole. A comparison between the simulated tetrahedral intermediate and a time-resolved crystallographic structure assigned as predominantly reflecting the tetrahedral intermediate suggests that the experimental structure may not precisely represent an optimal arrangement for catalysis in solution. Movement of loop residues 216-223 and P3 residue, seen both in the tetrahedral simulation and the experimental analysis, could be related to product release. Furthermore, an analysis of the geometric data obtained from the simulations and the pH 5 crystal structure of the acyl-enzyme suggests that since His-57 is protonated, in some aspects, this crystal structure resembles the tetrahedral intermediate.     

NMR structural refinement of a tandem G.A mismatched decamer d(CCAAGATTGG)2 via the hybrid matrix procedure

A complete relaxation matrix approach employing a matrix eigenvalue/eigenvector solution to the Bloch equations is used to evaluate the NMR solution structure of a tandemly positioned G.A double mismatch decamer oligodeoxyribonucleotide duplex, d(CCAAGATTGG)2. An iterative refinement method using a hybrid relaxation matrix combined with restrained molecular dynamics calculations is shown to provide structures having good agreement with the experimentally derived structures. Distances incorporated into the MD simulations have been calculated from the relaxation rate matrix evaluated from a hybrid NOESY volume matrix whose elements are obtained from the merging of experimental and calculated NOESY intensities. Starting from both A- and B-DNA and mismatch syn and anti models, it is possible to calculate structures that are in good atomic RMS agreement with each other (less than 1.6 A RMS) but differ from the reported crystal structure (greater than 3.6 A). Importantly, the hybrid matrix derived structures are in excellent agreement with the experimental solution conformation as determined by comparison of the 200-ms simulated and experimental NOESY spectra, while the crystallographic data provide spectra that are grossly different.     

Crystal structure of T-protein of the glycine cleavage system. Cofactor binding, insights into H-protein recognition, and molecular basis for understanding nonketotic hyperglycinemia

The glycine cleavage system catalyzes the oxidative decarboxylation of glycine in bacteria and in mitochondria of animals and plants. Its deficiency in human causes nonketotic hyperglycinemia, an inborn error of glycine metabolism. T-protein, one of the four components of the glycine cleavage system,is a tetrahydrofolate dependent aminomethyltransferase. It catalyzes the transfer of the methylene carbon unit to tetrahydrofolate from the methylamine group covalently attached to the lipoamide arm of H-protein. To gain insight into the T-protein function at the molecular level, we have determined the first crystal structure of T-protein from Thermotoga maritima by the multiwavelength anomalous diffraction method of x-ray crystallography and refined four structures: the apoform; the tetrahydrofolate complex; the folinic acid complex; and the lipoic acid complex. The overall fold of T-protein is similar to that of the C-terminal tetrahydrofolate-binding region (residues 421-830) of Arthrobacter globiformis dimethylglycine oxidase. Tetrahydrofolate (or folinic acid) is bound near the center of the tripartite T-protein. Lipoic acid is bound adjacent to the tetrahydrofolate binding pocket, thus defining the interaction surface for H-protein binding. A homology model of the human T-protein provides the structural framework for understanding the molecular mechanisms underlying the development of nonketotic hyperglycinemia due to missense mutations of the human T-protein.     

Substrate binding to mononuclear metallo-beta-lactamase from Bacillus cereus

Structure and dynamics of substrate binding (cefotaxime) to the catalytic pocket of the mononuclear zinc-beta-lactamase from Bacillus cereus are investigated by molecular dynamics simulations. The calculations, which are based on the hydrogen-bond pattern recently proposed by Dal Peraro et al. (J Biol Inorg Chem 2002; 7:704-712), are carried out for both the free and the complexed enzyme. In the resting state, active site pattern and temperature B-factors are in agreement with crystallographic data. In the complexed form, cefotaxime is accommodated into a stable orientation in the catalytic pocket within the nanosecond timescale, interacting with the enzyme zinc-bound hydroxide and the surrounding loops. The beta-lactam ring remains stable and very close to the hydroxide nucleophile agent, giving a stable representation of the productive enzyme-substrate complex.     

Molecular dynamics characterization of the active cavity of carboxypeptidase A and some of its inhibitor adducts.

Molecular dynamics (MD) calculations have been performed on carboxypeptidase A and on its adducts with inhibitors, such as d-phenylalanine (dPhe) and acetate. The catalytically essential zinc ion present in the protein was explicitly included in all the simulations. The simulation was carried out over a sphere of 15 A centered on the zinc ion. The crystallographic water molecules were explicitly taken into account; then the protein was solvated with a 18 A sphere of water molecules. MD calculations were carried out for 45-60 ps. There is no large deviation from the available X-ray structures of native and the dPhe adduct for the MD structures. Average MD structures were calculated starting from the X-ray structure of the dPhe adduct, and, from a structure obtained by docking the inhibitor in the native structure. Comparison between these two structures and with that of the native protein shows that some of the key variations produced by inhibitor binding are reproduced by MD calculations. Addition of acetate induces structural changes relevant for the understanding of the interaction network in the active cavity. The structural variations induced by different inhibitors are examined. The effects of these interactions on the catalytic mechanism and on the binding of substrate are discussed.