Mutagenesis Objective Search and Selection Tool (MOSST): an algorithm to predict structure-function related mutations in proteins

<p>Abstract</p> <p>Background</p> <p>Functionally relevant artificial or natural mutations are difficult to assess or predict if no structure-function information is available for a protein. This is especially important to correctly identify functionally significant non-synonymous single nucleotide polymorphisms (nsSNPs) or to design a site-directed mutagenesis strategy for a target protein. A new and powerful methodology is proposed to guide these two decision strategies, based only on conservation rules of physicochemical properties of amino acids extracted from a multiple alignment of a protein family where the target protein belongs, with no need of explicit structure-function relationships.</p> <p>Results</p> <p>A statistical analysis is performed over each amino acid position in the multiple protein alignment, based on different amino acid physical or chemical characteristics, including hydrophobicity, side-chain volume, charge and protein conformational parameters. The variances of each of these properties at each position are combined to obtain a global statistical indicator of the conservation degree of each property. Different types of physicochemical conservation are defined to characterize relevant and irrelevant positions. The differences between statistical variances are taken together as the basis of hypothesis tests at each position to search for functionally significant mutable sites and to identify specific mutagenesis targets. The outcome is used to statistically predict physicochemical consensus sequences based on different properties and to calculate the amino acid propensities at each position in a given protein. Hence, amino acid positions are identified that are putatively responsible for function, specificity, stability or binding interactions in a family of proteins. Once these key functional positions are identified, position-specific statistical distributions are applied to divide the 20 common protein amino acids in each position of the protein's primary sequence into a group of functionally non-disruptive amino acids and a second group of functionally deleterious amino acids.</p> <p>Conclusions</p> <p>With this approach, not only conserved amino acid positions in a protein family can be labeled as functionally relevant, but also non-conserved amino acid positions can be identified to have a physicochemically meaningful functional effect. These results become a discriminative tool in the selection and elaboration of rational mutagenesis strategies for the protein. They can also be used to predict if a given nsSNP, identified, for instance, in a genomic-scale analysis, can have a functional implication for a particular protein and which nsSNPs are most likely to be functionally silent for a protein. This analytical tool could be used to rapidly and automatically discard any irrelevant nsSNP and guide the research focus toward functionally significant mutations. Based on preliminary results and applications, this technique shows promising performance as a valuable bioinformatics tool to aid in the development of new protein variants and in the understanding of function-structure relationships in proteins.</p

Probing Molecular Mechanisms of the Hsp90 Chaperone: Biophysical Modeling Identifies Key Regulators of Functional Dynamics

Deciphering functional mechanisms of the Hsp90 chaperone machinery is an important objective in cancer biology aiming to facilitate discovery of targeted anti-cancer therapies. Despite significant advances in understanding structure and function of molecular chaperones, organizing molecular principles that control the relationship between conformational diversity and functional mechanisms of the Hsp90 activity lack a sufficient quantitative characterization. We combined molecular dynamics simulations, principal component analysis, the energy landscape model and structure-functional analysis of Hsp90 regulatory interactions to systematically investigate functional dynamics of the molecular chaperone. This approach has identified a network of conserved regions common to the Hsp90 chaperones that could play a universal role in coordinating functional dynamics, principal collective motions and allosteric signaling of Hsp90. We have found that these functional motifs may be utilized by the molecular chaperone machinery to act collectively as central regulators of Hsp90 dynamics and activity, including the inter-domain communications, control of ATP hydrolysis, and protein client binding. These findings have provided support to a long-standing assertion that allosteric regulation and catalysis may have emerged via common evolutionary routes. The interaction networks regulating functional motions of Hsp90 may be determined by the inherent structural architecture of the molecular chaperone. At the same time, the thermodynamics-based “conformational selection” of functional states is likely to be activated based on the nature of the binding partner. This mechanistic model of Hsp90 dynamics and function is consistent with the notion that allosteric networks orchestrating cooperative protein motions can be formed by evolutionary conserved and sparsely connected residue clusters. Hence, allosteric signaling through a small network of distantly connected residue clusters may be a rather general functional requirement encoded across molecular chaperones. The obtained insights may be useful in guiding discovery of allosteric Hsp90 inhibitors targeting protein interfaces with co-chaperones and protein binding clients

Interaction between contrasting rice genotypes and soil physical conditions induced by hydraulic stresses typical of alternate wetting and drying irrigation of soil

Background and aims: Alternate wetting and drying (AWD) saves water in paddy rice production but could influence soil physical conditions and root growth. This study investigated the interaction between contrasting rice genotypes, soil structure and mechanical impedance influenced by hydraulic stresses typical of AWD. Methods: Contrasting rice genotypes, IR64 and deeper- rooting Black Gora were grown in various soil conditions for 2 weeks. For the AWD treatments the soil was either maintained in a puddled state, equilibrated to −5 kPa (WET), or dried to −50 kPa and then rewetted at thewater potential of −5 kPa (DRY-WET). There was an additional manipulated macropore structure treatment, i.e. the soil was broken into aggregates, packed into cores and equilibrated to −5 kPa (REPACKED). A flooded treatment (puddled soil remained flooded until harvest) was set as a control (FLOODED). Soil bulk density, penetration resistance and X-ray Computed Tomography (CT) derived macropore structure were measured. Total root length, root surface area, root volume, average diameter, and tip number were determined by WinRhizo. Results: AWD induced formation of macropores and slightly increased soil mechanical impedance. The total root length of the AWD and REPACKED treatments were 1.7–2.2 and 3.5–4.2 times greater than that of the FLOODED treatment. There was no significant difference between WET and DRY-WET treatments. The differences between genotypes were minimal. Conclusions: AWD influenced soil physical properties and some root characteristics of rice seedlings, but drying soil initially to −50 kPa versus −5 kPa had no impact. Macropores formed intentionally from repacking caused a large change in root characteristics

Rapid and Intense Phosphate Desorption Kinetics When Saltwater Intrudes into Carbonate Rock

It is important to understand how phosphate sorption dynamics of coastal carbonate aquifers are affected by seawater intrusion, because many coastal aquifers are composed of carbonate rocks and subject to an increase in saltwater intrusion during relative sea-level rise. Twelve carbonate rock and unconsolidated sediment specimens were acquired from a test corehole spanning the full thickness of the Biscayne aquifer in southeastern Florida. All 12 samples exhibit low phosphorus content but variable contents of iron. Column leaching experiments were conducted with two carbonate aquifer samples, alternating between freshwater and saltwater flow. With the first influx of saltwater, phosphate concentration in leachate increased rapidly from a freshwater value of approximately 0.2 μM to peaks of between 0.8 and 1.6 μM. The phosphate concentration began to diminish as saltwater continued to flow, but sustained desorption continued for over 2 h. Overall, seawater drove sorption behavior much more than chemical composition for the aquifer rocks and sediment from the seven rock samples for which we did isotherm sorption experiments. Our results indicate that an immediate and intense pulse of phosphate desorption from carbonate rock and sediment with low phosphorus content occurs in response to an influx of seawater and that the duration of desorption will vary by layer within a single aquifer

Visioning the Future: Scenarios Modeling of the Florida Coastal Everglades

In this paper, we provide screening-level analysis of plausible Everglades ecosystem response by 2060 to sea level rise (0.50 m) interacting with macroclimate change (1.5 °C warming, 7% increase in evapotranspiration, and rainfall that either increases or decreases by 10%). We used these climate scenarios as input to the Ecological Landscape Model to simulate changes to seven interactive hydro-ecological metrics. Mangrove forest and other marine influences migrated up to 15 km inland in both scenarios, delineated by the saltwater front. Freshwater habitat area decreased by 25–30% under our two climate change scenarios and was largely replaced by mangroves and, in the increased rainfall scenario, open water as well. Significant mangroves drowned along northern Florida Bay in both climate change scenarios due to sea level rise. Increased rainfall of 10% provided significant benefits to the spatial and temporal salinity regime within the marine-influenced zone, providing a more gradual and natural adjustment for at-risk flora and fauna. However, increased rainfall also increased the risk of open water, due to water depths that inhibited mangrove establishment and reduced peat accumulation rates. We infer that ecological effects related to sea level rise may occur in the extreme front-edge of saltwater intrusion, that topography will control the incursion of this zone as sea level rises, and that differences in freshwater availability will have ecologically significant effects on ecosystem resilience through the temporal and spatial pattern of salinity changes