Angiotensin Analogs with Divergent Bias Stabilize Distinct Receptor Conformations.

"Biased" G protein-coupled receptor (GPCR) agonists preferentially activate pathways mediated by G proteins or ß-arrestins. Here, we use double electron-electron resonance spectroscopy to probe the changes that ligands induce in the conformational distribution of the angiotensin II type I receptor. Monitoring distances between 10 pairs of nitroxide labels distributed across the intracellular regions enabled mapping of four underlying sets of conformations. Ligands from different functional classes have distinct, characteristic effects on the conformational heterogeneity of the receptor. Compared to angiotensin II, the endogenous agonist, agonists with enhanced Gq coupling more strongly stabilize an "open" conformation with an accessible transducer-binding site. ß-arrestin-biased agonists deficient in Gq coupling do not stabilize this open conformation but instead favor two more occluded conformations. These data suggest a structural mechanism for biased ligand action at the angiotensin receptor that can be exploited to rationally design GPCR-targeting drugs with greater specificity of action.

Molecular dynamics-based identification of binding pathways and two distinct high-affinity sites for succinate in succinate receptor 1/GPR91.

SUCNR1 is an auto- and paracrine sensor of the metabolic stress signal succinate. Using unsupervised molecular dynamics (MD) simulations (170.400 ns) and mutagenesis across human, mouse, and rat SUCNR1, we characterize how a five-arginine motif around the extracellular pole of TM-VI determines the initial capture of succinate in the extracellular vestibule (ECV) to either stay or move down to the orthosteric site. Metadynamics demonstrate low-energy succinate binding in both sites, with an energy barrier corresponding to an intermediate stage during which succinate, with an associated water cluster, unlocks the hydrogen-bond-stabilized conformationally constrained extracellular loop (ECL)-2b. Importantly, simultaneous binding of two succinate molecules through either a "sequential" or "bypassing" mode is a frequent endpoint. The mono-carboxylate NF-56-EJ40 antagonist enters SUCNR1 between TM-I and -II and does not unlock ECL-2b. It is proposed that occupancy of both high-affinity sites is required for selective activation of SUCNR1 by high local succinate concentrations.

Structure of the human signal peptidase complex reveals the determinants for signal peptide cleavage.

The signal peptidase complex (SPC) is an essential membrane complex in the endoplasmic reticulum (ER), where it removes signal peptides (SPs) from a large variety of secretory pre-proteins with exquisite specificity. Although the determinants of this process have been established empirically, the molecular details of SP recognition and removal remain elusive. Here, we show that the human SPC exists in two functional paralogs with distinct proteolytic subunits. We determined the atomic structures of both paralogs using electron cryo-microscopy and structural proteomics. The active site is formed by a catalytic triad and abuts the ER membrane, where a transmembrane window collectively formed by all subunits locally thins the bilayer. Molecular dynamics simulations indicate that this unique architecture generates specificity for SPs based on the length of their hydrophobic segments.

Structure, mechanism, and inhibition of Hedgehog acyltransferase.

The Sonic Hedgehog (SHH) morphogen pathway is fundamental for embryonic development and stem cell maintenance and is implicated in various cancers. A key step in signaling is transfer of a palmitate group to the SHH N terminus, catalyzed by the multi-pass transmembrane enzyme Hedgehog acyltransferase (HHAT). We present the high-resolution cryo-EM structure of HHAT bound to substrate analog palmityl-coenzyme A and a SHH-mimetic megabody, revealing a heme group bound to HHAT that is essential for HHAT function. A structure of HHAT bound to potent small-molecule inhibitor IMP-1575 revealed conformational changes in the active site that occlude substrate binding. Our multidisciplinary analysis provides a detailed view of the mechanism by which HHAT adapts the membrane environment to transfer an acyl chain across the endoplasmic reticulum membrane. This structure of a membrane-bound O-acyltransferase (MBOAT) superfamily member provides a blueprint for other protein-substrate MBOATs and a template for future drug discovery.

Transport mechanism of DgoT, a bacterial homolog of SLC17 organic anion transporters.

The solute carrier 17 (SLC17) family contains anion transporters that accumulate neurotransmitters in secretory vesicles, remove carboxylated monosaccharides from lysosomes, or extrude organic anions from the kidneys and liver. We combined classical molecular dynamics simulations, Markov state modeling and hybrid first principles quantum mechanical/classical mechanical (QM/MM) simulations with experimental approaches to describe the transport mechanisms of a model bacterial protein, the D-galactonate transporter DgoT, at atomic resolution. We found that protonation of D46 and E133 precedes galactonate binding and that substrate binding induces closure of the extracellular gate, with the conserved R47 coupling substrate binding to transmembrane helix movement. After isomerization to an inward-facing conformation, deprotonation of E133 and subsequent proton transfer from D46 to E133 opens the intracellular gate and permits galactonate dissociation either in its unprotonated form or after proton transfer from E133. After release of the second proton, apo DgoT returns to the outward-facing conformation. Our results provide a framework to understand how various SLC17 transport functions with distinct transport stoichiometries can be attained through subtle variations in proton and substrate binding/unbinding.

Advances in computational methods for ligand binding kinetics.

Binding kinetic parameters can be correlated with drug efficacy, which in recent years led to the development of various computational methods for predicting binding kinetic rates and gaining insight into protein-drug binding paths and mechanisms. In this review, we introduce and compare computational methods recently developed and applied to two systems, trypsin-benzamidine and kinase-inhibitor complexes. Methods involving enhanced sampling in molecular dynamics simulations or machine learning can be used not only to predict kinetic rates, but also to reveal factors modulating the duration of residence times, selectivity, and drug resistance to mutations. Methods which require less computational time to make predictions are highlighted, and suggestions to reduce the error of computed kinetic rates are presented.

The Ebola virus VP40 matrix layer undergoes endosomal disassembly essential for membrane fusion.

Ebola viruses (EBOVs) assemble into filamentous virions, whose shape and stability are determined by the matrix viral protein 40 (VP40). Virus entry into host cells occurs via membrane fusion in late endosomes; however, the mechanism of how the remarkably long virions undergo uncoating, including virion disassembly and nucleocapsid release into the cytosol, remains unknown. Here, we investigate the structural architecture of EBOVs entering host cells and discover that the VP40 matrix disassembles prior to membrane fusion. We reveal that VP40 disassembly is caused by the weakening of VP40-lipid interactions driven by low endosomal pH that equilibrates passively across the viral envelope without a dedicated ion channel. We further show that viral membrane fusion depends on VP40 matrix integrity, and its disassembly reduces the energy barrier for fusion stalk formation. Thus, pH-driven structural remodeling of the VP40 matrix acts as a molecular switch coupling viral matrix uncoating to membrane fusion during EBOV entry.

Mechano-regulation of Peptide-MHC Class I Conformations Determines TCR Antigen Recognition.

TCRs recognize cognate pMHCs to initiate T cell signaling and adaptive immunity. Mechanical force strengthens TCR-pMHC interactions to elicit agonist-specific catch bonds to trigger TCR signaling, but the underlying dynamic structural mechanism is unclear. We combined steered molecular dynamics (SMD) simulation, single-molecule biophysical approaches, and functional assays to collectively demonstrate that mechanical force induces conformational changes in pMHCs to enhance pre-existing contacts and activates new interactions at the TCR-pMHC binding interface to resist bond dissociation under force, resulting in TCR-pMHC catch bonds and T cell activation. Intriguingly, cancer-associated somatic mutations in HLA-A2 that may restrict these conformational changes suppressed TCR-pMHC catch bonds. Structural analysis also indicated that HLA polymorphism might alter the equilibrium of these conformational changes. Our findings not only reveal critical roles of force-induced conformational changes in pMHCs for activating TCR-pMHC catch bonds but also have implications for T cell-based immunotherapy.

Molecular mechanism underlying SNARE-mediated membrane fusion enlightened by all-atom molecular dynamics simulations.

The SNAP receptor (SNARE) proteins syntaxin-1, SNAP-25, and synaptobrevin mediate neurotransmitter release by forming tight SNARE complexes that fuse synaptic vesicles with the plasma membranes in microseconds. Membrane fusion is generally explained by the action of proteins on macroscopic membrane properties such as curvature, elastic modulus, and tension, and a widespread model envisions that the SNARE motifs, juxtamembrane linkers, and C-terminal transmembrane regions of synaptobrevin and syntaxin-1 form continuous helices that act mechanically as semirigid rods, squeezing the membranes together as they assemble ("zipper") from the N to the C termini. However, the mechanism underlying fast SNARE-induced membrane fusion remains unknown. We have used all-atom molecular dynamics simulations to investigate this mechanism. Our results need to be interpreted with caution because of the limited number and length of the simulations, but they suggest a model of membrane fusion that has a natural physicochemical basis, emphasizes local molecular events over general membrane properties, and explains extensive experimental data. In this model, the central event that initiates fast (microsecond scale) membrane fusion occurs when the SNARE helices zipper into the juxtamembrane linkers which, together with the adjacent transmembrane regions, promote encounters of acyl chains from both bilayers at the polar interface. The resulting hydrophobic nucleus rapidly expands into stalk-like structures that gradually progress to form a fusion pore, aided by the SNARE transmembrane regions and without clearly discernible intermediates. The propensity of polyunsaturated lipids to participate in encounters that initiate fusion suggests that these lipids may be important for the high speed of neurotransmitter release.

Asymmetric allostery in estrogen receptor-α homodimers drives responses to the ensemble of estrogens in the hormonal milieu.

The estrogen receptor-α (ER) is thought to function only as a homodimer but responds to a variety of environmental, metazoan, and therapeutic estrogens at subsaturating doses, supporting binding mixtures of ligands as well as dimers that are only partially occupied. Here, we present a series of flexible ER ligands that bind to receptor dimers with individual ligand poses favoring distinct receptor conformations-receptor conformational heterodimers-mimicking the binding of two different ligands. Molecular dynamics simulations showed that the pairs of different ligand poses changed the correlated motion across the dimer interface to generate asymmetric communication between the dimer interface, the ligands, and the surface binding sites for epigenetic regulatory proteins. By examining the binding of the same ligand in crystal structures of ER in the agonist vs. antagonist conformers, we also showed that these allosteric signals are bidirectional. The receptor conformer can drive different ligand binding modes to support agonist vs. antagonist activity profiles, a revision of ligand binding theory that has focused on unidirectional signaling from the ligand to the coregulator binding site. We also observed differences in the allosteric signals between ligand and coregulator binding sites in the monomeric vs. dimeric receptor, and when bound by two different ligands, states that are physiologically relevant. Thus, ER conformational heterodimers integrate two different ligand-regulated activity profiles, representing different modes for ligand-dependent regulation of ER activity.

Identification of the potassium-binding site in serotonin transporter.

Clearance of serotonin (5-hydroxytryptamine, 5-HT) from the synaptic cleft after neuronal signaling is mediated by serotonin transporter (SERT), which couples this process to the movement of a Na+ ion down its chemical gradient. After release of 5-HT and Na+ into the cytoplasm, the transporter faces a rate-limiting challenge of resetting its conformation to be primed again for 5-HT and Na+ binding. Early studies of vesicles containing native SERT revealed that K+ gradients can provide an additional driving force, via K+ antiport. Moreover, under appropriate conditions, a H+ ion can replace K+. Intracellular K+ accelerates the resetting step. Structural studies of SERT have identified two binding sites for Na+ ions, but the K+ site remains enigmatic. Here, we show that K+ antiport can drive substrate accumulation into vesicles containing SERT extracted from a heterologous expression system, allowing us to study the residues responsible for K+ binding. To identify candidate binding residues, we examine many cation binding configurations using molecular dynamics simulations, predicting that K+ binds to the so-called Na2 site. Site-directed mutagenesis of residues in this site can eliminate the ability of both K+ and H+ to drive 5-HT accumulation into vesicles and, in patch clamp recordings, prevent the acceleration of turnover rates and the formation of a channel-like state by K+ or H+. In conclusion, the Na2 site plays a pivotal role in orchestrating the sequential binding of Na+ and then K+ (or H+) ions to facilitate 5-HT uptake in SERT.

Mechanism of proton-powered c-ring rotation in a mitochondrial ATP synthase.

Proton-powered c-ring rotation in mitochondrial ATP synthase is crucial to convert the transmembrane protonmotive force into torque to drive the synthesis of adenosine triphosphate (ATP). Capitalizing on recent cryo-EM structures, we aim at a structural and energetic understanding of how functional directional rotation is achieved. We performed multi-microsecond atomistic simulations to determine the free energy profiles along the c-ring rotation angle before and after the arrival of a new proton. Our results reveal that rotation proceeds by dynamic sliding of the ring over the a-subunit surface, during which interactions with conserved polar residues stabilize distinct intermediates. Ordered water chains line up for a Grotthuss-type proton transfer in one of these intermediates. After proton transfer, a high barrier prevents backward rotation and an overall drop in free energy favors forward rotation, ensuring the directionality of c-ring rotation required for the thermodynamically disfavored ATP synthesis. The essential arginine of the a-subunit stabilizes the rotated configuration through a salt bridge with the c-ring. Overall, we describe a complete mechanism for the rotation step of the ATP synthase rotor, thereby illuminating a process critical to all life at atomic resolution.

Sound-mediated nucleation and growth of amyloid fibrils.

Mechanical energy, specifically in the form of ultrasound, can induce pressure variations and temperature fluctuations when applied to an aqueous media. These conditions can both positively and negatively affect protein complexes, consequently altering their stability, folding patterns, and self-assembling behavior. Despite much scientific progress, our current understanding of the effects of ultrasound on the self-assembly of amyloidogenic proteins remains limited. In the present study, we demonstrate that when the amplitude of the delivered ultrasonic energy is sufficiently low, it can induce refolding of specific motifs in protein monomers, which is sufficient for primary nucleation; this has been revealed by MD. These ultrasound-induced structural changes are initiated by pressure perturbations and are accelerated by a temperature factor. Furthermore, the prolonged action of low-amplitude ultrasound enables the elongation of amyloid protein nanofibrils directly from natively folded monomeric lysozyme protein, in a controlled manner, until it reaches a critical length. Using solution X-ray scattering, we determined that nanofibrillar assemblies, formed either under the action of sound or from natively fibrillated lysozyme, share identical structural characteristics. Thus, these results provide insights into the effects of ultrasound on fibrillar protein self-assembly and lay the foundation for the potential use of sound energy in protein chemistry.

An α-ketoglutarate conformational switch controls iron accessibility, activation, and substrate selection of the human FTO protein.

The fat mass and obesity-associated fatso (FTO) protein is a member of the Alkb family of dioxygenases and catalyzes oxidative demethylation of N6-methyladenosine (m6A), N1-methyladenosine (m1A), 3-methylthymine (m3T), and 3-methyluracil (m3U) in single-stranded nucleic acids. It is well established that the catalytic activity of FTO proceeds via two coupled reactions. The first reaction involves decarboxylation of alpha-ketoglutarate (αKG) and formation of an oxyferryl species. In the second reaction, the oxyferryl intermediate oxidizes the methylated nucleic acid to reestablish Fe(II) and the canonical base. However, it remains unclear how binding of the nucleic acid activates the αKG decarboxylation reaction and why FTO demethylates different methyl modifications at different rates. Here, we investigate the interaction of FTO with 5-mer DNA oligos incorporating the m6A, m1A, or m3T modifications using solution NMR, molecular dynamics (MD) simulations, and enzymatic assays. We show that binding of the nucleic acid to FTO activates a two-state conformational equilibrium in the αKG cosubstrate that modulates the O2 accessibility of the Fe(II) catalyst. Notably, the substrates that provide better stabilization to the αKG conformation in which Fe(II) is exposed to O2 are demethylated more efficiently by FTO. These results indicate that i) binding of the methylated nucleic acid is required to expose the catalytic metal to O2 and activate the αKG decarboxylation reaction, and ii) the measured turnover of the demethylation reaction (which is an ensemble average over the entire sample) depends on the ability of the methylated base to favor the Fe(II) state accessible to O2.

Secondary structure determines electron transport in peptides.

Proteins play a key role in biological electron transport, but the structure-function relationships governing the electronic properties of peptides are not fully understood. Despite recent progress, understanding the link between peptide conformational flexibility, hierarchical structures, and electron transport pathways has been challenging. Here, we use single-molecule experiments, molecular dynamics (MD) simulations, nonequilibrium Green's function-density functional theory (NEGF-DFT), and unsupervised machine learning to understand the role of secondary structure on electron transport in peptides. Our results reveal a two-state molecular conductance behavior for peptides across several different amino acid sequences. MD simulations and Gaussian mixture modeling are used to show that this two-state molecular conductance behavior arises due to the conformational flexibility of peptide backbones, with a high-conductance state arising due to a more defined secondary structure (beta turn or 310 helices) and a low-conductance state occurring for extended peptide structures. These results highlight the importance of helical conformations on electron transport in peptides. Conformer selection for the peptide structures is rationalized using principal component analysis of intramolecular hydrogen bonding distances along peptide backbones. Molecular conformations from MD simulations are used to model charge transport in NEGF-DFT calculations, and the results are in reasonable qualitative agreement with experiments. Projected density of states calculations and molecular orbital visualizations are further used to understand the role of amino acid side chains on transport. Overall, our results show that secondary structure plays a key role in electron transport in peptides, which provides broad avenues for understanding the electronic properties of proteins.

Comprehensive classification of TP53 somatic missense variants based on their impact on p53 structural stability.

Somatic variation is a major type of genetic variation contributing to human diseases including cancer. Of the vast quantities of somatic variants identified, the functional impact of many somatic variants, in particular the missense variants, remains unclear. Lack of the functional information prevents the translation of rich variation data into clinical applications. We previously developed a method named Ramachandran Plot-Molecular Dynamics Simulations (RP-MDS), aiming to predict the function of germline missense variants based on their effects on protein structure stability, and successfully applied to predict the deleteriousness of unclassified germline missense variants in multiple cancer genes. We hypothesized that regardless of their different genetic origins, somatic missense variants and germline missense variants could have similar effects on the stability of their affected protein structure. As such, the RP-MDS method designed for germline missense variants should also be applicable to predict the function of somatic missense variants. In the current study, we tested our hypothesis by using the somatic missense variants in TP53 as a model. Of the 397 somatic missense variants analyzed, RP-MDS predicted that 195 (49.1%) variants were deleterious as they significantly disturbed p53 structure. The results were largely validated by using a p53-p21 promoter-green fluorescent protein (GFP) reporter gene assay. Our study demonstrated that deleterious somatic missense variants can be identified by referring to their effects on protein structural stability.

Water dissociation at the water-rutile TiO2(110) interface from ab initio-based deep neural network simulations.

The interaction of water with TiO2 surfaces is of crucial importance in various scientific fields and applications, from photocatalysis for hydrogen production and the photooxidation of organic pollutants to self-cleaning surfaces and bio-medical devices. In particular, the equilibrium fraction of water dissociation at the TiO2-water interface has a critical role in the surface chemistry of TiO2, but is difficult to determine both experimentally and computationally. Among TiO2 surfaces, rutile TiO2(110) is of special interest as the most abundant surface of TiO2's stable rutile phase. While surface-science studies have provided detailed information on the interaction of rutile TiO2(110) with gas-phase water, much less is known about the TiO2(110)-water interface, which is more relevant to many applications. In this work, we characterize the structure of the aqueous TiO2(110) interface using nanosecond timescale molecular dynamics simulations with ab initio-based deep neural network potentials that accurately describe water/TiO2(110) interactions over a wide range of water coverages. Simulations on TiO2(110) slab models of increasing thickness provide insight into the dynamic equilibrium between molecular and dissociated adsorbed water at the interface and allow us to obtain a reliable estimate of the equilibrium fraction of water dissociation. We find a dissociation fraction of 22 ± 6% with an associated average hydroxyl lifetime of 7.6 ± 1.8 ns. These quantities are both much larger than corresponding estimates for the aqueous anatase TiO2(101) interface, consistent with the higher water photooxidation activity that is observed for rutile relative to anatase.

Phospholamban inhibits the cardiac calcium pump by interrupting an allosteric activation pathway.

Phospholamban (PLB) is a transmembrane micropeptide that regulates the sarcoplasmic reticulum Ca2+-ATPase (SERCA) in cardiac muscle, but the physical mechanism of this regulation remains poorly understood. PLB reduces the Ca2+ sensitivity of active SERCA, increasing the Ca2+ concentration required for pump cycling. However, PLB does not decrease Ca2+ binding to SERCA when ATP is absent, suggesting PLB does not inhibit SERCA Ca2+ affinity. The prevailing explanation for these seemingly conflicting results is that PLB slows transitions in the SERCA enzymatic cycle associated with Ca2+ binding, altering transport Ca2+ dependence without actually affecting the equilibrium binding affinity of the Ca2+-coordinating sites. Here, we consider another hypothesis, that measurements of Ca2+ binding in the absence of ATP overlook important allosteric effects of nucleotide binding that increase SERCA Ca2+ binding affinity. We speculated that PLB inhibits SERCA by reversing this allostery. To test this, we used a fluorescent SERCA biosensor to quantify the Ca2+ affinity of non-cycling SERCA in the presence and absence of a non-hydrolyzable ATP-analog, AMPPCP. Nucleotide activation increased SERCA Ca2+ affinity, and this effect was reversed by co-expression of PLB. Interestingly, PLB had no effect on Ca2+ affinity in the absence of nucleotide. These results reconcile the previous conflicting observations from ATPase assays versus Ca2+ binding assays. Moreover, structural analysis of SERCA revealed a novel allosteric pathway connecting the ATP- and Ca2+-binding sites. We propose this pathway is disrupted by PLB binding. Thus, PLB reduces the equilibrium Ca2+ affinity of SERCA by interrupting allosteric activation of the pump by ATP.

Probing the role of protein conformational changes in the mechanism of prenylated-FMN-dependent phenazine-1-carboxylic acid decarboxylase.

Phenazine-1-carboxylic acid decarboxylase (PhdA) is a prenylated-FMN-dependent (prFMN) enzyme belonging to the UbiD family of decarboxylases. Many UbiD-like enzymes catalyze (de)carboxylation reactions on aromatic rings and conjugated double bonds and are potentially valuable industrial catalysts. We have investigated the mechanism of PhdA using a slow turnover substrate, 2,3-dimethylquinoxaline-5-carboxylic acid (DQCA). Detailed analysis of the pH dependence and solvent deuterium isotope effects associated with the reaction uncovered unusual kinetic behavior. At low substrate concentrations, a substantial inverse solvent isotope effect (SIE) is observed on Vmax/KM of ∼ 0.5 when reaction rates of DQCA in H2O and D2O are compared. Under the same conditions, a normal SIE of 4.15 is measured by internal competition for proton transfer to the product. These apparently contradictory results indicate that the SIE values report on different steps in the mechanism. A proton inventory analysis of the reaction under Vmax/KM and Vmax conditions points to a "medium effect" as the source of the inverse SIE. Molecular dynamics simulations of the effect of D2O on PhdA structure support that D2O reduces the conformational lability of the enzyme and results in a more compact structure, akin to the active, "closed" conformer observed in crystal structures of some UbiD-like enzymes. Consistent with the simulations, PhdA was found to be more stable in D2O and to bind DQCA more tightly, leading to the observed rate enhancement under Vmax/KM conditions.

Substitution Models of Protein Evolution with Selection on Enzymatic Activity.

Substitution models of evolution are necessary for diverse evolutionary analyses including phylogenetic tree and ancestral sequence reconstructions. At the protein level, empirical substitution models are traditionally used due to their simplicity, but they ignore the variability of substitution patterns among protein sites. Next, in order to improve the realism of the modeling of protein evolution, a series of structurally constrained substitution models were presented, but still they usually ignore constraints on the protein activity. Here, we present a substitution model of protein evolution with selection on both protein structure and enzymatic activity, and that can be applied to phylogenetics. In particular, the model considers the binding affinity of the enzyme-substrate complex as well as structural constraints that include the flexibility of structural flaps, hydrogen bonds, amino acids backbone radius of gyration, and solvent-accessible surface area that are quantified through molecular dynamics simulations. We applied the model to the HIV-1 protease and evaluated it by phylogenetic likelihood in comparison with the best-fitting empirical substitution model and a structurally constrained substitution model that ignores the enzymatic activity. We found that accounting for selection on the protein activity improves the fitting of the modeled functional regions with the real observations, especially in data with high molecular identity, which recommends considering constraints on the protein activity in the development of substitution models of evolution.