RESUMO
The so-called "hydrophobic gating" is widely discussed as a putative mechanism to control water and ion conduction via ion channels. This effect can occur in narrow areas of the channels pore lined by non-polar residues. In the closed state of the channel, such regions may spontaneously transit to a dehydrated state to block water and ions transport without full pore occlusion. In the open state, the hydrophobic gate is wide enough to provide sustainable hydration and conduction. Apparently, the transport through the open hydrophobic gate may by facilitated by some polar residues that assist polar/charged substances to overcome the energy barrier created by nonpolar environment. In this work, we investigated the behavior of Na+ ions and their hydration shells in the open pore of the rat TRPV1 ion channel by molecular dynamics simulations. We show that polar protein groups coordinate water molecules in such a way as to restore the hydration shell of ions in the hydrophobic gate that ensures ion transport through the gate in a fully hydrated state.
Assuntos
Ativação do Canal Iônico , Canais Iônicos , Animais , Ratos , Canais Iônicos/química , Canais Iônicos/metabolismo , Simulação de Dinâmica Molecular , Íons/metabolismo , Água/química , Canais de Cátion TRPVRESUMO
Atomistic aspects of the structural organization, dynamics, and functioning of hydrated lipid bilayers-model cell membranes-are primarily governed by the fine balance of intermolecular interactions between all constituents of these systems. Besides the hydrophobic effect, which shapes the overall skeleton of lipid membranes, a very important contribution to their behavior is made by hydrogen bonds (H-bonds) between lipid head groups. The latter determine crucial phenomena in cell membranes, such as dynamic ultrananodomain organization, hydration, and fine-tuning of microscopic physicochemical properties that allow the membrane to adapt quickly when binding/insertion external agents (proteins, etc.). The characteristics of such H-bonds (strength, spatial localization, etc.) dramatically depend on the local polarity properties of the lipid-water environment. In this work, we calculated free energies of H-bonded complexes between typical donor (NH3+, NH, OH) and acceptor (CâO, OH, COO-, COOH) groups of lipids in vacuo and in a set of explicit solvents with dielectric constants (ε) from 1 to 78.3, which mimic membrane environment at different depths. This was done using Monte Carlo simulations and an assessment of the corresponding potential of mean force profiles. The strongest H-bonded complexes were observed in the nonpolar environment, and their strength increased sharply with decreasing ε below 17. When ε changed, the largest free energy gain (>10.8 kcal/mol) was observed for pairs of acceptors CâO and O(H) with donor NH3+. The complexation of the same acceptors with NH donor in this range of ε values was rather less sensitive to the environmental polarity, by â¼1.5 kcal/mol. Dielectric-dependent interactions of polar lipid groups with water were evaluated as well. The results explain the delicate balance that determines the unique pattern of H-bonds for a particular lipid bilayer. Understanding the factors that regulate the propensity for H-bonding in lipid bilayers provides a fundamental basis for the rational design of new membrane nano objects with predefined properties.
Assuntos
Membrana Celular/química , Cicloexanos/química , Impedância Elétrica , Ligação de Hidrogênio , Metanol/química , Modelos Moleculares , Conformação Molecular , Termodinâmica , Água/químicaRESUMO
The interaction between the secondary structure elements is the key process, determining the spatial structure and activity of a membrane protein. Transmembrane (TM) helix-helix interaction is known to be especially important for the function of so-called type I or bitopic membrane proteins. In the present work, we present the approach to study the helix-helix interaction in the TM domains of membrane proteins in various lipid environment using solution NMR spectroscopy and phospholipid bicelles. The technique is based on the ability of bicelles to form particles with the size, depending on the lipid/detergent ratio. To implement the approach, we report the experimental parameters of "ideal bicelle" models for four kinds of zwitterionic phospholipids, which can be also used in other structural studies. We show that size of bicelles and type of the rim-forming detergent do not affect substantially the spatial structure and stability of the model TM dimer. On the other hand, the effect of bilayer thickness on the free energy of the dimer is dramatic, while the structure of the protein is unchanged in various lipids with fatty chains having a length from 12 to 18 carbon atoms. The obtained data is analyzed using the computer simulations to find the physical origin of the observed effects.
Assuntos
Simulação por Computador , Proteínas de Membrana/metabolismo , Ressonância Magnética Nuclear Biomolecular/métodos , Mapeamento de Interação de Proteínas , Sequência de Aminoácidos , Detergentes/química , Dimerização , Micelas , Simulação de Dinâmica Molecular , Fosfolipídeos/química , Estabilidade Proteica , Estrutura Secundária de ProteínaRESUMO
Human-secreted Ly-6/uPAR-related protein-2 (SLURP-2) regulates the growth and differentiation of epithelial cells. Previously, the auto/paracrine activity of SLURP-2 was considered to be mediated via its interaction with the α3ß2 subtype of the nicotinic acetylcholine receptors (nAChRs). Here, we describe the structure and pharmacology of a recombinant analogue of SLURP-2. Nuclear magnetic resonance spectroscopy revealed a 'three-finger' fold of SLURP-2 with a conserved ß-structural core and three protruding loops. Affinity purification using cortical extracts revealed that SLURP-2 could interact with the α3, α4, α5, α7, ß2, and ß4 nAChR subunits, revealing its broader pharmacological profile. SLURP-2 inhibits acetylcholine-evoked currents at α4ß2 and α3ß2-nAChRs (IC50 ~0.17 and >3 µM, respectively) expressed in Xenopus oocytes. In contrast, at α7-nAChRs, SLURP-2 significantly enhances acetylcholine-evoked currents at concentrations <1 µM but induces inhibition at higher concentrations. SLURP-2 allosterically interacts with human M1 and M3 muscarinic acetylcholine receptors (mAChRs) that are overexpressed in CHO cells. SLURP-2 was found to promote the proliferation of human oral keratinocytes via interactions with α3ß2-nAChRs, while it inhibited cell growth via α7-nAChRs. SLURP-2/mAChRs interactions are also probably involved in the control of keratinocyte growth. Computer modeling revealed possible SLURP-2 binding to the 'classical' orthosteric agonist/antagonist binding sites at α7 and α3ß2-nAChRs.
Assuntos
Potenciais Evocados/fisiologia , Proteínas Ligadas por GPI/metabolismo , Queratinócitos/metabolismo , Receptores Muscarínicos/metabolismo , Receptores Nicotínicos/metabolismo , Receptor Nicotínico de Acetilcolina alfa7/metabolismo , Proteínas Adaptadoras de Transdução de Sinal , Adulto , Animais , Sítios de Ligação/fisiologia , Células CHO , Linhagem Celular , Proliferação de Células/fisiologia , Simulação por Computador , Cricetulus , Epilepsia do Lobo Temporal/patologia , Feminino , Humanos , Pessoa de Meia-Idade , Ressonância Magnética Nuclear Biomolecular , Oócitos/metabolismo , Células PC12 , Ligação Proteica/fisiologia , Ratos , XenopusAssuntos
Algoritmos , Conotoxinas/química , Desenho de Fármacos , Antagonistas Nicotínicos/química , Receptor Nicotínico de Acetilcolina alfa7/química , Acetilcolina/farmacologia , Sequência de Aminoácidos , Animais , Aplysia , Conotoxinas/genética , Relação Dose-Resposta a Droga , Humanos , Ligantes , Lymnaea , Potenciais da Membrana/efeitos dos fármacos , Potenciais da Membrana/fisiologia , Modelos Moleculares , Dados de Sequência Molecular , Agonistas Nicotínicos/farmacologia , Antagonistas Nicotínicos/farmacologia , Oócitos , Xenopus laevis , Receptor Nicotínico de Acetilcolina alfa7/agonistas , Receptor Nicotínico de Acetilcolina alfa7/antagonistas & inibidores , Receptor Nicotínico de Acetilcolina alfa7/metabolismoRESUMO
An efficient computational approach is developed to quantify the free energy of a spontaneous association of the α-helices of proteins in the membrane environment. The approach is based on the numerical decomposition of the free energy profiles of the transmembrane (TM) helices into components corresponding to protein-protein, protein-lipid, and protein-water interactions. The method was tested for the TM segments of human glycophorin A (GpA) and two mutant forms, Gly83Ala and Thr87Val. It was shown that lipids make a significant negative contribution to the free energy of dimerization, while amino acid residues forming the interface of the helix-helix contact may be unfavorable in terms of free energy. The detailed balance between different energy contributions is highly dependent on the amino acid sequence of the TM protein segment. The results show the dominant role of the environment in the interaction of membrane proteins that is changing our notion of the driving force behind the spontaneous association of TM α-helices. Adequate estimation of the contribution of the water-lipid environment thus becomes an extremely urgent task for a rational design of new molecules targeting bitopic membrane proteins, including receptor tyrosine kinases.
RESUMO
Polysialic acid (PSA) is a natural anionic polymer typically occurring on the outer surface of cell membranes. PSA is involved in cell signaling and intermolecular interactions with proteins and peptides. The antimicrobial potential of peptides is usually evaluated in model membranes consisting of lipid bilayers but devoid of either PSA or its analogs. The goal of this work was to investigate the possible effect of PSA on the structure of melittin (Mlt) and latarcins Ltc1K, Ltc2a, and the activity of these peptides with respect to model membranes. These peptides are linear cationic ones derived from the venom of bee (Mlt) and spider (both latarcins). The length of each of the peptides is 26 amino acid residues, and they all have antimicrobial activity. However, they differ with respect to conformational mobility, hydrophobic characteristics, and overall charge. In this work, using circular dichroism spectroscopy, we show that the peptides adopt an α-helical conformation upon interaction with either PSA or phospholipid liposomes formed of either zwitterionic or anionic phospholipids or their mixtures. The extent of helicity depends on the amino acid sequence and properties of the medium. Based on small angle X-ray scattering data and the analysis of the fluorescence spectrum of the Trp residue in Mlt, we conclude that the peptide forms an oligomeric complex consisting of α-helical Mlt and several PSA molecules. Both latarcins, unlike Mlt, the most hydrophobic of the peptides, interact weakly with zwitterionic liposomes. However, they bind anionic liposomes or those composed of anionic/zwitterionic lipid mixtures. Latarcin Ltc1K forms associates on liposomes composed of zwitterionic/anionic lipid mixture. The structure of the peptide associates is either disordered or of ß-sheet conformation. In all other cases the studied peptides adopt predominately α-helical conformation. In addition, we demonstrate that PSA inhibits membranolytic activity of Mlt and latarcin Ltc1K. These data suggest that the peptides, due to their high conformational lability, can vary structural and amphiphilic properties in the presence of PSA. As a result, various scenarios of the interaction of the peptides with membranes, whose surface is abundant with anionic polysaccharides, can take place. This can account for difficulties in understanding the structure-functional relationships in interactions of linear cationic peptides with biological membranes.
Assuntos
Bicamadas Lipídicas/metabolismo , Peptídeos/metabolismo , Fosfolipídeos/química , Ácidos Siálicos/metabolismo , Sequência de Aminoácidos , Peptídeos Catiônicos Antimicrobianos/química , Peptídeos Catiônicos Antimicrobianos/metabolismo , Cátions/química , Dicroísmo Circular , Interações Hidrofóbicas e Hidrofílicas , Bicamadas Lipídicas/química , Lipossomos/química , Lipossomos/metabolismo , Meliteno/química , Meliteno/metabolismo , Dados de Sequência Molecular , Peptídeos/química , Ligação Proteica , Estrutura Secundária de Proteína , Ácidos Siálicos/química , Venenos de Aranha/química , Venenos de Aranha/metabolismoRESUMO
Natural polycationic membrane-active peptides typically lack disulfide bonds and exhibit fusion, cell-penetrating, antimicrobial activities. They are mostly unordered in solution, but adopt a helical structure, when bound to phospholipid membranes. Structurally different are cardiotoxins (or cytotoxins, CTs) from cobra venom. They are fully ß- structured molecules, characterized by the three-finger fold (TFF). Affinity of CTs to lipid bilayer was shown to depend on amino acid sequence in the tips of the three loops. In the present review, CT-membrane interactions are analyzed through the prism of data on binding of the toxins to phospholipid liposomes and detergent micelles, as well as their structural and computational studies in membrane mimicking environments. We assess different hydrophobicity scales to compare membrane partitioning of various CTs and their membrane effects. A comparison of hydrophobic/hydrophilic properties of CTs and linear polycationic peptides provides a key to their biological activity and creates a fundamental basis for rational design of new membrane-interacting compounds, including new promising drugs. For instance, from the viewpoint of the data obtained on model lipid membranes, cytotoxic activity of CTs against cancer cells is discussed.
Assuntos
Antineoplásicos/química , Antineoplásicos/farmacologia , Membrana Celular/efeitos dos fármacos , Proteínas Cardiotóxicas de Elapídeos/química , Proteínas Cardiotóxicas de Elapídeos/farmacologia , Sequência de Aminoácidos , Animais , Membrana Celular/metabolismo , Humanos , Lipossomos/metabolismo , Micelas , Modelos Moleculares , Dados de Sequência Molecular , Neoplasias/tratamento farmacológico , Neoplasias/metabolismo , Fosfolipídeos/metabolismoRESUMO
Specific interactions between transmembrane α-helices, to a large extent, determine the biological function of integral membrane proteins upon normal development and in pathological states of an organism. Various membrane-like media, partially those mimicking the conditions of multicomponent biological membranes, are used to study the structural and thermodynamic features that define the character of oligomerization of transmembrane helical segments. The choice of the composition of the membrane-mimicking medium is conducted in an effort to obtain a biologically relevant conformation of the protein complex and a sample that would be stable enough to allow to perform a series of long-term experiments with its use. In the present work, heteronuclear NMR spectroscopy and molecular dynamics simulations were used to demonstrate that the two most widely used media (detergent DPC micelles and lipid DMPC/DHPC bicelles) enable to perform structural studies of the specific interactions between transmembrane α-helices by the example of dimerizing the transmembrane domain of the bitopic protein glycophorin A. However, a number of peculiarities place lipid bicelles closer to natural lipid bilayers in terms of their physical properties.
RESUMO
Knowledge of 3D-structure of protein-ligand complex is a major prerequisite for understanding the functioning mechanism of cellular proteins and membrane receptors. This is also of a great help in rational drug design projects. In the present paper we briefly review the molecular docking approaches used to predict possible orientation of a ligand in the protein binding site. The recent trends to improve the accuracy and efficiency of docking algorithms are demonstrated with the results obtained in Laboratory of Biomolecular Modeling. Particular attention is paid to protein-ligand hydrophobic and stacking interactions responsible for molecular recognition of ligand fragments. Such type of interactions are not always adequately represented in scoring criteria of docking applications that leads to mismatch in 3D-structure complexes predictions. That is why further inquiry of methods to account for these interactions is now the area of active research.
Assuntos
Algoritmos , Modelos Moleculares , Nucleotídeos/química , Peptídeos/química , Proteínas/química , Software , Sítios de Ligação , Nucleotídeos/metabolismo , Peptídeos/metabolismo , Ligação Proteica , Proteínas/metabolismoRESUMO
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.
Assuntos
Antagonistas de Receptores Adrenérgicos beta 2 , Antagonistas Adrenérgicos beta/metabolismo , Interações Hidrofóbicas e Hidrofílicas , Modelos Moleculares , Receptores Adrenérgicos beta 2/metabolismo , Estrutura Molecular , Ligação ProteicaRESUMO
Membrane-spanning segments of numerous proteins (e.g. receptor tyrosine kinases) represent a novel class of pharmacologically important targets, whose activity can be modulated by specially designed artificial peptides, the so-called interceptors. Rational construction of such peptides requires understanding of the main factors driving peptide-peptide association in lipid membranes. Here we present a new method for rapid prediction of the spatial structure of transmembrane (TM) helix-helix complexes. It is based on computer simulations in membrane-like media and subsequent refinement/validation of the results using experimental studies of TM helix dimerization in a bacterial membrane by means of the ToxR system. The approach was applied to TM fragments of the ephrin receptor A1 (EphA1). A set of spatial structures of the dimer was proposed based on Monte Carlo simulations in an implicit membrane followed by molecular dynamics relaxation in an explicit lipid bilayer. The resulting models were employed for rational design of wild-type and mutant genetic constructions for ToxR assays. The computational and the experimental data are self-consistent and provide an unambiguous spatial model of the TM dimer of EphA1. The results of this work can be further used to develop new biologically active 'peptide interceptors' specifically targeting membrane domains of proteins.
Assuntos
Proteínas de Bactérias/química , Simulação por Computador , Proteínas de Ligação a DNA/química , Modelos Moleculares , Receptor EphA1/química , Fatores de Transcrição/química , Proteínas de Bactérias/genética , Proteínas de Ligação a DNA/genética , Simulação de Dinâmica Molecular , Método de Monte Carlo , Multimerização Proteica , Estrutura Secundária de Proteína , Receptor EphA1/metabolismo , Fatores de Transcrição/genética , beta-Galactosidase/química , beta-Galactosidase/metabolismoRESUMO
The cytotoxic properties of cytotoxins (CTs) from snake venom are mediated by their interaction with the cell membrane. The hydrophobic pattern containing the tips of loops I-III and flanked by polar residues is known to be a membrane-binding motif of CTs. However, this is not enough to explain the difference in activity among various CTs which are similar in sequence and in 3D structure. The mechanism of further CT-membrane interaction leading to pore formation and cell death still remains unknown. Published experimental data on the specific interaction between CT and low molecular weight anionic components (sulphatide) of the bilayer point to the existence of corresponding ligand binding sites on the surface of toxin molecules. In this work we study the membrane-lytic properties of CT I, CT II (Naja oxiana), and Ct 4 (Naja kaouthia), which belong to different structural and functional types (P- and S-type) of CTs, by measuring the intensity of a fluorescent dye, calcein released from liposomes containing a phosphatidylserine (PS) lipid as an anionic component. Using molecular docking simulations, we find and characterize three sites in CT molecules that can potentially bind the PS polar head. Based on the data obtained, we suggest a hypothesis that CTs can specifically interact with one or more of the anionic lipids (in particular, with PS) contained in the membrane, thus facilitating the interaction between CTs and the lipid bilayer of a cell membrane.
RESUMO
Knowledge of the spatial structure of proteins is a prerequisite for both awareness of their functional mechanisms and the framework for rational drug discovery and design. Meanwhile, direct structural determination is often hampered or impractical due to the complexity, expensiveness, and limited capabilities of experimental techniques. These issues are especially pronounced for integral membrane proteins. On numerous occasions, the theoretical prediction of protein structures may facilitate the process by exploiting physical or empirical principles. This paper surveys modern techniques for the prediction of the spatial structure of proteins using computer algorithms, and the main emphasis is placed on the most "complex" targets - membrane proteins (MPs). The first part of the review describes de novo methods based on empirical physical principles; in the second part, a comparative modeling philosophy, which accounts for the structure of related proteins, is described. Special focus is made regarding pharmacologically relevant classes of G-coupled receptors, receptor tyrosine ki-nases, and other MPs. Algorithms for the assessment of the models quality and potential fields of application of computer models are discussed.
Assuntos
Algoritmos , Proteínas de Membrana/química , Modelos Moleculares , Animais , Humanos , Conformação Proteica , Relação Estrutura-AtividadeRESUMO
Molecular docking is a powerful computational method that has been widely used in many biomolecular studies to predict geometry of a protein-ligand complex. However, while its conformational search algorithms are usually able to generate correct conformation of a ligand in the binding site, the scoring methods often fail to discriminate it among many false variants. We propose to treat this problem by applying more precise ligand-specific scoring filters to re-rank docking solutions. In this way specific features of interactions between protein and different types of compounds can be implicitly taken into account. New scoring functions were constructed including hydrogen bonds, hydrophobic and hydrophilic complementarity terms. These scoring functions also discriminate ligands by the size of the molecule, the total hydrophobicity, and the number of peptide bonds for peptide ligands. Weighting coefficients of the scoring functions were adjusted using a training set of 60 protein-ligand complexes. The proposed method was then tested on the results of docking obtained for an additional 70 complexes. In both cases the success rate was 5-8% better compared to the standard functions implemented in popular docking software.
Assuntos
Modelos Moleculares , Ligação de Hidrogênio , Interações Hidrofóbicas e Hidrofílicas , Ligantes , Ligação ProteicaRESUMO
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.
Assuntos
Proteínas/química , Análise de Sequência de Proteína/métodos , Algoritmos , Modelos Teóricos , Estrutura Secundária de Proteína , Estrutura Terciária de Proteína , Reprodutibilidade dos Testes , Rodopsina/químicaRESUMO
The determination of the intermediate state structures of the bacteriorhodopsin photocycle has lead to an unprecedented level of understanding of the catalytic process exerted by a membrane protein. However, the crystallographic structures of the intermediate states are only relevant if the working cycle is not impaired by the crystal lattice. Therefore, we applied visible and Fourier transform infrared spectroscopy (FTIR) microspectroscopy with microsecond time resolution to compare the photoreaction of a single bacteriorhodopsin crystal to that of bacteriorhodopsin residing in the native purple membrane. The analysis of the FTIR difference spectra of the resolved intermediate states reveals great similarity in structural changes taking place in the crystal and in PM. However, the kinetics of the photocycle are significantly altered in the three-dimensional crystal as compared to PM. Strikingly, the L state decay is accelerated in the crystal, whereas the M decay is delayed. The physical origin of this deviation and the implications for trapping of intermediate states are discussed. As a methodological advance, time-resolved step-scan FTIR spectroscopy on a single protein crystal is demonstrated for the first time which may be used in the future to gauge the functionality of other crystallized proteins with the molecular resolution of vibrational spectroscopy.
Assuntos
Bacteriorodopsinas/química , Bacteriorodopsinas/efeitos da radiação , Cristalografia/métodos , Espectroscopia de Infravermelho com Transformada de Fourier/métodos , Bacteriorodopsinas/ultraestrutura , Relação Dose-Resposta à Radiação , Cinética , Luz , Fotobiologia/métodos , Fotoquímica/métodos , Doses de Radiação , Fatores de TempoRESUMO
Sensory rhodopsin II (SRII) from Halobacterium salinarum is heterologously expressed in Escherichia coli with a yield of 3-4 mg of purified SRII per liter cell culture. UV/Vis absorption spectroscopy display bands characteristic for native SRII. The resonance Raman spectrum provides evidence for a strongly hydrogen-bonded Schiff base like in mammalian rhodopsin but unlike to the homologous pSRII from Natronobacterium pharaonis. Laser flash spectroscopy indicates that SRII in detergent as well as after reconstitution into polar lipids shows its typical photochemical properties with prolonged photocycle kinetics. The first functional heterologous expression of SRII from H. salinarum provides the basis for studies with its cognate transducer HtrII to investigate the molecular processes involved in phototransduction as well as in chemotransduction.
Assuntos
Escherichia coli/genética , Halobacterium salinarum/genética , Rodopsinas Sensoriais/genética , Rodopsinas Sensoriais/metabolismo , Eletroforese , Cinética , Proteínas Recombinantes/genética , Proteínas Recombinantes/isolamento & purificação , Proteínas Recombinantes/metabolismo , Rodopsinas Sensoriais/isolamento & purificação , Análise Espectral RamanRESUMO
For the development of a method for the prediction of single point mutations substantially affecting protein thermostability, we studied the effect of the E85R and R82E mutations on the thermostability of thioredoxins from Escherichia coli (Trx) and Bacillus acidocaldarius (BacTrx), respectively. The basic method of investigation was the molecular dynamics simulation of 3D protein models in a particular solvent at different temperatures (300 and 373 K). Some thermolabile regions in Trx, BacTrx, and their mutants were revealed by analyzing the temperature effect on the molecular dynamics of the protein molecule. The effect of single point mutations on the temperature changes of the protein conformation mobility in several thermolabile regions was found. The results of the calculations are in accord with the experimental data indicating that the mutation E85R increases Trx thermostability, whereas the mutation R82E decreases BacTrx thermostability. The thermostability of these proteins was revealed to depend on ionic interactions between the thermolabile regions. The single point mutations change the parameters of these interactions and make them more favorable in the E85R-Trx mutant and less favorable in the R82E-BacTrx mutant. The English version of the paper: Russian Journal of Bioorganic Chemistry, 2004, vol. 30, no. 5; see also http: // www.maik.ru.