Molecular Modeling of the HR2 and Transmembrane Domains of the SARS-CoV-2 S Protein in the Prefusion State.

SARS-CoV-2, the causative agent of COVID-19, remains the focus of research worldwide. SARS-CoV-2 entry into the cell starts with its S protein binding to the angiotensin-converting enzyme-2 (ACE2) expressed on the cell surface. The knowledge of the S protein's spatial structure is indispensable for understanding the molecular principles of its work. The S protein structure has been almost fully described using experimental approaches with the only exception for the protein's endodomain, the transmembrane domain, and the ectodomain parts adjacent to the latter. The paper reports molecular modelling of the S protein fragment corresponding to its coiled coil HR2 domain and fully palmitoylated transmembrane domain. Model stability in lipid bilayer was confirmed by all-atom and coarse-grained molecular dynamics simulations. It has been demonstrated that palmitoylation leads to a significant decrease in transmembrane domain mobility and local bilayer thickening, which may be relevant for protein trimerization.

Complexes of Peptide Blockers with Kv1.6 Pore Domain: Molecular Modeling and Studies with KcsA-Kv1.6 Channel.

Potassium voltage-gated Kv1.6 channel, which is distributed primarily in neurons of central and peripheral nervous systems, is of significant physiological importance. To date, several high-affinity Kv1.6-channel blockers are known, but the lack of selective ones among them hampers the studies of tissue localization and functioning of Kv1.6 channels. Here we present an approach to advanced understanding of interactions of peptide toxin blockers with a Kv1.6 pore. It combines molecular modeling studies and an application of a new bioengineering system based on a KcsA-Kv1.6 hybrid channel for the quantitative fluorescent analysis of blocker-channel interactions. Using this system we demonstrate that peptide toxins agitoxin 2, kaliotoxin1 and OSK1 have similar high affinity to the extracellular vestibule of the K+-conducting pore of Kv1.6, hetlaxin is a low-affinity ligand, whereas margatoxin and scyllatoxin do not bind to Kv1.6 pore. Binding of toxins to Kv1.6 pore has considerable inverse dependence on the ionic strength. Model structures of KcsA-Kv1.6 and Kv1.6 complexes with agitoxin 2, kaliotoxin 1 and OSK1 were obtained using homology modeling and molecular dynamics simulation. Interaction interfaces, which are formed by 15-19 toxin residues and 10 channel residues, are described and compared. Specific sites of Kv1.6 pore recognition are identified for targeting of peptide blockers. Analysis of interactions between agitoxin 2 derivatives with point mutations (S7K, S11G, L19S, R31G) and KcsA-Kv1.6 confirms reliability of the calculated complex structure.

Modeling of the Binding of Peptide Blockers to Voltage-Gated Potassium Channels: Approaches and Evidence.

Modeling of the structure of voltage-gated potassium (KV) channels bound to peptide blockers aims to identify the key amino acid residues dictating affinity and provide insights into the toxin-channel interface. Computational approaches open up possibilities for in silico rational design of selective blockers, new molecular tools to study the cellular distribution and functional roles of potassium channels. It is anticipated that optimized blockers will advance the development of drugs that reduce over activation of potassium channels and attenuate the associated malfunction. Starting with an overview of the recent advances in computational simulation strategies to predict the bound state orientations of peptide pore blockers relative to KV-channels, we go on to review algorithms for the analysis of intermolecular interactions, and then take a look at the results of their application.

Construction of TNF-binding proteins by grafting hypervariable regions of F10 antibody on human fibronectin domain scaffold.

Tumor necrosis factor (TNF) plays a key role in the pathogenesis of various diseases. To study the possibility of constructing TNF-binding proteins by grafting hypervariable regions of immunoglobulins (CDR), we have replaced amino acid sequences of loops from the tenth type III domain of human fibronectin ((10)Fn3) by amino acid sequences of CDR from the light and heavy chains of the anti-TNF antibody F10. The assessment of TNF-binding properties of the resulting proteins by ELISA has revealed the highest activity of Hd3 containing sequences CDR-H1 and CDR-H2 of the antibody F10 and of Hd2 containing sequences CDR-H1 and CDR-H3. The proteins constructed by us on the fibronectin domain scaffold specifically bound TNF during Western blotting and also weakened its cytotoxic effect on L929 line cells. The highest neutralizing activity was demonstrated by the proteins Hd2 and Hd3, which induced, respectively, 10- and 50-fold increase in the EC(50) of TNF.

Analysis of hydrophobic interactions of antagonists with the beta2-adrenergic receptor.

The adrenergic receptors mediate a wide variety of physiological responses, including vasodilatation and vasoconstriction, heart rate modulation, and others. Beta-adrenergic antagonists ('beta-blockers') thus constitute a widely used class of drugs in cardiovascular medicine as well as in management of anxiety, migraine, and glaucoma. The importance of the hydrophobic effect has been evidenced for a wide range of beta-blocker properties. To better understand the role of the hydrophobic effect in recognition of beta-blockers by their receptor, we carried out a molecular docking study combined with an original approach to estimate receptor-ligand hydrophobic interactions. The proposed method is based on automatic detection of molecular fragments in ligands and the analysis of their interactions with receptors separately. A series of beta-blockers, based on phenylethanolamines and phenoxypropanolamines, were docked to the beta2-adrenoceptor binding site in the crystal structure. Hydrophobic complementarity between the ligand and the receptor was calculated using the PLATINUM web-server (http://model.nmr.ru/platinum). Based on the analysis of the hydrophobic match for molecular fragments of beta-blockers, we have developed a new scoring function which efficiently predicts dissociation constant (pKd) with strong correlations (r(2) approximately 0.8) with experimental data.

A novel method for packing quality assessment of transmembrane alpha-helical domains in proteins.

Here we present a novel method for assessment of packing quality for transmembrane (TM) domains of alpha-helical membrane proteins (MPs), based on analysis of available high-resolution experimental structures of MPs. The presented concept of protein-membrane environment classes permits quantitative description of packing characteristics in terms of membrane accessibility and polarity of the nearest protein groups. We demonstrate that the method allows identification of native-like conformations among the large set of theoretical MP models. The developed "membrane scoring function" will be of use for optimization of TM domain packing in theoretical models of MPs, first of all G-protein coupled receptors.