Molecular Dynamics Simulations on the Mechanism of Transporting Methylamine and Ammonia by Ammonium Transporter AmtB

A protet-based, protonic charge transfer model of energy coupling in oxidative and photosynthetic phosphorylation

Textbooks of biochemistry will explain that the otherwise endergonic reactions of ATP synthesis can be driven by the exergonic reactions of respiratory electron transport, and that these two half-reactions are catalyzed by protein complexes embedded in the same, closed membrane. These views are correct. The textbooks also state that, according to the chemiosmotic coupling hypothesis, a (or the) kinetically and thermodynamically competent intermediate linking the two half-reactions is the electrochemical difference of protons that is in equilibrium with that between the two bulk phases that the coupling membrane serves to separate. This gradient consists of a membrane potential term Δψ and a pH gradient term ΔpH, and is known colloquially as the protonmotive force or pmf. Artificial imposition of a pmf can drive phosphorylation, but only if the pmf exceeds some 150-170mV; to achieve in vivo rates the imposed pmf must reach 200mV. The key question then is 'does the pmf generated by electron transport exceed 200mV, or even 170mV?' The possibly surprising answer, from a great many kinds of experiment and sources of evidence, including direct measurements with microelectrodes, indicates it that it does not. Observable pH changes driven by electron transport are real, and they control various processes; however, compensating ion movements restrict the Δψ component to low values. A protet-based model, that I outline here, can account for all the necessary observations, including all of those inconsistent with chemiosmotic coupling, and provides for a variety of testable hypotheses by which it might be refined.

Design and development of high affinity dual anticancer peptide-inhibitors against p53-MDM2/X interaction

AIMS

Inhibition of P53-MDM2/X interaction is known as an effective cancer therapy strategy. In this regard, pDI peptide was introduced previously with the potential of targeting MDM2. In this research, the large-scale peptide mutation screening was used to achieve the best sequence of pDI with the highest affinity for inhibition activity against MDM2/X.

MAIN METHODS

Three mutant peptides of pDI as dual inhibitor peptides including single mutations of pDIm/4W, pDIm/11M and double mutations of pDIdm/4W11M were presented with the high affinities to inhibit both MDM2/X. The selected mutants were then evaluated comprehensively to confirm their ability as potent MDM2/X inhibitors, using a theoretical simulation approach.

KEY FINDINGS

MD simulations analyses confirmed their dual inhibition potential against both MDM2/X interactions with p53 protein. The developed pDIm and mainly pDIdm peptides showed stable conformations over the simulation time with conserved secondary structure and effective interaction with MDM2/X by physical binding such as hydrogen bonding. Besides, umbrella sampling free energy calculation indicated higher binding energy, ΔG_binding, of pDIm-MDM2/X and pDIdm-MDM2/X compared to pDI-MDM2/X.

SIGNIFICANCE

The optimized and improved mutant pDI, pDIdm, with more effective ΔG_binding values of -30 and -25 kcal/mol to MDMX and MDM2, respectively, is recommended as a promising anticancer agent and suitable candidate for experimental evaluations.

Impact of Fixed Nitrogen Availability on Dehalococcoides mccartyi Reductive Dechlorination Activity

Biostimulation to promote reductive dechlorination is widely practiced, but the value of adding an exogenous nitrogen (N) source (e.g., NH₄⁺) during treatment is unclear. This study investigates the effect of NH₄⁺ availability on organohalide-respiring Dehalococcoides mccartyi (Dhc) growth and reductive dechlorination in enrichment cultures derived from groundwater (PW4) and river sediment (TC) impacted with chlorinated ethenes. In PW4 cultures, the addition of NH₄⁺ increased cis-1,2-dichloroethene (cDCE)-to-ethene dechlorination rates about 5-fold (20.6 ± 1.6 versus 3.8 ± 0.5 μM Cl^- d^-1), and the total number of Dhc 16S rRNA gene copies were about 43-fold higher in incubations with NH₄⁺ ((1.8 ± 0.9) × 10⁸ mL^-1) compared to incubations without NH₄⁺ ((4.1 ± 0.8) × 10⁷ mL^-1). In TC cultures, NH₄⁺ also stimulated cDCE-to-ethene dechlorination and Dhc growth. Quantitative polymerase chain reaction (qPCR) revealed that Cornell-type Dhc capable of N₂ fixation dominated PW4 cultures without NH₄⁺, but their relative abundance decreased in cultures with NH₄⁺ amendment (i.e., 99 versus 54% of total Dhc). Pinellas-type Dhc incapable of N₂ fixation were responsible for cDCE dechlorination in TC cultures, and diazotrophic community members met their fixed N requirement in the medium without NH₄⁺. Responses to NH₄⁺ were apparent at the community level, and N₂-fixing bacterial populations increased in incubations without NH₄⁺. Quantitative assessment of Dhc nitrogenase genes, transcripts, and proteomics data linked Cornell-type Dhc nifD and nifK expression with fixed N limitation. NH₄⁺ additions also demonstrated positive effects on Dhc in situ dechlorination activity in the vicinity of well PW4. These findings demonstrate that biostimulation with NH₄⁺ can enhance Dhc reductive dechlorination rates; however, a "do nothing" approach that relies on indigenous diazotrophs can achieve similar dechlorination end points and avoids the potential for stalled dechlorination due to inhibitory levels of NH₄⁺ or transformation products (i.e., nitrous oxide).

A pore-occluding phenylalanine gate prevents ion slippage through plant ammonium transporters

Throughout all kingdoms of life, highly conserved transport proteins mediate the passage of ammonium across membranes. These transporters share a high homology and a common pore structure. Whether NH₃, NH₄⁺ or NH₃ + H⁺ is the molecularly transported substrate, still remains unclear for distinct proteins. High-resolution protein structures of several ammonium transporters suggested two conserved pore domains, an external NH₄⁺ recruitment site and a pore-occluding twin phenylalanine gate, to take over a crucial role in substrate determination and selectivity. Here, we show that while the external recruitment site seems essential for AtAMT1;2 function, single mutants of the double phenylalanine gate were not reduced in their ammonium transport capacity. Despite an unchanged ammonium transport rate, a single mutant of the inner phenylalanine showed reduced N-isotope selection that was proposed to be associated with ammonium deprotonation during transport. Even though ammonium might pass the mutant AMT pore in the ionic form, the transporter still excluded potassium ions from being transported. Our results, highlight the importance of the twin phenylalanine gate in blocking uncontrolled ammonium ion flux.

Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation

Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.

Probing Difference in Binding Modes of Inhibitors to MDMX by Molecular Dynamics Simulations and Different Free Energy Methods

The p53-MDMX interaction has attracted extensive attention of anti-cancer drug development in recent years. This current work adopted molecular dynamics (MD) simulations and cross-correlation analysis to investigate conformation changes of MDMX caused by inhibitor bindings. The obtained information indicates that the binding cleft of MDMX undergoes a large conformational change and the dynamic behavior of residues obviously change by the presence of different structural inhibitors. Two different methods of binding free energy predictions were employed to carry out a comparable insight into binding mechanisms of four inhibitors PMI, pDI, WK23 and WW8 to MDMX. The data show that the main factor controlling the inhibitor bindings to MDMX arises from van der Waals interactions. The binding free energies were further divided into contribution of each residue and the derived information gives a conclusion that the hydrophobic interactions, such as CH-CH, CH-π and π-π interactions, are responsible for the inhibitor associations with MDMX.

A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation

Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.

Revealing origin of decrease in potency of darunavir and amprenavir against HIV-2 relative to HIV-1 protease by molecular dynamics simulations

Clinical inhibitors Darunavir (DRV) and Amprenavir (APV) are less effective on HIV-2 protease (PR2) than on HIV-1 protease (PR1). To identify molecular basis associated with the lower inhibition, molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations were performed to investigate the effectiveness of the PR1 inhibitors DRV and APV against PR1/PR2. The rank of predicted binding free energies agrees with the experimental determined one. Moreover, our results show that two inhibitors bind less strongly to PR2 than to PR1, again in agreement with the experimental findings. The decrease in binding free energies for PR2 relative to PR1 is found to arise from the reduction of the van der Waals interactions induced by the structural adjustment of the triple mutant V32I, I47V and V82I. This result is further supported by the difference between the van der Waals interactions of inhibitors with each residue in PR2 and in PR1. The results from the principle component analysis suggest that inhibitor binding tends to make the flaps of PR2 close and the one of PR1 open. We expect that this study can theoretically provide significant guidance and dynamics information for the design of potent dual inhibitors targeting PR1/PR2.

Binding modes of three inhibitors 8CA, F8A and I4A to A-FABP studied based on molecular dynamics simulation

Adipocyte fatty-acid binding protein (A-FABP) is an important target of drug designs treating some diseases related to lipid-mediated biology. Molecular dynamics (MD) simulations coupled with solvated interaction energy method (SIE) were carried out to study the binding modes of three inhibitors 8CA, F8A and I4A to A-FABP. The rank of our predicted binding affinities is in accordance with experimental data. The results show that the substitution in the position 5 of N-benzyl and the seven-membered ring of N-benzyl-indole carboxylic acids strengthen the I4A binding, while the substitution in the position 2 of N-benzyl weakens the F8A binding. Computational alanine scanning and dynamics analyses were performed and the results suggest that the polar interactions of the positively charged residue R126 with the three inhibitors provide a significant contribution to inhibitor bindings. This polar interaction induces the disappearance of the correlated motion of the C terminus of A-FABP relative to the N terminus and favors the stability of the binding complex. This study is helpful for the rational design of potent inhibitors within the fields of metabolic disease, inflammation and atherosclerosis.

Epinephrine depletion exacerbates the fasting-induced protein breakdown in fast-twitch skeletal muscles

The physiological role of epinephrine in the regulation of skeletal muscle protein metabolism under fasting is unknown. We examined the effects of plasma epinephrine depletion, induced by adrenodemedullation (ADMX), on muscle protein metabolism in fed and 2-day-fasted rats. In fed rats, ADMX for 10 days reduced muscle mass, the cross-sectional area of extensor digitorum longus (EDL) muscle fibers, and the phosphorylation levels of Akt. In addition, ADMX led to a compensatory increase in muscle sympathetic activity, as estimated by the rate of norepinephrine turnover; this increase was accompanied by high rates of muscle protein synthesis. In fasted rats, ADMX exacerbated fasting-induced proteolysis in EDL but did not affect the low rates of protein synthesis. Accordingly, ADMX activated lysosomal proteolysis and further increased the activity of the ubiquitin (Ub)-proteasome system (UPS). Moreover, expression of the atrophy-related Ub ligases atrogin-1 and MuRF1 and the autophagy-related genes LC3b and GABARAPl1 were upregulated in EDL muscles from ADMX-fasted rats compared with sham-fasted rats, and ADMX reduced cAMP levels and increased fasting-induced Akt dephosphorylation. Unlike that observed for EDL muscles, soleus muscle proteolysis and Akt phosphorylation levels were not affected by ADMX. In isolated EDL, epinephrine reduced the basal UPS activity and suppressed overall proteolysis and atrogin-1 and MuRF1 induction following fasting. These data suggest that epinephrine released from the adrenal medulla inhibits fasting-induced protein breakdown in fast-twitch skeletal muscles, and these antiproteolytic effects on the UPS and lysosomal system are apparently mediated through a cAMP-Akt-dependent pathway, which suppresses ubiquitination and autophagy.

Nitrogen assimilation in Escherichia coli: putting molecular data into a systems perspective

We present a comprehensive overview of the hierarchical network of intracellular processes revolving around central nitrogen metabolism in Escherichia coli. The hierarchy intertwines transport, metabolism, signaling leading to posttranslational modification, and transcription. The protein components of the network include an ammonium transporter (AmtB), a glutamine transporter (GlnHPQ), two ammonium assimilation pathways (glutamine synthetase [GS]-glutamate synthase [glutamine 2-oxoglutarate amidotransferase {GOGAT}] and glutamate dehydrogenase [GDH]), the two bifunctional enzymes adenylyl transferase/adenylyl-removing enzyme (ATase) and uridylyl transferase/uridylyl-removing enzyme (UTase), the two trimeric signal transduction proteins (GlnB and GlnK), the two-component regulatory system composed of the histidine protein kinase nitrogen regulator II (NRII) and the response nitrogen regulator I (NRI), three global transcriptional regulators called nitrogen assimilation control (Nac) protein, leucine-responsive regulatory protein (Lrp), and cyclic AMP (cAMP) receptor protein (Crp), the glutaminases, and the nitrogen-phosphotransferase system. First, the structural and molecular knowledge on these proteins is reviewed. Thereafter, the activities of the components as they engage together in transport, metabolism, signal transduction, and transcription and their regulation are discussed. Next, old and new molecular data and physiological data are put into a common perspective on integral cellular functioning, especially with the aim of resolving counterintuitive or paradoxical processes featured in nitrogen assimilation. Finally, we articulate what still remains to be discovered and what general lessons can be learned from the vast amounts of data that are available now.

A computational analysis of binding modes and conformation changes of MDM2 induced by p53 and inhibitor bindings

Molecular dynamics (MD) simulations followed by principal component analysis were performed to study the conformational change of MDM2 induced by p53 and two inhibitor (P4 and MI63a) bindings. The results show that the hydrophobic cleft of MDM2 is very flexible and adaptive to different structural binding partners. The cleft tends to become wider and more stable as MDM2 binds to the three binding partners, while unbound MDM2 shows a narrower and pretty flexible cleft, which agrees with recent experimental data and theoretical studies. It was also found that the binding of P4 and p53 stabilizes the motion of the loop L2 linking the helix α2 and β strand (β3), but the presence of MI63a makes the motion of L2 disordered. In addition, the binding free energies of the three partners to MDM2 were calculated using molecular mechanics generalized Born surface area to explain the binding modes of these three partners to MDM2. This study will be helpful not only for better understanding the functional, concerted motion of MDM2, but also for the rational design of potent anticancer drugs targeting the p53-MDM2 interaction.

Ammonium transport proteins with changes in one of the conserved pore histidines have different performance in ammonia and methylamine conduction

Two conserved histidine residues are located near the mid-point of the conduction channel of ammonium transport proteins. The role of these histidines in ammonia and methylamine transport was evaluated by using a combination of in vivo studies, molecular dynamics (MD) simulation, and potential of mean force (PMF) calculations. Our in vivo results showed that a single change of either of the conserved histidines to alanine leads to the failure to transport methylamine but still facilitates good growth on ammonia, whereas double histidine variants completely lose their ability to transport both methylamine and ammonia. Molecular dynamics simulations indicated the molecular basis of the in vivo observations. They clearly showed that a single histidine variant (H168A or H318A) of AmtB confines the rather hydrophobic methylamine more strongly than ammonia around the mutated sites, resulting in dysfunction in conducting the former but not the latter molecule. PMF calculations further revealed that the single histidine variants form a potential energy well of up to 6 kcal/mol for methylamine, impairing conduction of this substrate. Unlike the single histidine variants, the double histidine variant, H168A/H318A, of AmtB was found to lose its unidirectional property of transporting both ammonia and methylamine. This could be attributed to a greatly increased frequency of opening of the entrance gate formed by F215 and F107, in this variant compared to wild-type, with a resultant lowering of the energy barrier for substrate to return to the periplasm.

Thermodynamics of transport through the ammonium transporter Amt-1 investigated with free energy calculations

Amt-1 from Archaeoglobus fulgidus (AfAmt-1) belongs to the Amt/Rh family of ammonium/ammonia transporting membrane proteins. The transport mode and the precise microscopic permeation mechanism utilized by these proteins are intensely debated. Open questions concern the identity of the transported substrate (ammonia and/or ammonium) and whether the transport is passive or active. To address these questions, we studied the overall thermodynamics of the different transport modes as a function of the environmental conditions. Then, we investigated the thermodynamics of the underlying microscopic transport mechanisms with free energy calculations within a continuum electrostatics model. The formalism developed for this purpose is of general utility in the calculation of binding free energies for ligands with multiple protonation forms or other binding forms. The results of our calculations are compared to the available experimental and theoretical data on Amt/Rh proteins and discussed in light of the current knowledge on the physiological conditions experienced by microorganisms and plants. We found that microscopic models of electroneutral and electrogenic transport modes are in principle thermodynamically viable. However, only the electrogenic variants have a net thermodynamic driving force under the physiological conditions experienced by microorganisms and plants. Thus, the transport mechanism of AfAmt-1 is most likely electrogenic.