-
Fluctuation relations to calculate protein redox potentials from molecular dynamics simulations
Authors:
A. S. F. Oliveira,
J. Rubio,
C. E. M. Noble,
J. L. R. Anderson,
J. Anders,
A. J. Mulholland
Abstract:
The tunable design of protein redox potentials promises to open a range of applications in biotechnology and catalysis. Here we introduce a method to calculate redox potential changes by combining fluctuation relations with molecular dynamics simulations. It involves the simulation of reduced and oxidized states, followed by the instantaneous conversion between them. Energy differences introduced…
▽ More
The tunable design of protein redox potentials promises to open a range of applications in biotechnology and catalysis. Here we introduce a method to calculate redox potential changes by combining fluctuation relations with molecular dynamics simulations. It involves the simulation of reduced and oxidized states, followed by the instantaneous conversion between them. Energy differences introduced by the perturbations are obtained using the Kubo-Onsager approach. Using a detailed fluctuation relation coupled with Bayesian inference, these are post-processed into estimates for the redox potentials in an efficient manner. This new method, denoted MD+CB, is tested on a de novo four-helix bundle heme protein (the m4D2 `maquette') and five designed mutants, including some mutants characterized experimentally in this work. The MD+CB approach is found to perform reliably, giving redox potential shifts with reasonably good correlation (0.85) to the experimental values for the mutants. The MD+CB approach also compares well with redox potential shift predictions using a continuum electrostatic method. The estimation method employed within the MD+CB approach is straightforwardly transferable to standard equilibrium MD simulations, and holds promise for redox protein engineering and design applications.
△ Less
Submitted 11 January, 2024; v1 submitted 25 February, 2023;
originally announced February 2023.
-
Protein-Polymer Mixtures in the Colloid Limit: Aggregation, Sedimentation and Crystallization
Authors:
Rui Cheng,
Jingwen Li,
Ioatzin Ríos de Anda,
Thomas W. C. Taylor,
Malcolm A. Faers,
J. L. Ross Anderson,
Annela M. Seddon,
C. Patrick Royall
Abstract:
While proteins have been treated as particles with a spherically symmetric interaction, of course in reality the situation is rather more complex. A simple step towards higher complexity is to treat the proteins as non--spherical particles and that is the approach we pursue here. We investigate the phase behavior of enhanced green fluorescent protein (eGFP) under the addition of a non--adsorbing p…
▽ More
While proteins have been treated as particles with a spherically symmetric interaction, of course in reality the situation is rather more complex. A simple step towards higher complexity is to treat the proteins as non--spherical particles and that is the approach we pursue here. We investigate the phase behavior of enhanced green fluorescent protein (eGFP) under the addition of a non--adsorbing polymer, polyethylene glycol (PEG). From small angle x-ray scattering we infer that the eGFP undergoes dimerization and we treat the dimers as spherocylinders with aspect ratio $L/D-1 = 1.05$. Despite the complex nature of the proteins, we find that the phase behaviour is similar to that of hard spherocylinders with ideal polymer depletant, exhibiting aggregation and, in a small region of the phase diagram, crystallization. By comparing our measurements of the onset of aggregation with predictions for hard colloids and ideal polymers [S.V. Savenko and M. Dijkstra, J. Chem. Phys 124, 234902 (2006) and F. lo Verso et al., Phys. Rev. E 73, 061407 (2006)] we find good agreement, which suggests that the eGFP proteins are consistent with hard spherocylinders and ideal polymer.
△ Less
Submitted 16 June, 2021;
originally announced June 2021.
-
Massive viral replication and cytopathic effects in early COVID-19 pneumonia
Authors:
L. Schifanella,
J. L. Anderson,
M. Galli,
M. Corbellino,
A. Lai,
G. Wieking,
B. Grzywacz,
N. R. Klatt,
A. T. Haase,
T. W. Schacker
Abstract:
SARS-CoV-2 is the cause of COVID-19 acute respiratory illness that like its predecessors, MERS and SARS, can be severe and fatal 1-4. By April of 2020, COVID-19 infections had become a worldwide pandemic with nearly 3 million infections and over 200,000 deaths. The relative contributions of virus replication and cytopathic effects or immunopathological host responses to the severe and fatal outcom…
▽ More
SARS-CoV-2 is the cause of COVID-19 acute respiratory illness that like its predecessors, MERS and SARS, can be severe and fatal 1-4. By April of 2020, COVID-19 infections had become a worldwide pandemic with nearly 3 million infections and over 200,000 deaths. The relative contributions of virus replication and cytopathic effects or immunopathological host responses to the severe and fatal outcomes of COVID-19 lung infections have as yet to be determined. Here we show that SARS-CoV-2 replication and cytopathic effects in type II alveolar pneumocytes causes focal lung injury in an individual with no history of pulmonary symptoms. These findings point to the potential benefit of early effective antiviral treatment to prevent progression to severe and fatal COVID-19 pneumonia.
△ Less
Submitted 30 April, 2020;
originally announced May 2020.
-
Decorated Protein Networks: Functional Nanomaterials with Tunable Domain Size
Authors:
Ioatzin Rıos de Anda,
Angelique Coutable-Pennarun,
Chris Brasnett,
Stephen Whitelam,
Annela Seddon,
John Russo,
J. L. Ross Anderson,
C. Patrick Royall1
Abstract:
The implementation of natural and artificial proteins with designer properties and functionalities offers unparalleled opportunity for functional nanoarchitectures formed through self-assembly. However, to exploit the opportunities offered we require the ability to control protein assembly into the desired architecture while avoiding denaturation and therefore retaining protein functionality. Here…
▽ More
The implementation of natural and artificial proteins with designer properties and functionalities offers unparalleled opportunity for functional nanoarchitectures formed through self-assembly. However, to exploit the opportunities offered we require the ability to control protein assembly into the desired architecture while avoiding denaturation and therefore retaining protein functionality. Here we address this challenge with a model system of fluorescent proteins. Using techniques of self-assembly manipulation inspired by soft matter where interactions between components are controlled to yield the desired structure, we show that it is possible to assemble networks of proteins of one species which we can decorate with another, whose coverage we can tune. Consequently, the interfaces between domains of each component can also be tuned, with applications for example in energy transfer. Our model system of fluorescent proteins eGFP and mCherry retain their fluorescence throughout the assembly process, thus demonstrating that functionality is preserved.
△ Less
Submitted 13 November, 2019;
originally announced November 2019.
-
Ensemble Kalman filter data assimilation in a Babcock-Leighton solar dynamo model: an observation system simulation experiment for reconstructing meridional flow-speed
Authors:
Mausumi Dikpati,
Jeffrey L. Anderson,
Dhrubaditya Mitra
Abstract:
Accurate knowledge of time-variation in meridional flow-speed and profile is crucial for estimating a solar cycle's features, which are ultimately responsible for causing space climate variations. However, no consensus has been reached yet about the Sun's meridional circulation pattern observations and theories. By implementing an Ensemble Kalman Filter (EnKF) data assimilation in a Babcock-Leight…
▽ More
Accurate knowledge of time-variation in meridional flow-speed and profile is crucial for estimating a solar cycle's features, which are ultimately responsible for causing space climate variations. However, no consensus has been reached yet about the Sun's meridional circulation pattern observations and theories. By implementing an Ensemble Kalman Filter (EnKF) data assimilation in a Babcock-Leighton solar dynamo model using Data Assimilation Research Testbed (DART) framework, we find that the best reconstruction of time-variation in meridional flow-speed can be obtained when ten or more observations are used with an updating time of 15 days and a $\le 10\%$ observational error. Increasing ensemble-size from 16 to 160 improves reconstruction. Comparison of reconstructed flow-speed with "true-state" reveals that EnKF data assimilation is very powerful for reconstructing meridional flow-speeds and suggests that it can be implemented for reconstructing spatio-temporal patterns of meridional circulation.
△ Less
Submitted 21 August, 2014;
originally announced August 2014.
-
Evaluating potential for data assimilation in a flux-transport dynamo model by assessing sensitivity and response to meridional flow variation
Authors:
Mausumi Dikpati,
Jeffrey L. Anderson
Abstract:
We estimate here a flux-transport dynamo model's response time to changes in meridional flow speed. Time-variation in meridional flow primarily determines the shape of a cycle in this class of dynamo models. In order to simultaneously predict the shape, amplitude and timing of a solar cycle by implementing an Ensemble Kalman Filter in the framework of Data Assimilation Research Testbed (DART), it…
▽ More
We estimate here a flux-transport dynamo model's response time to changes in meridional flow speed. Time-variation in meridional flow primarily determines the shape of a cycle in this class of dynamo models. In order to simultaneously predict the shape, amplitude and timing of a solar cycle by implementing an Ensemble Kalman Filter in the framework of Data Assimilation Research Testbed (DART), it is important to know the model's sensitivity to flow variation. Guided by observations we consider a smooth increase or decrease in meridional flow speed for a specified time (a few months to a few years), after which the flow speed comes back to the steady speed, and implement that time-varying meridional flow at different phases of solar cycle. We find that the model's response time to change in flow speed peaks at four to six months if the flow change lasts for one year. The longer the changed flow lasts, the longer the model takes to respond. Magnetic diffusivity has no influence in model's response to flow variation as long as the dynamo operates in the advection-dominated regime. Experiments with more complex flow variations indicate that the shape and amplitude of flow-perturbation have no influence in the estimate of model's response time.
△ Less
Submitted 26 June, 2012;
originally announced June 2012.
-
Approximate Equations of Motion for Compact Spinning Bodies in General Relativity
Authors:
James L. Anderson
Abstract:
Approximate equations are derived for the motion of a gyroscope on the earth's gravitational field using the Einstein, Infeld, Hoffmann surface integral method. This method does not require a knowledge of the energy-momentum-stress tensor associated with the gyroscope and uses only its exterior field for its characterization. The resulting equations of motion differ from those of previous deriva…
▽ More
Approximate equations are derived for the motion of a gyroscope on the earth's gravitational field using the Einstein, Infeld, Hoffmann surface integral method. This method does not require a knowledge of the energy-momentum-stress tensor associated with the gyroscope and uses only its exterior field for its characterization. The resulting equations of motion differ from those of previous derivations.
△ Less
Submitted 21 November, 2005; v1 submitted 16 November, 2005;
originally announced November 2005.
-
Chaotic Behavior of Positronium in a Magnetic Field
Authors:
J. L. Anderson,
R. K. Murawski,
G. Schmidt
Abstract:
Classical motion of positronium embedded in a magnetic field is studied, and computation shows the emergence of chaotic orbits. Recent work investigating quantum behavior of this system predicts extremely long lifetimes. Chaos assisted tunneling, however may lead to significant shortening of the lifetime of this system.
Classical motion of positronium embedded in a magnetic field is studied, and computation shows the emergence of chaotic orbits. Recent work investigating quantum behavior of this system predicts extremely long lifetimes. Chaos assisted tunneling, however may lead to significant shortening of the lifetime of this system.
△ Less
Submitted 19 June, 2002;
originally announced June 2002.
-
Does General Relativity Require a Metric
Authors:
James L. Anderson
Abstract:
The nexus between the gravitational field and the space-time metric was an essential element in Einstein's development of General Relativity and led him to his discovery of the field equations for the gravitational field/metric. I will argue here that the metric is in fact an inessential element of this theory and can be dispensed with entirely. Its sole function in the theory was to describe th…
▽ More
The nexus between the gravitational field and the space-time metric was an essential element in Einstein's development of General Relativity and led him to his discovery of the field equations for the gravitational field/metric. I will argue here that the metric is in fact an inessential element of this theory and can be dispensed with entirely. Its sole function in the theory was to describe the space-time measurements made by ideal clocks and rods. However, the behavior of model clocks and measuring rods can be derived directly from the field equations of General Relativity using the Einstein-Infeld-Hoffmann (EIH) approiximation procedure. Therefore one does not need to introduce these ideal clocks and rods and hence has no need of a metric.
△ Less
Submitted 13 December, 1999;
originally announced December 1999.
-
Asymptotic conditions of motion for radiating charged particles
Authors:
James L. Anderson
Abstract:
Approximate asymptotic conditions on the motion of compact, electrically charged particles are derived within the framework of general relativity using the Einstein- Infeld-Hoffmann (EIH) surface integral method. While superficially similar to the Abraham-Lorentz and Lorentz-Dirac (ALD) equations of motion, these conditions differ from them in several fundamental ways. They are not equations of…
▽ More
Approximate asymptotic conditions on the motion of compact, electrically charged particles are derived within the framework of general relativity using the Einstein- Infeld-Hoffmann (EIH) surface integral method. While superficially similar to the Abraham-Lorentz and Lorentz-Dirac (ALD) equations of motion, these conditions differ from them in several fundamental ways. They are not equations of motion in the usual sense but rather a set of conditions which these motions must obey in the asymptotic future of an initial value surface. In addition to being asymptotic, these conditions of motion are approximate and apply, as do the original EIH equations, only to slowly moving systems. Also, they do not admit the run- away solutions of these other equations. As in the original EIH work, they are integrability conditions gotten from integrating the empty-space (i.e., source free) Einstein-Maxwell equations of general relativity over closed two-surfaces surrounding the sources of the fields governed by these equations. No additional ad hoc assumptions, such as the form of a force law or the introduction of inertial reaction terms, needed to derive the ALD equations are required for this purpose. Nor is there a need for any of the infinite mass renormalizations that are required in deriving these other equations.
△ Less
Submitted 21 September, 1997;
originally announced September 1997.
-
Thompson Scattering in an Expanding Universe
Authors:
James L. Anderson
Abstract:
The Thompson cross section for scattering of electromagnetic waves by a free electron in an expanding universe is derived here. The equations of motion of the electron are obtained using the Einstein-Infeld-Hoffmann (EIH) surface integral method. These integrals are evaluated approximately by perturbing off an Einstein-deSitter cosmological field. It is found that the Thompson cross section vari…
▽ More
The Thompson cross section for scattering of electromagnetic waves by a free electron in an expanding universe is derived here. The equations of motion of the electron are obtained using the Einstein-Infeld-Hoffmann (EIH) surface integral method. These integrals are evaluated approximately by perturbing off an Einstein-deSitter cosmological field. It is found that the Thompson cross section varies with time as the inverse square of the cosmic scale factor R(t).
△ Less
Submitted 14 September, 1997;
originally announced September 1997.