Use of the root contact concept,an empirical leaf conductance model and pressure-volume curves in simulating crop water relations

A simulation model “DanStress” was developed for studying the integrated effects of soil, crop and climatic conditions on water relations and water use of field grown cereal crops. The root zone was separated into 0.1 m deep layers of topsoil and subsoil. For each layer the water potential at the root surface was calculated by a single root model, and the uptake of water across the root was calculated by a root contact model. Crop transpiration was calculated by Monteith's combination equation for vapour flow. Crop conductance to water vapour transfer for use in Monteith's combination equation was scaled up from an empirical stomatal conductance model used on sunlit and shaded crop surfaces of different crop layers. In the model, transpirational water loss originates from root water uptake and changes in crop water storage. Crop water capacitance, used for describing the water storage, was derived from the slope of pressure-volume (PV) curves of the leaves. PV curves were also used for deriving crop water potential, osmotic potential, and turgor pressure. The model could simulate detailed diurnal soil-crop water relations during a 23-day-drying cycle with time steps of one hour. During the grain filling period in spring barley (Hordeum distichum L.), grown in a sandy soil in the field, measured and predicted values of leaf water and osmotic potential, RWC, and leaf stomatal conductance were compared. Good agreement was obtained between measured and predicted values at different soil water deficits and climatic conditions. In the field, measured and predicted volumetric soil water contents (θ) of topsoil and subsoil layers were also compared during a drying cycle. Predicted and measured θ-values as a function of soil water deficits were similar suggesting that the root contact model approach was valid. From the investigation we concluded: (I) a model, which takes the degree of contact between root surface and soil water into account, can be used in sandy soil for calculation of root water uptake, so that the root conductance during soil water depletion only varies by the degree of contact; (II) crop conductance, used for calculation of crop transpiration, can be scaled up from an empirical single leaf stomatal conductance model controlled by the level of leaf water potential and micrometeorological conditions; (III) PV curves are usable for describing crop water status including crop water storage.     

Responses of population structure and genomic diversity to climate change and fishing pressure in a pelagic fish

The responses of marine species to environmental changes and anthropogenic pressures (e.g., fishing) interact with ecological and evolutionary processes that are not well understood. Knowledge of changes in the distribution range and genetic diversity of species and their populations into the future is essential for the conservation and sustainable management of resources. Almaco jack (Seriola rivoliana) is a pelagic fish with high importance to fisheries and aquaculture in the Pacific Ocean. In this study, we assessed contemporary genomic diversity and structure in loci that are putatively under selection (outlier loci) and determined their potential functions. Using a combination of genotype–environment association, spatial distribution models, and demogenetic simulations, we modeled the effects of climate change (under three different RCP scenarios) and fishing pressure on the species' geographic distribution and genomic diversity and structure to 2050 and 2100. Our results show that most of the outlier loci identified were related to biological and metabolic processes that may be associated with temperature and salinity. The contemporary genomic structure showed three populations—two in the Eastern Pacific (Cabo San Lucas and Eastern Pacific) and one in the Central Pacific (Hawaii). Future projections suggest a loss of suitable habitat and potential range contractions for most scenarios, while fishing pressure decreased population connectivity. Our results suggest that future climate change scenarios and fishing pressure will affect the genomic structure and genotypic composition of S. rivoliana and lead to loss of genomic diversity in populations distributed in the eastern-central Pacific Ocean, which could have profound effects on fisheries that depend on this resource.     

Helix-capping interaction in λ cro protein: A free energy simulation analysis

The stability mutant Tyr-26 → Asp was studied in the Cro protein from bacteriophage λ using free energy molecular dynamics simulations. The mutant was calculated to be more stable than the wild type by 3.0 ± 1.7 kcal/mol/monomer, in reasonable agreement with experiment (1.4 kcal/mol/monomer). Moreover, the aspartic acid in the mutant was found to form a capping interation with the amino terminus of the third α-helix of Cro. The simulations were analyzed to understand better the source of the stability of this helix-capping interaction and to examine the results in light of previous explanations of stabilizing helix caps-namely, a model of local unsatisfied hydrogen bonds at the helix termini and the helix macro dipole model. Analysis of the simulations shows that the stabilizing effect of this charged helical cap is due both to favorable hydrogen bonds with backbone NH groups at the helix terminus and to favorable electrostatic interactions (but not hydrogen bonds) with their carbonyls (effectively the next row of local dipoles in the helix). However, electrostatic interactions are weak or negligible with backbone dipolar groups in the helix further away from the terminus. Moreover, the importance of other local electrostatic interactions with polar side chains near the helix terminus, which are neglected in most treatments of this effect, are shown to be important. Thus, the results support a model that is intermediate between the two previous explanations: both unsatisfied hydrogen bonds at the helix terminus and other, local preoriented dipolar groups stabilize the helix cap. These findings suggest that similar interactions with preoriented dipolar groups may be important for cooperativity in other charge–dipole interactions and may be employed to advantage for molecular design. © 1994 Wiley-Liss, Inc.     

Molecular dynamics simulations of hydrophobic collapse of ubiquitin.

Nine nonnative conformations of ubiquitin, generated during two different thermal denaturation trajectories, were simulated under nearly native conditions (62 degrees C). The simulations included all protein and solvent atoms explicitly, and simulation times ranged from 1-2.4 ns. The starting structures had alpha-carbon root-mean-square deviations (RMSDs) from the crystal structure of 4-12 A and radii of gyration as high as 1.3 times that of the native state. In all but one case, the protein collapsed when the temperature was lowered and sampled conformations as compact as those reached in a control simulation beginning from the crystal structure. In contrast, the protein did not collapse when simulated in a 60% methanol:water mixture. The behavior of the protein depended on the starting structure: during simulation of the most native-like starting structures (<5 A RMSD to the crystal structure) the RMSD decreased, the number of native hydrogen bonds increased, and the secondary and tertiary structure increased. Intermediate starting structures (5-10 A RMSD) collapsed to the radius of gyration of the control simulation, hydrophobic residues were preferentially buried, and the protein acquired some native contacts. However, the protein did not refold. The least native starting structures (10-12 A RMSD) did not collapse as completely as the more native-like structures; instead, they experienced large fluctuations in radius of gyration and went through cycles of expansion and collapse, with improved burial of hydrophobic residues in successive collapsed states.     

Molecular Modelling Studies on the Catalytic Mechanism of Candida Rugosa Lipase

Quantum chemical and molecular dynamics investigations have been performed on model systems for Candida rugosa lipase (CRL) to study mechanistic and conformational features of the catalytic hydrolysis. Based on X-ray data, a simplified model of the CRL substrate complex was created for the PM3 and ab initio calculations, including the amino acid residues both of the catalytic triad and the oxyanion hole.The energetic and structural properties of significant species along the pathway of the hydrolysis of the model substrate acetic acid methyl ester have been calculated. By modifications of the residues of the oxyanion hole as well as the catalytic triad, the influence of these parts of the active site on the pathway of the reaction was analysed in more detail.Moreover, molecular dynamics simulations have been carried out on CRL adducts with (±)-cis-4-acetamido-cyclopent-2-ene-1-carboxylic esters with different lengths of their alkyl chain and their absolute configuration as substrates. For the MD simulations using the AMBER program, all amino acid residues and water molecules with a cut-off radius less than 1500 pm from the substrate were taken into account. From the analysis of the trajectories and histograms for significant hydrogen bonds in the active site of the enzyme adducts, some hints were obtained for the enantiodifferentiation and the chain dependence of the esters in catalytic hydrolysis by CRL.Electronic Supplementary Material available.     

Scrutiny of the mechanism of small molecule inhibitor preventing conformational transition of amyloid-β42 monomer: insights from molecular dynamics simulations

Alzheimer’s disease (AD) is a progressive neurodegenerative disorder that is characterized by loss of intellectual functioning of brain and memory loss. According to amyloid cascade hypothesis, aggregation of amyloid-β₄₂ (Aβ₄₂) peptide can generate toxic oligomers and their accumulation in the brain is responsible for the onset of AD. In spite of carrying out a large number of experimental studies on inhibition of Aβ₄₂ aggregation by small molecules, the detailed inhibitory mechanism remains elusive. In the present study, comparable molecular dynamics (MD) simulations were performed to elucidate the inhibitory mechanism of a sulfonamide inhibitor C1 (2,5-dichloro-N-(4-piperidinophenyl)-3-thiophenesulfonamide), reported for its in vitro and in vivo anti-aggregation activity against Aβ₄₂. MD simulations reveal that C1 stabilizes native α-helix conformation of Aβ₄₂ by interacting with key residues in the central helix region (13–26) with hydrogen bonds and π–π interactions. C1 lowers the solvent-accessible surface area of the central hydrophobic core (CHC), KLVFF (16–20), that confirms burial of hydrophobic residues leading to the dominance of helical conformation in the CHC region. The binding free energy analysis with MM–PBSA demonstrates that Ala2, Phe4, Tyr10, Gln15, Lys16, Leu17, Val18, Phe19, Phe20, Glu22, and Met35 contribute maximum to binding free energy (?43.1 kcal/mol) between C1 and Aβ₄₂ monomer. Overall, MD simulations reveal that C1 inhibits Aβ₄₂ aggregation by stabilizing native helical conformation and inhibiting the formation of aggregation-prone β-sheet conformation. The present results will shed light on the underlying inhibitory mechanism of small molecules that show potential in vitro anti-aggregation activity against Aβ₄₂.     

3D-QSAR,molecular docking,and molecular dynamic simulations for prediction of new Hsp90 inhibitors based on isoxazole scaffold

Heat shock protein 90(Hsp90), as a molecular chaperone, play a crucial role in folding and proper function of many proteins. Hsp90 inhibitors containing isoxazole scaffold are currently being used in the treatment of cancer as tumor suppressers. Here in the present studies, new compounds based on isoxazole scaffold were predicted using a combination of molecular modeling techniques including three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamic (MD) simulations. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were also done. The steric and electrostatic contour map of CoMFA and CoMSIA were created. Hydrophobic, hydrogen bond donor and acceptor of CoMSIA model also were generated, and new compounds were predicted by CoMFA and CoMSIA contour maps. To investigate the binding modes of the predicted compounds in the active site of Hsp90, a molecular docking simulation was carried out. MD simulations were also conducted to evaluate the obtained results on the best predicted compound and the best reported Hsp90 inhibitors in the 3D-QSAR model. Findings indicate that the predicted ligands were stable in the active site of Hsp90.     

Zika virus NS5 protein potential inhibitors: an enhanced in silico approach in drug discovery

The re-emerging Zika virus (ZIKV) is an arthropod-borne virus that has been described to have explosive potential as a worldwide pandemic. The initial transmission of the virus was through a mosquito vector, however, evolving modes of transmission has allowed the spread of the disease over continents. The virus has already been linked to irreversible chronic central nervous system conditions. The concerns of the scientific and clinical community are the consequences of Zika viral mutations, thus suggesting the urgent need for viral inhibitors. There have been large strides in vaccine development against the virus but there are still no FDA approved drugs available. Rapid rational drug design and discovery research is fundamental in the production of potent inhibitors against the virus that will not just mask the virus, but destroy it completely. In silico drug design allows for this prompt screening of potential leads, thus decreasing the consumption of precious time and resources. This study demonstrates an optimized and proven screening technique in the discovery of two potential small molecule inhibitors of ZIKV Methyltransferase and RNA dependent RNA polymerase. This in silico ‘per-residue energy decomposition pharmacophore’ virtual screening approach will be critical in aiding scientists in the discovery of not only effective inhibitors of Zika viral targets, but also a wide range of anti-viral agents.     

Structure,stability and water permeation of ([D-Leu-L-Lys-(D-Gln-L-Ala)3]) cyclic peptide nanotube: a molecular dynamics study

The structural stability of 8 × ([D-Leu-L-Lys-(D-Gln-L-Ala)₃]) cyclic peptide nanotube (CPN) in water and different phospholipid bilayers were explored by 100 ns independent molecular dynamics (MD) simulations. The role of non-bonded interaction energy between the side and main chains of cyclic peptide rings in different membrane environments assessed, wherein the repulsive electrostatic interaction energy between neighbouring cyclic peptide rings was found adequate to break hydrogen bond energy thereby to crumple CPN. Further, the water permeation across the CPN channel was studied in four types of phospholipid bilayers- DMPG (1,2-Dimyristoyl-sn-glycero-3-phosphorylglycerol), DMPS (1,2-Dimyristoyl-sn-glycero-3-phosphoserine), POPC (1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and POPE (1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine) from MD simulations. DMPS membrane shows higher non-bonded interaction energies (?1913.06 kJ/mol of electrostatic interaction energy and ?994.13 kJ/mol of van der Waals interaction energy) with CPN due to the presence of polar molecules in lipid structure. Thusly, the non-bonded interaction energies were essential towards the stability of CPN than hydrogen bonds between the nearby cyclic peptides. The result also reveals the role of side chains, hydrogen bonds and non-bonded interaction energies in an aqueous environment. The diffusion coefficient of water obtained from means square deviation calculation shows similar coefficients irrespective of the lipid surroundings. However, the permeation coefficients demonstrate water flow in the channel relies upon the environment.     

Curcumin derivatives and Aβ-fibrillar aggregates: An interactions’ study for diagnostic/therapeutic purposes in neurodegenerative diseases

Several neurodegenerative diseases, like Alzheimer’s (AD), are characterized by amyloid fibrillar deposition of misfolded proteins, and this feature can be exploited for both diagnosis and therapy design. In this paper, structural modifications of curcumin scaffold were examined in order to improve its bioavailability and stability in physiological conditions, as well as its ability to interfere with β-amyloid fibrils and aggregates. The acid-base behaviour of curcumin derivatives, their pharmacokinetic stability in physiological conditions, and in vitro ability to interfere with Aβ fibrils at different incubation time were investigated. The mechanisms governing these phenomena have been studied at atomic level by means of molecular docking and dynamic simulations. Finally, biological activity of selected curcuminoids has been investigated in vitro to evaluate their safety and efficiency in oxidative stress protection on hippocampal HT-22 mouse cells.Two aromatic rings, π-conjugated structure and H-donor/acceptor substituents on the aromatic rings showed to be the sine qua non structural features to provide interaction and disaggregation activity even at very low incubation time (2h). Computational simulations proved that upon binding the ligands modify the conformational dynamics and/or interact with the amyloidogenic region of the protofibril facilitating disaggregation. Significantly, in vitro results on hippocampal cells pointed out protection against glutamate toxicity and safety when administered at low concentrations (1?μM). On the overall, in view of its higher stability in physiological conditions with respect to curcumin, of his rapid binding to fibrillar aggregates and strong depolymerizing activity, phtalimmide derivative K2F21 appeared a good candidate for both AD diagnostic and therapeutic purposes.