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.

Entropy drives the ligand recognition in G-protein-coupled receptor subtypes.

Achieving ligand subtype selectivity within highly homologous subtypes of G-protein-coupled receptor (GPCR) is critical yet challenging for GPCR drug discovery, primarily due to the unclear mechanism underlying ligand subtype selectivity, which hampers the rational design of subtype-selective ligands. Herein, we disclose an unusual molecular mechanism of entropy-driven ligand recognition in cannabinoid (CB) receptor subtypes, revealed through atomic-level molecular dynamics simulations, cryoelectron microscopy structure, and mutagenesis experiments. This mechanism is attributed to the distinct conformational dynamics of the receptor's orthosteric pocket, leading to variations in ligand binding entropy and consequently, differential binding affinities, which culminate in specific ligand recognition. We experimentally validated this mechanism and leveraged it to design ligands with enhanced or ablated subtype selectivity. One such ligand demonstrated favorable pharmacokinetic properties and significant efficacy in rodent inflammatory analgesic models. More importantly, it is precisely due to the high subtype selectivity obtained based on this mechanism that this ligand does not show addictive properties in animal models. Our findings elucidate the unconventional role of entropy in CB receptor subtype selectivity and suggest a strategy for rational design of ligands to achieve entropy-driven subtype selectivity for many pharmaceutically important GPCRs.

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.

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.

The accomplices: Heparan sulfates and N-glycans foster SARS-CoV-2 spike:ACE2 receptor binding and virus priming.

Although it is well established that the SARS-CoV-2 spike glycoprotein binds to the host cell ACE2 receptor to initiate infection, far less is known about the tissue tropism and host cell susceptibility to the virus. Differential expression across different cell types of heparan sulfate (HS) proteoglycans, with variably sulfated glycosaminoglycans (GAGs), and their synergistic interactions with host and viral N-glycans may contribute to tissue tropism and host cell susceptibility. Nevertheless, their contribution remains unclear since HS and N-glycans evade experimental characterization. We, therefore, carried out microsecond-long all-atom molecular dynamics simulations, followed by random acceleration molecular dynamics simulations, of the fully glycosylated spike:ACE2 complex with and without highly sulfated GAG chains bound. By considering the model GAGs as surrogates for the highly sulfated HS expressed in lung cells, we identified key cell entry mechanisms of spike SARS-CoV-2. We find that HS promotes structural and energetic stabilization of the active conformation of the spike receptor-binding domain (RBD) and reorientation of ACE2 toward the N-terminal domain in the same spike subunit as the RBD. Spike and ACE2 N-glycans exert synergistic effects, promoting better packing, strengthening the protein:protein interaction, and prolonging the residence time of the complex. ACE2 and HS binding trigger rearrangement of the S2' functional protease cleavage site through allosteric interdomain communication. These results thus show that HS has a multifaceted role in facilitating SARS-CoV-2 infection, and they provide a mechanistic basis for the development of GAG derivatives with anti-SARS-CoV-2 potential.

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.

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.

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.

Arg92Leu-cTnT Alters the cTnC-cTnI Interface Disrupting PKA-Mediated Relaxation.

BACKGROUND: Impaired left ventricular relaxation, high filling pressures, and dysregulation of Ca2+ homeostasis are common findings contributing to diastolic dysfunction in hypertrophic cardiomyopathy (HCM). Studies have shown that impaired relaxation is an early observation in the sarcomere-gene-positive preclinical HCM cohort, which suggests the potential involvement of myofilament regulators in relaxation. A molecular-level understanding of mechanism(s) at the level of the myofilament is lacking. We hypothesized that mutation-specific, allosterically mediated, changes to the cTnC (cardiac troponin C)-cTnI (cardiac troponin I) interface can account for the development of early-onset diastolic dysfunction via decreased PKA accessibility to cTnI. METHODS: HCM mutations R92L-cTnT (cardiac troponin T; Arg92Leu) and Δ160E-cTnT (Glu160 deletion) were studied in vivo, in vitro, and in silico via 2-dimensional echocardiography, Western blotting, ex vivo hemodynamics, stopped-flow kinetics, time-resolved fluorescence resonance energy transfer, and molecular dynamics simulations. RESULTS: The HCM-causative mutations R92L-cTnT and Δ160E-cTnT result in different time-of-onset diastolic dysfunction. R92L-cTnT demonstrated early-onset diastolic dysfunction accompanied by a localized decrease in phosphorylation of cTnI. Constitutive phosphorylation of cTnI (cTnI-D23D24) was sufficient to recover diastolic function to non-Tg levels only for R92L-cTnT. Mutation-specific changes in Ca2+ dissociation rates associated with R92L-cTnT reconstituted with cTnI-D23D24 led us to investigate potential involvement of structural changes in the cTnC-cTnI interface as an explanation for these observations. We probed the interface via time-resolved fluorescence resonance energy transfer revealing a repositioning of the N-terminus of cTnI, closer to cTnC, and concomitant decreases in distance distributions at sites flanking the PKA consensus sequence. Implementing time-resolved fluorescence resonance energy transfer distances as constraints into our atomistic model identified additional electrostatic interactions at the consensus sequence. CONCLUSIONS: These data show that the early diastolic dysfunction observed in a subset of HCM is attributable to allosterically mediated structural changes at the cTnC-cTnI interface that impair accessibility of PKA, thereby blunting ß-adrenergic responsiveness and identifying a potential molecular target for therapeutic intervention.

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.

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.

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.

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.

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.

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.