Exploring the Levinthal limit in protein folding.

According to the thermodynamic hypothesis, the native state of proteins is uniquely defined by their amino acid sequence. On the other hand, according to Levinthal, the native state is just a local minimum of the free energy and a given amino acid sequence, in the same thermodynamic conditions, can assume many, very different structures that are as thermodynamically stable as the native state. This is the Levinthal limit explored in this work. Using computer simulations, we compare the interactions that stabilize the native state of four different proteins with those that stabilize three non-native states of each protein and find that the nature of the interactions is very similar for all such 16 conformers. Furthermore, an enhancement of the degree of fluctuation of the non-native conformers can be explained by an insufficient relaxation to their local free energy minimum. These results favor Levinthal's hypothesis that protein folding is a kinetic non-equilibrium process.

Extended secondary structures in proteins.

Super secondary structures of proteins have been systematically searched and classified, but not enough attention has been devoted to such large edifices beyond the basic identification of secondary structures. The objective of the present study is to show that the association of secondary structures that share some of their backbone residues is a commonplace in globular proteins, and that such deeper fusion of secondary structures, namely extended secondary structures (ESSs), helps stabilize the original secondary structures and the resulting tertiary structures. For statistical purposes, a set of 163 proteins from the protein databank was randomly selected and a few specific cases are structurally analyzed and characterized in more detail. The results point that about 30%of the residues from each protein, on average, participate in ESS. Alternatively, for the specific cases considered,our results were based on the secondary structures produced after extensive Molecular Dynamics simulation of a protein­aqueous solvent system. Based on the very small width of the time distribution of the root mean squared deviations, between the ESS taken along the simulation and the ESS from the mean structure of the protein, for each ESS, we conclude that the ESSs significantly increase the conformational stability by forming very stable aggregates.The ubiquity and specificity of the ESS suggest that the role they play in the structure of proteins, including the domains formation, deserves to be thoroughly investigated.

Extensive structural change of the envelope protein of dengue virus induced by a tuned ionic strength: conformational and energetic analyses.

The Dengue has become a global public health threat, with over 100 million infections annually; to date there is no specific vaccine or any antiviral drug. The structures of the envelope (E) proteins of the four known serotype of the dengue virus (DENV) are already known, but there are insufficient molecular details of their structural behavior in solution in the distinct environmental conditions in which the DENVs are submitted, from the digestive tract of the mosquito up to its replication inside the host cell. Such detailed knowledge becomes important because of the multifunctional character of the E protein: it mediates the early events in cell entry, via receptor endocytosis and, as a class II protein, participates determinately in the process of membrane fusion. The proposed infection mechanism asserts that once in the endosome, at low pH, the E homodimers dissociate and insert into the endosomal lipid membrane, after an extensive conformational change, mainly on the relative arrangement of its three domains. In this work we employ all-atom explicit solvent Molecular Dynamics simulations to specify the thermodynamic conditions in that the E proteins are induced to experience extensive structural changes, such as during the process of reducing pH. We study the structural behavior of the E protein monomer at acid pH solution of distinct ionic strength. Extensive simulations are carried out with all the histidine residues in its full protonated form at four distinct ionic strengths. The results are analyzed in detail from structural and energetic perspectives, and the virtual protein movements are described by means of the principal component analyses. As the main result, we found that at acid pH and physiological ionic strength, the E protein suffers a major structural change; for lower or higher ionic strengths, the crystal structure is essentially maintained along of all extensive simulations. On the other hand, at basic pH, when all histidine residues are in the unprotonated form, the protein structure is very stable for ionic strengths ranging from 0 to 225 mM. Therefore, our findings support the hypothesis that the histidines constitute the hot points that induce configurational changes of E protein in acid pH, and give extra motivation to the development of new ideas for antivirus compound design.

Characterization of temperature dependent and substrate-binding cleft movements in Bacillus circulans family 11 xylanase: a molecular dynamics investigation.

BACKGROUND: Xylanases (EC 3.2.1.8) hydrolyze xylan, one of the most abundant plant polysaccharides found in nature, and have many potential applications in biotechnology. METHODS: Molecular dynamics simulations were used to investigate the effects of temperature between 298 to 338 K and xylobiose binding on residues located in the substrate-binding cleft of the family 11 xylanase from Bacillus circulans (BcX). RESULTS: In the absence of xylobiose the BcX exhibits temperature dependent movement of the thumb region which adopts an open conformation exposing the active site at the optimum catalytic temperature (328 K). In the presence of substrate, the thumb region restricts access to the active site at all temperatures, and this conformation is maintained by substrate/protein hydrogen bonds involving active site residues, including hydrogen bonds between Tyr69 and the 2' hydroxyl group of the substrate. Substrate access to the active site is regulated by temperature dependent motions that are restricted to the thumb region, and the BcX/substrate complex is stabilized by extensive intermolecular hydrogen bonding with residues in the active site. GENERAL SIGNIFICANCE: These results call for a revision of both the "hinge-bending" model for the activity of group 11 xylanases, and the role of Tyr69 in the catalytic mechanism.

Use of molecular dynamics data in biochemistry courses: An amphipathy scale to determine protein α-helix transmembrane segments.

The aim of this manuscript is to explain the application of an amphipathy scale obtained from molecular dynamics simulations and to demonstrate how it can be useful in the protein structure field. It is shown that this scale is easy to be used with the advantage of revealing domains of transmembrane α-helix of proteins without the need of knowing anything besides the protein primary structure. In addition, it allows the students to correlate concepts of protein structure and function, energy minimization, molecular dynamics simulations, and protein location.

Study of the interaction of human defensins with cell membrane models: relationships between structure and biological activity.

The HNP-1, HNP-2, and HNP-3 defensins are human antimicrobial peptides produced in response to microbial invasion. Their properties are distinct, with a more potent action for HNP-3. In this study, the relationship between their structural dissimilarities and their different microbial actions was evaluated by molecular dynamics simulation. Structural determinants related to their intra- and intermolecular interactions were defined for each HNP using a simplified membrane model consisting of a water/n-hexane interface. The hydrophobic portion of the HNPs promotes their diffusion to the interface with a concomitant, slight change in the structure induced by the intermolecular electrostatic interactions between the HPN molecules and the interface. As a consequence, different orientations are probably adopted by the HNPs at the interface, which may explain their different actions. The interaction of HNP-1 and HNP-2 with the surfaces was also studied using Langmuir monolayers as a biomimetic system. It was found that peptides adsorb rapidly at n-hexane/water interfaces as well as at phospholipid Langmuir monolayers but not at the air/liquid interface. This reveals that the presence of an organic phase is required for the exposure of the hydrophobic groups of the peptides. In addition, adsorption kinetics and surface pressure-area isotherms for Langmuir monolayers suggested that the lipid-peptide interaction is strongly influenced by the monolayer electrical charge and packing, depending also on the HPN structure. This study supports a model in which defensins, acting in a dimeric form, are able to disrupt membranes. The model also shows that the individual structures of the HNPs are responsible for their different actions on microbes.

Correlation of temperature induced conformation change with optimum catalytic activity in the recombinant G/11 xylanase A from Bacillus subtilis strain 168 (1A1).

The 1.7A resolution crystal structure of recombinant family G/11 beta-1,4-xylanase (rXynA) from Bacillus subtilis 1A1 shows a jellyroll fold in which two curved beta-sheets form the active-site and substrate-binding cleft. The onset of thermal denaturation of rXynA occurs at 328 K, in excellent agreement with the optimum catalytic temperature. Molecular dynamics simulations at temperatures of 298-328 K demonstrate that below the optimum temperature the thumb loop and palm domain adopt a closed conformation. However, at 328 K these two domains separate facilitating substrate access to the active-site pocket, thereby accounting for the optimum catalytic temperature of the rXynA.

Liquid correlation across the walls in a slit pore: effect on the wetting and drying transition.

The liquid structure next to the walls of a slit pore, immersed in a model simple liquid, is studied through a liquid theory and grand canonical Monte Carlo simulations. A liquid correlation across slit walls, of finite width, is found. This correlation modifies the structure and capillary partial wetting and drying transitions of the nonhomogeneous fluid, when close to its liquid-vapor coexistence curve.

The pH dependence of flavivirus envelope protein structure: insights from molecular dynamics simulations.

The flavivirus membrane fusion is triggered by the acid pH of the endosomes after virus endocytosis. The proposed mechanism involves changes in the protonation state of conserved histidine residues of the E protein present in the viral surface that undergoes a series of structural rearrangements that result in the fusion between the endosome and viral bilayers. We studied the pH dependence of E protein rearrangements of dengue virus type 2, used as a model, in the pH range experimented by the virus along the fusion process. We employed a low computational cost scheme to explore the behavior of the E protein by molecular dynamics (MD) simulations of complete systems that include the protein, the solvent, and ions. The procedure alternates cyclically the update of the ionization states of the protein residues with common MD steps applied to the new ionization configuration. Important pH-dependent protein structure rearrangements consistent with the changes of the protonation states of conserved histidine residues were observed. The involvement of other conserved residues in the flavivirus in the rearrangements was also identified. The results show interesting correlations with a proposed model for the fusion mechanism, as well as the experimentally identified key residues, contributing to a better understanding of the structural changes in protein E that lead to the fusion process.

New pockets in dengue virus 2 surface identified by molecular dynamics simulation.

One of the factors limiting the search of new compounds based on the structure of target proteins involved in diseases is the limited amount of target structural information. Great advances in the search for lead compounds could be achieved to find new cavities in protein structures that are generated using well established computational chemistry tools. In the case of dengue, the discovery of pockets in the crystallographic structure of the E protein has contributed to the search for lead compounds aimed at interfering in conformational transitions involved in the pH-dependent fusion process. This is a complex mechanism triggered by the acid pH of the endosomes that leads to the initial changes in the E protein assembly at the virus surface. In the present work, an arrangement of three ectodomain portions of the E protein present on the surface of the mature dengue virus was studied through long all-atom molecular dynamics simulations with explicit solvent. In order to identify new pockets and to evaluate the influence of the acid pH on these pockets, the physiological neutral pH conditions and the acid pH of the endosomes that trigger the fusion process were modeled. Several pockets presenting pH-dependent characteristics were found in the contact regions between the chains. Pockets at the protein-protein interfaces induced by a monomer in another monomer were also found. Some of the pockets are good candidates for the design of lead compounds that could interfere in the rearrangements in E proteins along the fusion process contributing to the development of specific inhibitors of the dengue disease.

Effect of the thermostat in the molecular dynamics simulation on the folding of the model protein chignolin.

Molecular dynamics simulations of the model protein chignolin with explicit solvent were carried out, in order to analyze the influence of the Berendsen thermostat on the evolution and folding of the peptide. The dependence of the peptide behavior on temperature was tested with the commonly employed thermostat scheme consisting of one thermostat for the protein and another for the solvent. The thermostat coupling time of the protein was increased to infinity, when the protein is not in direct contact with the thermal bath, a situation known as minimally invasive thermostat. In agreement with other works, it was observed that only in the last situation the instantaneous temperature of the model protein obeys a canonical distribution. As for the folding studies, it was shown that, in the applications of the commonly utilized thermostat schemes, the systems are trapped in local minima regions from which it has difficulty escaping. With the minimally invasive thermostat the time that the protein needs to fold was reduced by two to three times. These results show that the obstacles to the evolution of the extended peptide to the folded structure can be overcome when the temperature of the peptide is not directly controlled.

Structure of the Ebola fusion peptide in a membrane-mimetic environment and the interaction with lipid rafts.

The fusion peptide EBO(16) (GAAIGLAWIPYFGPAA) comprises the fusion domain of an internal sequence located in the envelope fusion glycoprotein (GP2) of the Ebola virus. This region interacts with the cellular membrane of the host and leads to membrane fusion. To gain insight into the mechanism of the peptide-membrane interaction and fusion, insertion of the peptide was modeled by experiments in which the tryptophan fluorescence and (1)H NMR were monitored in the presence of sodium dodecyl sulfate micelles or in the presence of detergent-resistant membrane fractions. In the presence of SDS micelles, EBO(16) undergoes a random coil-helix transition, showing a tendency to self-associate. The three-dimensional structure displays a 3(10)-helix in the central part of molecule, similar to the fusion peptides of many known membrane fusion proteins. Our results also reveal that EBO(16) can interact with detergent-resistant membrane fractions and strongly suggest that Trp-8 and Phe-12 are important for structure maintenance within the membrane bilayer. Replacement of tryptophan 8 with alanine (W8A) resulted in dramatic loss of helical structure, proving the importance of the aromatic ring in stabilizing the helix. Molecular dynamics studies of the interaction between the peptide and the target membrane also corroborated the crucial participation of these aromatic residues. The aromatic-aromatic interaction may provide a mechanism for the free energy coupling between random coil-helical transition and membrane anchoring. Our data shed light on the structural "domains" of fusion peptides and provide a clue for the development of a drug that might block the early steps of viral infection.

LMProt: an efficient algorithm for Monte Carlo sampling of protein conformational space.

A new and efficient Monte Carlo algorithm for sampling protein configurations in the continuous space is presented; the efficiency of this algorithm, named Local Moves for Proteins (LMProt), was compared to other alternative algorithms. For this purpose, we used an intrachain interaction energy function that is proportional to the root mean square deviation (rmsd) with respect to alpha-carbons from native structures of real proteins. For phantom chains, the LMProt method is approximately 10(4) and 20 times faster than the algorithms Thrashing (no local moves) and Sevenfold Way (local moves), respectively. Additionally, the LMProt was tested for real chains (excluded-volume all-atoms model); proteins 5NLL (138 residues) and 1BFF (129 residues) were used to determine the folding success xi as a function of the number eta of residues involved in the chain movements, and as a function of the maximum amplitude of atomic displacement delta r(max). Our results indicate that multiple local moves associated with relative chain flexibility, controlled by appropriate adjustments for eta and delta r(max), are essential for configurational search efficiency.