Three-dimensional (3D) structure prediction of the American and African oil-palms β-ketoacyl-[ACP] synthase-II protein by comparative modelling

Background: The fatty-acid profile of the vegetable oils determines its properties and nutritional value. Palm-oil obtained from the African oil-palm [Elaeis guineensis Jacq. (Tenera)] contains 44% palmitic acid (C_16:0), but, palm-oil obtained from the American oilpalm [Elaeis oleifera] contains only 25% C16:0. In part, the b-ketoacyl-[ACP] synthase II (KASII) [EC: 2.3.1.179] protein is responsible for the high level of C_16:0 in palm-oil derived from the African oil-palm. To understand more about E. guineensis KASII (EgKASII) and E. oleifera KASII (EoKASII) proteins, it is essential to know its structures. Hence, this study was undertaken. Objective: The objective of this study was to predict three-dimensional (3D) structure of EgKASII and EoKASII proteins using molecular modelling tools. Materials and Methods: The amino-acid sequences for KASII proteins were retrieved from the protein database of National Center for Biotechnology Information (NCBI), USA. The 3D structures were predicted for both proteins using homology modelling and ab-initio technique approach of protein structure prediction. The molecular dynamics (MD) simulation was performed to refine the predicted structures. The predicted structure models were evaluated and root mean square deviation (RMSD) and root mean square fluctuation (RMSF) values were calculated. Results: The homology modelling showed that EgKASII and EoKASII proteins are 78% and 74% similar with Streptococcus pneumonia KASII and Brucella melitensis KASII, respectively. The EgKASII and EoKASII structures predicted by using ab-initio technique approach shows 6% and 9% deviation to its structures predicted by homology modelling, respectively. The structure refinement and validation confirmed that the predicted structures are accurate. Conclusion: The 3D structures for EgKASII and EoKASII proteins were predicted. However, further research is essential to understand the interaction of EgKASII and EoKASII proteins with its substrates.     

Atomic-level insights into metabolite recognition and specificity of the SAM-II riboswitch

Although S-adenosylhomocysteine (SAH), a metabolic by-product of S-adenosylmethionine (SAM), differs from SAM only by a single methyl group and an overall positive charge, SAH binds the SAM-II riboswitch with more than 1000-fold less affinity than SAM. Using atomistic molecular dynamics simulations, we investigated the molecular basis of such high selectivity in ligand recognition by SAM-II riboswitch. The biosynthesis of SAM exclusively generates the (S,S) stereoisomer, and (S,S)-SAM can spontaneously convert to the (R,S) form. We, therefore, also examined the effects of (R,S)-SAM binding to SAM-II and its potential biological function. We find that the unfavorable loss in entropy in SAM-II binding is greater for (S,S)- and (R,S)-SAM than SAH, which is compensated by stabilizing electrostatic interactions with the riboswitch. The positively charged sulfonium moiety on SAM acts as the crucial anchor point responsible for the formation of key ionic interactions as it fits favorably in the negatively charged binding pocket. In contrast, SAH, with its lone pair of electrons on the sulfur, experiences repulsion in the binding pocket of SAM-II and is enthalpically destabilized. In the presence of SAH, similar to the unbound riboswitch, the pseudoknot structure of SAM-II is not completely formed, thus exposing the Shine-Dalgarno sequence. Unlike SAM, this may further facilitate ribosomal assembly and translation initiation. Our analysis of the conformational ensemble sampled by SAM-II in the absence of ligands and when bound to SAM or SAH reveals that ligand binding follows a combination of conformational selection and induced-fit mechanisms.     

Stochastic Gating and Drug-Ribosome Interactions

Gentamicin is a potent antibiotic that is used in combination therapy for inhalation anthrax disease. The drug is also often used in therapy for methicillin-resistant Staphylococcusaureus. Gentamicin works by flipping a conformational switch on the ribosome, disrupting the reading head (i.e., 16S ribosomal decoding bases 1492-1493) used for decoding messenger RNA. We use explicit solvent all-atom molecular simulation to study the thermodynamics of the ribosomal decoding site and its interaction with gentamicin. The replica exchange molecular dynamics simulations used an aggregate sampling of 15 μs when summed over all replicas, allowing us to explicitly calculate the free-energy landscape, including a rigorous treatment of enthalpic and entropic effects. Here, we show that the decoding bases flip on a timescale faster than that of gentamicin binding, supporting a stochastic gating mechanism for antibiotic binding, rather than an induced-fit model where the bases only flip in the presence of a ligand. The study also allows us to explore the nonspecific binding landscape near the binding site and reveals that, rather than a two-state bound/unbound scenario, drug dissociation entails shuttling between many metastable local minima in the free-energy landscape. Special care is dedicated to validation of the obtained results, both by direct comparison to experiment and by estimation of simulation convergence.     

Hydrophobic cooperativity as a mechanism for amyloid nucleation

The kinetics of amyloid fibril formation are in most cases explained by classical nucleation theory, yet the mechanisms behind nucleation are not well understood. We show using molecular dynamics simulations that the hydrophobic cooperativity in the self-association of the model amyloidogenic peptide STVIYE is sufficient to allow for nucleation-dependent polymerization with a pentamer critical nucleus. The role of electrostatics was also investigated. Novel considerations of the electrostatic solvation energy using the Born-Onsager equation are put forth to rationalize the aggregation of charged peptides and provide new insight into the energetic differences between parallel and antiparallel beta-sheets. Together these results help explain the influence of molecular charge in the class of fibril-forming hexapeptides recently designed by Serrano and collaborators.     

Alteration of the C-Terminal Ligand Specificity of the Erbin PDZ Domain by Allosteric Mutational Effects

Modulation of protein binding specificity is important for basic biology and for applied science. Here we explore how binding specificity is conveyed in PDZ (postsynaptic density protein-95/discs large/zonula occludens-1) domains, small interaction modules that recognize various proteins by binding to an extended C terminus. Our goal was to engineer variants of the Erbin PDZ domain with altered specificity for the most C-terminal position (position 0) where a Val is strongly preferred by the wild-type domain. We constructed a library of PDZ domains by randomizing residues in direct contact with position 0 and in a loop that is close to but does not contact position 0. We used phage display to select for PDZ variants that bind to 19 peptide ligands differing only at position 0. To verify that each obtained PDZ domain exhibited the correct binding specificity, we selected peptide ligands for each domain. Despite intensive efforts, we were only able to evolve Erbin PDZ domain variants with selectivity for the aliphatic C-terminal side chains Val, Ile and Leu. Interestingly, many PDZ domains with these three distinct specificities contained identical amino acids at positions that directly contact position 0 but differed in the loop that does not contact position 0. Computational modeling of the selected PDZ domains shows how slight conformational changes in the loop region propagate to the binding site and result in different binding specificities. Our results demonstrate that second-sphere residues could be crucial in determining protein binding specificity.     

High temperature unfolding of Bacillus anthracis amidase-03 by molecular dynamics simulations

The stability of amidase-03 structure (a cell wall hydrolase protein) from Bacillus anthracis was studied using classical molecular dynamics (MD) simulation. This protein (GenBank accession number: {"type":"entrez-protein","attrs":{"text":"NP_844822","term_id":"30262445","term_text":"NP_844822"}}NP_844822) contains an amidase-03 domain which is known to exhibit the catalytic activity of N-acetylmuramoyl-L-alanine amidase (digesting MurNAc-Lalanine linkage of bacterial cell wall). The amidase-03 enzyme has stability at high temperature due to the core formed by the combination of several secondary structure elements made of β-sheets. We used root-mean-square-displacement (RMSD) of the simulated structure from its initial state to demonstrate the unfolding of the enzyme using its secondary structural elements. Results show that amidase-03 unfolds in transition state ensemble (TSE). The data suggests that α-helices unfold before β-sheets from the core during simulation.     

DNA condensation by TmHU studied by optical tweezers,AFM and molecular dynamics simulations

The compaction of DNA by the HU protein from Thermotoga maritima (TmHU) is analysed on a single-molecule level by the usage of an optical tweezers-assisted force clamp. The condensation reaction is investigated at forces between 2 and 40 pN applied to the ends of the DNA as well as in dependence on the TmHU concentration. At 2 and 5 pN, the DNA compaction down to 30% of the initial end-to-end distance takes place in two regimes. Increasing the force changes the progression of the reaction until almost nothing is observed at 40 pN. Based on the results of steered molecular dynamics simulations, the first regime of the length reduction is assigned to a primary level of DNA compaction by TmHU. The second one is supposed to correspond to the formation of higher levels of structural organisation. These findings are supported by results obtained by atomic force microscopy.     

Donor-strand exchange in chaperone-assisted pilus assembly revealed in atomic detail by molecular dynamics

Adhesive multi-subunit fibres are assembled on the surface of many pathogenic bacteria via the chaperone-usher pathway. In the periplasm, a chaperone donates a β-strand to a pilus subunit to complement its incomplete immunoglobulin-like fold. At the outer membrane, this is replaced with a β-strand formed from the N-terminal extension (Nte) of an incoming pilus subunit by a donor-strand exchange (DSE) mechanism. This reaction has previously been shown to proceed via a concerted mechanism, in which the Nte interacts with the chaperone:subunit complex before the chaperone has been displaced, forming a ternary intermediate. Thereafter, the pilus and chaperone β-strands have been postulated to undergo a strand swap by a ‘zip-in-zip-out’ mechanism, whereby the chaperone strand zips out, residue by residue, as the Nte simultaneously zips in, although direct experimental evidence for a zippering mechanism is still lacking. Here, molecular dynamics simulations have been used to probe the DSE mechanism during formation of the Saf pilus from Salmonella enterica at the atomic level, allowing the direct investigation of the zip-in-zip-out hypothesis. The simulations provide an explanation of how the incoming Nte is able to dock and initiate DSE due to inherent dynamic fluctuations within the chaperone:subunit complex. In the simulations, the chaperone donor strand was seen to unbind from the pilus subunit, residue by residue, in direct support of the zip-in-zip-out hypothesis. In addition, an interaction of a residue towards the N-terminus of the Nte with a specific binding pocket (P*) on the adjacent pilus subunit was seen to stabilise the DSE product against unbinding, which also proceeded in the simulations by a zippering mechanism. Together, the study provides an in-depth picture of DSE, including the first atomistic insights into the molecular events occurring during the zip-in-zip-out mechanism.     

K+/Na+ selectivity in K channels and valinomycin: over-coordination versus cavity-size constraints

Potassium channels and valinomycin molecules share the exquisite ability to select K(+) over Na(+). Highly selective K channels maintain a special local environment around their binding sites devoid of competing hydrogen bond donor groups, which enables spontaneous transfer of K(+) from states of low coordinations in water into states of over-coordination by eight carbonyl ligands. In such a phase-activated state, electrostatic interactions from these 8-fold binding sites, constrained to maintain high coordinations, result in K(+)/Na(+) selectivity with no need for a specific cavity size. Under such conditions, however, direct coordination from five or six carbonyl ligands does not result in selectivity. Yet, valinomycin molecules achieve selectivity by providing only six carbonyl ligands. Does valinomycin use additional coordinating ligands from the solvent or does it have special structural features not present in K channels? Quantum chemical investigations undertaken here demonstrate that valinomycin selectivity is due to cavity size constraints that physically prevent it from collapsing onto the smaller sodium ion. Valinomycin enforces these constraints by using a combination of intramolecular hydrogen bonds and other structural features, including its specific ring size and the spacing between its connected ligands. Results of these investigations provide a consistent explanation for the experimental data available for the ion-complexation properties of valinomycin in solvents of varying polarity. Together, investigations of these two systems reveal how nature, despite being popular for its parsimony in recycling functional motifs, can use different combinations of phase, coordination number, cavity size, and rigidity (constraints) to achieve K(+)/Na(+) selectivity.     

Evidence for proton shuffling in a thioredoxin-like protein during catalysis

Proteins of the thioredoxin (Trx) superfamily catalyze disulfide-bond formation, reduction and isomerization in substrate proteins both in prokaryotic and in eukaryotic cells. All members of the Trx family with thiol-disulfide oxidoreductase activity contain the characteristic Cys-X-X-Cys motif in their active site. Here, using Poisson-Boltzmann-based protonation-state calculations based on 100-ns molecular dynamics simulations, we investigate the catalytic mechanism of DsbL, the most oxidizing Trx-like protein known to date. We observed several correlated transitions in the protonation states of the buried active-site cysteine and a neighboring lysine coupled to the exposure of the active-site thiolate. These results support the view of an internal proton shuffling mechanism during oxidation crucial for the uptake of two electrons from the substrate protein. Intramolecular disulfide-bond formation is probably steered by the conformational switch facilitating interaction with the active-site thiolate. A consistent catalytic mechanism for DsbL, probably conferrable to other proteins of the same class, is presented. Our results suggest a functional role of hydration entropy of active-site groups.     

Engineered Myosin VI Motors Reveal Minimal Structural Determinants of Directionality and Processivity

Myosins have diverse mechanical properties reflecting a range of cellular roles. A major challenge is to understand the structural basis for generating novel functions from a common motor core. Myosin VI (M6) is specialized for processive motion toward the (−) end of actin filaments. We have used engineered M6 motors to test and refine the “redirected power stroke” model for (−) end directionality and to explore poorly understood structural requirements for processive stepping. Guided by crystal structures and molecular modeling, we fused artificial lever arms to the catalytic head of M6 at several positions, retaining varying amounts of native structure. We found that an 18-residue α-helical insert is sufficient to reverse the directionality of the motor, with no requirement for any calmodulin light chains. Further, we observed robust processive stepping of motors with artificial lever arms, demonstrating that processivity can arise without optimizing lever arm composition or mechanics.     

Molecular Dynamics Simulations and Structure-Guided Mutagenesis Provide Insight into the Architecture of the Catalytic Core of the Ectoine Hydroxylase

Many bacteria amass compatible solutes to fend-off the detrimental effects of high osmolarity on cellular physiology and water content. These solutes also function as stabilizers of macromolecules, a property for which they are referred to as chemical chaperones. The tetrahydropyrimidine ectoine is such a compatible solute and is widely synthesized by members of the Bacteria. Many ectoine producers also synthesize the stress protectant 5-hydroxyectoine from the precursor ectoine, a process that is catalyzed by the ectoine hydroxylase (EctD). The EctD enzyme is a member of the non-heme-containing iron(II) and 2-oxoglutarate-dependent dioxygenase superfamily. A crystal structure of the EctD protein from the moderate halophile Virgibacillus salexigens has previously been reported and revealed the coordination of the iron catalyst, but it lacked the substrate ectoine and the co-substrate 2-oxoglutarate. Here we used this crystal structure as a template to assess the likely positioning of the ectoine and 2-oxoglutarate ligands within the active site by structural comparison, molecular dynamics simulations, and site-directed mutagenesis. Collectively, these approaches suggest the positioning of the iron, ectoine, and 2-oxoglutarate ligands in close proximity to each other and with a spatial orientation that will allow the region-selective and stereo-specific hydroxylation of (4S)-ectoine to (4S,5S)-5-hydroxyectoine. Our study thus provides a view into the catalytic core of the ectoine hydroxylase and suggests an intricate network of interactions between the three ligands and evolutionarily highly conserved residues in members of the EctD protein family.     

Structural Insights into Substrate Recognition and Product Expulsion in CTX-M Enzymes

β-Lactamase-mediated resistance to β-lactam antibiotics poses a major threat to our antibiotic armamentarium. Among β-lactamases, a significant threat comes from enzymes that hydrolyze extended-spectrum cephalosporins such as cefotaxime. Among the enzymes that exhibit this phenotype, the CTX-M family is found worldwide. These enzymes have a small active site, which makes it difficult to explain how they hydrolyze the bulky extended-spectrum cephalosporins into the binding site. We investigated noncovalent substrate recognition and product release in CTX-M enzymes using steered molecular dynamics simulation and X-ray diffraction. An arginine residue located far from the binding site favors the capture and tracking of substrates during entrance into the catalytic pocket. We show that the accommodation of extended-spectrum cephalosporins by CTX-M enzymes induced subtle changes in the active site and established a high density of electrostatic interactions. Interestingly, the product of the catalytic reaction initiates its own release because of steric hindrances and electrostatic repulsions. This suggests that there exists a general mechanism for product release for all members of the β-lactamase family and probably for most carboxypeptidases.     

Breeding system diversification and evolution in American Poa supersect. Homalopoa (Poaceae: Poeae: Poinae)

Background and Aims Poa subgenus Poa supersect. Homalopoa has diversified extensively in the Americas. Over half of the species in the supersection are diclinous; most of these are from the New World, while a few are from South-East Asia. Diclinism in Homalopoa can be divided into three main types: gynomonoecism, gynodioecism and dioecism. Here the sampling of species of New World Homalopoa is expanded to date its origin and diversification in North and South America and examine the evolution and origin of the breeding system diversity.Methods A total of 124 specimens were included in the matrix, of which 89 are species of Poa supersect. Homalopoa sections Acutifoliae, Anthochloa, Brizoides, Dasypoa, Dioicopoa, Dissanthelium, Homalopoa sensu lato (s.l.), Madropoa and Tovarochloa, and the informal Punapoa group. Bayesian and parsimony analyses were conducted on the data sets based on four markers: the nuclear ribosomal internal tanscribed spacer (ITS) and external transcribed spacer (ETS), and plastid trnT-L and trnL-F. Dating analyses were performed on a reduced Poa matrix and enlarged Poaceae outgroup to utilize fossils as calibration points. A relaxed Bayesian molecular clock method was used.Key Results Hermaphroditism appears to be pleisiomorphic in the monophyletic Poa supersect. Homalopoa, which is suggested to have originated in Eurasia 8·4–4·2 million years ago (Mya). The ancestor of Poa supersect. Homalopoa radiated throughout the New World in the Late Miocene–Early Pliocene, with major lineages originating during the Pliocene to Pleistocene (5–2 Mya). Breeding systems are linked to geographic areas, showing an evolutionary pattern associated with different habitats. At least three major pathways from hermaphroditism to diclinism are inferred in New World Homalopoa: two leading to dioecism, one via gynodioecism in South America and another directly from hermaphroditism in North America, a result that needs to be checked with a broader sampling of diclinous species in North America. A third pathway leads from hermaphroditism to gynomonoecism in Andean species of South America, with strictly pistillate species evolving in the highest altitudes.Conclusions Divergence dating provides a temporal context to the evolution of breeding systems in New World Poa supersect. Homalopoa. The results are consistent with the infrageneric classification in part; monophyletic sections are confirmed, it is proposed to reclassify species of sect. Acutifoliae, Dasypoa and Homalopoa s.l. and it is acknowledged that revision of the infrageneric taxonomy of the gynomonoecious species is needed.     

An experimental study of the population and evolutionary dynamics of Vibrio cholerae O1 and the bacteriophage JSF4

Studies of Vibrio cholerae in the environment and infected patients suggest that the waning of cholera outbreaks is associated with rise in the density of lytic bacteriophage. In accordance with mathematical models, there are seemingly realistic conditions where phage predation could be responsible for declines in the incidence of cholera. Here, we present the results of experiments with the El Tor strain of V. cholerae (N16961) and a naturally occurring lytic phage (JSF4), exploring the validity of the main premise of this model: that phage predation limits the density of V. cholerae populations. At one level, the results of our experiments are inconsistent with this hypothesis. JSF4-resistant V. cholerae evolve within a short time following their confrontation with these viruses and their populations become limited by resources rather than phage predation. At a larger scale, however, the results of our experiments are not inconsistent with the hypothesis that bacteriophage modulate outbreaks of cholera. We postulate that the resistant bacteria that evolved play an insignificant role in the ecology or pathogenicity of V. cholerae. Relative to the phage-sensitive cells from whence they are derived, the evolved JSF4-resistant V. cholerae have fitness costs and other characters that are likely to impair their ability to compete with the sensitive cells in their natural habitat and may be avirulent in human hosts. The results of this in vitro study make predictions that can be tested in natural populations of V. cholerae and cholera-infected patients.     

Structural and thermodynamic effects of post-translational modifications in mutant and wild type Cu, Zn superoxide dismutase

Aggregation of Cu,Zn superoxide dismutase (SOD1) is implicated in amyotrophic lateral sclerosis. Glutathionylation and phosphorylation of SOD1 is omnipresent in the human body, even in healthy individuals, and has been shown to increase SOD1 dimer dissociation, which is the first step on the pathway toward SOD1 aggregation. We found that post-translational modification of SOD1, especially glutathionylation, promotes dimer dissociation. We discovered an intermediate state in the pathway to dissociation, a conformational change that involves a “loosening” of the β-barrels and a loss or shift of dimer interface interactions. In modified SOD1, this intermediate state is stabilized as compared to unmodified SOD1. The presence of post-translational modifications could explain the environmental factors involved in the speed of disease progression. Because post-translational modifications such as glutathionylation are often induced by oxidative stress, post-translational modification of SOD1 could be a factor in the occurrence of sporadic cases of amyotrophic lateral sclerosis, which represent 90% of all cases of the disease.     

Computational design and experimental testing of the fastest-folding β-sheet protein

One of the most important and elusive goals of molecular biology is the formulation of a detailed, atomic-level understanding of the process of protein folding. Fast-folding proteins with low free-energy barriers have proved to be particularly productive objects of investigation in this context, but the design of fast-folding proteins was previously driven largely by experiment. Dramatic advances in the attainable length of molecular dynamics simulations have allowed us to characterize in atomic-level detail the folding mechanism of the fast-folding all-β WW domain FiP35. In the work reported here, we applied the biophysical insights gained from these studies to computationally design an even faster-folding variant of FiP35 containing only naturally occurring amino acids. The increased stability and high folding rate predicted by our simulations were subsequently validated by temperature-jump experiments. The experimentally measured folding time was 4.3 μs at 80 °C—about three times faster than the fastest previously known protein with β-sheet content and in good agreement with our prediction. These results provide a compelling demonstration of the potential utility of very long molecular dynamics simulations in redesigning proteins well beyond their evolved stability and folding speed.     

Mechanism of DNA Recognition by the Restriction Enzyme EcoRV

EcoRV, a restriction enzyme in Escherichia coli, destroys invading foreign DNA by cleaving it at the center step of a GATATC sequence. In the EcoRV-cognate DNA crystallographic complex, a sharp kink of 50° has been found at the center base-pair step (TA). Here, we examine the interplay between the intrinsic propensity of the cognate sequence to kink and the induction by the enzyme by performing all-atom molecular dynamics simulations of EcoRV unbound and interacting with three DNA sequences: the cognate sequence, GATATC (TA); the non-cognate sequence, GAATTC (AT); and with the cognate sequence methylated on the first adenine GA_CH3TATC (TA-CH₃). In the unbound EcoRV, the cleft between the two C-terminal subdomains is found to be open. Binding to AT narrows the cleft and forms a partially bound state. However, the intrinsic bending propensity of AT is insufficient to allow tight binding. In contrast, the cognate TA sequence is easier to bend, allowing specific, high-occupancy hydrogen bonds to form in the complex. The absence of cleavage for this methylated sequence is found to arise from the loss of specific hydrogen bonds between the first adenine of the recognition sequence and Asn185. On the basis of the results, we suggest a three-step recognition mechanism. In the first step, EcoRV, in an open conformation, binds to the DNA at a random sequence and slides along it. In the second step, when the two outer base pairs, GAxxTC, are recognized, the R loops of the protein become more ordered, forming strong hydrogen-bonding interactions, resulting in a partially bound EcoRV-DNA complex. In the third step, the flexibility of the center base pair is probed, and in the case of the full cognate sequence the DNA bends, the complex strengthens and the protein and DNA interact more closely, allowing cleavage.     

Effects of Pratylenchus zeae and Quinisulcius acutus Alone and in Combination on Sorghum

Host-parasite relationships of Pratylenchus zeae and Quinisulcius acutus, alone or in combination, were studied on sorghum in the greenhouse and laboratory. Q. acutus at 1,000 or 5,000 nematodes per 15-cm-d pot and P. zeae at 500 nematodes per pot significantly suppressed plant height and fresh and oven dry shoot and root weights. A mixture of 1,000 Q. acutus and 500 P. zeae per pot resulted in greatest suppression of growth. Roots of plants inoculated with Q. acutus alone were reduced in number and size and showed lesions and discoloration. Reproduction of this nematode 42 days after inoculation was much greater in treatments of 100 or 1,000 than 5,000 nematodes. The population density of the two species at 6 weeks after inoculation was significantly less when combined than for each species alone. When the two species were combined, reproduction of P. zeae was greater than that of Q. acutus, but the final populations per gram of root weight were the same. Q. acutus fed ectoparasitically on epidermal cells of sorghum roots in the zone of elongation and differentiation when observed under in vitro conditions.