Hydrophobicity of transmembrane proteins: spatially profiling the distribution

Main Factors Shaping Amino Acid Usage Across Evolution

The standard genetic code determines that in most species, including viruses, there are 20 amino acids that are coded by 61 codons, while the other three codons are stop triplets. Considering the whole proteome each species features its own amino acid frequencies, given the slow rate of change, closely related species display similar GC content and amino acids usage. In contrast, distantly related species display different amino acid frequencies. Furthermore, within certain multicellular species, as mammals, intragenomic differences in the usage of amino acids are evident. In this communication, we shall summarize some of the most prominent and well-established factors that determine the differences found in the amino acid usage, both across evolution and intragenomically.

Molecular Characterisation of a Novel and Highly Divergent Passerine Adenovirus 1

Wild birds harbour a large number of adenoviruses that remain uncharacterised with respect to their genomic organisation, diversity, and evolution within complex ecosystems. Here, we present the first complete genome sequence of an atadenovirus from a passerine bird that is tentatively named Passerine adenovirus 1 (PaAdV-1). The PaAdV-1 genome is 39,664 bp in length, which was the longest atadenovirus to be sequenced, to the best of our knowledge, and contained 42 putative genes. Its genome organisation was characteristic of the members of genus Atadenovirus; however, the novel PaAdV-1 genome was highly divergent and showed the highest sequence similarity with psittacine adenovirus-3 (55.58%). Importantly, PaAdV-1 complete genome was deemed to contain 17 predicted novel genes that were not present in any other adenoviruses sequenced to date, with several of these predicted novel genes encoding proteins that harbour transmembrane helices. Subsequent analysis of the novel PaAdV-1 genome positioned phylogenetically to a distinct sub-clade with all others sequenced atadenoviruses and did not show any obvious close evolutionary relationship. This study concluded that the PaAdV-1 complete genome described here is not closely related to any other adenovirus isolated from avian or other natural host species and that it should be considered a separate species.

Cumulative hydropathic topology of a voltage-gated sodium channel at atomic resolution

Voltage-gated sodium channels (NavChs) are biological pores that control the flow of sodium ions through the cell membrane. In humans, mutations in genes encoding NavChs can disrupt physiological cellular activity thus leading to a wide spectrum of diseases. Here, we present a topological connection between the functional architecture of a NavAb bacterial channel and accumulation of atomic hydropathicity around its pore. This connection is established via a scaling analysis methodology that elucidates how intrachannel hydropathic density variations translate into hydropathic dipole field configurations along the pore. Our findings suggest the existence of a nonrandom cumulative hydropathic topology that is organized parallel to the membrane surface so that pore's stability, as well as, gating behavior are guaranteed. Given the biophysical significance of the hydropathic effect, our study seeks to provide a computational framework for studying cumulative hydropathic topological properties of NavChs and pore-forming proteins in general.

Ensembles of Hydrophobicity Scales as Potent Classifiers for Chimeric Virus-Like Particle Solubility - An Amino Acid Sequence-Based Machine Learning Approach

Virus-like particles (VLPs) are protein-based nanoscale structures that show high potential as immunotherapeutics or cargo delivery vehicles. Chimeric VLPs are decorated with foreign peptides resulting in structures that confer immune responses against the displayed epitope. However, insertion of foreign sequences often results in insoluble proteins, calling for methods capable of assessing a VLP candidate's solubility in silico. The prediction of VLP solubility requires a model that can identify critical hydrophobicity-related parameters, distinguishing between VLP-forming aggregation and aggregation leading to insoluble virus protein clusters. Therefore, we developed and implemented a soft ensemble vote classifier (sEVC) framework based on chimeric hepatitis B core antigen (HBcAg) amino acid sequences and 91 publicly available hydrophobicity scales. Based on each hydrophobicity scale, an individual decision tree was induced as classifier in the sEVC. An embedded feature selection algorithm and stratified sampling proved beneficial for model construction. With a learning experiment, model performance in the space of model training set size and number of included classifiers in the sEVC was explored. Additionally, seven models were created from training data of 24-384 chimeric HBcAg constructs, which were validated by 100-fold Monte Carlo cross-validation. The models predicted external test sets of 184-544 chimeric HBcAg constructs. Best models showed a Matthew's correlation coefficient of >0.6 on the validation and the external test set. Feature selection was evaluated for classifiers with best and worst performance in the chimeric HBcAg VLP solubility scenario. Analysis of the associated hydrophobicity scales allowed for retrieval of biological information related to the mechanistic backgrounds of VLP solubility, suggesting a special role of arginine for VLP assembly and solubility. In the future, the developed sEVC could further be applied to hydrophobicity-related problems in other domains, such as monoclonal antibodies.

Modeling Membrane Curvature Generation due to Membrane⁻Protein Interactions

To alter and adjust the shape of the plasma membrane, cells harness various mechanisms of curvature generation. Many of these curvature generation mechanisms rely on the interactions between peripheral membrane proteins, integral membrane proteins, and lipids in the bilayer membrane. Mathematical and computational modeling of membrane curvature generation has provided great insights into the physics underlying these processes. However, one of the challenges in modeling these processes is identifying the suitable constitutive relationships that describe the membrane free energy including protein distribution and curvature generation capability. Here, we review some of the commonly used continuum elastic membrane models that have been developed for this purpose and discuss their applications. Finally, we address some fundamental challenges that future theoretical methods need to overcome to push the boundaries of current model applications.

Hydrophobic Effect from Conjugated Chemicals or Drugs on In Vivo Biodistribution of RNA Nanoparticles

Liver or other organ accumulation of drugs is one of the major problems that leads to toxicity and side effects in therapy using chemicals or other macromolecules. It has been shown that specially designed RNA nanoparticles can specifically target cancer cells, silence oncogenic genes, and stop cancer growth with little or no accumulation in the liver or other vital organs. It is well known that physical properties of nanoparticles such as size, shape, and surface chemistry affect biodistribution and pharmacokinetic profiles in vivo. This study examined how the hydrophobicity of chemicals conjugated to RNA nanoparticles affect in vivo biodistribution. Weaker organ accumulation was observed for hydrophobic chemicals after they were conjugated to RNA nanoparticles, revealing RNA's ability to solubilize hydrophobic chemicals. It was found that different chemicals conjugated to RNA nanoparticles resulted in the alteration of RNA hydrophobicity. Stronger hydrophobicity induced by chemical conjugates resulted in higher accumulation of RNA nanoparticles in vital organs in mice. This study provides new insights for handling drug insolubility, therapeutic toxicity, and organ clearance in drug development.

Interaction and conformational dynamics of membrane-spanning protein helices

Within 1 or 2 decades, the reputation of membrane-spanning alpha-helices has changed dramatically. Once mostly regarded as dull membrane anchors, transmembrane domains are now recognized as major instigators of protein-protein interaction. These interactions may be of exquisite specificity in mediating assembly of stable membrane protein complexes from cognate subunits. Further, they can be reversible and regulatable by external factors to allow for dynamic changes of protein conformation in biological function. Finally, these helices are increasingly regarded as dynamic domains. These domains can move relative to each other in different functional protein conformations. In addition, small-scale backbone fluctuations may affect their function and their impact on surrounding lipid shells. Elucidating the ways by which these intricate structural features are encoded by the amino acid sequences will be a fascinating subject of research for years to come.

Localization of ligand binding site in proteins identified in silico

Knowledge-based models for protein folding assume that the early-stage structural form of a polypeptide is determined by the backbone conformation, followed by hydrophobic collapse. Side chain-side chain interactions, mostly of hydrophobic character, lead to the formation of the hydrophobic core, which seems to stabilize the structure of the protein in its natural environment. The fuzzy-oil-drop model is employed to represent the idealized hydrophobicity distribution in the protein molecule. Comparing it with the one empirically observed in the protein molecule reveals that they are not in agreement. It is shown in this study that the irregularity of hydrophobic distributions is aim-oriented. The character and strength of these irregularities in the organization of the hydrophobic core point to the specificity of a particular protein's structure/function. When the location of these irregularities is determined versus the idealized fuzzy-oil-drop, function-related areas in the protein molecule can be identified. The presented model can also be used to identify ways in which protein-protein complexes can possibly be created. Active sites can be predicted for any protein structure according to the presented model with the free prediction server at http://www.bioinformatics.cm-uj.krakow.pl/activesite. The implication based on the model presented in this work suggests the necessity of active presence of ligand during the protein folding process simulation.

Hydrophobic collapse in late-stage folding (in silico) of bovine pancreatic trypsin inhibitor

Hydrophobic collapse is commonly considered as a process of significance for protein folding. Many models have been applied for description of this event. A model introducing an external force field mimicking the hydrophobic environment and simultaneously the driving force for the folding process is presented in this paper. Bovine pancreatic trypsin inhibitor (BPTI) was taken as a test protein. An early-stage folding (in silico) model presented elsewhere was used to create the starting structure for hydrophobic collapse. The resulting structure was energy-refined using molecular dynamics simulation in an explicit solvent. The similarity versus the crystal structure of BPTI is estimated using visual analysis, residue-residue contacts, phi, psi angle distributions, RMSD, accessible solvent area, radii of gyration and hydrodynamic radii. A program allowing creation of early-stage folding structural forms to be created for any protein is available from http://bioinformatics.cm-uj.krakow.pl/earlystage. The program for late-stage folding simulation is available on request.

Energetics and stability of transmembrane helix packing: A replica-exchange simulation with a knowledge-based membrane potential

The energetics and stability of the packing of transmembrane helices were investigated by Monte Carlo simulations with the replica-exchange method. The helices were modeled with a united atom representation, and the CHARMM19 force field was employed. Based on known experimental structures of membrane proteins, an implicit knowledge-based potential was developed to describe the helix-membrane interactions at the residue level, whose validity was tested through prediction of the orientations when single helices were inserted into a membrane. Two systems were studied in this article, namely the glycophorin A dimer, and helices A and B of Bacteriorhodopsin. For the glycophorin A dimer, the most stable structure (0.5 A away from the experimental structure) is mainly stabilized by the favorable helix-helix interactions, and has the most population regardless of the helix-membrane interaction. However, for helices A and B of Bacteriorhodopsin, it was found that the packing determined by helix-helix interactions is nonspecific, and a native-like structure (0.2 A from the experimental one) can be identified from several structural analogs as the most stable one only after applying the membrane potential. Our results suggest that the contribution from the helix-membrane interaction could be critical in the correct packing of transmembrane helices in the membrane.

The Mighty Arginine, the Stable Quaternary Amines, the Powerful Aromatics, and the Aggressive Phosphate: Their Role in the Noncovalent Minuet

In the age of proteomics, the role of certain amino acid residues and some post-translational modifications in noncovalent complex formation are gaining in importance, as the understanding of interactions between biological molecules, is at the heart of the structure function relationship puzzle. In this work, mass spectrometry is used to highlight ammonium- or guanidinium-aromatic interactions through Cation-pi bonds and ammonium- or guanidinium-phosphate interactions through salt bridge formation. Such interactions are crucial factors in certain ligand-receptor interactions and receptor-receptor interactions. In addition, the ability of phosphorylated residues and phosphorylated lipids to form noncovalent complexes with guanidinium and quaternary ammonium (mostly through Coulombic interactions) is demonstrated, and could explain the stability of certain membrane embedded protein, or a possible role for phosphorylation in protein-protein interactions. Dougherty's work demonstrates cation-pi interactions in intra-protein interactions and folding, the present work explores inter-peptide interactions, i.e., the formation of noncovalent complexes between peptides' epitopes containing adjacent aromatic residues and ones containing adjacent Arg as a model to better understand the role of cation-pi complexes in protein-protein interaction. Complexes of peptides containing aromatic residues with quaternary amines as well as the interaction of aromatic compounds, with the guanidinium group of Arg are also investigated. Considering that an inordinate number of therapeutic compounds contain aromatic rings and quaternary amines, the above-described interactions could possibly be of great importance in better understanding their mechanism of action.