Theoretical and Molecular Docking Study of Ketoconazole on Heptakis(2,3,6‐tri‐O‐methyl)‐β‐cyclodextrin as Chiral Selector

A molecular docking study, using molecular mechanics calculations with AutoDock and semi‐empirical PM3 calculations, was used to predict the enantiodiscrimination of heptakis(2,3,6‐tri‐O‐methyl)‐β‐cyclodextrin (TMβCD) and ketoconazole (KTZ) enantiomers. A Density Functional Theory (DFT) single‐point calculation at the level of B3LYP/6‐311G (d,p) was performed for the PM3‐optimized complexes to obtain more accurate binding energy and the electronic structures of the complexes. The difference in energies of the inclusion complexes between the KTZ enantiomers and TMβCD is probably a measure of chiral discrimination, which results in the separation of the enantiomers as observed in the experimental studies. Chirality 28:209–214, 2016. © 2015 Wiley Periodicals, Inc.     

Phylogenetic congruence between subtropical trees and their associated fungi

Recent studies have detected phylogenetic signals in pathogen–host networks for both soil‐borne and leaf‐infecting fungi, suggesting that pathogenic fungi may track or coevolve with their preferred hosts. However, a phylogenetically concordant relationship between multiple hosts and multiple fungi in has rarely been investigated. Using next‐generation high‐throughput DNA sequencing techniques, we analyzed fungal taxa associated with diseased leaves, rotten seeds, and infected seedlings of subtropical trees. We compared the topologies of the phylogenetic trees of the soil and foliar fungi based on the internal transcribed spacer (ITS) region with the phylogeny of host tree species based on matK, rbcL, atpB, and 5.8S genes. We identified 37 foliar and 103 soil pathogenic fungi belonging to the Ascomycota and Basidiomycota phyla and detected significantly nonrandom host–fungus combinations, which clustered on both the fungus phylogeny and the host phylogeny. The explicit evidence of congruent phylogenies between tree hosts and their potential fungal pathogens suggests either diffuse coevolution among the plant–fungal interaction networks or that the distribution of fungal species tracked spatially associated hosts with phylogenetically conserved traits and habitat preferences. Phylogenetic conservatism in plant–fungal interactions within a local community promotes host and parasite specificity, which is integral to the important role of fungi in promoting species coexistence and maintaining biodiversity of forest communities.     

1,3‐propanediol binds deep inside the channel to inhibit water permeation through aquaporins

Aquaporins and aquaglyceroporins (AQPs) are membrane channel proteins responsible for transport of water and for transport of glycerol in addition to water across the cell membrane, respectively. They are expressed throughout the human body and also in other forms of life. Inhibitors of human AQPs have been sought for therapeutic treatment for various medical conditions including hypertension, refractory edema, neurotoxic brain edema, and so forth. Conducting all‐atom molecular dynamics simulations, we computed the binding affinity of acetazolamide to human AQP4 that agrees closely with in vitro experiments. Using this validated computational method, we found that 1,3‐propanediol (PDO) binds deep inside the AQP4 channel to inhibit that particular aquaporin efficaciously. Furthermore, we used the same method to compute the affinities of PDO binding to four other AQPs and one aquaglyceroporin whose atomic coordinates are available from the protein data bank (PDB). For bovine AQP1, human AQP2, AQP4, AQP5, and Plasmodium falciparum PfAQP whose structures were resolved with high resolution, we obtained definitive predictions on the PDO dissociation constant. For human AQP1 whose PDB coordinates are less accurate, we estimated the dissociation constant with a rather large error bar. Taking into account the fact that PDO is generally recognized as safe by the US FDA, we predict that PDO can be an effective diuretic which directly modulates water flow through the protein channels. It should be free from the serious side effects associated with other diuretics that change the hydro‐homeostasis indirectly by altering the osmotic gradients.     

Protein side chain conformation predictions with an MMGBSA energy function

The prediction of protein side chain conformations from backbone coordinates is an important task in structural biology, with applications in structure prediction and protein design. It is a difficult problem due to its combinatorial nature. We study the performance of an “MMGBSA” energy function, implemented in our protein design program Proteus, which combines molecular mechanics terms, a Generalized Born and Surface Area (GBSA) solvent model, with approximations that make the model pairwise additive. Proteus is not a competitor to specialized side chain prediction programs due to its cost, but it allows protein design applications, where side chain prediction is an important step and MMGBSA an effective energy model. We predict the side chain conformations for 18 proteins. The side chains are first predicted individually, with the rest of the protein in its crystallographic conformation. Next, all side chains are predicted together. The contributions of individual energy terms are evaluated and various parameterizations are compared. We find that the GB and SA terms, with an appropriate choice of the dielectric constant and surface energy coefficients, are beneficial for single side chain predictions. For the prediction of all side chains, however, errors due to the pairwise additive approximation overcome the improvement brought by these terms. We also show the crucial contribution of side chain minimization to alleviate the rigid rotamer approximation. Even without GB and SA terms, we obtain accuracies comparable to SCWRL4, a specialized side chain prediction program. In particular, we obtain a better RMSD than SCWRL4 for core residues (at a higher cost), despite our simpler rotamer library. Proteins 2016; 84:803–819. © 2016 Wiley Periodicals, Inc.     

A multi‐resolution model to capture both global fluctuations of an enzyme and molecular recognition in the ligand‐binding site

In multi‐resolution simulations, different system components are simultaneously modeled at different levels of resolution, these being smoothly coupled together. In the case of enzyme systems, computationally expensive atomistic detail is needed in the active site to capture the chemistry of ligand binding. Global properties of the rest of the protein also play an essential role, determining the structure and fluctuations of the binding site; however, these can be modeled on a coarser level. Similarly, in the most computationally efficient scheme only the solvent hydrating the active site requires atomistic detail. We present a methodology to couple atomistic and coarse‐grained protein models, while solvating the atomistic part of the protein in atomistic water. This allows a free choice of which protein and solvent degrees of freedom to include atomistically. This multi‐resolution methodology can successfully model stable ligand binding, and we further confirm its validity by exploring the reproduction of system properties relevant to enzymatic function. In addition to a computational speedup, such an approach can allow the identification of the essential degrees of freedom playing a role in a given process, potentially yielding new insights into biomolecular function. Proteins 2016; 84:1902–1913. © 2016 Wiley Periodicals, Inc.     

Morphology versus DNA barcoding: two sides of the same coin. A case study of Ceutorhynchus erysimi and C. contractus identification

Genotyping of 2 well‐known weevil species from the genus Ceutorhynchus (Coleoptera: Curculionidae) distributed in west Palearctic, C. erysimi and C. contractus, revealed phenotype versus genotype inconsistencies in a set of 56 specimens (25 C. erysimi and 31 C. contractus) collected from 25 locations in Serbia and Montenegro. An analysis of the mitochondrial cytochrome oxidase subunit I gene (COI), widely used as a barcoding region, and a nuclear gene, elongation factor‐1α (EF‐1α), revealed stable genetic divergence among these species. The average uncorrected pairwise distances for the COI and EF‐1α genes were 3.8%, and 1.3%, respectively, indicating 2 genetically well‐segregated species. However, the genetic data were not congruent with the phenotypic characteristics of the studied specimens. In the first place, C. erysimi genotypes were attached to specimens with phenotypic characteristics of C. contractus. Species‐specific PCR‐RFLP assays for the barcoding gene COI were applied for the molecular identification of 101 additional specimens of both morphospecies (33 C. erysimi and 68 C. contractus) and were found to confirm this incongruity. The discrepancy between the genetic and morphological data raises the question of the accuracy of using a barcoding approach, as it may result in misleading conclusions about the taxonomic position of the studied organism. Additionally, the typological species concept shows considerable weakness when genetic data are not supported with phenotypic characteristics as in case of asymmetric introgression, which may cause certain problems, especially in applied studies such as biological control programs in which the biological properties of the studied organisms are the main focus.     

A Novel Topology of Proline-rich Transmembrane Protein 2 (PRRT2): HINTS FOR AN INTRACELLULAR FUNCTION AT THE SYNAPSE*

Proline-rich transmembrane protein 2 (PRRT2) has been identified as the single causative gene for a group of paroxysmal syndromes of infancy, including epilepsy, paroxysmal movement disorders, and migraine. On the basis of topology predictions, PRRT2 has been assigned to the recently characterized family of Dispanins, whose members share the two-transmembrane domain topology with a large N terminus and short C terminus oriented toward the outside of the cell. Because PRRT2 plays a role at the synapse, it is important to confirm the exact orientation of its N and C termini with respect to the plasma membrane to get clues regarding its possible function. Using a combination of different experimental approaches, including live immunolabeling, immunogold electron microscopy, surface biotinylation and computational modeling, we demonstrate a novel topology for this protein. PRRT2 is a type II transmembrane protein in which only the second hydrophobic segment spans the plasma membrane, whereas the first one is associated with the internal surface of the membrane and forms a helix-loop-helix structure without crossing it. Most importantly, the large proline-rich N-terminal domain is not exposed to the extracellular space but is localized intracellularly, and only the short C terminus is extracellular (N_cyt/C_exo topology). Accordingly, we show that PRRT2 interacts with the Src homology 3 domain-bearing protein Intersectin 1, an intracellular protein involved in synaptic vesicle cycling. These findings will contribute to the clarification of the role of PRRT2 at the synapse and the understanding of pathogenic mechanisms on the basis of PRRT2-related neurological disorders.     

First Report of Black Leaf Spot Caused by Alternaria alternata on Ramie in China

In 2012 and 2013, black leaf spot disease was observed on ramie plants in Hunan and Hubei Provinces, China. In the field, the symptoms of this disease included dark green to black big spots on leaves, often resulting in upwardly curled leaf margins. The pathogen isolates were identified as Alternaria alternata (Fr.) Keissler on the basis of morphology and sequence similarity of 99–100% to the published data for internal transcribed spacer (ITS), and glyceraldehyde‐3‐phosphate dehydrogenase (gdp). To our knowledge, this is the first report of Alternaria leaf spot of ramie in China.     

Nonnative Electrostatic Interactions Can Modulate Protein Folding: Molecular Dynamics with a Grain of Salt

In recent years, a growing number of protein folding studies have focused on the unfolded state, which is now recognized as playing a major role in the folding process. Some of these studies show that interactions occurring in the unfolded state can significantly affect the stability and kinetics of the protein folding reaction. In this study, we modeled the effect of electrostatic interactions, both native and nonnative, on the folding of three protein systems that underwent selective charge neutralization or reversal or complete charge suppression. In the case of the N-terminal L9 protein domain, our results directly attribute the increase in thermodynamic stability to destabilization of the unfolded ensemble, reaffirming the experimental observations. These results provide a deeper structural insight into the ensemble of the unfolded state and predict a new mutation site for increased protein stability. In the second case, charge reversal mutations of RNase Sa affected protein stability, with the destabilizing mutations being less destabilizing at higher salt concentrations, indicating the formation of charge-charge interactions in the unfolded state. In the N-terminal L9 and RNase Sa systems, changes in electrostatic interactions in the unfolded state that cause an increase in free energy had an overall compaction effect that suggests a decrease in entropy. In the third case, in which we compared the β-lactalbumin and hen egg-white lysozyme protein homologues, we successfully eliminated differences between the folding kinetics of the two systems by suppressing electrostatic interactions, supporting previously reported findings. Our coarse-grained molecular dynamics study not only reproduces experimentally reported findings but also provides a detailed molecular understanding of the elusive unfolded-state ensemble and how charge-charge interactions can modulate the biophysical characteristics of folding.