bims-lances 2024-02-25 papers

Biochemistry. 2024 Feb 20.

Lipid-Dependent Activation of the Orphan G Protein-Coupled Receptor, GPR3.

Isabella C Russell, Xin Zhang, Fabian Bumbak, Samantha M McNeill, Tracy M Josephs, Michael G Leeming, George Christopoulos, Hariprasad Venugopal, Maria M Flocco, Patrick M Sexton, Denise Wootten, Matthew J Belousoff.

The class A orphan G protein-coupled receptor (GPCR), GPR3, has been implicated in a variety of conditions, including Alzheimer's and premature ovarian failure. GPR3 constitutively couples with Gαs, resulting in the production of cAMP in cells. While tool compounds and several putative endogenous ligands have emerged for the receptor, its endogenous ligand, if it exists, remains a mystery. As novel potential drug targets, the structures of orphan GPCRs have been of increasing interest, revealing distinct modes of activation, including autoactivation, presence of constitutively activating mutations, or via cryptic ligands. Here, we present a cryo-electron microscopy (cryo-EM) structure of the orphan GPCR, GPR3 in complex with DNGαs and Gβ1γ2. The structure revealed clear density for a lipid-like ligand that bound within an extended hydrophobic groove, suggesting that the observed "constitutive activity" was likely due to activation via a lipid that may be ubiquitously present. Analysis of conformational variance within the cryo-EM data set revealed twisting motions of the GPR3 transmembrane helices that appeared coordinated with changes in the lipid-like density. We propose a mechanism for the binding of a lipid to its putative orthosteric binding pocket linked to the GPR3 dynamics.

DOI: https://doi.org/10.1021/acs.biochem.3c00647

Cell. 2024 Feb 13. pii: S0092-8674(24)00099-0. [Epub ahead of print]

Structural mechanisms of α7 nicotinic receptor allosteric modulation and activation.

Sean M Burke, Mariia Avstrikova, Colleen M Noviello, Nuriya Mukhtasimova, Jean-Pierre Changeux, Ganesh A Thakur, Steven M Sine, Marco Cecchini, Ryan E Hibbs.

The α7 nicotinic acetylcholine receptor is a pentameric ligand-gated ion channel that plays an important role in cholinergic signaling throughout the nervous system. Its unique physiological characteristics and implications in neurological disorders and inflammation make it a promising but challenging therapeutic target. Positive allosteric modulators overcome limitations of traditional α7 agonists, but their potentiation mechanisms remain unclear. Here, we present high-resolution structures of α7-modulator complexes, revealing partially overlapping binding sites but varying conformational states. Structure-guided functional and computational tests suggest that differences in modulator activity arise from the stable rotation of a channel gating residue out of the pore. We extend the study using a time-resolved cryoelectron microscopy (cryo-EM) approach to reveal asymmetric state transitions for this homomeric channel and also find that a modulator with allosteric agonist activity exploits a distinct channel-gating mechanism. These results define mechanisms of α7 allosteric modulation and activation with implications across the pentameric receptor superfamily.

Keywords: cryo-EM; electrophysiology; ligand-gated ion channel; molecular dynamics; neurotransmitter receptor; nicotinic acetylcholine receptor; positive allosteric modulator; structural biology; structure-guided drug design

DOI: https://doi.org/10.1016/j.cell.2024.01.032

Inorg Chem. 2024 Feb 20.

Multiscale Study on the Intramolecular C-S Bond Formation Catalyzed by P450 Monooxygenase CxnD Involved in the Biosynthesis of Chuangxinmycin: The Critical Roles of Noncrystal Water Molecule and Conformational Change.

Xinyi Li, Yongjun Liu.

Cytochrome P450 monooxygenase CxnD catalyzes intramolecular C-S bond formation in the biosynthesis of chuangxinmycin, which is representative of the synthesis of sulfur-containing natural heterocyclic compounds. The intramolecular cyclization usually requires the activation of two reaction sites and a large conformational change; thus, illuminating its detailed reaction mechanism remains challengeable. Here, the reaction pathway of CxnD-catalyzed C-S bond formation was clarified by a series of calculations, including Gaussian accelerated molecular dynamics simulations and quantum mechanical-molecular mechanical calculations. Our results revealed that the C-S formation follows a diradical coupling mechanism. CxnD first employs Cpd I to abstract the hydrogen atom from the imino group of the indole ring, and then, the resulted Cpd II further extracts another hydrogen atom from the thiol group of the side chain to afford a diradical intermediate, in which a noncrystal water molecule entering into the active site after the formation of Cpd I was proved to play an indispensable role. Moreover, the diradical intermediate cannot directly perform the coupling reaction. It should first undergo a series of conformational changes leading to the proximity of two reaction sites. It is the flexibility of the active site of the enzyme and the side chain of the substrate that makes the diradical coupling to be successful.

DOI: https://doi.org/10.1021/acs.inorgchem.3c03748

Struct Dyn. 2024 Jan;11(1): 014702

Conformational dynamics of adenylate kinase in crystals.

Junhyung Kim, Sojin Moon, Tod D Romo, Yifei Yang, Euiyoung Bae, George N Phillips.

Adenylate kinase is a ubiquitous enzyme in living systems and undergoes dramatic conformational changes during its catalytic cycle. For these reasons, it is widely studied by genetic, biochemical, and biophysical methods, both experimental and theoretical. We have determined the basic crystal structures of three differently liganded states of adenylate kinase from Methanotorrus igneus, a hyperthermophilic organism whose adenylate kinase is a homotrimeric oligomer. The multiple copies of each protomer in the asymmetric unit of the crystal provide a unique opportunity to study the variation in the structure and were further analyzed using advanced crystallographic refinement methods and analysis tools to reveal conformational heterogeneity and, thus, implied dynamic behaviors in the catalytic cycle.

DOI: https://doi.org/10.1063/4.0000205

J Biomol Struct Dyn. 2024 Feb 23. 1-16

The conformational dynamics of Hepatitis C Virus E2 glycoprotein with the increasing number of N-glycosylation unraveled by molecular dynamics simulations.

Sayan Poddar, Rajarshi Roy, Parimal Kar.

The Hepatitis C Virus (HCV), responsible for causing hepatitis and a significant contributor to liver disorders, presents a challenge for treatment due to its high genetic variability. Despite efforts, there is still no effective medication available for this virus. One of the promising targets for drug development involves targeting glycoprotein E2. However, our understanding of the dynamic behavior of E2 and its associated glycans remains limited. In this study, we investigated the dynamic characteristics of E2 with varying degrees of glycosylation using all-atom molecular dynamics simulations. We also explored glycan's interactions with the protein and among themselves. An overall increase in correlation between the vital protein regions was observed with an increase in glycan number. The protein dynamics is followed by the analysis of glycan dynamics, where the flexibility of the individual glycans was analyzed in their free and bound state, which revealed a decrease in their fluctuation in some cases. Furthermore, we generated the free energy landscape of individual N-glycan linkages in both free and bound states and observed both increases and decreases in flexibility, which can be attributed to the formation and breakage of hydrogen bonds with amino acids. Finally, we found that for a high glycosylation system, glycans interact with glycoprotein and form hydrogen bonds among themselves. Moreover, the hydrogen bond profiles of a given glycan can vary when influenced by other glycans. In summary, our study provides valuable insights into the dynamics of the core region of HCV E2 glycoprotein and its associated glycans.Communicated by Ramaswamy H. Sarma.

Keywords: E2 glycoprotein; Hepatitis C virus (HCV); N-glycans; conformational dynamics; molecular dynamics

DOI: https://doi.org/10.1080/07391102.2024.2319679

J Chem Inf Model. 2024 Feb 19.

AlphaFold2-Enabled Atomistic Modeling of Structure, Conformational Ensembles, and Binding Energetics of the SARS-CoV-2 Omicron BA.2.86 Spike Protein with ACE2 Host Receptor and Antibodies: Compensatory Functional Effects of Binding Hotspots in Modulating Mechanisms of Receptor Binding and Immune Escape.

Nishank Raisinghani, Mohammed Alshahrani, Grace Gupta, Sian Xiao, Peng Tao, Gennady Verkhivker.

The latest wave of SARS-CoV-2 Omicron variants displayed a growth advantage and increased viral fitness through convergent evolution of functional hotspots that work synchronously to balance fitness requirements for productive receptor binding and efficient immune evasion. In this study, we combined AlphaFold2-based structural modeling approaches with atomistic simulations and mutational profiling of binding energetics and stability for prediction and comprehensive analysis of the structure, dynamics, and binding of the SARS-CoV-2 Omicron BA.2.86 spike variant with ACE2 host receptor and distinct classes of antibodies. We adapted several AlphaFold2 approaches to predict both the structure and conformational ensembles of the Omicron BA.2.86 spike protein in the complex with the host receptor. The results showed that the AlphaFold2-predicted structural ensemble of the BA.2.86 spike protein complex with ACE2 can accurately capture the main conformational states of the Omicron variant. Complementary to AlphaFold2 structural predictions, microsecond molecular dynamics simulations reveal the details of the conformational landscape and produced equilibrium ensembles of the BA.2.86 structures that are used to perform mutational scanning of spike residues and characterize structural stability and binding energy hotspots. The ensemble-based mutational profiling of the receptor binding domain residues in the BA.2 and BA.2.86 spike complexes with ACE2 revealed a group of conserved hydrophobic hotspots and critical variant-specific contributions of the BA.2.86 convergent mutational hotspots R403K, F486P, and R493Q. To examine the immune evasion properties of BA.2.86 in atomistic detail, we performed structure-based mutational profiling of the spike protein binding interfaces with distinct classes of antibodies that displayed significantly reduced neutralization against the BA.2.86 variant. The results revealed the molecular basis of compensatory functional effects of the binding hotspots, showing that BA.2.86 lineage may have evolved to outcompete other Omicron subvariants by improving immune evasion while preserving binding affinity with ACE2 via through a compensatory effect of R493Q and F486P convergent mutational hotspots. This study demonstrated that an integrative approach combining AlphaFold2 predictions with complementary atomistic molecular dynamics simulations and robust ensemble-based mutational profiling of spike residues can enable accurate and comprehensive characterization of structure, dynamics, and binding mechanisms of newly emerging Omicron variants.

DOI: https://doi.org/10.1021/acs.jcim.3c01857

Bioinformatics. 2024 Feb 23. pii: btae106. [Epub ahead of print]

Optimizing representations for integrative structural modeling using bayesian model selection.

Shreyas Arvindekar, Aditi S Pathak, Kartik Majila, Shruthi Viswanath.

MOTIVATION: Integrative structural modeling combines data from experiments, physical principles, statistics of previous structures, and prior models to obtain structures of macromolecular assemblies that are challenging to characterize experimentally. The choice of model representation is a key decision in integrative modeling, as it dictates the accuracy of scoring, efficiency of sampling, and resolution of analysis. But currently, the choice is usually made ad hoc, manually.
RESULTS: Here, we report NestOR (Nested Sampling for Optimizing Representation), a fully automated, statistically rigorous method based on Bayesian model selection to identify the optimal coarse-grained representation for a given integrative modeling setup. Given an integrative modeling setup, it determines the optimal representations from given candidate representations based on their model evidence and sampling efficiency. The performance of NestOR was evaluated on a benchmark of four macromolecular assemblies.
AVAILABILITY: NestOR is implemented in the Integrative Modeling Platform (https://integrativemodeling.org) and is available at https://github.com/isblab/nestor.Data for the benchmark is at https://www.doi.org/10.5281/zenodo.10360718.
SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

DOI: https://doi.org/10.1093/bioinformatics/btae106