Integrating Multiple Accelerated Molecular Dynamics To Improve Accuracy of Free Energy Calculations

Mutations influence the conformational dynamics of the GDP/KRAS complex

Mutations near allosteric sites can have a significant impact on the function of KRAS. Three specific mutations, K104Q, G12D/K104Q, and G12D/G75A, which are located near allosteric positions, were selected to investigate the molecular mechanisms behind mutation-induced influences on the activity of KRAS. Gaussian accelerated molecular dynamics (GaMD) simulations followed by the principal component analysis (PCA) were performed to improve the sampling of conformational states. The results revealed that these mutations significantly alter the structural flexibility, correlated motions, and dynamic behavior of the switch regions that are essential for KRAS binding to effectors or regulators. Furthermore, the mutations have a significant impact on the hydrogen bonding interactions between GDP and the switch regions, as well as on the electrostatic interactions of magnesium ions (Mg2+) with these regions. Our results verified that these mutations strongly influence the binding of KRAS to its effectors or regulators and allosterically regulate the activity. We believe that this work can provide valuable theoretical insights into a deeper understanding of KRAS function.Communicated by Ramaswamy H. Sarma.

Roles of Accelerated Molecular Dynamics Simulations in Predictions of Binding Kinetic Parameters

Rational predictions on binding kinetics parameters of drugs to targets play significant roles in future drug designs. Full conformational samplings of targets are requisite for accurate predictions of binding kinetic parameters. In this review, we mainly focus on the applications of enhanced sampling technologies in calculations of binding kinetics parameters and residence time of drugs. The methods involved in molecular dynamics simulations are applied to not only probe conformational changes of targets but also reveal calculations of residence time that is significant for drug efficiency. For this review, special attention are paid to accelerated molecular dynamics (aMD) and Gaussian aMD (GaMD) simulations that have been adopted to predict the association or disassociation rate constant. We also expect that this review can provide useful information for future drug design.

Decoding the Identification Mechanism of an SAM-III Riboswitch on Ligands through Multiple Independent Gaussian-Accelerated Molecular Dynamics Simulations

S-Adenosyl-l-methionine (SAM)-responsive riboswitches play a central role in the regulation of bacterial gene expression at the level of transcription attenuation or translation inhibition. In this study, multiple independent Gaussian-accelerated molecular dynamics simulations were performed to decipher the identification mechanisms of SAM-III (S_MK) on ligands SAM, SAH, and EEM. The results reveal that ligand binding highly affects the structural flexibility, internal dynamics, and conformational changes of SAM-III. The dynamic analysis shows that helices P3 and P4 as well as two junctions J23 and J24 of SAM-III are highly susceptible to ligand binding. Analyses of free energy landscapes suggest that ligand binding induces different free energy profiles of SAM-III, which leads to the difference in identification sites of SAM-III on ligands. The information on ligand-nucleotide interactions not only uncovers that the π-π, cation-π, and hydrogen bonding interactions drive identification of SAM-III on the three ligands but also reveals that different electrostatic properties of SAM, SAH, and EEM alter the active sites of SAM-III. Meanwhile, the results also verify that the adenine group of SAM, SAH, and EEM is well recognized by conserved nucleotides G7, A29, U37, A38, and G48. We expect that this study can provide useful information for understanding the applications of SAM-III in chemical, synthetic RNA biology, and biomedical fields.

Enhancing sampling with free-energy calculations

In recent years, considerable progress has been made to enhance sampling and help address biological questions, including, but not limited to conformational transitions in biomolecules and protein-ligand reversible binding, hitherto intractable by brute-force computer simulations. Many of these advances result from the development of a palette of methods aimed at exploring rare events through reliable free-energy calculations. The advent of new, often conceptually related methods has also rendered difficult the choice of the best suited option for a given problem. Here, we focus on geometrical transformations and algorithms designed to enhance sampling along adequately chosen progress variables, tracing their theoretical foundations, and showing how they are connected and can be blended together for improved performance.

Enhanced sampling without borders: on global biasing functions and how to reweight them

Molecular dynamics (MD) simulations are a powerful tool to follow the time evolution of biomolecular motions in atomistic resolution. However, the high computational demand of these simulations limits the timescales of motions that can be observed. To resolve this issue, so called enhanced sampling techniques are developed, which extend conventional MD algorithms to speed up the simulation process. Here, we focus on techniques that apply global biasing functions. We provide a broad overview of established enhanced sampling methods and promising new advances. As the ultimate goal is to retrieve unbiased information from biased ensembles, we also discuss benefits and limitations of common reweighting schemes. In addition to concisely summarizing critical assumptions and implications, we highlight the general application opportunities as well as uncertainties of global enhanced sampling.

Q61 mutant-mediated dynamics changes of the GTP-KRAS complex probed by Gaussian accelerated molecular dynamics and free energy landscapes

Understanding the molecular mechanism of the GTP-KRAS binding is significant for improving the target roles of KRAS in cancer treatment. In this work, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations were applied to decode the effect of Q61A, Q61H and Q61L on the activity of KRAS. Dynamics analyses based on MR-GaMD trajectory reveal that motion modes and dynamics behavior of the switch domain in KRAS are heavily affected by the three Q61 mutants. Information of free energy landscapes (FELs) shows that Q61A, Q61H and Q61L induce structural disorder of the switch domain and disturb the activity of KRAS. Analysis of the interaction network uncovers that the decrease in the stability of hydrogen bonding interactions (HBIs) of GTP with residues V29 and D30 induced by Q61A, Q61H and Q61L is responsible for the structural disorder of the switch-I and that in the occupancy of the hydrogen bond between GTP and residue G60 leads to the structural disorder of the switch-II. Thus, the high disorder of the switch domain caused by three current Q61 mutants produces a significant effect on binding of KRAS to its effectors. This work is expected to provide useful information for further understanding function and target roles of KRAS in anti-cancer drug development.

Computer-aided design of molecular machines: techniques, paradigms and difficulties

With their development in the past decade, molecular machines, which achieve specific tasks by responding to external stimuli, have gradually come to be regarded as powerful tools for a wide range of applications, rather than interesting molecular toys. This conceptual change in turn motivates scientists to design molecular machines with complex architectures. Due to the lack of general principles bridging the functions and the chemical structures of molecular machines, experience-based design becomes difficult with the increase of size and complexity of the architectures. Computer-aided molecular-machine design, therefore, has attracted widespread attention on account of its ability to model and investigate complex molecular architectures without too much time and expense required for synthetic experiments. Using leading-edge numerical-simulation techniques, the mechanisms underlying achieving tasks through response to external stimuli of a large number of existing molecular machines have been successfully explored. Based on the experience of studying existing molecular machines, generalized methodologies of predicting the properties and working principles of molecular candidates have been established, paving the way for de novo computer-aided design of molecular machines. In this perspective, we introduce cutting-edge techniques that have been applied for investigating and designing molecular machines. We show paradigms of computer-aided design of molecular machines, which can serve as guidelines for the investigation of new supramolecular architectures. Moreover, we discuss the limitations and possible future developments of current techniques and methodologies in the field of computer-aided design of molecular machines.

Advanced Sampling Methods for Multiscale Simulation of Disordered Proteins and Dynamic Interactions

Intrinsically disordered proteins (IDPs) are highly prevalent and play important roles in biology and human diseases. It is now also recognized that many IDPs remain dynamic even in specific complexes and functional assemblies. Computer simulations are essential for deriving a molecular description of the disordered protein ensembles and dynamic interactions for a mechanistic understanding of IDPs in biology, diseases, and therapeutics. Here, we provide an in-depth review of recent advances in the multi-scale simulation of disordered protein states, with a particular emphasis on the development and application of advanced sampling techniques for studying IDPs. These techniques are critical for adequate sampling of the manifold functionally relevant conformational spaces of IDPs. Together with dramatically improved protein force fields, these advanced simulation approaches have achieved substantial success and demonstrated significant promise towards the quantitative and predictive modeling of IDPs and their dynamic interactions. We will also discuss important challenges remaining in the atomistic simulation of larger systems and how various coarse-grained approaches may help to bridge the remaining gaps in the accessible time- and length-scales of IDP simulations.

Overcoming Free-Energy Barriers with a Seamless Combination of a Biasing Force and a Collective Variable-Independent Boost Potential

Amid collective-variable (CV)-based importance-sampling algorithms, a hybrid of the extended adaptive biasing force and the well-tempered metadynamics algorithms (WTM-eABF) has proven particularly cost-effective for exploring the rugged free-energy landscapes that underlie biological processes. However, as an inherently CV-based algorithm, this hybrid scheme does not explicitly accelerate sampling in the space orthogonal to the chosen CVs, thereby limiting its efficiency and accuracy, most notably in those cases where the slow degrees of freedom of the process at hand are not accounted for in the model transition coordinate. Here, inspired by Gaussian-accelerated molecular dynamics (GaMD), we introduce the same CV-independent harmonic boost potential into WTM-eABF, yielding a hybrid algorithm coined GaWTM-eABF. This algorithm leans on WTM-eABF to explore the transition coordinate with a GaMD-mollified potential and recovers the unbiased free-energy landscape through thermodynamic integration followed by proper reweighting. As illustrated in our numerical tests, GaWTM-eABF effectively overcomes the free-energy barriers in orthogonal space and correctly recovers the unbiased potential of mean force (PMF). Furthermore, applying both GaWTM-eABF and WTM-eABF to two biologically relevant processes, namely, the reversible folding of (i) deca-alanine and (ii) chignolin, our results indicate that GaWTM-eABF reduces the uncertainty in the PMF calculation and converges appreciably faster than WTM-eABF. Obviating the need of multiple-copy strategies, GaWTM-eABF is a robust, computationally efficient algorithm to surmount the free-energy barriers in orthogonal space and maps with utmost fidelity the free-energy landscape along selections of CVs. Moreover, our strategy that combines WTM-eABF with GaMD can be easily extended to other biasing-force algorithms.

Binding mechanism of inhibitors to p38α MAP kinase deciphered by using multiple replica Gaussian accelerated molecular dynamics and calculations of binding free energies

The p38α MAP Kinase has been an important target of drug design for treatment of inflammatory diseases and cancers. This work applies multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations and the molecular mechanics generalized Born surface area (MM-GBSA) method to probe the binding mechanism of inhibitors L51, R24 and 1AU to p38α. Dynamics analyses show that inhibitor bindings exert significant effect on conformational changes of the active helix α2 and the conserved DFG loop. The rank of binding free energies calculated with MM-GBSA not only agrees well with that determined by the experimental IC50 values but also suggests that mutual compensation between the enthalpy and entropy changes can improve binding of inhibitors to p38α. The analyses of free energy landscapes indicate that the L51, R24 and 1AU bound p38α display a DFG-out conformation. The residue-based free energy decomposition method is used to evaluate contributions of separate residues to the inhibitor-p38α binding and the results imply that residues V30, V38, L74, L75, I84, T106, H107, L108, M109, L167, F169 and D168 can be utilized as efficient targets of potent inhibitors toward p38α.

Inhibition mechanism and hot-spot prediction of nine potential drugs for SARS-CoV-2 Mpro by large-scale molecular dynamic simulations combined with accurate binding free energy calculations

Coronavirus disease 2019 (COVID-19), which is caused by a new coronavirus known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is spreading around the world. However, a universally effective treatment regimen has not been developed to date. The main protease (M^pro), a key enzyme of SARS-CoV-2, plays a crucial role in the replication and transcription of this virus in cells and has become the ideal target for rational antiviral drug design. In this study, we performed molecular dynamics simulations three times for these complexes of M^pro (monomeric and dimeric) and nine potential drugs that have a certain effect on the treatment of COVID-19 to explore their binding mechanism. In addition, a total of 12 methods for calculating binding free energy were employed to determine the optimal drug. Ritonavir, Arbidol, and Chloroquine consistently showed an outstanding binding ability to monomeric M^pro under various methods. Ritonavir, Arbidol, and Saquinavir presented the best performance when binding to a dimer, which was independent of the protonated state of Hie41 (protonated at N_ε) and Hid41 (protonated at N_δ), and these findings suggest that Chloroquine may not effectively inhibit the activity of dimeric M^proin vivo. Furthermore, three common hot-spot residues of Met165, Hie41, and Gln189 of monomeric M^pro systems dominated the binding of Ritonavir, Arbidol, and Chloroquine. In dimeric M^pro, Gln189, Met165, and Met49 contributed significantly to binding with Ritonavir, Arbidol, and Saquinavir; therefore, Gln189 and Met165 might serve as the focus in the discovery and development of anti-COVID-19 drugs. In addition, the van der Waals interaction played a significant role in the binding process, and the benzene ring of the drugs showed an apparent inhibitory effect on the normal function of M^pro. The binding cavity had great flexibility to accommodate these different drugs. The results would be notably helpful for enabling a detailed understanding of the binding mechanisms for these important drug-M^pro interactions and provide valuable guidance for the design of potent inhibitors.

Quality of force fields and sampling methods in simulating pepX peptides: a case study for intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are a group of proteins that lack well-defined structures under native conditions and carry out crucial physiological functions in various biochemical pathways. Due to the heterogeneous nature of IDPs, molecular dynamics simulations have been extensively adopted to investigate the conformational ensembles and dynamic properties of these proteins. However, their accuracy remains limited by the development of force fields and sampling algorithms. Here, we evaluated the quality of both force fields and enhanced sampling algorithms based on five short pepX peptides. Our results show that the more extended conformational ensembles sampled by the AMOEBA polarizable force field present a higher ability to reproduce experimental NMR observables than AMBER and CHARMM classical force fields. Moreover, a better agreement with experiments is achieved in the simulation of IaMD (integrated accelerated molecular dynamics) than in aMD (accelerated molecular dynamics). The results together indicate that the combination of AMOEBA force field and IaMD enhanced sampling might be a better choice for simulating IDPs. This work may provide important clues for developments and applications of force fields and enhanced sampling methods in future simulations of IDPs.

Adaptively Iterative Multiscale Switching Simulation Strategy and Applications to Protein Folding and Structure Prediction

Structure prediction is an important means to quickly understand new protein functions. However, the prediction of effects of proteins that have no detectable templates is still to be improved. Molecular dynamics simulation is supposed to be the primary research tool for structure predictions, but it still has limitations of huge computational cost in all-atom (AA) models and rough accuracy in coarse-grained (CG) models. We propose a universal multiscale simulation strategy named AIMS in which simulations can iteratively switch among multiple resolutions in order to adaptively trade off AA accuracy and CG high-efficiency. AIMS follows the idea of CG-guided enhanced sampling so that final results always keep AA accuracy. We successfully achieve four ab initio and four data-assisted protein structure predictions using AIMS. The prediction result is an ensemble rather than a structure and provides special insights on folding metastable states. AIMS is estimated to achieve a computational speed about 40 times faster than that of conventional AA simulations.

Conformational Changes of Glutamine 5'-Phosphoribosylpyrophosphate Amidotransferase for Two Substrates Analogue Binding: Insight from Conventional Molecular Dynamics and Accelerated Molecular Dynamics Simulations

Glutamine 5′-phosphoribosylpyrophosphate amidotransferase (GPATase) catalyzes the synthesis of phosphoribosylamine, pyrophosphate, and glutamate from phosphoribosylpyrophosphate, as well as glutamine at two sites (i.e., glutaminase and phosphoribosylpyrophosphate sites), through a 20 Å NH₃ channel. In this study, conventional molecular dynamics (cMD) simulations and enhanced sampling accelerated molecular dynamics (aMD) simulations were integrated to characterize the mechanism for coordination catalysis at two separate active sites in the enzyme. Results of cMD simulations illustrated the mechanism by which two substrate analogues, namely, DON and cPRPP, affect the structural stability of GPATase from the perspective of dynamic behavior. aMD simulations obtained several key findings. First, a comparison of protein conformational changes in the complexes of GPATase–DON and GPATase–DON–cPRPP showed that binding cPRPP to the PRTase flexible loop (K326 to L350) substantially effected the formation of the R73-DON salt bridge. Moreover, only the PRTase flexible loop in the GPATase–DON–cPRPP complex could remain closed and had sufficient space for cPRPP binding, indicating that binding of DON to the glutamine loop had an impact on the PRTase flexible loop. Finally, both DON and cPRPP tightly bonded to the two domains, thereby inducing the glutamine loop and the PRTase flexible loop to move close to each other. This movement facilitated the transfer of NH3 via the NH3 channel. These theoretical results are useful to the ongoing research on efficient inhibitors related to GPATase.

Unveiling conformational dynamics changes of H-Ras induced by mutations based on accelerated molecular dynamics

Uncovering molecular basis with regard to the conformational change of two switches I and II in the GppNHp (GNP)-bound H-Ras is highly significant for the understanding of Ras signaling. For this purpose, accelerated molecular dynamics (aMD) simulations and principal component (PC) analysis are integrated to probe the effect of mutations G12V, T35S and Q61K on conformational transformation between two switches of the GNP-bound H-Ras. The RMSF and cross-correlation analyses suggest that three mutations exert a vital effect on the flexibility and internal dynamics of two switches in the GNP-bound H-Ras. The results stemming from PC analysis indicate that two switches in the GNP-bound WT H-Ras tend to form a closed state in most conformations, while those in the GNP-bound mutated H-Ras display transformation between different states. This conclusion is further supported by free energy landscapes constructed by using the distances of residues 12 away from 35 and 35 away from 61 as reaction coordinates and different experimental studies. Interaction scanning is performed on aMD trajectories and the information shows that conformational transformations of two switches I and II induced by mutations extremely affect the GNP-residue interactions. Meanwhile, the scanning results also signify that residues G15, A18, F28, K117, A146 and K147 form stable contacts with GNP, while residues D30, E31, Y32, D33, P34 and E62 in two switches I and II produce unstable contacts with GNP. This study not only reveals dynamic behavior changes of two switches in H-Ras induced by mutations, but also unveils general principles and mechanisms with regard to functional conformational changes of H-Ras.

Making NSCLC Crystal Clear: How Kinase Structures Revolutionized Lung Cancer Treatment

The parallel advances of different scientific fields provide a contemporary scenario where collaboration is not a differential, but actually a requirement. In this context, crystallography has had a major contribution on the medical sciences, providing a “face” for targets of diseases that previously were known solely by name or sequence. Worldwide, cancer still leads the number of annual deaths, with 9.6 million associated deaths, with a major contribution from lung cancer and its 1.7 million deaths. Since the relationship between cancer and kinases was unraveled, these proteins have been extensively explored and became associated with drugs that later attained blockbuster status. Crystallographic structures of kinases related to lung cancer and their developed and marketed drugs provided insight on their conformation in the absence or presence of small molecules. Notwithstanding, these structures were also of service once the initially highly successful drugs started to lose their effectiveness in the emergence of mutations. This review focuses on a subclassification of lung cancer, non-small cell lung cancer (NSCLC), and major oncogenic driver mutations in kinases, and how crystallographic structures can be used, not only to provide awareness of the function and inhibition of these mutations, but also how these structures can be used in further computational studies aiming at addressing these novel mutations in the field of personalized medicine.

Arbitrary Resolution with Two Bead Types Coarse-Grained Strategy and Applications to Protein Recognition

Molecular recognition is a fundamental step in essentially any biological process. However, the kinetic processes during association and dissociation are difficult to be efficiently sampled by direct all-atom molecular dynamics simulations because of the large spatial and temporal scales. Here we propose an arbitrary resolution with two bead types (ART) coarse-grained (CG) strategy that is adept in molecular recognition. ART is a universal user-customized CG strategy that can generate a system-specific CG force field anytime and be applied to any system with an arbitrary CG resolution according to research requirements. ART CG simulations can be very efficiently performed with implicit solvation in prevalent simulation packages and provide interfaces for any enhanced sampling method. We used three applications, HLA-HIV epitope recognition, barnase-barstar association, and trimeric TRAF2 self-assembly, to validate the feasibility of the ART CG strategy, its advantages in protein recognition, and its high performance in simulations. Regular CG simulations can successfully achieve valid protein recognitions without any prior bound structure.

Enhanced sampling and free energy calculations for protein simulations

Molecular dynamics simulation is a powerful computational technique to study biomolecular systems, which complements experiments by providing insights into the structural dynamics relevant to biological functions at atomic scale. It can also be used to calculate the free energy landscapes of the conformational transitions to better understand the functions of the biomolecules. However, the sampling of biomolecular configurations is limited by the free energy barriers that need to be overcome, leading to considerable gaps between the timescales reached by MD simulation and those governing biological processes. To address this issue, many enhanced sampling methodologies have been developed to increase the sampling efficiency of molecular dynamics simulations and free energy calculations. Usually, enhanced sampling algorithms can be classified into methods based on collective variables (CV-based) and approaches which do not require predefined CVs (CV-free). In this chapter, the theoretical basis of free energy estimation is briefly reviewed first, followed by the reviews of the most common CV-based and CV-free methods including the presentation of some examples and recent developments. Finally, the combination of different enhanced sampling methods is discussed.

Large-Scale Biomolecular Conformational Transitions Explored by a Combined Elastic Network Model and Enhanced Sampling Molecular Dynamics

Biomolecules often undergo large-scale conformational transitions when carrying out their functions. However, it is still challenging for conventional molecular dynamics simulations to provide adequate structural dynamics information to interpret associated mechanisms. Here, we present a combined elastic network model and enhanced sampling-based strategy (iterANM-IaMD) by adopting iterANM to construct initial conformation space and enhanced sampling IaMD to explore the free energy landscape along specific large-scale conformational transitions. We applied this strategy to three functionally and structurally distinct proteins (adenylate kinase, calmodulin, and p38α kinase), which undergo striking conformational change upon ligand binding. The simulation results for both free and ligand-bound proteins show qualitative and quantitative agreement with existing studies, suggesting iterANM-IaMD as an accurate and efficient tool to investigate structural dynamics involved in complicated biological processes. Our work also provides insights into the relationship between the dynamics and functionality of biomolecules.

Improving the performance of the MM/PBSA and MM/GBSA methods in recognizing the native structure of the Bcl-2 family using the interaction entropy method

Correct discrimination of native structure plays an important role in drug design. IE method significantly improves the performance of MM/PB(GB)SA method in discriminating native and decoy structures in protein–ligand/protein systems of Bcl-2 family.

Taming Rugged Free Energy Landscapes Using an Average Force

The observation of complex structural transitions in biological and abiological molecular objects within time scales amenable to molecular dynamics (MD) simulations is often hampered by significant free energy barriers associated with entangled movements. Importance-sampling algorithms, a powerful class of numerical schemes for the investigation of rare events, have been widely used to extend simulations beyond the time scale common to MD. However, probing processes spanning milliseconds through microsecond molecular simulations still constitutes in practice a daunting challenge because of the difficulty of taming the ruggedness of multidimensional free energy surfaces by means of naive transition coordinates. To address this limitation, in recent years we have elaborated importance-sampling methods relying on an adaptive biasing force (ABF). In this Account, we review recent developments of algorithms aimed at mapping rugged free energy landscapes that correspond to complex processes of physical, chemical, and biological relevance. Through these developments, we have broadened the spectrum of applications of the popular ABF algorithm while improving its computational efficiency, notably for multidimensional free energy calculations. One major algorithmic advance, coined meta-eABF, merges the key features of metadynamics and an extended Lagrangian variant of ABF (eABF) by simultaneously shaving the barriers and flooding the valleys of the free energy landscape, and it possesses a convergence rate up to 5-fold greater than those of other importance-sampling algorithms. Through faster convergence and enhanced ergodic properties, meta-eABF represents a significant step forward in the simulation of millisecond-time-scale events. Here we introduce extensions of the algorithm, notably its well-tempered and replica-exchange variants, which further boost the sampling efficiency while gaining in numerical stability, thus allowing quantum-mechanical/molecular-mechanical free energy calculations to be performed at a lower cost. As a paradigm to bridge microsecond simulations to millisecond events by means of free energy calculations, we have applied the ABF family of algorithms to decompose complex movements in molecular objects of biological and abiological nature. We show here how water lubricates the shuttling of an amide-based rotaxane by altering the mechanism that underlies the concerted translation and isomerization of the macrocycle. Introducing novel collective variables in a computational workflow for the rigorous determination of standard binding free energies, we predict with utmost accuracy the thermodynamics of protein-ligand reversible association. Because of their simplicity, versatility, and robust mathematical foundations, the algorithms of the ABF family represent an appealing option for the theoretical investigation of a broad range of problems relevant to physics, chemistry, and biology.

Structural and functional insights into the tetrameric photosystem I from heterocyst-forming cyanobacteria

Two large protein-cofactor complexes, photosystem I and photosystem II, are the central components of photosynthesis in the thylakoid membranes. Here, we report the 2.37-Å structure of a tetrameric photosystem I complex from a heterocyst-forming cyanobacterium Anabaena sp. PCC 7120. Four photosystem I monomers, organized in a dimer of dimer, form two distinct interfaces that are largely mediated by specifically orientated polar lipids, such as sulfoquinovosyl diacylglycerol. The structure depicts a more closely connected network of chlorophylls across monomer interfaces than those seen in trimeric PSI from thermophilic cyanobacteria, possibly allowing a more efficient energy transfer between monomers. Our physiological data also revealed a functional link of photosystem I oligomerization to cyclic electron flow and thylakoid membrane organization.

α-Ketoglutarate-Activated NF-κB Signaling Promotes Compensatory Glucose Uptake and Brain Tumor Development

The rapid proliferation of cancer cells and dysregulated vasculature within the tumor leads to limited nutrient accessibility. Cancer cells often rewire their metabolic pathways for adaption to nutrient stress, and the underlying mechanism remains largely unknown. Glutamate dehydrogenase 1 (GDH1) is a key enzyme in glutaminolysis that converts glutamate to α-ketoglutarate (α-KG). Here, we show that, under low glucose, GDH1 is phosphorylated at serine (S) 384 and interacts with RelA and IKKβ. GDH1-produced α-KG directly binds to and activates IKKβ and nuclear factor κB (NF-κB) signaling, which promotes glucose uptake and tumor cell survival by upregulating GLUT1, thereby accelerating gliomagenesis. In addition, GDH1 S384 phosphorylation correlates with the malignancy and prognosis of human glioblastoma. Our finding reveals a unique role of α-KG to directly regulate signal pathway, uncovers a distinct mechanism of metabolite-mediated NF-κB activation, and also establishes the critical role of α-KG-activated NF-κB in brain tumor development.

Enhanced QM/MM sampling for free energy calculation of chemical reactions: A case study of double proton transfer

Free energy calculations for chemical reactions with a steep energy barrier require well defined reaction coordinates (RCs). However, when multiple parallel channels exist along selected RC, the application of conventional enhanced samplings is difficult to generate correct sampling within limited simulation time and thus cannot give correct prediction about the favorable pathways, the relative stability of multiple products or intermediates. Here, we implement the selective integrated tempering sampling (SITS) method with quantum mechanical and molecular mechanical (QM/MM) potential to investigate the chemical reactions in solution. The combined SITS-QM/MM scheme is used to identify possible reaction paths, intermediate and product states, and the free energy profiles for the different reaction paths. Two double proton transfer reactions were studied to validate the implemented method and simulation protocol, from which the independent and correlated proton transfer processes are identified in two representative systems, respectively. This protocol can be generalized to various kinds of chemical reactions for both academic studies and industry applications, such as in exploration and optimization of potential reactions in DNA encoded compound library and halogen or deuterium substitution of the hit discovery and lead optimization stages of drug design via providing a better understanding of the reaction mechanism along the designed chemical reaction pathways.

Structural insight into precursor tRNA processing by yeast ribonuclease P

Ribonuclease P (RNase P) is a universal ribozyme responsible for processing the 5'-leader of pre-transfer RNA (pre-tRNA). Here, we report the 3.5-angstrom cryo-electron microscopy structures of Saccharomyces cerevisiae RNase P alone and in complex with pre-tRNA^Phe The protein components form a hook-shaped architecture that wraps around the RNA and stabilizes RNase P into a "measuring device" with two fixed anchors that recognize the L-shaped pre-tRNA. A universally conserved uridine nucleobase and phosphate backbone in the catalytic center together with the scissile phosphate and the O3' leaving group of pre-tRNA jointly coordinate two catalytic magnesium ions. Binding of pre-tRNA induces a conformational change in the catalytic center that is required for catalysis. Moreover, simulation analysis suggests a two-metal-ion S_N2 reaction pathway of pre-tRNA cleavage. These results not only reveal the architecture of yeast RNase P but also provide a molecular basis of how the 5'-leader of pre-tRNA is processed by eukaryotic RNase P.