Large-scale movement of functional domains facilitates aminoacylation by human mitochondrial phenylalanyl-tRNA synthetase

Structural studies suggest rearrangement of the RNA-binding and catalytic domains of human mitochondrial PheRS (mtPheRS) is required for aminoacylation. Crosslinking the catalytic and RNA-binding domains resulted in a “closed” form of mtPheRS that still catalyzed ATP-dependent Phe activation, but was no longer able to transfer Phe to tRNA and complete the aminoacylation reaction. SAXS experiments indicated the presence of both the closed and open forms of mtPheRS in solution. Together, these results indicate that conformational flexibility of the two functional modules in mtPheRS is essential for its phenylalanylation activity. This is consistent with the evolution of the aminoacyl-tRNA synthetases as modular enzymes consisting of separate domains that display independent activities.     

Structural modelling of the complex of leucyl-tRNA synthetase and mis-aminoacylated tRNA

To assure fidelity of translation, class Ia aminoacyl-tRNA synthetases (aaRSs) edit mis-aminoacylated tRNAs. Mis-attached amino acids and structural water molecules are not included simultaneously in the current crystal structures of the aaRS•tRNA complexes, where the 3′-ends (adenine 76; A76) are bound to the editing sites. A structural model of the completely solvated leucyl-tRNA synthetase complexed with valyl-tRNA^Leu was constructed by exploiting molecular dynamics simulations modified for the present modelling. The results showed that the ribose conformation of A76 is distinct from those observed in the above-mentioned crystal structures, which could be derived from structural constraints in a sandwiched manner induced by the mis-attached valine and tRNA^Leu.     

Computational studies reveal phosphorylation-dependent changes in the unstructured R domain of CFTR

The cystic fibrosis transmembrane conductance regulator (CFTR) is a cAMP-dependent chloride channel that is mutated in cystic fibrosis, an inherited disease of high morbidity and mortality. The phosphorylation of its ∼ 200 amino acid R domain by protein kinase A is obligatory for channel gating under normal conditions. The R domain contains more than ten PKA phosphorylation sites. No individual site is essential but phosphorylation of increasing numbers of sites enables progressively greater channel activity. In spite of numerous studies of the role of the R domain in CFTR regulation, its mechanism of action remains largely unknown. This is because neither its structure nor its interactions with other parts of CFTR have been completely elucidated. Studies have shown that the R domain lacks well-defined secondary structural elements and is an intrinsically disordered region of the channel protein. Here, we have analyzed the disorder pattern and employed computational methods to explore low-energy conformations of the R domain. The specific disorder and secondary structure patterns detected suggest the presence of molecular recognition elements (MoREs) that may mediate phosphorylation-regulated intra- and inter-domain interactions. Simulations were performed to generate an ensemble of accessible R domain conformations. Although the calculated structures may represent more compact conformers than occur in vivo, their secondary structure propensities are consistent with predictions and published experimental data. Equilibrium simulations of a mimic of a phosphorylated R domain showed that it exhibited an increased radius of gyration. In one possible interpretation of these findings, by changing its size, the globally unstructured R domain may act as an entropic spring to perturb the packing of membrane-spanning sequences that constitute the ion permeability pathway and thereby activate channel gating.     

Crystal structure of tryptophanyl-tRNA synthetase complexed with adenosine-5' tetraphosphate: evidence for distributed use of catalytic binding energy in amino acid activation by class I aminoacyl-tRNA synthetases

Tryptophanyl-tRNA synthetase (TrpRS) is a functionally dimeric ligase, which specifically couples hydrolysis of ATP to AMP and pyrophosphate to the formation of an ester bond between tryptophan and the cognate tRNA. TrpRS from Bacillus stearothermophilus binds the ATP analogue, adenosine-5' tetraphosphate (AQP) competitively with ATP during pyrophosphate exchange. Estimates of binding affinity from this competitive inhibition and from isothermal titration calorimetry show that AQP binds 200 times more tightly than ATP both under conditions of induced-fit, where binding is coupled to an unfavorable conformational change, and under exchange conditions, where there is no conformational change. These binding data provide an indirect experimental measurement of +3.0 kcal/mol for the conformational free energy change associated with induced-fit assembly of the active site. Thermodynamic parameters derived from the calorimetry reveal very modest enthalpic changes, consistent with binding driven largely by a favorable entropy change. The 2.5 A structure of the TrpRS:AQP complex, determined de novo by X-ray crystallography, resembles that of the previously described, pre-transition state TrpRS:ATP complexes. The anticodon-binding domain untwists relative to the Rossmann-fold domain by 20% of the way toward the orientation observed for the Products complex. An unexpected tetraphosphate conformation allows the gamma and deltad phosphate groups to occupy positions equivalent to those occupied by the beta and gamma phosphates of ATP. The beta-phosphate effects a 1.11 A extension that relocates the alpha-phosphate toward the tryptophan carboxylate while the PPi mimic moves deeper into the KMSKS loop. This configuration improves interactions between enzyme and nucleotide significantly and uniformly in the adenosine and PPi binding subsites. A new hydrogen bond forms between S194 from the class I KMSKS signature sequence and the PPi mimic. These complementary thermodynamic and structural data are all consistent with the conclusion that the tetraphosphate mimics a transition-state in which the KMSKS loop develops increasingly tight bonds to the PPi leaving group, weakening linkage to the Palpha as it is relocated by an energetically favorable domain movement. Consistent with extensive mutational data on Tyrosyl-tRNA synthetase, this aspect of the mechanism develops high transition-state affinity for the adenosine and pyrophosphate moieties, which move significantly, relative to one another, during the catalytic step.     

Identification of the nucleophilic factors and the productive complex for the editing reaction by leucyl-tRNA synthetase

To ensure fidelity of translation, several aminoacyl-tRNA synthetases (aaRSs) possess editing capability to hydrolyse mis-aminoacylated tRNAs. In this report, based on our previously-modelled structure of leucyl-tRNA synthetase (LeuRS) complexed with valyl-tRNA^Leu, further structural modelling has been performed along with molecular dynamics simulations. This enabled the identification of the nucleophile, which is different from that suggested by the crystal structure of the LeuRS • Nva2AA complex. Our results revealed that the 3′ hydroxyl group of A76 acts as a “gate” to regulate the accessibility of the nucleophile; thus, the opening of the gate leads to the productive complex for the reaction.     

Exit Strategies for Charged tRNA from GluRS

For several class I aminoacyl-tRNA synthetases (aaRSs), the rate-determining step in aminoacylation is the dissociation of charged tRNA from the enzyme. In this study, the following factors affecting the release of the charged tRNA from aaRSs are computationally explored: the protonation states of amino acids and substrates present in the active site, and the presence and the absence of AMP and elongation factor Tu.Through molecular modeling, internal pK_a calculations, and molecular dynamics simulations, distinct, mechanistically relevant post-transfer states with charged tRNA bound to glutamyl-tRNA synthetase from Thermus thermophilus (Glu-tRNA^Glu) are considered. The behavior of these nonequilibrium states is characterized as a function of time using dynamical network analysis, local energetics, and changes in free energies to estimate transitions that occur during the release of the tRNA. The hundreds of nanoseconds of simulation time reveal system characteristics that are consistent with recent experimental studies.Energetic and network results support the previously proposed mechanism in which the transfer of amino acid to tRNA is accompanied by the protonation of AMP to H-AMP. Subsequent migration of proton to water reduces the stability of the complex and loosens the interface both in the presence and in the absence of AMP. The subsequent undocking of AMP or tRNA then proceeds along thermodynamically competitive pathways. Release of the tRNA acceptor stem is further accelerated by the deprotonation of the α-ammonium group on the charging amino acid. The proposed general base is Glu41, a residue binding the α-ammonium group that is conserved in both structure and sequence across nearly all class I aaRSs. This universal handle is predicted through pK_a calculations to be part of a proton relay system for destabilizing the bound charging amino acid following aminoacylation. Addition of elongation factor Tu to the aaRS·tRNA complex stimulates the dissociation of the tRNA core and the tRNA acceptor stem.     

Molecular dynamics-solvated interaction energy studies of protein-protein interactions: the MP1-p14 scaffolding complex

Using the MP1-p14 scaffolding complex from the mitogen-activated protein kinase signaling pathway as model system, we explored a structure-based computational protocol to probe and characterize binding affinity hot spots at protein-protein interfaces. Hot spots are located by virtual alanine-scanning consensus predictions over three different energy functions and two different single-structure representations of the complex. Refined binding affinity predictions for select hot-spot mutations are carried out by applying first-principle methods such as the molecular mechanics generalized Born surface area (MM-GBSA) and solvated interaction energy (SIE) to the molecular dynamics (MD) trajectories for mutated and wild-type complexes. Here, predicted hot-spot residues were actually mutated to alanine, and crystal structures of the mutated complexes were determined. Two mutated MP1-p14 complexes were investigated, the p14(Y56A)-mutated complex and the MP1(L63A,L65A)-mutated complex. Alternative ways to generate MD ensembles for mutant complexes, not relying on crystal structures for mutated complexes, were also investigated. The SIE function, fitted on protein-ligand binding affinities, gave absolute binding affinity predictions in excellent agreement with experiment and outperformed standard MM-GBSA predictions when tested on the MD ensembles of Ras-Raf and Ras-RalGDS protein-protein complexes. For wild-type and mutant MP1-p14 complexes, SIE predictions of relative binding affinities were supported by a yeast two-hybrid assay that provided semiquantitative relative interaction strengths. Results on the MP1-mutated complex suggested that SIE predictions deteriorate if mutant MD ensembles are approximated by just mutating the wild-type MD trajectory. The SIE data on the p14-mutated complex indicated feasibility for generating mutant MD ensembles from mutated wild-type crystal structure, despite local structural differences observed upon mutation. For energetic considerations, this would circumvent costly needs to produce and crystallize mutated complexes. The sensitized protein-protein interface afforded by the p14(Y56A) mutation identified here has practical applications in screening-based discovery of first-generation small-molecule hits for further development into specific modulators of the mitogen-activated protein kinase signaling pathway.     

Molecular modeling studies on CNG channel from bovine retinal rod: a structural model of the cyclic nucleotide-binding domain

A dimeric model of the cyclic nucleotide-binding domain of the all-alpha homomeric cyclic nucleotide-gated channel from bovine retinal rod is constructed. The model, based on the structure of the fairly homologous catabolite gene activator protein (Weber and Steitz, J Mol Biol 1987;198:311-326), is obtained by use of comparative modeling and molecular dynamics simulations. Our model provides a structural basis for the experimentally measured difference in activity between cAMP and cGMP, as well as the different solvent accessibilities of GLY597 in the complex with cGMP, with cAMP and in the protein in free state. In addition, it provides support for the rearrangement of the domain C helix on ligand binding and releasing proposed by Matulef et al. (Neuron 1999;24:443-452).     

Molecular dynamics simulation of highly charged proteins: comparison of the particle-particle particle-mesh and reaction field methods for the calculation of electrostatic interactions

Molecular dynamics (MD) simulations of the activation domain of porcine procarboxypeptidase B (ADBp) were performed to examine the effect of using the particle-particle particle-mesh (P3M) or the reaction field (RF) method for calculating electrostatic interactions in simulations of highly charged proteins. Several structural, thermodynamic, and dynamic observables were derived from the MD trajectories, including estimated entropies and solvation free energies and essential dynamics (ED). The P3M method leads to slightly higher atomic positional fluctuations and deviations from the crystallographic structure, along with somewhat lower values of the total energy and solvation free energy. However, the ED analysis of the system leads to nearly identical results for both simulations. Because of the strong similarity between the results, both methods appear well suited for the simulation of highly charged globular proteins in explicit solvent. However, the lower computational demand of the RF method in the present implementation represents a clear advantage over the P3M method.     

The role of multiple hydrogen-bonding groups in specific alcohol binding sites in proteins: insights from structural studies of LUSH

It is now generally accepted that many of the physiological effects of alcohol consumption are a direct result of binding to specific sites in neuronal proteins such as ion channels or other components of neuronal signaling cascades. Binding to these targets generally occurs in water-filled pockets and leads to alterations in protein structure and dynamics. However, the precise interactions required to confer alcohol sensitivity to a particular protein remain undefined.Using information from the previously solved crystal structures of the Drosophila melanogaster protein LUSH in complexes with short-chain alcohols, we have designed and tested the effects of specific amino acid substitutions on alcohol binding. The effects of these substitutions, specifically S52A, T57S, and T57A, were examined using a combination of molecular dynamics, X-ray crystallography, fluorescence spectroscopy, and thermal unfolding. These studies reveal that the binding of ethanol is highly sensitive to small changes in the composition of the alcohol binding site. We find that T57 is the most critical residue for binding alcohols; the T57A substitution completely abolishes binding, while the T57S substitution differentially affects ethanol binding compared to longer-chain alcohols. The additional requirement for a potential hydrogen-bond acceptor at position 52 suggests that both the presence of multiple hydrogen-bonding groups and the identity of the hydrogen-bonding residues are critical for defining an ethanol binding site. These results provide new insights into the detailed chemistry of alcohol's interactions with proteins.     

Stability and ATP binding of the nucleotide-binding domain of the Wilson disease protein: effect of the common H1069Q mutation

Perturbation of the human copper-transporter Wilson disease protein (ATP7B) causes intracellular copper accumulation and severe pathology, known as Wilson disease (WD). Several WD mutations are clustered within the nucleotide-binding subdomain (N-domain), including the most common mutation, H1069Q. To gain insight into the biophysical behavior of the N-domain under normal and disease conditions, we have characterized wild-type and H1069Q recombinant N-domains in vitro and in silico. The mutant has only twofold lower ATP affinity compared to that of the wild-type N-domain. Both proteins unfold in an apparent two-state reaction at 20 °C and ATP stabilizes the folded state. The thermal unfolding reactions are irreversible and, for the same scan rate, the wild-type protein is more resistant to perturbation than the mutant. For both proteins, ATP increases the activation barrier towards thermal denaturation. Molecular dynamics simulations identify specific differences in both ATP orientation and protein structure that can explain the absence of catalytic activity for the mutant N-domain. Taken together, our results provide biophysical characteristics that may be general to N-domains in other P_1B-ATPases as well as identify changes that may be responsible for the H1069Q WD phenotype in vivo.     

Structural insights into the regulation and the recognition of histone marks by the SET domain of NSD1

The development of epigenetic therapies fuels cancer hope. DNA-methylation inhibitors, histone-deacetylase and histone-methyltransferase (HMTase) inhibitors are being developed as the utilization of epigenetic targets is emerging as an effective and valuable approach to chemotherapy as well as chemoprevention of cancer. The nuclear receptor binding SET domain (NSD) protein is a family of three HMTases, NSD1, NSD2/MMSET/WHSC1, and NSD3/WHSC1L1 that are critical in maintaining the chromatin integrity. A growing number of studies have reported alterations or amplifications of NSD1, NSD2, or NSD3 in numerous carcinogenic events. Reducing NSDs activity through specific lysine-HMTase inhibitors appears promising to help suppressing cancer growth. However, little is known about the NSD pathways and our understanding of the histone lysine-HMTase mechanism is partial. To shed some light on both the recognition and the regulation of epigenetic marks by the SET domain of the NSD family, we investigate the structural mechanisms of the docking of the histone-H4 tail on the SET domain of NSD1. Our finding exposes a key regulatory and recognition mechanism driven by the flexibility of a loop at the interface of the SET and postSET region. Finally, we prospect the special value of this regulatory region for developing specific and selective NSD inhibitors for the epigenetic therapy of cancers.     

Enhanced plasminogen activation by staphylokinase in the presence of streptokinase beta/betagamma domains: plasminogen kringles play a role

Presence of isolated beta or betagamma domains of streptokinase (SK) increased the catalytic activity of staphylokinase (SAK)-plasmin (Pm) complex up to 60%. In contrast, fusion of SK beta or betagamma domains with the C-terminal end of SAK drastically reduced the catalytic activity of the activator complex. The enhancement effect mediated by beta or betagamma domain on Pg activator activity of SAK-Pm complex was reduced greatly (45%) in the presence of isolated kringles of Pg, whereas, kringles did not change cofactor activity of SAK fusion proteins (carrying beta or betagamma domains) significantly. When catalytic activity of SAK-microPm (catalytic domain of Pm lacking kringle domains) complex was examined in the presence of isolated beta and betagamma domains, no enhancement effect on Pg activation was observed, whereas, enzyme complex formed between microplasmin and SAK fusion proteins (SAKbeta and SAKbetagamma) displayed 50-70% reduction in their catalytic activity. The present study, thus, suggests that the exogenously present beta and betagamma interact with Pg/Pm via kringle domains and elevate catalytic activity of SAK-Pm activator complex resulting in enhanced substrate Pg activation. Fusion of beta or betagamma domains with SAK might alter these intermolecular interactions resulting in attenuated functional activity of SAK.     

The role of short-range Cys171-Cys178 disulfide bond in maintaining cutinase active site integrity: A molecular dynamics simulation

Understanding structural determinants in enzyme active site integrity can provide a good knowledge to design efficient novel catalytic machineries. Fusarium solani pisi cutinase with classic triad Ser-His-Asp is a promising enzyme to scrutinize these structural determinants. We performed two MD simulations: one, with the native structure, and the other with the broken Cys171-Cys178 disulfide bond. This disulfide bond stabilizes a turn in active site on which catalytic Asp175 is located. Functionally important H-bonds and atomic fluctuations in catalytic pocket have been changed. We proposed that this disulfide bond within active site can be considered as an important determinant of cutinase active site structural integrity.     

Magnesium-cationic dummy atom molecules enhance representation of DNA polymerase beta in molecular dynamics simulations: improved accuracy in studies of structural features and mutational effects

Human DNA polymerase beta (pol beta) fills gaps in DNA as part of base excision DNA repair. Due to its small size it is a convenient model enzyme for other DNA polymerases. Its active site contains two Mg(2+) ions, of which one binds an incoming dNTP and one catalyzes its condensation with the DNA primer strand. Simulating such binuclear metalloenzymes accurately but computationally efficiently is a challenging task. Here, we present a magnesium-cationic dummy atom approach that can easily be implemented in molecular mechanical force fields such as the ENZYMIX or the AMBER force fields. All properties investigated here, namely, structure and energetics of both Michaelis complexes and transition state (TS) complexes were represented more accurately using the magnesium-cationic dummy atom model than using the traditional one-atom representation for Mg(2+) ions. The improved agreement between calculated free energies of binding of TS models to different pol beta variants and the experimentally determined activation free energies indicates that this model will be useful in studying mutational effects on catalytic efficiency and fidelity of DNA polymerases. The model should also have broad applicability to the modeling of other magnesium-containing proteins.     

Rotational and hinge dynamics of discoidal high density lipoproteins probed by interchain disulfide bond formation

To develop a detailed double belt model for discoidal HDL, we previously scored inter-helical salt bridges between all possible registries of two stacked antiparallel amphipathic helical rings of apolipoprotein (apo) A-I. The top score was the antiparallel apposition of helix 5 with 5 followed closely by appositions of helix 5 with 4 and helix 5 with 6. The rationale for the current study is that, for each of the optimal scores, a pair of identical residues can be identified in juxtaposition directly on the contact edge between the two antiparallel helical belts of apoA-I. Further, these residues are always in the ‘9th position’ in one of the eighteen 11-mer repeats that make up the lipid-associating domain of apoA-I. To illustrate our terminology, 129j (LL5/5) refers to the juxtaposition of the Cα atoms of G129 (in a ‘9th position’) in the pairwise helix 5 domains. We reasoned that if identical residues in the double belt juxtapositions were mutated to a cysteine and kept under reducing conditions during disc formation, we would have a precise method for determining registration in discoidal HDL by formation of a disulfide-linked apoA-I homodimer. Using this approach, we conclude that 129j (LL5/5) is the major rotamer orientation for double belt HDL and propose that the small ubiquitous gap between the pairwise helix 5 portions of the double belt in larger HDL discoidal particles is significantly dynamic to hinge off the disc edge under certain conditions, e.g., in smaller particles or perhaps following binding of the enzyme LCAT. This article is part of a Special Issue entitled Advances in High Density Lipoprotein Formation and Metabolism: A Tribute to John F. Oram (1945-2010).     

Quinol oxidation by c-type cytochromes: structural characterization of the menaquinol binding site of NrfHA

Membrane-bound cytochrome c quinol dehydrogenases play a crucial role in bacterial respiration by oxidizing menaquinol and transferring electrons to various periplasmic oxidoreductases. In this work, the menaquinol oxidation site of NrfH was characterized by the determination of the X-ray structure of Desulfovibrio vulgaris NrfHA nitrite reductase complex bound to 2-heptyl-4-hydroxyquinoline-N-oxide, which is shown to act as a competitive inhibitor of NrfH quinol oxidation activity. The structure, at 2.8-Å resolution, reveals that the inhibitor binds close to NrfH heme 1, where it establishes polar contacts with two essential residues: Asp89, the residue occupying the heme distal ligand position, and Lys82, a strictly conserved residue. The menaquinol binding cavity is largely polar and has a wide opening to the protein surface. Coarse-grained molecular dynamics simulations suggest that the quinol binding site of NrfH and several other respiratory enzymes lie in the head group region of the membrane, which probably facilitates proton transfer to the periplasm. Although NrfH is not a multi-span membrane protein, its quinol binding site has several characteristics similar to those of quinone binding sites previously described. The data presented here provide the first characterization of the quinol binding site of the cytochrome c quinol dehydrogenase family.     

Theoretical mimicry of biomembranes

The study of membrane proteins requires a proper consideration of the specific environment provided by the biomembrane. The compositional complexity of this environment poses great challenges to all experimental and theoretical approaches. In this article a rather simple theoretical concept is discussed for its ability to mimic the biomembrane. The biomembrane is approximated by three mimicry solvents forming individual continuum layers of characteristic physical properties. Several specific structural problems are studied with a focus on the biological significance of such an approach. Our results support the general perception that the biomembrane is crucial for correct positioning and embedding of its constituents. The described model provides a semi-quantitative tool of potential interest to many problems in structural membrane biology.