How GPCR Phosphorylation Patterns Orchestrate Arrestin-Mediated Signaling.

Binding of arrestin to phosphorylated G-protein-coupled receptors (GPCRs) controls many aspects of cell signaling. The number and arrangement of phosphates may vary substantially for a given GPCR, and different phosphorylation patterns trigger different arrestin-mediated effects. Here, we determine how GPCR phosphorylation influences arrestin behavior by using atomic-level simulations and site-directed spectroscopy to reveal the effects of phosphorylation patterns on arrestin binding and conformation. We find that patterns favoring binding differ from those favoring activation-associated conformational change. Both binding and conformation depend more on arrangement of phosphates than on their total number, with phosphorylation at different positions sometimes exerting opposite effects. Phosphorylation patterns selectively favor a wide variety of arrestin conformations, differently affecting arrestin sites implicated in scaffolding distinct signaling proteins. We also reveal molecular mechanisms of these phenomena. Our work reveals the structural basis for the long-standing "barcode" hypothesis and has important implications for design of functionally selective GPCR-targeted drugs.

Mechanobiology of T Cell Activation: To Catch a Bond.

T cell activation is a critical event in the adaptive immune response, indispensable for cell-mediated and humoral immunity as well as for immune regulation. Recent years have witnessed an emerging trend emphasizing the essential role that physical force and mechanical properties play at the T cell interface. In this review, we integrate current knowledge of T cell antigen recognition and the different models of T cell activation from the perspective of mechanobiology, focusing on the interaction between the T cell receptor (TCR) and the peptide-major histocompatibility complex (pMHC) antigen. We address the shortcomings of TCR affinity alone in explaining T cell functional outcomes and the rising status of force-regulated TCR bond lifetimes, most notably the TCR catch bond. Ultimately, T cell activation and the ensuing physiological responses result from mechanical interaction between TCRs and the pMHC.

Atoms to Phenotypes: Molecular Design Principles of Cellular Energy Metabolism.

We report a 100-million atom-scale model of an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium, that reveals the cascade of energy conversion steps culminating in the generation of ATP from sunlight. Molecular dynamics simulations of this vesicle elucidate how the integral membrane complexes influence local curvature to tune photoexcitation of pigments. Brownian dynamics of small molecules within the chromatophore probe the mechanisms of directional charge transport under various pH and salinity conditions. Reproducing phenotypic properties from atomistic details, a kinetic model evinces that low-light adaptations of the bacterium emerge as a spontaneous outcome of optimizing the balance between the chromatophore's structural integrity and robust energy conversion. Parallels are drawn with the more universal mitochondrial bioenergetic machinery, from whence molecular-scale insights into the mechanism of cellular aging are inferred. Together, our integrative method and spectroscopic experiments pave the way to first-principles modeling of whole living cells.

Structural Basis of Membrane Protein Chaperoning through the Mitochondrial Intermembrane Space.

The exchange of metabolites between the mitochondrial matrix and the cytosol depends on ß-barrel channels in the outer membrane and α-helical carrier proteins in the inner membrane. The essential translocase of the inner membrane (TIM) chaperones escort these proteins through the intermembrane space, but the structural and mechanistic details remain elusive. We have used an integrated structural biology approach to reveal the functional principle of TIM chaperones. Multiple clamp-like binding sites hold the mitochondrial membrane proteins in a translocation-competent elongated form, thus mimicking characteristics of co-translational membrane insertion. The bound preprotein undergoes conformational dynamics within the chaperone binding clefts, pointing to a multitude of dynamic local binding events. Mutations in these binding sites cause cell death or growth defects associated with impairment of carrier and ß-barrel protein biogenesis. Our work reveals how a single mitochondrial "transfer-chaperone" system is able to guide α-helical and ß-barrel membrane proteins in a "nascent chain-like" conformation through a ribosome-free compartment.

Whole-Cell Models and Simulations in Molecular Detail.

Comprehensive data about the composition and structure of cellular components have enabled the construction of quantitative whole-cell models. While kinetic network-type models have been established, it is also becoming possible to build physical, molecular-level models of cellular environments. This review outlines challenges in constructing and simulating such models and discusses near- and long-term opportunities for developing physical whole-cell models that can connect molecular structure with biological function.

A universal allosteric mechanism for G protein activation.

G proteins play a central role in signal transduction and pharmacology. Signaling is initiated by cell-surface receptors, which promote guanosine triphosphate (GTP) binding and dissociation of Gα from the Gßγ subunits. Structural studies have revealed the molecular basis of subunit association with receptors, RGS proteins, and downstream effectors. In contrast, the mechanism of subunit dissociation is poorly understood. We use cell signaling assays, molecular dynamics (MD) simulations, and biochemistry and structural analyses to identify a conserved network of amino acids that dictates subunit release. In the presence of the terminal phosphate of GTP, a glycine forms a polar network with an arginine and glutamate, putting torsional strain on the subunit binding interface. This "G-R-E motif" secures GTP and, through an allosteric link, discharges the Gßγ dimer. Replacement of network residues prevents subunit dissociation regardless of agonist or GTP binding. These findings reveal the molecular basis of the final committed step of G protein activation.

Sum rule comparison of narrowband and broadband sum frequency generation spectra and comparison with theory.

Earlier sum frequency generation (SFG) experiments involve one infrared and one visible laser, and a measurement of the intensity of the response, yielding data on the surface sensitive properties of the sample. Recently, both the real and imaginary components of the susceptibility were measured in two different sets of experiments. In one set, a broadband infrared laser was used, permitting observations at very short times, while in another set the infrared laser was narrowband, permitting higher spectral resolution. The differences in the spectrum obtained by the two will be most evident in studying narrow absorption bands, e.g., the band due to dangling OH bonds at a water interface. The direct comparisons in the integrated amplitude (sum rule) of the imaginary part of the dangling OH bond region differ by a factor of 3. Due to variations in experimental setup and data processing, corrections were made for the quartz reference, Fresnel factors, and the incident visible laser wavelength. After the corrections, the agreement differs now by the factors of 1.1 within broadband and narrowband groups and the two groups now differ by a factor of 1.5. The 1.5 factor may arise from the extra heating of the more powerful broadband laser system on the water surface. The convolution from the narrowband SFG spectrum to the broadband SFG spectrum is also investigated and it does not affect the sum rule. Theory and narrowband experiments are compared using the sum rule and agree to a factor of 1.3 with no adjustable parameters.

Homologous mutations in human ß, embryonic, and perinatal muscle myosins have divergent effects on molecular power generation.

Mutations at a highly conserved homologous residue in three closely related muscle myosins cause three distinct diseases involving muscle defects: R671C in ß-cardiac myosin causes hypertrophic cardiomyopathy, R672C and R672H in embryonic skeletal myosin cause Freeman-Sheldon syndrome, and R674Q in perinatal skeletal myosin causes trismus-pseudocamptodactyly syndrome. It is not known whether their effects at the molecular level are similar to one another or correlate with disease phenotype and severity. To this end, we investigated the effects of the homologous mutations on key factors of molecular power production using recombinantly expressed human ß, embryonic, and perinatal myosin subfragment-1. We found large effects in the developmental myosins but minimal effects in ß myosin, and magnitude of changes correlated partially with clinical severity. The mutations in the developmental myosins dramatically decreased the step size and load-sensitive actin-detachment rate of single molecules measured by optical tweezers, in addition to decreasing overall enzymatic (ATPase) cycle rate. In contrast, the only measured effect of R671C in ß myosin was a larger step size. Our measurements of step size and bound times predicted velocities consistent with those measured in an in vitro motility assay. Finally, molecular dynamics simulations predicted that the arginine to cysteine mutation in embryonic, but not ß, myosin may reduce pre-powerstroke lever arm priming and ADP pocket opening, providing a possible structural mechanism consistent with the experimental observations. This paper presents direct comparisons of homologous mutations in several different myosin isoforms, whose divergent functional effects are a testament to myosin's highly allosteric nature.

C-type inactivation and proton modulation mechanisms of the TASK3 channel.

The TWIK-related acid-sensitive K+ channel 3 (TASK3) belongs to the two-pore domain (K2P) potassium channel family, which regulates cell excitability by mediating a constitutive "leak" potassium efflux in the nervous system. Extracellular acidification inhibits TASK3 channel, but the molecular mechanism by which channel inactivation is coupled to pH decrease remains unclear. Here, we report the cryo-electron microscopy structures of human TASK3 at neutral and acidic pH. Structural comparison revealed selectivity filter (SF) rearrangements upon acidification, characteristic of C-type inactivation, but with a unique structural basis. The extracellular mouth of the SF was prominently dilated and simultaneously blocked by a hydrophobic gate. His98 protonation shifted the conformational equilibrium between the conductive and C-type inactivated SF toward the latter by engaging a cation-π interaction with Trp78, consistent with molecular dynamics simulations and electrophysiological experiments. Our work illustrated how TASK3 is gated in response to extracellular pH change and implies how physiological stimuli might directly modulate the C-type gating of K2P channels.

Discovery of lirafugratinib (RLY-4008), a highly selective irreversible small-molecule inhibitor of FGFR2.

Fibroblast growth factor receptor (FGFR) kinase inhibitors have been shown to be effective in the treatment of intrahepatic cholangiocarcinoma and other advanced solid tumors harboring FGFR2 alterations, but the toxicity of these drugs frequently leads to dose reduction or interruption of treatment such that maximum efficacy cannot be achieved. The most common adverse effects are hyperphosphatemia caused by FGFR1 inhibition and diarrhea due to FGFR4 inhibition, as current therapies are not selective among the FGFRs. Designing selective inhibitors has proved difficult with conventional approaches because the orthosteric sites of FGFR family members are observed to be highly similar in X-ray structures. In this study, aided by analysis of protein dynamics, we designed a selective, covalent FGFR2 inhibitor. In a key initial step, analysis of long-timescale molecular dynamics simulations of the FGFR1 and FGFR2 kinase domains allowed us to identify differential motion in their P-loops, which are located adjacent to the orthosteric site. Using this insight, we were able to design orthosteric binders that selectively and covalently engage the P-loop of FGFR2. Our drug discovery efforts culminated in the development of lirafugratinib (RLY-4008), a covalent inhibitor of FGFR2 that shows substantial selectivity over FGFR1 (~250-fold) and FGFR4 (~5,000-fold) in vitro, causes tumor regression in multiple FGFR2-altered human xenograft models, and was recently demonstrated to be efficacious in the clinic at doses that do not induce clinically significant hyperphosphatemia or diarrhea.

Approaching the catalytic mechanism of protein lysine methyltransferases by biochemical and simulation techniques.

Protein lysine methyltransferases (PKMTs) transfer up to three methyl groups to the side chains of lysine residues in proteins and fulfill important regulatory functions by controlling protein stability, localization and protein/protein interactions. The methylation reactions are highly regulated, and aberrant methylation of proteins is associated with several types of diseases including neurologic disorders, cardiovascular diseases, and various types of cancer. This review describes novel insights into the catalytic machinery of various PKMTs achieved by the combined application of biochemical experiments and simulation approaches during the last years, focusing on clinically relevant and well-studied enzymes of this group like DOT1L, SMYD1-3, SET7/9, G9a/GLP, SETD2, SUV420H2, NSD1/2, different MLLs and EZH2. Biochemical experiments have unraveled many mechanistic features of PKMTs concerning their substrate and product specificity, processivity and the effects of somatic mutations observed in PKMTs in cancer cells. Structural data additionally provided information about the substrate recognition, enzyme-substrate complex formation, and allowed for simulations of the substrate peptide interaction and mechanism of PKMTs with atomistic resolution by molecular dynamics and hybrid quantum mechanics/molecular mechanics methods. These simulation technologies uncovered important mechanistic details of the PKMT reaction mechanism including the processes responsible for the deprotonation of the target lysine residue, essential conformational changes of the PKMT upon substrate binding, but also rationalized regulatory principles like PKMT autoinhibition. Further developments are discussed that could bring us closer to a mechanistic understanding of catalysis of this important class of enzymes in the near future. The results described here illustrate the power of the investigation of enzyme mechanisms by the combined application of biochemical experiments and simulation technologies.

Mechanistic View of hnRNPA2 Low-Complexity Domain Structure, Interactions, and Phase Separation Altered by Mutation and Arginine Methylation.

hnRNPA2, a component of RNA-processing membraneless organelles, forms inclusions when mutated in a syndrome characterized by the degeneration of neurons (bearing features of amyotrophic lateral sclerosis [ALS] and frontotemporal dementia), muscle, and bone. Here we provide a unified structural view of hnRNPA2 self-assembly, aggregation, and interaction and the distinct effects of small chemical changes-disease mutations and arginine methylation-on these assemblies. The hnRNPA2 low-complexity (LC) domain is compact and intrinsically disordered as a monomer, retaining predominant disorder in a liquid-liquid phase-separated form. Disease mutations D290V and P298L induce aggregation by enhancing and extending, respectively, the aggregation-prone region. Co-aggregating in disease inclusions, hnRNPA2 LC directly interacts with and induces phase separation of TDP-43. Conversely, arginine methylation reduces hnRNPA2 phase separation, disrupting arginine-mediated contacts. These results highlight the mechanistic role of specific LC domain interactions and modifications conserved across many hnRNP family members but altered by aggregation-causing pathological mutations.

Universality in RNA and DNA deformations induced by salt, temperature change, stretching force, and protein binding.

Nucleic acid deformations play important roles in many biological processes. The physical understanding of nucleic acid deformation by environmental stimuli is limited due to the challenge in the precise measurement of RNA and DNA deformations and the complexity of interactions in RNA and DNA. Magnetic tweezers experiments provide an excellent opportunity to precisely measure DNA and RNA twist changes induced by environmental stimuli. In this work, we applied magnetic tweezers to measure double-stranded RNA twist changes induced by salt and temperature changes. We observed RNA unwinds when lowering salt concentration, or increasing temperature. Our molecular dynamics simulations revealed the mechanism: lowering salt concentration or increasing temperature enlarges RNA major groove width, which causes twist decrease through twist-groove coupling. Combining these results with previous results, we found some universality in RNA and DNA deformations induced by three different stimuli: salt change, temperature, and stretching force. For RNA, these stimuli first modify the major groove width, which is transduced into twist change through twist-groove coupling. For DNA, these stimuli first modify diameter, which is transduced into twist change through twist-diameter coupling. Twist-groove coupling and twist-diameter coupling appear to be utilized by protein binding to reduce DNA and RNA deformation energy cost upon protein binding.

Coenzyme Q10 trapping in mitochondrial complex I underlies Leber's hereditary optic neuropathy.

How does a single amino acid mutation occurring in the blinding disease, Leber's hereditary optic neuropathy (LHON), impair electron shuttling in mitochondria? We investigated changes induced by the m.3460 G>A mutation in mitochondrial protein ND1 using the tools of Molecular Dynamics and Free Energy Perturbation simulations, with the goal of determining the mechanism by which this mutation affects mitochondrial function. A recent analysis suggested that the mutation's replacement of alanine A52 with a threonine perturbs the stability of a region where binding of the electron shuttling protein, Coenzyme Q10, occurs. We found two functionally opposing changes involving the role of Coenzyme Q10. The first showed that quantum electron transfer from the terminal Fe/S complex, N2, to the Coenzyme Q10 headgroup, docked in its binding pocket, is enhanced. However, this positive adjustment is overshadowed by our finding that the mobility of Coenzyme Q10 in its oxidized and reduced states, entering and exiting its binding pocket, is disrupted by the mutation in a manner that leads to conditions promoting the generation of reactive oxygen species. An increase in reactive oxygen species caused by the LHON mutation has been proposed to be responsible for this optic neuropathy.

Generation of antibodies to an extracellular region of the transporters Glut1/Glut4 by immunization with a designed antigen.

Monoclonal antibodies are one of the fastest growing class of drugs. Nevertheless, relatively few biologics target multispanning membrane proteins because of technical challenges. To target relatively small extracellular regions of multiple membrane-spanning proteins, synthetic peptides, which are composed of amino acids corresponding to an extracellular region of a membrane protein, are often utilized in antibody discovery. However, antibodies to these peptides often do not recognize parental membrane proteins. In this study, we designed fusion proteins in which an extracellular helix of the membrane protein glucose transporter 1 (Glut1) was grafted onto the scaffold protein Adhiron. In the initial design, the grafted fragment did not form a helical conformation. Molecular dynamics simulations of full-length Glut1 suggested the importance of intramolecular interactions formed by surrounding residues in the formation of the helical conformation. A fusion protein designed to maintain such intramolecular interactions did form the desired helical conformation in the grafted region. We then immunized an alpaca with the designed fusion protein and obtained VHH (variable region of heavy-chain antibodies) using the phage display method. The binding of these VHH antibodies to the recombinant Glut1 protein was evaluated by surface plasmon resonance, and their binding to Glut1 on the cell membrane was further validated by flow cytometry. Furthermore, we also succeeded in the generation of a VHH against another integral membrane protein, glucose transporter 4 (Glut4) with the same strategy. These illustrates that our combined biochemical and computational approach can be applied to designing other novel fusion proteins for generating site-specific antibodies.

Shear transformations in metallic glasses without excessive and predefinable defects.

Plastic flow in metallic glasses (MGs) is known to be mediated by shear transformations (STs), which have been hypothesized to preferentially initiate from identifiable local "defect" regions with loose atomic packing. Here we show that the above idea is incorrect, i.e., STs do not arise from signature structural defects that can be recognized a priori. This conclusion is reached via a realistic MG model obtained by combining molecular dynamics (MD) and Monte Carlo simulations, achieving liquid solidification at an effective cooling rate as slow as 500 K/s to approach that typical in experiments for producing bulk MGs. At shear stresses before global yielding, only about 2% of the total atoms participate in STs, each event involving typically ~10 atoms. These observations rectify the excessive content of "liquid-like regions" retained from unrealistically fast quench in MD-produced glass models. Our findings also shed light on the indeterministic aspect of the ST sites/zones, which emerge with varying spatial locations and distribution depending on specific mechanical loading conditions.

The dielectric function profile across the water interface through surface-specific vibrational spectroscopy and simulations.

The dielectric properties of interfacial water on subnanometer length scales govern chemical reactions, carrier transfer, and ion transport at interfaces. Yet, the nature of the interfacial dielectric function has remained under debate as it is challenging to access the interfacial dielectric with subnanometer resolution. Here we use the vibrational response of interfacial water molecules probed using surface-specific sum-frequency generation (SFG) spectra to obtain exquisite depth resolution. Different responses originate from water molecules at different depths and report back on the local interfacial dielectric environment via their spectral amplitudes. From experimental and simulated SFG spectra at the air/water interface, we find that the interfacial dielectric constant changes drastically across an ∼1 Šthin interfacial water region. The strong gradient of the interfacial dielectric constant leads, at charged planar interfaces, to the formation of an electric triple layer that goes beyond the standard double-layer model.

In situ optical spectroscopy of crystallization: One crystal nucleation at a time.

While crystallization is a ubiquitous and an important process, the microscopic picture of crystal nucleation is yet to be established. Recent studies suggest that the nucleation process can be more complex than the view offered by the classical nucleation theory. Here, we implement single crystal nucleation spectroscopy (SCNS) by combining Raman microspectroscopy and optical trapping induced crystallization to spectroscopically investigate one crystal nucleation at a time. Raman spectral evolution during a single glycine crystal nucleation from water, measured by SCNS and analyzed by a nonsupervised spectral decomposition technique, uncovered the Raman spectrum of prenucleation aggregates and their critical role as an intermediate species in the dynamics. The agreement between the spectral feature of prenucleation aggregates and our simulation suggests that their structural order emerges through the dynamic formation of linear hydrogen-bonded networks. The present work provides a strong impetus for accelerating the investigation of crystal nucleation by optical spectroscopy.

Single-Molecule Discrimination of Saccharides Using Carbon Nitride Nanopores.

Structural complexity brings a huge challenge to the analysis of sugar chains. As a single-molecule sensor, nanopores have the potential to provide fingerprint information on saccharides. Traditionally, direct single-molecule saccharide detection with nanopores is hampered by their small size and weak affinity. Here, a carbon nitride nanopore device is developed to discern two types of trisaccharide molecules (LeApN and SLeCpN) with minor structural differences. The resolution of LeApN and SLeCpN in the mixture reaches 0.98, which has never been achieved in solid-state nanopores so far. Monosaccharide (GlcNAcpN) and disaccharide (LacNAcpN) can also be discriminated using this system, indicating that the versatile carbon nitride nanopores possess a monosaccharide-level resolution. This study demonstrates that the carbon nitride nanopores have the potential for conducting structure analysis on single-molecule saccharides.

Controllable Pseudoelasticity in Metallic Nanocrystals by Grain Boundary Engineering.

Materials with pseudoelasticity can recover from large strains exceeding their elastic limits during unloading, making them promising damage-tolerant building blocks for advanced nanodevices. Nevertheless, a practical approach to realize controllable pseudoelastic behavior at nanoscale remains challenging. Here, we proposed a grain boundary (GB) engineering approach to endow metallic nanocrystals with a controllable pseudoelasticity. Both in situ nanomechanical testing and atomistic simulations demonstrate that such controllable pseudoelasticity is governed by the extension and contraction of an inherent stacking fault array at the GB. By precisely tuning GB misorientation and inclination, our simulation results reveal that metallic nanocrystals can exhibit tailored pseudoelastic performance across a broad spectrum of GBs in different face-centered cubic metals. These findings enrich our understanding of the intrinsic pseudoelasticity of GBs and provide a GB engineering approach toward metallic materials with reversible deformability.