G protein-coupled receptors (GPCRs) are ubiquitous and essential in modulating virtually all physiological processes. These receptors share a similar structural design consisting of the seven-transmembrane alpha-helical segments. The active conformations of the receptors are stabilized by an agonist and couple to structurally highly conserved heterotrimeric G proteins. One of the most important unanswered questions is how GPCRs couple to their cognate G proteins. Phototransduction represents an excellent model system for understanding G protein signaling, owing to the high expression of rhodopsin in rod photoreceptors and the multidisciplinary experimental approaches used to study this GPCR. Here, we describe how a G protein (transducin) docks on to an oligomeric GPCR (rhodopsin), revealing structural details of this critical interface in the signal transduction process. This conceptual model takes into account recent structural information on the receptor and G protein, as well as oligomeric states of GPCRs.

No abstract

No abstract

Vasoactive intestinal peptide receptor 1 (VPAC1) is a member of a secretin-like subfamily of G protein-coupled receptors. Its endogenous neuropeptide (VIP), secreted by neurons and immune cells, modulates various physiological functions such as exocrine and endocrine secretions, immune response, smooth muscles relaxation, vasodilation, and fetal development. As a drug target, VPAC1 has been selected for therapy of inflammatory diseases but drug discovery is still hampered by lack of its crystal structure. In this study we presented the homology model of this receptor constructed with the well-known web service GPCRM. The VPAC1 model is composed of extracellular and transmembrane domains that form a complex with an endogenous hormone VIP. Using the homology model of VPAC1 the mechanism of action of potential drug candidates for VPAC1 was described. Only two series of small-molecule antagonists of confirmed biological activity for VPAC1 have been described thus far. Molecular docking and a series of molecular dynamics simulations were performed to elucidate their binding to VPAC1 and resulting antagonist effect. The presented work provides the basis for the possible binding mode of VPAC1 antagonists and determinants of their molecular recognition in the context of other class B GPCRs. Until the crystal structure of VPAC1 will be released, the presented homology model of VPAC1 can serve as a scaffold for drug discovery studies and is available from the author upon request.

Recently, molecular covalent docking has been extensively developed to design new classes of inhibitors that form chemical bonds with their biological targets. This strategy for the design of such inhibitors, in particular boron-based inhibitors, holds great promise for the vast family of β-lactamases produced, inter alia, by Gram-negative antibiotic-resistant bacteria. However, the description of covalent docking processes requires a quantum-mechanical approach, and so far, only a few studies of this type have been presented. This study accurately describes the covalent docking process between two model inhibitors - representing two large families of inhibitors based on boronic-acid and bicyclic boronate scaffolds, and three β-lactamases which belong to the A, C, and D classes. Molecular fragments containing boron can be converted from a neutral, trigonal, planar state with sp2 hybridization to the anionic, tetrahedral sp3 state in a process sometimes referred to as morphing. This study applies multi-scale modeling methods, in particular, the hybrid QM/MM approach which has predictive power reaching well beyond conventional molecular modeling. Time-dependent QM/MM simulations indicated several structural changes and geometric preferences, ultimately leading to covalent docking processes. With current computing technologies, this approach is not computationally expensive, can be used in standard molecular modeling and molecular design works, and can effectively support experimental research which should allow for a detailed understanding of complex processes important to molecular medicine. In particular, it can support the rational design of covalent boron-based inhibitors for β-lactamases as well as for many other enzyme systems of clinical relevance, including SARS-CoV-2 proteins.

No abstract

Charcot-Marie-Tooth disease type 2A (CMT2A) is an autosomal dominant axonal peripheral neuropathy caused by mutations in the mitofusin 2 gene (MFN2). Mitofusin 2 is a GTPase protein present in the outer mitochondrial membrane and responsible for regulation of mitochondrial network architecture via the fusion of mitochondria. As that fusion process is known to be strongly dependent on the GTPase activity of mitofusin 2, it is postulated that the MFN2 mutation within the GTPase domain may lead to impaired GTPase activity, and in turn to mitochondrial dysfunction. The work described here has therefore sought to verify the effects of MFN2 mutation within its GTPase domain on mitochondrial and endoplasmic reticulum morphology, as well as the mtDNA content in a cultured primary fibroblast obtained from a CMT2A patient harboring a de novo Arg274Trp mutation. In fact, all the parameters studied were affected significantly by the presence of the mutant MFN2 protein. However, using the stable model for mitofusin 2 obtained by us, we were next able to determine that the Arg274Trp mutation does not impact directly upon GTP binding. Such results were also confirmed for GTP-hydrolysis activity of MFN2 protein in patient fibroblast. We therefore suggest that the biological malfunctions observable with the disease are not consequences of impaired GTPase activity, but rather reflect an impaired contribution of the GTPase domain to other MFN2 activities involving that region, for example protein-protein interactions.

Mechanical unfolding of single bacteriorhodopsins from a membrane bilayer is studied using molecular dynamics simulations. The initial conformation of the lipid membrane is determined through all-atom simulations and then its coarse-grained representation is used in the studies of stretching. A Go-like model with a realistic contact map and with Lennard-Jones contact interactions is applied to model the protein-membrane system. The model qualitatively reproduces the experimentally observed differences between force-extension patterns obtained on bacteriorhodopsin at different temperatures and predicts a lack of symmetry in the choice of the terminus to pull by. It also illustrates the decisive role of the interactions of the protein with the membrane in determining the force pattern and thus the stability of transmembrane proteins.