Fast charging of energy-dense lithium-ion batteries.

Lithium-ion batteries with nickel-rich layered oxide cathodes and graphite anodes have reached specific energies of 250-300 Wh kg-1 (refs. 1,2), and it is now possible to build a 90 kWh electric vehicle (EV) pack with a 300-mile cruise range. Unfortunately, using such massive batteries to alleviate range anxiety is ineffective for mainstream EV adoption owing to the limited raw resource supply and prohibitively high cost. Ten-minute fast charging enables downsizing of EV batteries for both affordability and sustainability, without causing range anxiety. However, fast charging of energy-dense batteries (more than 250 Wh kg-1 or higher than 4 mAh cm-2) remains a great challenge3,4. Here we combine a material-agnostic approach based on asymmetric temperature modulation with a thermally stable dual-salt electrolyte to achieve charging of a 265 Wh kg-1 battery to 75% (or 70%) state of charge in 12 (or 11) minutes for more than 900 (or 2,000) cycles. This is equivalent to a half million mile range in which every charge is a fast charge. Further, we build a digital twin of such a battery pack to assess its cooling and safety and demonstrate that thermally modulated 4C charging only requires air convection. This offers a compact and intrinsically safe route to cell-to-pack development. The rapid thermal modulation method to yield highly active electrochemical interfaces only during fast charging has important potential to realize both stability and fast charging of next-generation materials, including anodes like silicon and lithium metal.

Computationally Discovered Potentiating Role of Glycans on NMDA Receptors.

N-methyl-D-aspartate receptors (NMDARs) are glycoproteins in the brain central to learning and memory. The effects of glycosylation on the structure and dynamics of NMDARs are largely unknown. In this work, we use extensive molecular dynamics simulations of GluN1 and GluN2B ligand binding domains (LBDs) of NMDARs to investigate these effects. Our simulations predict that intra-domain interactions involving the glycan attached to residue GluN1-N440 stabilize closed-clamshell conformations of the GluN1 LBD. The glycan on GluN2B-N688 shows a similar, though weaker, effect. Based on these results, and assuming the transferability of the results of LBD simulations to the full receptor, we predict that glycans at GluN1-N440 might play a potentiator role in NMDARs. To validate this prediction, we perform electrophysiological analysis of full-length NMDARs with a glycosylation-preventing GluN1-N440Q mutation, and demonstrate an increase in the glycine EC50 value. Overall, our results suggest an intramolecular potentiating role of glycans on NMDA receptors.

The pathway of ligand entry from the membrane bilayer to a lipid G protein-coupled receptor.

The binding process through the membrane bilayer of lipid-like ligands to a protein target is an important but poorly explored recognition process at the atomic level. In this work we succeeded in resolving the binding of the lipid inhibitor ML056 to the sphingosine-1-phosphate receptor 1 (S1P1R) using unbiased molecular dynamics simulations with an aggregate sampling of over 800 µs. The binding pathway is a multi-stage process consisting of the ligand diffusing in the bilayer leaflet to contact a "membrane vestibule" at the top of TM 7, subsequently moving from this lipid-facing vestibule to the orthosteric binding cavity through a channel formed by TMs 1 and 7 and the N-terminal of the receptor. Unfolding of the N-terminal alpha-helix increases the volume of the channel upon ligand entry, helping to reach the crystallographic pose that also corresponds to the predicted favorable pose. The relaxation timescales of the binding process show that the binding of the ligand to the "membrane vestibule" is the rate-limiting step in the multi microseconds timescale. We comment on the significance and parallels of the binding process in the context of other binding studies.

Spontaneous inward opening of the dopamine transporter is triggered by PIP2-regulated dynamics of the N-terminus.

We present the dynamic mechanism of concerted motions in a full-length molecular model of the human dopamine transporter (hDAT), a member of the neurotransmitter/sodium symporter (NSS) family, involved in state-to-state transitions underlying function. The findings result from an analysis of unbiased atomistic molecular dynamics simulation trajectories (totaling >14 µs) of the hDAT molecule immersed in lipid membrane environments with or without phosphatidylinositol 4,5-biphosphate (PIP2) lipids. The N-terminal region of hDAT (N-term) is shown to have an essential mechanistic role in correlated rearrangements of specific structural motifs relevant to state-to-state transitions in the hDAT. The mechanism involves PIP2-mediated electrostatic interactions between the N-term and the intracellular loops of the transporter molecule. Quantitative analyses of collective motions in the trajectories reveal that these interactions correlate with the inward-opening dynamics of hDAT and are allosterically coupled to the known functional sites of the transporter. The observed large-scale motions are enabled by specific reconfiguration of the network of ionic interactions at the intracellular end of the protein. The isomerization to the inward-facing state in hDAT is accompanied by concomitant movements in the extracellular vestibule and results in the release of an Na(+) ion from the Na2 site and destabilization of the substrate dopamine in the primary substrate binding S1 site. The dynamic mechanism emerging from the findings highlights the involvement of the PIP2-regulated interactions between the N-term and the intracellular loop 4 in the functionally relevant conformational transitions that are also similar to those found to underlie state-to-state transitions in the leucine transporter (LeuT), a prototypical bacterial homologue of the NSS.

Progress in studying intrinsically disordered proteins with atomistic simulations.

Intrinsically disordered proteins are increasingly the focus of biological research since their significance was acknowledged over a decade ago. Due to their importance in biomolecular interactions, they are found to play key roles in many diseases such as cancers and amyloidoses. However, because they lack stable structure they pose a challenge for many experimental methods that are traditionally used to study proteins. Atomistic molecular dynamics simulations can help get around many of the problems faced by such methods provided appropriate timescales are sampled and underlying empirical force fields are applicable. This review presents recent works that highlight the power and potential of atomistic simulations to transform the investigatory pipeline by providing critical insights into the behavior and interactions of intrinsically disordered proteins.

High throughput molecular dynamics for drug discovery.

Molecular dynamics simulations hold the promise to be an important tool for biological research and drug discovery. Historically, however, there were several obstacles for it to become a practical research tool. Limitations in computer hardware had previously made it difficult to simulate for long enough to see interesting biological processes. Recent improvements in hardware and algorithms have largely removed this issue, leaving data analysis as the main obstacle. Advances in Markov state modeling appear to be on the way to remove this obstacle. We outline these advances here and discuss numerous recent studies that demonstrate that molecular dynamics simulations will start to be an important tool for pharmaceutical research.

Kinetic modulation of a disordered protein domain by phosphorylation.

Phosphorylation is a major post-translational mechanism of regulation that frequently targets disordered protein domains, but it remains unclear how phosphorylation modulates disordered states of proteins. Here we determine the kinetics and energetics of a disordered protein domain the kinase-inducible domain (KID) of the transcription factor CREB and that of its phosphorylated form pKID, using high-throughput molecular dynamic simulations. We identify the presence of a metastable, partially ordered state with a 60-fold slowdown in conformational kinetics that arises due to phosphorylation, kinetically stabilizing residues known to participate in an early binding intermediate. We show that this effect is only partially reconstituted by mutation to glutamate, indicating that the phosphate is uniquely required for the long-lived state to arise. This mechanism of kinetic modulation could be important for regulation beyond conformational equilibrium shifts.

Membrane lipids are key modulators of the endocannabinoid-hydrolase FAAH.

Lipid composition is expected to play an important role in modulating membrane enzyme activity, in particular if the substrates are themselves lipid molecules. A paradigmatic case is FAAH (fatty acid amide hydrolase), an enzyme critical in terminating endocannabinoid signalling and an important therapeutic target. In the present study, using a combined experimental and computational approach, we show that membrane lipids modulate the structure, subcellular localization and activity of FAAH. We report that the FAAH dimer is stabilized by the lipid bilayer and shows a higher membrane-binding affinity and enzymatic activity within membranes containing both cholesterol and the natural FAAH substrate AEA (anandamide). Additionally, co-localization of cholesterol, AEA and FAAH in mouse neuroblastoma cells suggests a mechanism through which cholesterol increases the substrate accessibility of FAAH.