SWEET13 transport of sucrose, but not gibberellin, restores male fertility in Arabidopsis sweet13;14.

SWEET sucrose transporters play important roles in the allocation of sucrose in plants. Some SWEETs were shown to also mediate transport of the plant growth regulator gibberellin (GA). The close physiological relationship between sucrose and GA raised the questions of whether there is a functional connection and whether one or both of the substrates are physiologically relevant. To dissect these two activities, molecular dynamics were used to map the binding sites of sucrose and GA in the pore of SWEET13 and predicted binding interactions that might be selective for sucrose or GA. Transport assays confirmed these predictions. In transport assays, the N76Q mutant had 7x higher relative GA3 activity, and the S142N mutant only transported sucrose. The impaired pollen viability and germination in sweet13;14 double mutants were complemented by the sucrose-selective SWEET13S142N, but not by the SWEET13N76Q mutant, indicating that sucrose is the physiologically relevant substrate and that GA transport capacity is dispensable in the context of male fertility. Therefore, GA supplementation to counter male sterility may act indirectly via stimulating sucrose supply in male sterile mutants. These findings are also relevant in the context of the role of SWEETs in pathogen susceptibility.

Micellar Composition Affects Lipid Accretion Kinetics in Molecular Dynamics Simulations: Support for Lipid Network Reproduction.

Mixed lipid micelles were proposed to facilitate life through their documented growth dynamics and catalytic properties. Our previous research predicted that micellar self-reproduction involves catalyzed accretion of lipid molecules by the residing lipids, leading to compositional homeostasis. Here, we employ atomistic Molecular Dynamics simulations, beginning with 54 lipid monomers, tracking an entire course of micellar accretion. This was done to examine the self-assembly of variegated lipid clusters, allowing us to measure entry and exit rates of monomeric lipids into pre-micelles with different compositions and sizes. We observe considerable rate-modifications that depend on the assembly composition and scrutinize the underlying mechanisms as well as the energy contributions. Lastly, we describe the measured potential for compositional homeostasis in our simulated mixed micelles. This affirms the basis for micellar self-reproduction, with implications for the study of the origin of life.

How does binding of agonist ligands control intrinsic molecular dynamics in human NMDA receptors?

NMDA-type glutamate receptors (NMDAR) are ligand-gated ion channels that contribute to excitatory neurotransmission in the central nervous system. NMDAR dysfunction has been found to be involved in various neurological disorders. Recent crystallographic and EM studies have shown the static structure of different states of the non-human NMDARs. Here we describe a model of a human NMDA receptor (hNMDAR) and its molecular dynamics (MD) before and after the binding of agonist ligands, glutamate and glycine. It is shown that the binding of ligands promotes a global reduction in molecular flexibility that produces a more tightly packed conformation than the unbound hNMDAR, and a higher cooperative regularity of moving. The ligand-induced synchronization of motion, identified on all structural levels of the modular hNMDA receptor is apparently a fundamental factor in channel gating. Although the time scale of the MD simulations (300 ns) was not sufficient to observe the complete gating event, the obtained data has shown the ligand-induced stabilization of hNMDAR that conforms the "going to be open state". We propose a mechanistic dynamic model of the ligand-dependent gating mechanism in the hNMDA receptor. At the binding of the ligands, the differently twisted conformations of the highly flexible receptor are stabilized in unique conformation with a linear molecular axis, which is a condition that is optimal for pore development. By searching the receptor surface, we have identified three new pockets, which are different from the pockets described in the literature as the potential and known positive allosteric modulator binding sites. A successful docking of two NMDAR modulators to their binding sites validates the model of a human NMDA receptor as a biological relevant target.

Probing the stereospecificity of tyrosyl- and glutaminyl-tRNA synthetase with molecular dynamics.

The stereospecificity of aminoacyl-tRNA synthetases helps exclude d-amino acids from protein synthesis and could perhaps be engineered to allow controlled d-amino acylation of tRNA. We use molecular dynamics simulations to probe the stereospecificity of the class I tyrosyl- and glutaminyl-tRNA synthetases (TyrRS, GlnRS), including wildtype enzymes and three point mutants suggested by three different protein design methods. l/d binding free energy differences are obtained by alchemically and reversibly transforming the ligand from L to D in simulations of the protein-ligand complex. The D81Q mutation in Escherichia coli TyrRS is homologous to the D81R mutant shown earlier to have inverted stereospecificity. D81Q is predicted to lead to a rotated ligand backbone and an increased, not a decreased l-Tyr preference. The E36Q mutation in Methanococcus jannaschii TyrRS has a predicted l/d binding free energy difference ΔΔG of just 0.5±0.9kcal/mol, compared to 3.1±0.8kcal/mol for the wildtype enzyme (favoring l-Tyr). The ligand ammonium position is preserved in the d-Tyr complex, while the carboxylate is shifted. Wildtype GlnRS has a similar preference for l-glutaminyl adenylate; the R260Q mutant has an increased preference, even though Arg260 makes a large contribution to the wildtype ΔΔG value.

Protein:Ligand binding free energies: A stringent test for computational protein design.

A computational protein design method is extended to allow Monte Carlo simulations where two ligands are titrated into a protein binding pocket, yielding binding free energy differences. These provide a stringent test of the physical model, including the energy surface and sidechain rotamer definition. As a test, we consider tyrosyl-tRNA synthetase (TyrRS), which has been extensively redesigned experimentally. We consider its specificity for its substrate l-tyrosine (l-Tyr), compared to the analogs d-Tyr, p-acetyl-, and p-azido-phenylalanine (ac-Phe, az-Phe). We simulate l- and d-Tyr binding to TyrRS and six mutants, and compare the structures and binding free energies to a more rigorous "MD/GBSA" procedure: molecular dynamics with explicit solvent for structures and a Generalized Born + Surface Area model for binding free energies. Next, we consider l-Tyr, ac- and az-Phe binding to six other TyrRS variants. The titration results are sensitive to the precise rotamer definition, which involves a short energy minimization for each sidechain pair to help relax bad contacts induced by the discrete rotamer set. However, when designed mutant structures are rescored with a standard GBSA energy model, results agree well with the more rigorous MD/GBSA. As a third test, we redesign three amino acid positions in the substrate coordination sphere, with either l-Tyr or d-Tyr as the ligand. For two, we obtain good agreement with experiment, recovering the wildtype residue when l-Tyr is the ligand and a d-Tyr specific mutant when d-Tyr is the ligand. For the third, we recover His with either ligand, instead of wildtype Gln.

Redesigning the stereospecificity of tyrosyl-tRNA synthetase.

D-Amino acids are largely excluded from protein synthesis, yet they are of great interest in biotechnology. Unnatural amino acids have been introduced into proteins using engineered aminoacyl-tRNA synthetases (aaRSs), and this strategy might be applicable to D-amino acids. Several aaRSs can aminoacylate their tRNA with a D-amino acid; of these, tyrosyl-tRNA synthetase (TyrRS) has the weakest stereospecificity. We use computational protein design to suggest active site mutations in Escherichia coli TyrRS that could increase its D-Tyr binding further, relative to L-Tyr. The mutations selected all modify one or more sidechain charges in the Tyr binding pocket. We test their effect by probing the aminoacyl-adenylation reaction through pyrophosphate exchange experiments. We also perform extensive alchemical free energy simulations to obtain L-Tyr/D-Tyr binding free energy differences. Agreement with experiment is good, validating the structural models and detailed thermodynamic predictions the simulations provide. The TyrRS stereospecificity proves hard to engineer through charge-altering mutations in the first and second coordination shells of the Tyr ammonium group. Of six mutants tested, two are active towards D-Tyr; one of these has an inverted stereospecificity, with a large preference for D-Tyr. However, its activity is low. Evidently, the TyrRS stereospecificity is robust towards charge rearrangements near the ligand. Future design may have to consider more distant and/or electrically neutral target mutations, and possibly design for binding of the transition state, whose structure however can only be modeled.

An allosteric signaling pathway of human 3-phosphoglycerate kinase from force distribution analysis.

3-Phosphogycerate kinase (PGK) is a two domain enzyme, which transfers a phosphate group between its two substrates, 1,3-bisphosphoglycerate bound to the N-domain and ADP bound to the C-domain. Indispensable for the phosphoryl transfer reaction is a large conformational change from an inactive open to an active closed conformation via a hinge motion that should bring substrates into close proximity. The allosteric pathway resulting in the active closed conformation has only been partially uncovered. Using Molecular Dynamics simulations combined with Force Distribution Analysis (FDA), we describe an allosteric pathway, which connects the substrate binding sites to the interdomain hinge region. Glu192 of alpha-helix 7 and Gly394 of loop L14 act as hinge points, at which these two secondary structure elements straighten, thereby moving the substrate-binding domains towards each other. The long-range allosteric pathway regulating hPGK catalytic activity, which is partially validated and can be further tested by mutagenesis, highlights the virtue of monitoring internal forces to reveal signal propagation, even if only minor conformational distortions, such as helix bending, initiate the large functional rearrangement of the macromolecule.

Importance of aspartate residues in balancing the flexibility and fine-tuning the catalysis of human 3-phosphoglycerate kinase.

The exact role of the metal ion, usually Mg(2+), in the catalysis of human 3-phosphoglycerate kinase, a well-studied two-domain enzyme, has not been clarified. Here we have prepared single and double alanine mutants of the potential metal-binding residues, D374 and D218. While all mutations weaken the catalytic interactions with Mg(2+), they surprisingly strengthen binding of both MgADP and MgATP, and the effects are even more pronounced for ADP and ATP. Thermodynamic parameters of binding indicate an increase in the binding entropy as a reason for the strengthening. In agreement with the experimental results, computer-simulated annealing calculations for the complexes of these mutants have supported the mobility of the nucleotide phosphates and, as a consequence, formation of their new interaction(s) within the active site. A similar type of mobility is suggested to be a characteristic feature of the nucleotide site of the wild-type enzyme, too, both in its inactive open conformation and in the active closed conformation. This mobility of the nucleotide phosphates that is regulated by the aspartate side chains of D218 and D374 through the complexing Mg(2+) is suggested to be essential in enzyme function.

Ligand chirality effects on the dynamics of human 3-phosphoglycerate kinase: comparison between D- and L-nucleotides.

l-nucleoside analogues are now largely used as antiviral drugs for the treatment of viral infections like HBV, HCV and HIV. However, in order to be fully active, they need to be phosphorylated by cellular or viral kinases. Human 3-phosphogycerate kinase (hPGK) was shown to catalyze the last step of activation of l-enantiomers and thus constitutes an attractive target for theoretical predictions of its phosphorylation efficiency. Molecular dynamics simulations were carried out with four different nucleotides (d-/l-ADP and d-/l-CDP) in complex with hPGK and 1,3-bisphospho-d-glycerate (bPG). The binding affinities of CDPs (both enantiomers) for hPGK were found very weak while d- and l-ADP were better substrates. Interestingly, the binding affinity of the bPG substrate was found to be lower in presence of d-ADP than l-ADP which indicates a potential antagonistic effect on one substrate to the other. A detailed analysis of the simulations unravels important dynamic conditions for efficient phosphorylation. Indeed, as previously described for the natural substrate, the hinge bending motion of the domains upon substrates binding should be more correlated and directional. Interestingly, the unforeseen finding was the larger dynamics freedoms observed for the substrates that was favored by the protein atoms flexibility around the nucleobase binding site.

Substrate binding modifies the hinge bending characteristics of human 3-phosphoglycerate kinase: a molecular dynamics study.

3-Phosphogycerate kinase (PGK) is a two domain enzyme, with a binding site of the 1,3-bisphosphoglycerate on the N-domain and of the ADP on the C-domain. To transfer a phosphate group the enzyme has to undergo a hinge bending motion from open to closed conformation to bring the substrates to close proximity. Molecular dynamics simulation was used to elucidate the effect of ligand binding onto the domain motions of this enzyme. The simulation results of the apo form indicate a hinge bending motion in the ns timescale while the time period of the hinge bending motion of the complex form is clearly over the 20 ns simulation time. The apo form exhibits several hinge points that contribute to the hinge bending motion while upon binding the ligands, the hinge bending becomes strictly restrained with one dominant hinge point in the vicinity of the substrates. At the same time, ligand binding results in an enhanced correlation of internal domain motions.