Search

Ions are involved in multiple biological processes and may exist bound to biomolecules or may be associated with their surface. Although the presence of ions in nucleic acids has traditionally gained more interest, ion-protein interactions, often with a marked dependency on pH, are beginning to gather attention. Here we present a detailed analysis on the binding and distribution of ions around β-lactoglobulin using a constant-pH MD (CpHMD) method, at a pH range 3-8, and compare it with the more traditional Poisson-Boltzmann (PB) model and the existing experimental data. Most analyses used ion concentration maps built around the protein, obtained from either the CpHMD simulations or PB calculations. The requirements of approximate charge neutrality and ionic strength equal to bulk, imposed on the MD box, imply that the absolute value of the ion excess should be half the protein charge, which is in agreement with experimental observation on other proteins (

Biometric measurements in the context of myopia are fundamental to detect eyes at risk of developing myopia and during the follow-up of patients with myopia control treatment. Thus, the accuracy of biometers has high clinical relevance.The Myopia Master is a new biometer based on partial coherence interferometry especially dedicated to the follow-up of myopic patients. This study aims to assess the repeatability of the Myopia Master and evaluate its agreement with a swept-source optical coherence interferometry biometer (IOL Master 700).This cross-sectional prospective study assessed the biometric parameters of two groups of myopes (age range: 8-16 years old), spectacle corrected (n = 60) and orthokeratology contact lens wearers (n = 60). One senior optometrist performed two consecutive measurements per instrument, which included axial length, mean keratometry and horizontal visible iris diameter (HVID). The repeatability of each device and the agreement between devices were assessed by the dispersion of the measurement differences, for AL, mean keratometry, corneal astigmatism and HVID.The two biometers measured approximately the same value in both measurements. Test-retest repeatability tended to be lower than clinical significant thresholds, in particular, for AL and mean keratometry. Corneal-related parameters tended to have lower repeatability in the orthokeratology group, especially mean keratometry. The agreement between instruments revealed statistically significant differences between devices with the SS-OCT measuring longer eyes, steeper corneas and larger HVID.In a paediatric population, the Myopia Master showed clinically acceptable repeatability levels, but the IOL Master 700 demonstrated superior repeatability. Eyes treated with orthokeratology may compromise the repeatability of the corneal-related parameters. The Myopia Master and the IOL Master 700 are repeatable devices appropriate for monitoring myopia progression, but the differences observed do not allow their use interchangeably.

Protein-protein association is often mediated by electrostatic interactions and modulated by pH. However, experimental and computational studies have often overlooked the effect of association on the protonation state of the protein. In this work, we present a methodological approach based on constant-pH molecular dynamics (MD), which aims to provide a detailed description of a pH-dependent protein-protein association, and apply it to the dimerization of β-lactoglobulin (BLG). A selection of analyses is performed using the data generated by constant-pH MD simulations of monomeric and dimeric forms of bovine BLG, in the pH range 3-8. First, we estimate free energies of dimerization using a computationally inexpensive approach based on the Wyman-Tanford linkage theory, calculated in a new way through the use of thermodynamically based splines. The individual free energy contribution of each titratable site is also calculated, allowing for identification of relevant residues. Second, the correlations between the proton occupancies of pairs of sites are calculated (using the Pearson coefficient), and extensive networks of correlated sites are observed at acidic pH values, sometimes involving distant pairs. In general, strongly correlated sites are also slow proton exchangers and contribute significantly to the pH-dependency of the dimerization free energy. Third, we use ionic density as a fingerprint of protein charge distribution and observe electrostatic complementarity between the monomer faces that form the dimer interface, more markedly at the isoionic point (where maximum dimerization occurs) than at other pH values, which might contribute to guide the association. Finally, the pH-dependent dimerization modes are inspected using PCA, among other analyses, and two states are identified: a

Under strong sunlight, plants avoid photooxidation by quenching the excess absorbed energy. Quenching is triggered by PsbS, a membrane protein that is activated and deactivated by the light-dependent pH changes in the thylakoid lumen. The mechanism of action of this protein is unknown, but it was suggested that several glutamates act as pH sensors. However, the pK a of glutamate is several pH units below the physiological values in the lumen. Thus, how can PsbS sense the pH of the lumen, and how does it respond to it? By applying a nonstandard molecular dynamics method that treats pH explicitly, we show that the lumen-exposed glutamates of PsbS have strongly shifted pK a values and that such shifts are crucial for the pH sensitivity in physiological conditions. We also demonstrate that protonation drives a systematic unfolding of a region key for protein-protein interactions, indicating that PsbS response to pH is a functional conformational switch.

Bridge falsework systems are one of the most common temporary structures used in the construction industry to support formwork during the construction of cast-in-place concrete bridges. This paper presents the results of advanced numerical studies of selected structural systems using a novel joint finite element and information gathered from an extensive experimental campaign of various types of joints commonly found in Cuplok falsework structures. Different hazardous scenarios identified as being critical to the structural performance of bridge falsework systems are analysed, such as ground settlements, bracing configurations and falsework systems using steel beam girders. For each considered hazardous scenario, structural behaviour and resistance are analysed and discussed. From the results, relevant practical information was obtained that can be used to reduce the risks associated with bridge falsework systems. For example, it was found that even a small value of isolated differential ground settlements could reduce the resistance of the system by more than 10%. In addition, for falsework systems using steel beam girders, it was found that inappropriate bracing of falsework towers could lead to a 50% reduction in system resistance.

© The Author(s) 2020. Open Access. This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/., The influenza virus fusion process, whereby the virus fuses its envelope with the host endosome membrane to release the genetic material, takes place in the acidic late endosome environment. Acidification triggers a large conformational change in the fusion protein, hemagglutinin (HA), which enables the insertion of the N-terminal region of the HA2 subunit, known as the fusion peptide, into the membrane of the host endosome. However, the mechanism by which pH modulates the molecular properties of the fusion peptide remains unclear. To answer this question, we performed the first constant-pH molecular dynamics simulations of the influenza fusion peptide in a membrane, extending for 40 µs of aggregated time. The simulations were combined with spectroscopic data, which showed that the peptide is twofold more active in promoting lipid mixing of model membranes at pH 5 than at pH 7.4. The realistic treatment of protonation introduced by the constant-pH molecular dynamics simulations revealed that low pH stabilizes a vertical membrane-spanning conformation and leads to more frequent contacts between the fusion peptide and the lipid headgroups, which may explain the increase in activity. The study also revealed that the N-terminal region is determinant for the peptide's effect on the membrane., This work was financially supported by FCT—Fundação para a Ciência e a Tecnologia, Portugal, through projects PTDC/QUI-BIQ/114774/2009, PTDC/CCI-BIO/28200/2017 and Pest-OE/EQB/LA0004/2011. This work was also financially supported by Project LISBOA-01-0145-FEDER-007660 (Microbiologia Molecular, Estrutural e Celular) funded by FEDER funds through COMPETE2020—Programa Operacional Competitividade e Internacionalização (POCI) and by national funds through FCT—Fundação para a Ciência e a Tecnologia. DL was supported by FCT post-doc fellowship SFRH/BPD/92537/2013.

Bridge falsework systems are frequently used during the construction of cast-in-place concrete bridges. These structures have a significant impact on the cost, construction rate and safety of the supported permanent structures. However, in recent years a high number of accidents involving bridge falsework systems have occurred worldwide. This paper concerns experimental tests on joints of falsework Cuplok systems. In this paper, the results of 192 tests carried out during the experimental programme are presented, consisting of five different experiments on three different types of joints. Relevant new results are obtained concerning the bending behaviour and resistance of spigot joints (used to connect two vertical elements) and forkhead joints (used to connect the falsework system to the formwork system), relating to which no published research has been found to date. Based on the results, joint models are evaluated and improved empirical analytical models are developed. Finally, probabilistic models for the most relevant variables controlling the behaviour and resistance of the Cuplok joints (used to connect horizontal elements to vertical elements) are presented that can be used in the probabilistic analysis of falsework Cuplok systems.

The Columbia River (CR) estuary is characterized by high river discharge and strong tides that generate high velocity flows and sharp density gradients. Its dynamics strongly affects the coastal ocean circulation. Tidal straining in turn modulates the stratification in the estuary. Simulating the hydrodynamics of the CR estuary and plume therefore requires a multi-scale model as both shelf and estuarine circulations are coupled. Such a model has to keep numerical dissipation as low as possible in order to correctly represent the plume propagation and the salinity intrusion in the estuary. Here, we show that the 3D baroclinic discontinuous Galerkin finite element model SLIM 3D is able to reproduce the main features of the CR estuary-to-ocean continuum. We introduce new vertical discretization and mode splitting that allow us to model a region characterized by complex bathymetry and sharp density and velocity gradients. Our model takes into account the major forcings, i.e. tides, surface wind stress and river discharge, on a single multi-scale grid. The simulation period covers the end of spring-early summer of 2006, a period of high river flow and strong changes in the wind regime. SLIM 3D is validated with in-situ data on the shelf and at multiple locations in the estuary and compared with an operational implementation of SELFE. The model skill in the estuary and on the shelf indicate that SLIM 3D is able to reproduce the key processes driving the river plume dynamics, such as the occurrence of bidirectional plumes or reversals of the inner shelf coastal currents.

This paper outlines strategies that would advance coastal ocean modeling, analysis and prediction as a complement to the observing and data management activities of the coastal components of the U.S. Integrated Ocean Observing System (IOOS��) and the Global Ocean Observing System (GOOS). The views presented are the consensus of a group of U.S. based researchers with a cross-section of coastal oceanography and ocean modeling expertise and community representation drawn from Regional and U.S. Federal partners in IOOS. Priorities for research and development are suggested that would enhance the value of IOOS observations through model-based synthesis, deliver better model-based information products, and assist the design, evaluation and operation of the observing system itself. The proposed priorities are: model coupling, data assimilation, nearshore processes, cyberinfrastructure and model skill assessment, modeling for observing system design, evaluation and operation, ensemble prediction, and fast predictors. Approaches are suggested to accomplish substantial progress in a 3-8 year timeframe. In addition, the group proposes steps to promote collaboration between research and operations groups in Regional Associations, U.S. Federal Agencies, and the international ocean research community in general that would foster coordination on scientific and technical issues, and strengthen federal-academic partnerships benefiting IOOS stakeholders and end users.

Bridge falsework systems are one of the most common temporary structures used to support the formwork during the construction of concrete bridges. In this paper, the development of a risk analysis ...

A sediment model coupled to the hydrodynamic model SELFE is validated against a benchmark combining a set of idealized tests and an application to a field-data rich energetic estuary. After sensitivity studies, model results for the idealized tests largely agree with previously reported results from other models in addition to analytical, semi-analytical, or laboratory results. Results of suspended sediment in an open channel test with fixed bottom are sensitive to turbulence closure and treatment for hydrodynamic bottom boundary. Results for the migration of a trench are very sensitive to critical stress and erosion rate, but largely insensitive to turbulence closure. The model is able to qualitatively represent sediment dynamics associated with estuarine turbidity maxima in an idealized estuary. Applied to the Columbia River estuary, the model qualitatively captures sediment dynamics observed by fixed stations and shipborne profiles. Representation of the vertical structure of suspended sediment degrades when stratification is underpredicted. Across all tests, skill metrics of suspended sediments lag those of hydrodynamics even when qualitatively representing dynamics. The benchmark is fully documented in an openly available repository to encourage unambiguous comparisons against other models.

Solution pH is a physicochemical property that has a key role in cellular regulation, and its impact at the molecular level is often difficult to study by experimental methods. In this context, several theoretical methods were developed to study pH effects in macromolecules. The stochastic titration constant-pH molecular dynamics method (CpHMD) was developed by coupling molecular sampling methods, which are appropriate to study the conformational ensemble of biomolecules, with continuum electrostatics approaches, which properly describe pH-dependent protonation states. However, in difficult cases, the protonation sampling can be too slow for the commonly accessible computational times. In this work, we combined a pH replica exchange scheme with this CpHMD method and explored several optimization strategies and possible limitations.

The present work aims to achieve an additional insight into the protein partitioning behavior in aqueous two phase systems (ATPSs), together with a study on the viability of a semi-empirical model based on continuum electrostatics to predict the protein partition characteristics. The partitioning behaviors of 14 globular proteins, with different properties, were explored in three polymer/polymer ATPSs. By the Collander equation, a linear correlation between protein partitioning coefficients in all systems was observed. Using the semi-empirical model it was possible to predict the partitioning behavior of proteins. The electrostatic energy depends on the protein size and ATPSs characteristics and varies in agreement with the difference in phase dielectric constants. Linear correlation of nonpolar energy, and the solvent accessible surface area was observed. Polymer structure and concentration have a significant influence on model viability. A good qualitative prediction of preferred phase for studied proteins was observed.

Dendrimers are a family of ramified synthetic molecules. pH effects and electrostatic interactions are known to be crucial players to explain the conformational and functional behaviors observed in these systems. Nonetheless, to date, no computational study involving these systems has explicitly addressed the protonation equilibrium taking place at different pH values for dendrimers containing multiple ionizable sites. Herein, we present the results of constant-pH molecular dynamics simulations performed at several pH values for four peptide dendrimers of different generations (from one to four) composed of the same type of amino acids: histidines, serines, and diaminopropionic acid. These dendrimers are known to catalyze the hydrolysis of pyrene sulfonate esters. Constant-pH MD simulations in the presence of substrate molecules at the optimum pH for catalysis are also reported. The results show that first and second generation dendrimers are almost structurally unresponsive to pH variations. For third and fourth generation dendrimers, pH plays a structuring role, with markedly different behaviors being observed when passing from acidic to neutral pH. Protonation-conformation coupling effects influence several intramolecular interactions, which, in turn, modulate the shape and structure at the different pH values. The atypical and highly pH-dependent protonation profiles of some histidine residues are also investigated. The interactions between dendrimers and substrates restrict the conformational space available to the dendrimers and enforce conformational homogeneity. This structuring effect is a consequence of the dendrimer-substrate interactions which occur through stabilizing hydrogen bonds and ion pairs between the substrate sulfonate groups and the dendrimer residues. Our results provide original fundamental data contributing to the development of novel pH-modulated dendritic systems and the improvement of the existing ones.

Ocean currents or temperature may substantially influence migration behavior in many marine species. However, high-resolution data on animal movement in the marine environment are scarce; therefore, analysts and managers must typically rely on unvalidated assumptions regarding movement, behavior, and habitat use. We used a spatially explicit, individual-based model of early marine migration with two stocks of yearling Chinook salmon to quantify the influence of external forces on estimates of swim speed, consumption, and growth. Model results suggest that salmon behaviorally compensate for changes in the strength and direction of ocean currents. These compensations can result in salmon swimming several times farther than their net movement (straight-line distance) would indicate. However, the magnitude of discrepancy between compensated and straight-line distances varied between oceanographic models. Nevertheless, estimates of relative swim speed among fish groups were less sensitive to the choice of model than estimates of absolute individual swim speed. By comparing groups of fish, this tool can be applied to management questions, such as how experiences and behavior may differ between groups of hatchery fish released early vs. later in the season. By taking into account the experiences and behavior of individual fish, as well as the influence of physical ocean processes, our approach helps illuminate the “black box” of juvenile salmon behavior in the early marine phase of the life cycle.

Cytochrome c oxidases (CcOs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes reduce dioxygen (O(2)) to water and, simultaneously, generate a transmembrane electrochemical proton gradient. Despite their importance in the aerobic metabolism and the large amount of structural and biochemical data available for the A1-type CcO family, there is still no consensually accepted description of the molecular mechanisms operating in this protein. A substantial number of questions about the CcO's working mechanism remain to be answered, including how the protonation behavior of some key residues is modulated during a reduction cycle and how is the conformation of the protein affected by protonation. The main objective of this work was to study the protonation-conformation coupling in CcOs and identify the molecular factors that control the protonation state of some key residues. In order to directly capture the interplay between protonation and conformational effects, we have performed constant-pH MD simulations of an A1-type CcO inserted into a lipid bilayer in two redox states (oxidized and reduced) at physiological pH. From the simulations, we were able to identify several groups with unusual titration behavior that are highly dependent on the protein redox state, including the A-propionate from heme a and the D-propionate from heme a3, two key groups possibly involved in proton pumping. The protonation state of these two groups is heavily influenced by subtle conformational changes in the protein (notably of R481(I) and R482(I)) and by small changes in the hydrogen bond network.

We propose a novel method for constructing a pseudo-synoptic view of estuarine features from non-synoptic observations captured by mobile platforms. The model-aided Lagrangian interpretation (MALI) method is based on relocating observations to a common reference moment in time along three-dimensional Lagrangian trajectories derived from a numerical model of estuarine circulation. The method relies on the model skill to capture large-scale circulation features, and on high-resolution in situ observations to characterize small-scale hydrographic structure. We demonstrate our technique by applying MALI to autonomous underwater vehicle observations in the Columbia River estuary, with the aid of a validated unstructured-grid finite-element numerical simulation. The method can be readily adapted to a broader range of environments, observational platforms, and model-data combinations.

Juvenile Chinook salmon population dynamics in the Columbia River estuary are influenced by physical processes, hatchery practices, and behavioral decision-making. To better understand how environmental forcing and swimming behavior influence estuarine migration and travel times, we developed an individual-based model (IBM) that uses 3-D outputs from a hydrodynamic model to simulate Lagrangian transport as well as swimming and bioenergetics sub-models to simulate active swimming and growth. Simulations were run in 2010 during the migration seasons for yearling and subyearling Chinook salmon. For both life history types, alternative behaviors were simulated, from random walks to behaviors that optimized efficient system migration for yearling Chinook salmon and growth for subyearling Chinook salmon. Simulation results compared well against observed data on travel times and common migration pathways; the simulated travel times for both yearling and subyearling Chinook salmon were within several hours of the observed travel times. In general, residence times and pathways were largely driven by river discharge and the phase of the tide. During periods of greater river discharge, simulated estuarine residence times were reduced and variability across individuals was minimal. The timing of estuarine exit was closely tied to the phase of the tide, with most simulated individuals exiting the system during the ebb phase. While travel times were largely driven by flow velocities, swimming behavior was likewise important. Simulated yearling Chinook salmon behaviors that optimized movement with surrounding flows resulted in reduced estuarine residence times when compared to passive and random walk behaviors. Similarly, simulated subyearling Chinook salmon behaviors that optimized growth directed individuals to shallow peripheral habitats, resulting in longer residence times and higher growth rates. Even if potentially important factors such as predator avoidance were not included, this IBM provides an informative tool to model migration pathways, growth, and residence times of juvenile salmon in an estuarine environment and could be used to inform management decisions by evaluating various scenarios.

Structural robustness has become an important design variable. However, based on the existing definitions of structural robustness it is often difficult to analyse and evaluate structural robustness, and sometimes not efficient since they mix structural robustness with several other structural variables. This paper concerns the development of a new structural robustness definition, and structural robustness and structural fragility indices. The basis for the development of the new indices is the analysis of the damage energy of structural systems for a given hazard scenario and involves a criterion to define an "unavoidable collapse" state. Illustrative examples are given detailing the steps and calculations needed to obtain values for both the structural robustness and the structural fragility indices. Finally, this paper presents the main advantages of the newly proposed definition and indices for the structural risk analysis over existing traditional methods.

Current constant-pH molecular dynamics (CpHMD) simulations provide a proper treatment of pH effects on the structure and dynamics of soluble biomolecules like peptides and proteins. However, addressing such effects on lipid membrane assemblies has remained problematic until now, despite the important role played by lipid ionization at physiological pH in a plethora of biological processes. Modeling (de)protonation events in these systems requires a proper consideration of the physicochemical features of the membrane environment, including a sound treatment of solution ions. Here, we apply our recent CpHMD-L method to the study of pH effects on a 25% DMPA/DMPC bilayer membrane model, closely reproducing the correct lipid phases of this system, namely, gel-fluid coexistence at pH 4 and a fluid phase at pH 7. A significant transition is observed for the membrane ionization and mechanical properties at physiological pH, providing a molecular basis for the well-established role of phosphatidic acid (PA) as a key player in the regulation of many cellular events. Also, as reported experimentally, we observed pH-induced PA-PA lipid aggregation at acidic pH. By including the titration of anionic phospholipids, the current methodology makes possible to simulate lipid bilayers with increased realism. To the best of our knowledge, this is the first simulation study dealing with a continuous phospholipid bilayer with pH titration of all constituent lipids.

The Columbia River has a dynamic and fast flushing estuary impacted by strong advection and mixing of riverine and oceanic waters, and high but variable loads of suspended particulate matter. Transient, but recurring water and nutrient fluxes from end-members impart strong spatial and temporal gradients, contributing to microbiological hotspots that play important ecological roles in the estuary. Investigations of corresponding microbiota require precisely timed samples that are contextualized by physical and biogeochemical data. To accomplish this, we embedded a robotic microbial sampler (Environmental Sample Processor, ESP) within the operations of an interdisciplinary observation and prediction system (Science and Technology University Research Network, SATURN; www.stccmop.org/saturn). Autonomous, adaptively sampled water collection by the ESP was implemented based on environmental conditions assessed from SATURN physical and biogeochemical sensors. Water was pumped from multiple depths to sensors and the ESP on dry land. If water met user-defined parameters, ESP sampling was automatically initiated. This strategy was tested during three deployments in summer 2013, during which operational tools for analysis and visualization were used to formulate well-constrained mission plans by providing estimates of the intensity and timing of oxygen-depleted ocean water intrusion and estuarine turbidity maxima. This allowed the effective characterization of the impact of these events on selected estuarine microbiota.

Pulmonary surfactant protein C (SP-C) is a small peptide with two covalently linked fatty acyl chains that plays a crucial role in the formation and stabilization of the pulmonary surfactant reservoirs during the compression and expansion steps of the respiratory cycle. Although its function is known to be tightly related to its highly hydrophobic character and key interactions maintained with specific lipid components, much is left to understand about its molecular mechanism of action. Also, although it adopts a mainly helical structure while associated with the membrane, factors as pH variation and deacylation have been shown to affect its stability and function. In this work, the conformational behavior of both the acylated and deacylated SP-C isoforms was studied in a DPPC bilayer under different pH conditions using constant-pH molecular dynamics simulations. Our findings show that both protein isoforms are remarkably stable over the studied pH range, even though the acylated isoform exhibits a labile helix-turn-helix motif rarely observed in the other isoform. We estimate similar tilt angles for the two isoforms over the studied pH range, with a generally higher degree of internalization of the basic N-terminal residues in the deacylated case, and observe and discuss some protonation-conformation coupling effects. Both isoforms establish contacts with the surrounding lipid molecules (preferentially with the sn-2 ester bonds) and have a local effect on the conformational behavior of the surrounding lipid molecules, the latter being more pronounced for acylated SP-C.

Oleic acid is a simple molecule with an aliphatic chain and a carboxylic group whose ionization and, consequently, intermolecular interactions are strongly dependent on the solution pH. The titration curve of these molecules was already obtained using different experimental methods, which have shown the lipid bilayer assemblies to be stable between pH 7.0 and 9.0. In this work, we take advantage of our recent implementations of periodic boundary conditions in Poisson-Boltzmann calculations and ionic strength treatment in simulations of charged lipid bilayers, and we studied the ionization dependent behavior of an oleic acid bilayer using a new extension of the stochastic titration constant-pH MD method. With this new approach, we obtained titration curves that are in good agreement with the experimental data. Also, we were able to estimate the slope of the titration curve from charge fluctuations, which is an important test of thermodynamic consistency for the sampling in a constant-pH MD method. The simulations were performed for ionizations up to 50%, because an experimentally observed macroscopic transition to micelles occurs above this value. As previously seen for a binary mixture of a zwitterionic and an anionic lipid, we were able to reproduce experimental results with simulation boxes usually far from neutrality. This observation further supports the idea that a charged membrane strongly influences the ion distribution in its vicinity and that neutrality is achieved significantly far from the bilayer surface. The good results obtained with this extension of the stochastic titration constant-pH MD method strongly supports its usefulness to sample the coupling between configuration and protonation in these types of biophysical systems. This method stands now as a powerful tool to study more realistic lipid bilayers where pH can influence both the lipids and the solutes interacting with them.