3D interaction homology: Hydropathic interaction environments of serine and cysteine are strikingly different and their roles adapt in membrane proteins

Characterization of myoinhibitory peptide signaling system and its implication in larval metamorphosis and spawning behavior in Pacific abalone

Myoinhibitory peptides (MIPs) affect various physiological functions, including juvenile hormone signaling, muscle contraction, larval development, and reproduction in invertebrates. Although MIPs are ligands for MIP and/or sex peptide receptors (MIP/SPRs) in diverse arthropods and model organisms belonging to Lophotrochozoa, the MIP signaling system has not yet been fully investigated in mollusks. In this study, we identified the MIP signaling system in the Pacific abalone Haliotis discus hannai (Hdh). Similar to the invertebrate MIPs, a total of eight paracopies of MIPs (named Hdh-MIP1 to Hdh-MIP8), harboring a WX5-7Wamide motif, except for Hdh-MIP2, were found in the Hdh-MIP precursor. Furthermore, we characterized a functional Hdh-MIPR, which responded to the Hdh-MIPs, except for Hdh-MIP2, possibly linked with the PKC/Ca2+ and PKA/cAMP signaling pathways. Hdh-MIPs delayed larval metamorphosis but increased the spawning behavior. These results suggest that the Hdh-MIP signaling system provides insights into the unique function of MIP in invertebrates.

Deciphering Molecular Mechanisms Involved in the Modulation of Human Aquaporins' Water Permeability by Zinc Cations: A Molecular Dynamics Approach

Aquaporins (AQPs) constitute a wide family of water channels implicated in all kind of physiological processes. Zinc is the second most abundant trace element in the human body and a few studies have highlighted regulation of AQP0 and AQP4 by zinc. In the present work, we addressed the putative regulation of AQPs by zinc cations in silico through molecular dynamics simulations of human AQP0, AQP2, AQP4, and AQP5. Our results align with other scales of study and several in vitro techniques, hence strengthening the reliability of this regulation by zinc. We also described two distinct putative molecular mechanisms associated with the increase or decrease in AQPs' water permeability after zinc binding. In association with other studies, our work will help deciphering the interaction networks existing between zinc and channel proteins.

Structure-Function Characterisation of Eop1 Effectors from the Erwinia-Pantoea Clade Reveals They May Acetylate Their Defence Target through a Catalytic Dyad

The YopJ group of acetylating effectors from phytopathogens of the genera Pseudomonas and Ralstonia have been widely studied to understand how they modify and suppress their host defence targets. In contrast, studies on a related group of effectors, the Eop1 group, lag far behind. Members of the Eop1 group are widely present in the Erwinia-Pantoea clade of Gram-negative bacteria, which contains phytopathogens, non-pathogens and potential biocontrol agents, implying that they may play an important role in agroecological or pathological adaptations. The lack of research in this group of YopJ effectors has left a significant knowledge gap in their functioning and role. For the first time, we perform a comparative analysis combining AlphaFold modelling, in planta transient expressions and targeted mutational analyses of the Eop1 group effectors from the Erwinia-Pantoea clade, to help elucidate their likely activity and mechanism(s). This integrated study revealed several new findings, including putative binding sites for inositol hexakisphosphate and acetyl coenzyme A and newly postulated target-binding domains, and raises questions about whether these effectors function through a catalytic triad mechanism. The results imply that some Eop1s may use a catalytic dyad acetylation mechanism that we found could be promoted by the electronegative environment around the active site.

Merging cultures and disciplines to create a drug discovery ecosystem at Virginia commonwealth university: Medicinal chemistry, structural biology, molecular and behavioral pharmacology and computational chemistry

The Department of Medicinal Chemistry, together with the Institute for Structural Biology, Drug Discovery and Development, at Virginia Commonwealth University (VCU) has evolved, organically with quite a bit of bootstrapping, into a unique drug discovery ecosystem in response to the environment and culture of the university and the wider research enterprise. Each faculty member that joined the department and/or institute added a layer of expertise, technology and most importantly, innovation, that fertilized numerous collaborations within the University and with outside partners. Despite moderate institutional support with respect to a typical drug discovery enterprise, the VCU drug discovery ecosystem has built and maintained an impressive array of facilities and instrumentation for drug synthesis, drug characterization, biomolecular structural analysis and biophysical analysis, and pharmacological studies. Altogether, this ecosystem has had major impacts on numerous therapeutic areas, such as neurology, psychiatry, drugs of abuse, cancer, sickle cell disease, coagulopathy, inflammation, aging disorders and others. Novel tools and strategies for drug discovery, design and development have been developed at VCU in the last five decades; e.g., fundamental rational structure-activity relationship (SAR)-based drug design, structure-based drug design, orthosteric and allosteric drug design, design of multi-functional agents towards polypharmacy outcomes, principles on designing glycosaminoglycans as drugs, and computational tools and algorithms for quantitative SAR (QSAR) and understanding the roles of water and the hydrophobic effect.

A force-sensitive mutation reveals a spindle assembly checkpoint-independent role for dynein in anaphase progression

The cytoplasmic dynein-1 (dynein) motor organizes cells by shaping microtubule networks and moving a large variety of cargoes along them. However, dynein's diverse roles complicate in vivo studies of its functions significantly. To address this issue, we have used gene editing to generate a series of missense mutations in Drosophila Dynein heavy chain (Dhc). We find that mutations associated with human neurological disease cause a range of defects in larval and adult flies, including impaired cargo trafficking in neurons. We also describe a novel mutation in the microtubule-binding domain (MTBD) of Dhc that, remarkably, causes metaphase arrest of mitotic spindles in the embryo but does not impair other dynein-dependent processes. We demonstrate that the mitotic arrest is independent of dynein's well-established roles in silencing the spindle assembly checkpoint. In vitro reconstitution and optical trapping assays reveal that the mutation only impairs the performance of dynein under load. In silico all-atom molecular dynamics simulations show that this effect correlates with increased flexibility of the MTBD, as well as an altered orientation of the stalk domain, with respect to the microtubule. Collectively, our data point to a novel role of dynein in anaphase progression that depends on the motor operating in a specific load regime. More broadly, our work illustrates how cytoskeletal transport processes can be dissected in vivo by manipulating mechanical properties of motors.

HINT, a code for understanding the interaction between biomolecules: a tribute to Donald J. Abraham

A long-lasting goal of computational biochemists, medicinal chemists, and structural biologists has been the development of tools capable of deciphering the molecule-molecule interaction code that produces a rich variety of complex biomolecular assemblies comprised of the many different simple and biological molecules of life: water, small metabolites, cofactors, substrates, proteins, DNAs, and RNAs. Software applications that can mimic the interactions amongst all of these species, taking account of the laws of thermodynamics, would help gain information for understanding qualitatively and quantitatively key determinants contributing to the energetics of the bimolecular recognition process. This, in turn, would allow the design of novel compounds that might bind at the intermolecular interface by either preventing or reinforcing the recognition. HINT, hydropathic interaction, was a model and software code developed from a deceptively simple idea of Donald Abraham with the close collaboration with Glen Kellogg at Virginia Commonwealth University. HINT is based on a function that scores atom-atom interaction using LogP, the partition coefficient of any molecule between two phases; here, the solvents are water that mimics the cytoplasm milieu and octanol that mimics the protein internal hydropathic environment. This review summarizes the results of the extensive and successful collaboration between Abraham and Kellogg at VCU and the group at the University of Parma for testing HINT in a variety of different biomolecular interactions, from proteins with ligands to proteins with DNA.

3D interaction homology: The hydrophobic residues alanine, isoleucine, leucine, proline and valine play different structural roles in soluble and membrane proteins

The aliphatic hydrophobic amino acid residues—alanine, isoleucine, leucine, proline and valine—are among the most common found in proteins. Their structural role in proteins is seemingly obvious: engage in hydrophobic interactions to stabilize secondary, and to a lesser extent, tertiary and quaternary structure. However, favorable hydrophobic interactions involving the sidechains of these residue types are generally less significant than the unfavorable set arising from interactions with polar atoms. Importantly, the constellation of interactions between residue sidechains and their environments can be recorded as three-dimensional maps that, in turn, can be clustered. The clustered average map sets compose a library of interaction profiles encoding interaction strengths, interaction types and the optimal 3D position for the interacting partners. This library is backbone angle-dependent and suggests solvent and lipid accessibility for each unique interaction profile. In this work, in addition to analysis of soluble proteins, a large set of membrane proteins that contained optimized artificial lipids were evaluated by parsing the structures into three distinct components: soluble extramembrane domain, lipid facing transmembrane domain, core transmembrane domain. The aliphatic residues were extracted from each of these sets and passed through our calculation protocol. Notable observations include: the roles of aliphatic residues in soluble proteins and in the membrane protein’s soluble domains are nearly identical, although the latter are slightly more solvent accessible; by comparing maps calculated with sidechain-lipid interactions to maps ignoring those interactions, the potential extent of residue-lipid and residue-interactions can be assessed and likely exploited in structure prediction and modeling; amongst these residue types, the levels of lipid engagement show isoleucine as the most engaged, while the other residues are largely interacting with neighboring helical residues.

Structural elucidation of ETHR-A and ETHR-B from Plutella xylostella and insight into non-conservative mutations in transmembrane helix-6

The development of Diamondback moth (DBM) depends on the ecdysis triggering hormone receptor (ETHR); a neuronal membrane G-protein coupled receptor (GPCR) connected to the metamorphosis cascade. Lepidopteran insect DBM is an infamous pest of cruciferous plants. This study examined the full-length coding sequences (CDS) of PxETHR-A and PxETHR-B from the DBM genome. The three-dimensional (3 D) models of both receptors were generated in an inactive state. The behaviour and stability of receptors were examined using molecular dynamics simulations in a lipid membrane system for 300 ns and established a GPCR family-based view. Secondary interactions within receptors were studied to know more about factors contributing to their stability. Multiple sequence alignment revealed conserved features of insect ETHRs those compared to the GPCR family proteins. These features were helpful during the evaluation of the molecular models of both receptors. Side-chain orientation of conserved residues, non-conserved and conserved hydrogen-bond networks (HBN) and hydrophobic clusters were examined in the structures of both receptors. The non-conserved residues L6.35, T6.39, C/S6.43, and L6.48, are present in a conserved position on the transmembrane helix-6 (TM6) of ETHRs. In TM6, PxETHR-A and PxETHR-B differ at positions C/S6.43 and Y/F6.51, both being part of the HBN.Communicated by Ramaswamy H. Sarma.

Systematized analysis of secondary structure dependence of key structural features of residues in soluble and membrane-bound proteins

Knowledge of three-dimensional protein structure is integral to most modern drug discovery efforts. Recent advancements have highlighted new techniques for 3D protein structure determination and, where structural data cannot be collected experimentally, prediction of protein structure. We have undertaken a major effort to use existing protein structures to collect, characterize, and catalogue the inter-atomic interactions that define and compose 3D structure by mapping hydropathic interaction environments as maps in 3D space. This work has been performed on a residue-by-residue basis, where we have seen evidence for relationships between environment character, residue solvent-accessible surface areas and their secondary structures. In this graphical review, we apply principles from our earlier studies and expand the scope to all common amino acid residue types in both soluble and membrane proteins. Key to this analysis is parsing the Ramachandran plot to an 8-by-8 chessboard to define secondary structure bins. Our analysis yielded a number of quantitative discoveries: 1) increased fraction of hydrophobic residues (alanine, isoleucine, leucine, phenylalanine and valine) in membrane proteins compared to their fractions in soluble proteins; 2) less burial coupled with significant increases in favorable hydrophobic interactions for hydrophobic residues in membrane proteins compared to soluble proteins; and 3) higher burial and more favorable polar interactions for polar residues now preferring the interior of membrane proteins. These observations and the supporting data should provide benchmarks for current studies of protein residues in different environments and may be able to guide future protein structure prediction efforts.

3D Interaction Homology: Computational Titration of Aspartic Acid, Glutamic Acid and Histidine Can Create pH-Tunable Hydropathic Environment Maps

Aspartic acid, glutamic acid and histidine are ionizable residues occupying various protein environments and perform many different functions in structures. Their roles are tied to their acid/base equilibria, solvent exposure, and backbone conformations. We propose that the number of unique environments for ASP, GLU and HIS is quite limited. We generated maps of these residue's environments using a hydropathic scoring function to record the type and magnitude of interactions for each residue in a 2703-protein structural dataset. These maps are backbone-dependent and suggest the existence of new structural motifs for each residue type. Additionally, we developed an algorithm for tuning these maps to any pH, a potentially useful element for protein design and structure building. Here, we elucidate the complex interplay between secondary structure, relative solvent accessibility, and residue ionization states: the degree of protonation for ionizable residues increases with solvent accessibility, which in turn is notably dependent on backbone structure.