Molecular dynamics of water in the neighborhood of aquaporins.

Present knowledge obtained by molecular dynamics (MD) simulation studies regarding the dynamics of water, both in the vicinity of biological membranes and within the proteinaceous water channels, also known as aquaporins (AQPs), is reviewed. A brief general summary of the water models most extensively employed in MD simulations (SPC, SPC/E, TIP3P, TIP4P), indicating their most relevant pros and cons, is likewise provided. Structural considerations of water are also discussed, based on different order parameters, which can be extracted from MD simulations as well as from experiments. Secondly, the behaviour of water in the neighbourhood of membranes by means of molecular dynamics simulations is addressed. Consequently, the comparison with previous experimental evidence is pointed out. In living cells, water is transported across the plasma membrane through the lipid bilayer and the aforementioned AQPs, which motivates this review to focus mostly on MD simulation studies of water within AQPs. Relevant contributions explaining peculiar properties of these channels are discussed, such as selectivity and gating. Water models used in these studies are also summarised. Finally, based on the information presented here, further MD studies are encouraged.

A new analogue of islet neogenesis associated protein with higher structural and plasma stability.

Islet Neogenesis Associated Protein pentadecapeptide (INGAP-PP) increases ß-cell mass and function in experimental animals. A short clinical trial also yielded promising results. However, HTD4010, a new peptide derived from INGAP-PP, was developed in order to optimize its specific effects by minimizing its side effects. To study and compare the tertiary structure, stability dynamics, and plasma stability of HTD4010, an INGAP-PP analogue. Both peptides were pre-incubated in human, rat and mouse plasma at 37 °C, and their presence was identified and quantified by high performance liquid chromatography at different time-points. GROMACS 2019 package and the Gromos 54A7 force field were used to evaluate overall correlated motion of the peptide molecule during molecular dynamics simulation by essential dynamics. HTD4010 exhibited significantly larger plasma stability than INGAP-PP, and its structural stability was almost 3.36-fold higher than INGAP-PP. These results suggest that HTD4010 may facilitate longer tissue interaction, thereby developing higher potential biological effects. If so, HTD4010 may become a promising therapeutic agent to treat people with diabetes. Communicated by Ramaswamy H. Sarma.

The behavior of the hydrophobic effect under pressure and protein denaturation.

It is well known that proteins denature under high pressure. The mechanism that underlies such a process is still not clearly understood, however, giving way to controversial interpretations. Using molecular dynamics simulation on systems that may be regarded experimentally as limiting examples of the effect of high pressure on globular proteins, such as lysozyme and apomyoglobin, we have effectively reproduced such similarities and differences in behavior as are interpreted from experiment. From the analysis of such data, we explain the experimental evidence at hand through the effect of pressure on the change of water structure, and hence the weakening of the hydrophobic effect that is known to be the main driving force in protein folding.

Phylogenomics of tomato chloroplasts using assembly and alignment-free method.

Phylogenetics and population genetics are central disciplines in evolutionary biology. Both are based on the comparison of single DNA sequences, or a concatenation of a number of these. However, with the advent of next-generation DNA sequencing technologies, the approaches that consider large genomic data sets are of growing importance for the elucidation of evolutionary relationships among species. Among these approaches, the assembly and alignment-free methods which allow an efficient distance computation and phylogeny reconstruction are of great importance. However, it is not yet clear under what quality conditions and abundance of genomic data such methods are able to infer phylogenies accurately. In the present study we assess the method originally proposed by Fan et al. for whole genome data, in the elucidation of Tomatoes' chloroplast phylogenetics using short read sequences. We find that this assembly and alignment-free method is capable of reproducing previous results under conditions of high coverage, given that low frequency k-mers (i.e. error prone data) are effectively filtered out. Finally, we present a complete chloroplast phylogeny for the best data quality candidates of the recently published 360 tomato genomes.

Integument CYP genes of the largest genome-wide cytochrome P450 expansions in triatomines participate in detoxification in deltamethrin-resistant Triatoma infestans.

Insect resistance to chemical insecticides is attributed to a combination of different mechanisms, such as metabolic resistance, knockdown resistance, and the cuticular resistance or penetration factor. The insect integument offers an efficient barrier against contact insecticides and its role as penetration factor has been previously reported; however, there is no information about its potential function in the metabolic resistance. Cytochrome P450 genes (CYP) are highly expressed in the fat body of several insects and thus play a key role in their metabolic resistance. Here, we describe new members that belong to the highly genome-wide expanded CYP3093A and CYP4EM subfamilies in the Chagas disease vectors Rhodnius prolixus and Triatoma infestans. We modeled the docking of deltamethrin in their active site and detected differences in some amino acids between both species that are critical for a correct interaction with the substrate. We also knocked down the two constitutively most expressed genes in the integument of resistant T. infestans nymphs (CYP3093A11 and CYP4EM10) in order to find clues on their participation in deltamethrin resistance. This is the first report on the role of the insect integument in detoxification events; although these two CYP genes do not fully explain the resistance observed in T. infestans.

Pressure denaturation of apomyoglobin: a molecular dynamics simulation study.

The effect of pressure on the structure and mobility of Sperm Wale Apomyoglobin was studied by Molecular Dynamics computer simulation at 1 bar and 3 kbar (1 atm=1.01325 bar=101.325 kPa). The results are in good agreement with the available experimental data, allowing further analysis of other features of the effect of pressure on the protein solution. From the analysis of Secondary Structures (SS) along the trajectories it is observed that alpha-helixes are favoured under pressure at the expense of bends, turns and 3-helixes. The studies of mobility show that although the general mobility is restricted under pressure this is not true for some particular residues. The studies of tertiary structure show important conformational changes. The evolution of the Solvent Accessed Surface (SAS) with pressure shows a notorious increase due almost completely to a biased raise in the hydrophobic area exposed, which consequently shows that the hydrophobic interaction is considerably weaker under high hydrostatic pressure conditions.

Effect of pressure on the conformation of proteins. A molecular dynamics simulation of lysozyme.

The effect of pressure on the structure and mobility of lysozyme was studied by molecular dynamics computer simulation at 1 and 3 kbar (1 atm = 1.01325 bar = 101.325 kPa). The results have good agreement with the available experimental data, allowing the analysis of other features of the effect of pressure on the protein solution. The studies of mobility show that although the general mobility is restricted under pressure this is not true for some particular residues. From the analysis of secondary structure along the trajectories it is observed that the conformation under pressure is more stable, suggesting that pressure acts as a 'conformer selector' on the protein. The difference in solvent-accessed surface (SAS) with pressure shows a clear inversion of the hydrophilic/hydrophobic SAS ratio, which consequently shows that the hydrophobic interaction is considerably weaker under high hydrostatic pressure conditions.

A molecular dynamics approach to ligand-receptor interaction in the aspirin-human serum albumin complex.

In this work, we present a study of the interaction between human serum albumin (HSA) and acetylsalicylic acid (ASA, C(9)H(8)O(4)) by molecular dynamics simulations (MD). Starting from an experimentally resolved structure of the complex, we performed the extraction of the ligand by means of the application of an external force. After stabilization of the system, we quantified the force used to remove the ASA from its specific site of binding to HSA and calculated the mechanical nonequilibrium external work done during this process. We obtain a reasonable value for the upper boundary of the Gibbs free energy difference (an equilibrium thermodynamic potential) between the complexed and noncomplexed states. To achieve this goal, we used the finite sampling estimator of the average work, calculated from the Jarzynski Equality. To evaluate the effect of the solvent, we calculated the so-called "viscous work," that is, the work done to move the aspirin in the same trajectory through the solvent in absence of the protein, so as to assess the relevance of its contribution to the total work. The results are in good agreement with the available experimental data for the albumin affinity constant for aspirin, obtained through quenching fluorescence methods.

Islet neogenesis associated protein (ingap): structural and dynamical properties of its active pentadecapeptide.

We have studied the structural and dynamical properties of the biologically active pentadecapeptide of the islet neogenesis associated protein (INGAP-PP) and of two other pentadecapeptides with the same amino acid composition but randomly scrambled primary sequences, using molecular dynamic simulations. Our data demonstrates that whilst the peptides with scrambled sequences show no definite prevalent structure in solution, INGAP-PP maintains a notably stable tertiary fold, namely, a conformer with a central beta-sheet and closed C-terminal. Such structure resembles the one corresponding to the amino acid sequence of human pancreatitis associated protein-1 (PAP-1), which presents 85% sequence homology with INGAP. These results could reasonably explain why the two scrambled sequences tested showed no biological activity, while INGAP-PP significantly increases beta-cells function and mass both in vitro and in vivo conditions. The capability of INGAP-PP to temporarily adopt other closely related conformations offers also a plausible explanation for the 50 fold experimental difference in potency between the active pentadecapeptide and the whole protein. They also suggest that the C-terminal region of INGAP-PP may plausibly be the locus for its interaction with the cell receptor. Consequently, the knowledge gathered through our data can help to obtain more potent INGAP-PP analogs, suitable for the prevention and treatment of diabetes.

Studying the unfolding kinetics of proteins under pressure using long molecular dynamic simulation runs.

The usefulness of computational methods such as molecular dynamics simulation has been extensively established for studying systems in equilibrium. Nevertheless, its application to complex non-equilibrium biological processes such as protein unfolding has been generally regarded as producing results which cannot be interpreted straightforwardly. In the present study, we present results for the kinetics of unfolding of apomyoglobin, based on the analysis of long simulation runs of this protein in solution at 3 kbar (1 atm = 1.01325, bar = 101,325 Pa). We hereby demonstrate that the analysis of the data collected within a simulated time span of 0.18 mus suffices for producing results, which coincide remarkably with the available unfolding kinetics experimental data. This not only validates molecular dynamics simulation as a valuable alternative for studying non-equilibrium processes, but also enables a detailed analysis of the actual structural mechanism which underlies the unfolding process of proteins under elusive denaturing conditions such as high pressure.

Clustering of Lennard-Jones particles in water: temperature and pressure effects.

While the hydrophobic effect is, for many systems, one of the most relevant interactions, it may be said that in the case of biological systems this effect becomes of determinant importance. Although the matter has been analyzed extensively, certain aspects are yet to be elucidated. Hence, the study on the behavior of the hydrophobic effect with temperature, and particularly with pressure deserves further investigation; model systems may help us in the task. We have analyzed the behavior of Lennard-Jones particles in water by means of molecular dynamics simulation under different conditions of size, concentration, temperature, and pressure. Following the formation of particle aggregates we can observe an increase of the hydrophobic effect with temperature and a strong weakening of the effect at high pressures. The results agree with the experimental evidence and show the ability of molecular dynamics simulation to account for the behavior of nonpolar substances under different conditions, provided that the intermolecular interactions used are adequate.